Publications

2019
Basu SS, McMinn MH, Giménez-Cassina Lopéz B, Regan MS, Randall EC, Clark AR, Cox CR, Agar NYR. Metal Oxide Laser Ionization Mass Spectrometry Imaging (MOLI MSI) Using Cerium(IV) Oxide. Anal Chem. 2019;91 (10) :6800-6807.Abstract
Matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI MSI) is a powerful technique for spatially resolved metabolomics. A variation on MALDI, termed metal oxide laser ionization (MOLI), capitalizes on the unique property of cerium(IV) oxide (CeO) to induce laser-catalyzed fatty acyl cleavage from lipids and has been utilized for bacterial identification. In this study, we present the development and utilization of CeO as an MSI catalyst. The method was developed using a MALDI TOF instrument in negative ion mode, equipped with a high frequency laser. Instrument parameters for MOLI MS fatty acid catalysis with CeO were optimized with phospholipid standards and fatty acid catalysis was confirmed using lipid extracts from reference bacterial strains, and sample preparation was optimized using mouse brain tissue. MOLI MSI was applied to the imaging of normal mouse brain revealing differentiable fatty acyl pools in myelinated and nonmyelinated regions. Similarly, MOLI MSI showed distinct fatty acyl composition in tumor regions of a patient derived xenograft mouse model of glioblastoma. To assess the potential of MOLI MSI to detect pathogens directly from tissue, a pseudoinfection model was prepared by spotting Escherichia coli lipid extracts on mouse brain tissue sections and imaged by MOLI MSI. The spotted regions were molecularly resolved from the supporting mouse brain tissue by the diagnostic odd-chained fatty acids and reflected control bacterial MOLI MS signatures. We describe MOLI MSI for the first time and highlight its potential for spatially resolved fatty acyl analysis, characterization of fatty acyl composition in tumors, and its potential for pathogen detection directly from tissue.
Panda A, Obmann VC, Lo W-C, Margevicius S, Jiang Y, Schluchter M, Patel IJ, Nakamoto D, Badve C, Griswold MA, et al. MR Fingerprinting and ADC Mapping for Characterization of Lesions in the Transition Zone of the Prostate Gland. Radiology. 2019;292 (3) :685-694.Abstract
BackgroundPreliminary studies have shown that MR fingerprinting-based relaxometry combined with apparent diffusion coefficient (ADC) mapping can be used to differentiate normal peripheral zone from prostate cancer and prostatitis. The utility of relaxometry and ADC mapping for the transition zone (TZ) is unknown.PurposeTo evaluate the utility of MR fingerprinting combined with ADC mapping for characterizing TZ lesions.Materials and MethodsTZ lesions that were suspicious for cancer in men who underwent MRI with T2-weighted imaging and ADC mapping ( values, 50-1400 sec/mm), MR fingerprinting with steady-state free precession, and targeted biopsy (60 in-gantry and 15 cognitive targeting) between September 2014 and August 2018 in a single university hospital were retrospectively analyzed. Two radiologists blinded to Prostate Imaging Reporting and Data System (PI-RADS) scores and pathologic diagnosis drew regions of interest on cancer-suspicious lesions and contralateral visually normal TZs (NTZs) on MR fingerprinting and ADC maps. Linear mixed models compared two-reader means of T1, T2, and ADC. Generalized estimating equations logistic regression analysis was used to evaluate both MR fingerprinting and ADC in differentiating NTZ, cancers and noncancers, clinically significant (Gleason score ≥ 7) cancers from clinically insignificant lesions (noncancers and Gleason 6 cancers), and characterizing PI-RADS version 2 category 3 lesions.ResultsIn 67 men (mean age, 66 years ± 8 [standard deviation]) with 75 lesions, targeted biopsy revealed 37 cancers (six PI-RADS category 3 cancers and 31 PI-RADS category 4 or 5 cancers) and 38 noncancers (31 PI-RADS category 3 lesions and seven PI-RADS category 4 or 5 lesions). The T1, T2, and ADC of NTZ (1800 msec ± 150, 65 msec ± 22, and [1.13 ± 0.19] × 10 mm/sec, respectively) were higher than those in cancers (1450 msec ± 110, 36 msec ± 11, and [0.57 ± 0.13] × 10 mm/sec, respectively; < .001 for all). The T1, T2, and ADC in cancers were lower than those in noncancers (1620 msec ± 120, 47 msec ± 16, and [0.82 ± 0.13] × 10 mm/sec, respectively; = .001 for T1 and ADC and = .03 for T2). The area under the receiver operating characteristic curve (AUC) for T1 plus ADC was 0.94 for separation. T1 and ADC in clinically significant cancers (1440 msec ± 140 and [0.58 ± 0.14] × 10 mm/sec, respectively) were lower than those in clinically insignificant lesions (1580 msec ± 120 and [0.75 ± 0.17] × 10 mm/sec, respectively; = .001 for all). The AUC for T1 plus ADC was 0.81 for separation. Within PI-RADS category 3 lesions, T1 and ADC of cancers (1430 msec ± 220 and [0.60 ± 0.17] × 10 mm/sec, respectively) were lower than those of noncancers (1630 msec ± 120 and [0.81 ± 0.13] × 10 mm/sec, respectively; = .006 for T1 and = .004 for ADC). The AUC for T1 was 0.79 for differentiating category 3 lesions.ConclusionMR fingerprinting-based relaxometry combined with apparent diffusion coefficient mapping may improve transition zone lesion characterization.© RSNA, 2019
Cheng C-C, Preiswerk F, Madore B. Multi-pathway Multi-echo Acquisition and Neural Contrast Translation to Generate a Variety of Quantitative and Qualitative Image Contrasts. Magn Reson Med. 2019.Abstract
PURPOSE: Clinical exams typically involve acquiring many different image contrasts to help discriminate healthy from diseased states. Ideally, 3D quantitative maps of all of the main MR parameters would be obtained for improved tissue characterization. Using data from a 7-min whole-brain multi-pathway multi-echo (MPME) scan, we aimed to synthesize several 3D quantitative maps (T and T ) and qualitative contrasts (MPRAGE, FLAIR, T -weighted, T -weighted, and proton density [PD]-weighted). The ability of MPME acquisitions to capture large amounts of information in a relatively short amount of time suggests it may help reduce the duration of neuro MR exams. METHODS: Eight healthy volunteers were imaged at 3.0T using a 3D isotropic (1.2 mm) MPME sequence. Spin-echo, MPRAGE, and FLAIR scans were performed for training and validation. MPME signals were interpreted through neural networks for predictions of different quantitative and qualitative contrasts. Predictions were compared to reference values at voxel and region-of-interest levels. RESULTS: Mean absolute errors (MAEs) for T and T maps were 216 ms and 11 ms, respectively. In ROIs containing white matter (WM) and thalamus tissues, the mean T /T predicted values were 899/62 ms and 1139/58 ms, consistent with reference values of 850/66 ms and 1126/58 ms, respectively. For qualitative contrasts, signals were normalized to those of WM, and MAEs for MPRAGE, FLAIR, T -weighted, T -weighted, and PD-weighted contrasts were 0.14, 0.15, 0.13, 0.16, and 0.05, respectively. CONCLUSIONS: Using an MPME sequence and neural-network contrast translation, whole-brain results were obtained with a variety of quantitative and qualitative contrast in ~6.8 min.
Zhou H, Zhang T, Jagadeesan J. Re-weighting and 1-Point RANSAC-Based P nP Solution to Handle Outliers. IEEE Trans Pattern Anal Mach Intell. 2019;41 (12) :3022-33.Abstract
The ability to handle outliers is essential for performing the perspective- n-point (P nP) approach in practical applications, but conventional RANSAC+P3P or P4P methods have high time complexities. We propose a fast P nP solution named R1PP nP to handle outliers by utilizing a soft re-weighting mechanism and the 1-point RANSAC scheme. We first present a P nP algorithm, which serves as the core of R1PP nP, for solving the P nP problem in outlier-free situations. The core algorithm is an optimal process minimizing an objective function conducted with a random control point. Then, to reduce the impact of outliers, we propose a reprojection error-based re-weighting method and integrate it into the core algorithm. Finally, we employ the 1-point RANSAC scheme to try different control points. Experiments with synthetic and real-world data demonstrate that R1PP nP is faster than RANSAC+P3P or P4P methods especially when the percentage of outliers is large, and is accurate. Besides, comparisons with outlier-free synthetic data show that R1PP nP is among the most accurate and fast P nP solutions, which usually serve as the final refinement step of RANSAC+P3P or P4P. Compared with REPP nP, which is the state-of-the-art P nP algorithm with an explicit outliers-handling mechanism, R1PP nP is slower but does not suffer from the percentage of outliers limitation as REPP nP.
Guenette JP, Ben-Shlomo N, Jayender J, Seethamraju RT, Kimbrell V, Tran N-A, Huang RY, Kim CJ, Kass JI, Corrales CE, et al. MR Imaging of the Extracranial Facial Nerve with the CISS Sequence. AJNR Am J Neuroradiol. 2019.Abstract
BACKGROUND AND PURPOSE: MR imaging is not routinely used to image the extracranial facial nerve. The purpose of this study was to determine the extent to which this nerve can be visualized with a CISS sequence and to determine the feasibility of using that sequence for locating the nerve relative to tumor. MATERIALS AND METHODS: Thirty-two facial nerves in 16 healthy subjects and 4 facial nerves in 4 subjects with parotid gland tumors were imaged with an axial CISS sequence protocol that included 0.8-mm isotropic voxels on a 3T MR imaging system with a 64-channel head/neck coil. Four observers independently segmented the 32 healthy subject nerves. Segmentations were compared by calculating average Hausdorff distance values and Dice similarity coefficients. RESULTS: The primary bifurcation of the extracranial facial nerve into the superior temporofacial and inferior cervicofacial trunks was visible on all 128 segmentations. The mean of the average Hausdorff distances was 1.2 mm (range, 0.3-4.6 mm). Dice coefficients ranged from 0.40 to 0.82. The relative position of the facial nerve to the tumor could be inferred in all 4 tumor cases. CONCLUSIONS: The facial nerve can be seen on CISS images from the stylomastoid foramen to the temporofacial and cervicofacial trunks, proximal to the parotid plexus. Use of a CISS protocol is feasible in the clinical setting to determine the location of the facial nerve relative to tumor.
Luo J, Sedghi A, Popuri K, Cobzas D, Zhang M, Preiswerk F, Toews M, Golby A, Sugiyama M, Wells WIIIM, et al. On the Applicability of Registration Uncertainty, in MICCAI 2019. Vol LNCS 11765. Shenzhen, China: Springer ; 2019 :410-9.Abstract
Estimating the uncertainty in (probabilistic) image registration enables, e.g., surgeons to assess the operative risk based on the trustworthiness of the registered image data. If surgeons receive inaccurately calculated registration uncertainty and misplace unwarranted confidence in the alignment solutions, severe consequences may result. For probabilistic image registration (PIR), the predominant way to quantify the registration uncertainty is using summary statistics of the distribution of transformation parameters. The majority of existing research focuses on trying out different summary statistics as well as means to exploit them. Distinctively, in this paper, we study two rarely examined topics: (1) whether those summary statistics of the transformation distribution most informatively represent the registration uncertainty; (2) Does utilizing the registration uncertainty always be beneficial. We show that there are two types of uncertainties: the transformation uncertainty, Ut, and label uncertainty Ul. The conventional way of using Ut to quantify Ul is inappropriate and can be misleading. By a real data experiment, we also share a potentially critical finding that making use of the registration uncertainty may not always be an improvement.
Luo MICCAI 2019
Miller K, Joldes GR, Bourantas G, Warfield SK, Hyde DE, Kikinis R, Wittek A. Biomechanical Modeling and Computer Simulation of the Brain during Neurosurgery. Int J Numer Method Biomed Eng. 2019 :e3250.Abstract
Computational biomechanics of the brain for neurosurgery is an emerging area of research recently gaining in importance and practical applications. This review paper presents the contributions of the Intelligent Systems for Medicine Laboratory and its collaborators to this field, discussing the modeling approaches adopted and the methods developed for obtaining the numerical solutions. We adopt a physics-based modeling approach and describe the brain deformation in mechanical terms (such as displacements, strains, and stresses), which can be computed using a biomechanical model, by solving a continuum mechanics problem. We present our modeling approaches related to geometry creation, boundary conditions, loading, and material properties. From the point of view of solution methods, we advocate the use of fully nonlinear modeling approaches, capable of capturing very large deformations and nonlinear material behavior. We discuss finite element and meshless domain discretization, the use of the total Lagrangian formulation of continuum mechanics, and explicit time integration for solving both time-accurate and steady-state problems. We present the methods developed for handling contacts and for warping 3D medical images using the results of our simulations. We present two examples to showcase these methods: brain shift estimation for image registration and brain deformation computation for neuronavigation in epilepsy treatment.
Fan G, Liu H, Wu Z, Li Y, Feng C, Wang D, Luo J, Wells WM, He S. Deep Learning-Based Automatic Segmentation of Lumbosacral Nerves on CT for Spinal Intervention: A Translational Study. AJNR Am J Neuroradiol. 2019;40 (6) :1074-81.Abstract
BACKGROUND AND PURPOSE: 3D reconstruction of a targeted area ("safe" triangle and Kambin triangle) may benefit the viability assessment of transforaminal epidural steroid injection, especially at the L5/S1 level. However, manual segmentation of lumbosacral nerves for 3D reconstruction is time-consuming. The aim of this study was to investigate the feasibility of deep learning-based segmentation of lumbosacral nerves on CT and the reconstruction of the safe triangle and Kambin triangle. MATERIALS AND METHODS: A total of 50 cases of spinal CT were manually labeled for lumbosacral nerves and bones using Slicer 4.8. The ratio of training/validation/testing was 32:8:10. A 3D U-Net was adopted to build the model SPINECT for automatic segmentations of lumbosacral structures. The Dice score, pixel accuracy, and Intersection over Union were computed to assess the segmentation performance of SPINECT. The areas of Kambin and safe triangles were measured to validate the 3D reconstruction. RESULTS: The results revealed successful segmentation of lumbosacral bone and nerve on CT. The average pixel accuracy for bone was 0.940, and for nerve, 0.918. The average Intersection over Union for bone was 0.897 and for nerve, 0.827. The Dice score for bone was 0.945, and for nerve, it was 0.905. There were no significant differences in the quantified Kambin triangle or safe triangle between manually segmented images and automatically segmented images ( > .05). CONCLUSIONS: Deep learning-based automatic segmentation of lumbosacral structures (nerves and bone) on routine CT is feasible, and SPINECT-based 3D reconstruction of safe and Kambin triangles is also validated.
Lemaire J-J, De Salles A, Coll G, El Ouadih Y, Chaix R, Coste J, Durif F, Makris N, Kikinis R. MRI Atlas of the Human Deep Brain. Front Neurol. 2019;10 :851.Abstract
Mastering detailed anatomy of the human deep brain in clinical neurosciences is challenging. Although numerous pioneering works have gathered a large dataset of structural and topographic information, it is still difficult to transfer this knowledge into practice, even with advanced magnetic resonance imaging techniques. Thus, classical histological atlases continue to be used to identify structures for stereotactic targeting in functional neurosurgery. Physicians mainly use these atlases as a template co-registered with the patient's brain. However, it is possible to directly identify stereotactic targets on MRI scans, enabling personalized targeting. In order to help clinicians directly identify deep brain structures relevant to present and future medical applications, we built a volumetric MRI atlas of the deep brain (MDBA) on a large scale (infra millimetric). Twelve hypothalamic, 39 subthalamic, 36 telencephalic, and 32 thalamic structures were identified, contoured, and labeled. Nineteen coronal, 18 axial, and 15 sagittal MRI plates were created. Although primarily designed for direct labeling, the anatomic space was also subdivided in twelfths of AC-PC distance, leading to proportional scaling in the coronal, axial, and sagittal planes. This extensive work is now available to clinicians and neuroscientists, offering another representation of the human deep brain ([https://hal.archives-ouvertes.fr/] [hal-02116633]). The atlas may also be used by computer scientists who are interested in deciphering the topography of this complex region.
Patel NA, Li G, Shang W, Wartenberg M, Heffter T, Burdette EC, Iordachita I, Tokuda J, Hata N, Tempany CM, et al. System Integration and Preliminary Clinical Evaluation of a Robotic System for MRI-Guided Transperineal Prostate Biopsy. J Med Robot Res. 2019;4 (2).Abstract
This paper presents the development, preclinical evaluation, and preliminary clinical study of a robotic system for targeted transperineal prostate biopsy under direct interventional magnetic resonance imaging (MRI) guidance. The clinically integrated robotic system is developed based on a modular design approach, comprised of surgical navigation application, robot control software, MRI robot controller hardware, and robotic needle placement manipulator. The system provides enabling technologies for MRI-guided procedures. It can be easily transported and setup for supporting the clinical workflow of interventional procedures, and the system is readily extensible and reconfigurable to other clinical applications. Preclinical evaluation of the system is performed with phantom studies in a 3 Tesla MRI scanner, rehearsing the proposed clinical workflow, and demonstrating an in-plane targeting error of 1.5mm. The robotic system has been approved by the institutional review board (IRB) for clinical trials. A preliminary clinical study is conducted with the patient consent, demonstrating the targeting errors at two biopsy target sites to be 4.0 and 3.7, which is sufficient to target a clinically significant tumor foci. First-in-human trials to evaluate the system's effectiveness and accuracy for MR image-guide prostate biopsy are underway.
Machado I, Toews M, George E, Unadkat P, Essayed W, Luo J, Teodoro P, Carvalho H, Martins J, Golland P, et al. Deformable MRI-Ultrasound Registration using Correlation-based Attribute Matching for Brain Shift Correction: Accuracy and Generality in Multi-site Data. Neuroimage. 2019;202 :116094.Abstract
Intraoperative tissue deformation, known as brain shift, decreases the benefit of using preoperative images to guide neurosurgery. Non-rigid registration of preoperative magnetic resonance (MR) to intraoperative ultrasound (US) has been proposed as a means to compensate for brain shift. We focus on the initial registration from MR to predurotomy US. We present a method that builds on previous work to address the need for accuracy and generality of MR-iUS registration algorithms in multi-site clinical data. To improve accuracy of registration, we use high-dimensional texture attributes instead of image intensities and propose to replace the standard difference-based attribute matching with correlation-based attribute matching. We also present a strategy that deals explicitly with the large field-of-view mismatch between MR and iUS images. We optimize key parameters across independent MR-iUS brain tumor datasets acquired at three different institutions, with a total of 43 tumor patients and 758 corresponding landmarks to validate the registration algorithm. Despite differences in imaging protocols, patient demographics and landmark distributions, our algorithm was able to reduce landmark errors prior to registration in three data sets (5.37 ± 4.27, 4.18 ± 1.97 and 6.18 ± 3.38 mm, respectively) to a consistently low level (2.28 ± 0.71, 2.08 ± 0.37 and 2.24 ± 0.78 mm, respectively). Our algorithm is compared to 15 other algorithms that have been previously tested on MR-iUS registration and it is competitive with the state-of-the-art on multiple datasets. We show that our algorithm has one of the lowest errors in all datasets (accuracy), and this is achieved while sticking to a fixed set of parameters for multi-site data (generality). In contrast, other algorithms/tools of similar performance need per-dataset parameter tuning (high accuracy but lower generality), and those that stick to fixed parameters have larger errors or inconsistent performance (generality but not the top accuracy). We further characterized landmark errors according to brain regions and tumor types, a topic so far missing in the literature. We found that landmark errors were higher in high-grade than low-grade glioma patients, and higher in tumor regions than in other brain regions.
Frisken S, Luo M, Juvekar P, Bunevicius A, Machado I, Unadkat P, Bertotti MM, Toews M, Wells WM, Miga MI, et al. A Comparison of Thin-Plate Spline Deformation and Finite Element Modeling to Compensate for Brain Shift during Tumor Resection. Int J Comput Assist Radiol Surg. 2019.Abstract
PURPOSE: Brain shift during tumor resection can progressively invalidate the accuracy of neuronavigation systems and affect neurosurgeons' ability to achieve optimal resections. This paper compares two methods that have been presented in the literature to compensate for brain shift: a thin-plate spline deformation model and a finite element method (FEM). For this comparison, both methods are driven by identical sparse data. Specifically, both methods are driven by displacements between automatically detected and matched feature points from intraoperative 3D ultrasound (iUS). Both methods have been shown to be fast enough for intraoperative brain shift correction (Machado et al. in Int J Comput Assist Radiol Surg 13(10):1525-1538, 2018; Luo et al. in J Med Imaging (Bellingham) 4(3):035003, 2017). However, the spline method requires no preprocessing and ignores physical properties of the brain while the FEM method requires significant preprocessing and incorporates patient-specific physical and geometric constraints. The goal of this work was to explore the relative merits of these methods on recent clinical data. METHODS: Data acquired during 19 sequential tumor resections in Brigham and Women's Hospital's Advanced Multi-modal Image-Guided Operating Suite between December 2017 and October 2018 were considered for this retrospective study. Of these, 15 cases and a total of 24 iUS to iUS image pairs met inclusion requirements. Automatic feature detection (Machado et al. in Int J Comput Assist Radiol Surg 13(10):1525-1538, 2018) was used to detect and match features in each pair of iUS images. Displacements between matched features were then used to drive both the spline model and the FEM method to compensate for brain shift between image acquisitions. The accuracies of the resultant deformation models were measured by comparing the displacements of manually identified landmarks before and after deformation. RESULTS: The mean initial subcortical registration error between preoperative MRI and the first iUS image averaged 5.3 ± 0.75 mm. The mean subcortical brain shift, measured using displacements between manually identified landmarks in pairs of iUS images, was 2.5 ± 1.3 mm. Our results showed that FEM was able to reduce subcortical registration error by a small but statistically significant amount (from 2.46 to 2.02 mm). A large variability in the results of the spline method prevented us from demonstrating either a statistically significant reduction in subcortical registration error after applying the spline method or a statistically significant difference between the results of the two methods. CONCLUSIONS: In this study, we observed less subcortical brain shift than has previously been reported in the literature (Frisken et al., in: Miller (ed) Biomechanics of the brain, Springer, Cham, 2019). This may be due to the fact that we separated out the initial misregistration between preoperative MRI and the first iUS image from our brain shift measurements or it may be due to modern neurosurgical practices designed to reduce brain shift, including reduced craniotomy sizes and better control of intracranial pressure with the use of mannitol and other medications. It appears that the FEM method and its use of geometric and biomechanical constraints provided more consistent brain shift correction and better correction farther from the driving feature displacements than the simple spline model. The spline-based method was simpler and tended to give better results for small deformations. However, large variability in the spline results and relatively small brain shift prevented this study from demonstrating a statistically significant difference between the results of the two methods.
Canalini L, Klein J, Miller D, Kikinis R. Segmentation-based Registration of Ultrasound Volumes for Glioma Resection in Image-guided Neurosurgery. Int J Comput Assist Radiol Surg. 2019;14 (10) :1697-1713.Abstract
PURPOSE: In image-guided surgery for glioma removal, neurosurgeons usually plan the resection on images acquired before surgery and use them for guidance during the subsequent intervention. However, after the surgical procedure has begun, the preplanning images become unreliable due to the brain shift phenomenon, caused by modifications of anatomical structures and imprecisions in the neuronavigation system. To obtain an updated view of the resection cavity, a solution is to collect intraoperative data, which can be additionally acquired at different stages of the procedure in order to provide a better understanding of the resection. A spatial mapping between structures identified in subsequent acquisitions would be beneficial. We propose here a fully automated segmentation-based registration method to register ultrasound (US) volumes acquired at multiple stages of neurosurgery. METHODS: We chose to segment sulci and falx cerebri in US volumes, which remain visible during resection. To automatically segment these elements, first we trained a convolutional neural network on manually annotated structures in volumes acquired before the opening of the dura mater and then we applied it to segment corresponding structures in different surgical phases. Finally, the obtained masks are used to register US volumes acquired at multiple resection stages. RESULTS: Our method reduces the mean target registration error (mTRE) between volumes acquired before the opening of the dura mater and during resection from 3.49 mm (± 1.55 mm) to 1.36 mm (± 0.61 mm). Moreover, the mTRE between volumes acquired before opening the dura mater and at the end of the resection is reduced from 3.54 mm (± 1.75 mm) to 2.05 mm (± 1.12 mm). CONCLUSION: The segmented structures demonstrated to be good candidates to register US volumes acquired at different neurosurgical phases. Therefore, our solution can compensate brain shift in neurosurgical procedures involving intraoperative US data.
Zaffino P, Pernelle G, Mastmeyer A, Mehrtash A, Zhang H, Kikinis R, Kapur T, Spadea MF. Fully Automatic Catheter Segmentation in MRI with 3D Convolutional Neural Networks: Application to MRI-guided Gynecologic Brachytherapy. Phys Med Biol. 2019;64 (16) :165008.Abstract
External-beam radiotherapy followed by High Dose Rate (HDR) brachytherapy is the standard-of-care for treating gynecologic cancers. The enhanced soft-tissue contrast provided by Magnetic Resonance Imaging (MRI) makes it a valuable imaging modality for diagnosing and treating these cancers. However, in contrast to Computed Tomography (CT) imaging, the appearance of the brachytherapy catheters, through which radiation sources are inserted to reach the cancerous tissue later on, is often variable across images. This paper reports, for the first time, a new deep-learning-based method for fully automatic segmentation of multiple closely spaced brachytherapy catheters in intraoperative MRI. &#13; Represented in the data are 50 gynecologic cancer patients treated by MRI-guided HDR brachytherapy. For each patient, a single intraoperative MRI was used. 826 catheters in the images were manually segmented by an expert radiation physicist who is also a trained radiation oncologist. The number of catheters in a patient ranged between 10 and 35. A deep 3-dimensional Convolutional Neural Network (CNN) model was developed and trained. In order to make the learning process more robust, the network was trained 5 times, each time using a different combination of shown patients. Finally, each test case was processed by the 5 networks and the final segmentation was generated by voting on the obtained 5 candidate segmentations. 4-fold validation was executed and all the patients were segmented.&#13; An average distance error of 2.0±3.4 mm was achieved. False positive and false negative catheters were 6.7% and 1.5% respectively. Average Dice score was equal to 0.60±0.17. The algorithm is available for use in the open source software platform 3D Slicer allowing for wide scale testing and research discussion. In conclusion, to the best of our knowledge, fully automatic segmentation of multiple closely spaced catheters from intraoperative MR images was achieved for the first time in gynecological brachytherapy.
Schwier M, van Griethuysen J, Vangel MG, Pieper S, Peled S, Tempany C, Aerts HJWL, Kikinis R, Fennessy FM, Fedorov A. Repeatability of Multiparametric Prostate MRI Radiomics Features. Sci Rep. 2019;9 (1) :9441.Abstract
In this study we assessed the repeatability of radiomics features on small prostate tumors using test-retest Multiparametric Magnetic Resonance Imaging (mpMRI). The premise of radiomics is that quantitative image-based features can serve as biomarkers for detecting and characterizing disease. For such biomarkers to be useful, repeatability is a basic requirement, meaning its value must remain stable between two scans, if the conditions remain stable. We investigated repeatability of radiomics features under various preprocessing and extraction configurations including various image normalization schemes, different image pre-filtering, and different bin widths for image discretization. Although we found many radiomics features and preprocessing combinations with high repeatability (Intraclass Correlation Coefficient > 0.85), our results indicate that overall the repeatability is highly sensitive to the processing parameters. Neither image normalization, using a variety of approaches, nor the use of pre-filtering options resulted in consistent improvements in repeatability. We urge caution when interpreting radiomics features and advise paying close attention to the processing configuration details of reported results. Furthermore, we advocate reporting all processing details in radiomics studies and strongly recommend the use of open source implementations.
Aryal M, Papademetriou I, Zhang Y-Z, Power C, McDannold N, Porter T. MRI Monitoring and Quantification of Ultrasound-Mediated Delivery of Liposomes Dually Labeled with Gadolinium and Fluorophore through the Blood-Brain Barrier. Ultrasound Med Biol. 2019;45 (7) :1733-42.Abstract
Magnetic resonance image-guided focused ultrasound has emerged as a viable non-invasive technique for the treatment of central nervous system-related diseases/disorders. Application of mechanical and thermal effects associated with focused transcranial ultrasound has been studied extensively in pre-clinical models, which has paved the way for clinical trials. However, in vivo treatment evaluation techniques on drug delivery application via blood-brain barrier opening has not been fully explored. Current treatment evaluation techniques via magnetic resonance imaging are hindered by systemic toxicity resulting from free gadolinium delivery. Here we propose a novel treatment evaluation strategy to overcome limitations by (i) synthesizing liposomes that are dually labeled with gadolinium, a magnetic resonance imaging (MRI) contrast agent, and rhodamine, a fluorophore; (ii) applying a focused ultrasound (FUS)-mediated BBB opening technique to deliver the liposomes across vascular barriers, achieving local gadolinium enhancement while reducing systemic and unwanted regional toxic effects associated with free gadolinium; and (iii) utilizing the MRI modality to confirm the delivery as it is already included in the FUS treatment in clinic. Liposomes were secondarily labeled with a fluorescent marker to confirm results obtained by MRI quantification postmortem. Two different sizes, 77.5 nm (group A) and 140 nm (group B), of gadolinium- and fluorescence-labeled liposomes were fabricated using thin-film hydration followed by extrusion methods and determined their stability up to 6 h under physiologic conditions. Gadolinium signal was detected on contrast-enhanced T1-weighted MRI 5 h after the delivery of liposomes via the BBB opening approach with an ultrasound pulse of 0.42 MPa (estimate in water) combined with microbubbles. MRI contrast was enhanced significantly in sonicated regions compared with non-sonicated regions of the brain. This was due to the accumulation of labeled liposomes, which was confirmed by detection of rhodamine fluorescence in histologic sections. The relative increase in MRI signal intensity was greater for smaller liposomes (mean diameter = 77.5 nm) than larger liposomes (mean diameter = 140 nm), which suggested a greater accumulation of the smaller liposomes in the brain after ultrasound-mediated opening of the BBB. Our findings suggest that the dual-labeled nanocarrier platform can be established, the FUS-mediated BBB opening approach can be used to deliver it through vascular barriers and MRI can be used to evaluate the extent of nanocarrier delivery.
Basu SS, Regan MS, Randall EC, Abdelmoula WM, Clark AR, Gimenez-Cassina Lopez B, Cornett DS, Haase A, Santagata S, Agar NYR. Rapid MALDI Mass Spectrometry Imaging for Surgical Pathology. NPJ Precis Oncol. 2019;3 :17.Abstract
Matrix assisted laser desorption ionization mass spectrometry imaging (MALDI MSI) is an emerging analytical technique, which generates spatially resolved proteomic and metabolomic images from tissue specimens. Conventional MALDI MSI processing and data acquisition can take over 30 min, limiting its clinical utility for intraoperative diagnostics. We present a rapid MALDI MSI method, completed under 5 min, including sample preparation and analysis, providing a workflow compatible with the clinical frozen section procedure.
Abdelmoula WM, Regan MS, Lopez BGC, Randall EC, Lawler S, Mladek AC, Nowicki MO, Marin BM, Agar JN, Swanson KR, et al. Automatic 3D Nonlinear Registration of Mass Spectrometry Imaging and Magnetic Resonance Imaging Data. Anal Chem. 2019;91 (9) :6206-16.Abstract
Multimodal integration between mass spectrometry imaging (MSI) and radiology-established modalities such as magnetic resonance imaging (MRI) would allow the investigations of key questions in complex biological systems such as the central nervous system. Such integration would provide complementary multiscale data to bridge the gap between molecular and anatomical phenotypes, potentially revealing new insights into molecular mechanisms underlying anatomical pathologies presented on MRI. Automatic coregistration between 3D MSI/MRI is a computationally challenging process due to dimensional complexity, MSI data sparsity, lack of direct spatial-correspondences, and nonlinear tissue deformation. Here, we present a new computational approach based on stochastic neighbor embedding to nonlinearly align 3D MSI to MRI data, identify and reconstruct biologically relevant molecular patterns in 3D, and fuse the MSI datacube to the MRI space. We demonstrate our method using multimodal high-spectral resolution matrix-assisted laser desorption ionization (MALDI) 9.4 T MSI and 7 T in vivo MRI data, acquired from a patient-derived, xenograft mouse brain model of glioblastoma following administration of the EGFR inhibitor drug of Erlotinib. Results show the distribution of some identified molecular ions of the EGFR inhibitor erlotinib, a phosphatidylcholine lipid, and cholesterol, which were reconstructed in 3D and mapped to the MRI space. The registration quality was evaluated on two normal mouse brains using the Dice coefficient for the regions of brainstem, hippocampus, and cortex. The method is generic and can therefore be applied to hyperspectral images from different mass spectrometers and integrated with other established in vivo imaging modalities such as computed tomography (CT) and positron emission tomography (PET).
Herz C, MacNeil K, Behringer PA, Tokuda J, Mehrtash A, Mousavi P, Kikinis R, Fennessy FM, Tempany CM, Tuncali K, et al. Open Source Platform for Transperineal In-bore MRI-guided Targeted Prostate Biopsy. IEEE Trans Biomed Eng. 2019.Abstract
OBJECTIVE: Accurate biopsy sampling of the suspected lesions is critical for the diagnosis and clinical management of prostate cancer (PCa). Transperineal in-bore MRI-guided prostate biopsy (tpMRgBx) is a targeted biopsy technique that was shown to be safe, efficient and accurate. Our goal was to develop an open source software platform to support evaluation, refinement and translation of this biopsy approach. METHODS: We developed SliceTracker, a 3D Slicer extension to support tpMRgBx. We followed modular design of the implementation to enable customization of the interface, and interchange of image segmentation and registration components to assess their effect on the processing time, precision and accuracy of the biopsy needle placement. The platform and supporting documentation were developed to enable the use of software by an operator with minimal technical training to facilitate translation. Retrospective evaluation studied registration accuracy, effect of the prostate segmentation approach, and re-identification time of biopsy targets. Prospective evaluation focused on the total procedure time and biopsy targeting error (BTE). RESULTS: Evaluation utilized data from 73 retrospective and 10 prospective tpMRgBx cases. Mean Landmark Registration Error (LRE) for retrospective evaluation was 1.88 ±2.63 mm and was not sensitive to the approach used for prostate gland segmentation. Prospectively, we observed target re-identification time of 4.60 ±2.40 min, and BTE of 2.40 ±0.98 mm. CONCLUSION: SliceTracker is modular and extensible open source platform for supporting image processing aspects of the tpMRgBx procedure. It has been successfully utilized to support clinical research procedures at our site.
Gao Y, Takagi K, Kato T, Shono N, Hata N. Continuum Robot with Follow the Leader Motion for Endoscopic Third Ventriculostomy and Tumor Biopsy. IEEE Trans Biomed Eng. 2019.Abstract
[Background] In a combined endoscopic third ventriculostomy (ETV) and endoscopic tumor biopsy (ETB) procedure, an optimal tool trajectory is mandatory to minimize trauma to surrounding cerebral tissue. [Objective] This paper presents wire-driven multi-section robot with push-pull driving wire. The robot is tested to attain follow-the-leader (FTL) motion to place surgical instruments through narrow passages while minimizing the trauma to tissues. [Methods] A wire-driven continuum robot with six sub-sections was developed and its kinematic model to achieve FTL motion was proposed. An accuracy test to assess the robot's ability to attain FTL motion along a set of elementary curved trajectory was performed. We also used hydrocephalus ventricular model created from human subject data to generate five ETV/ETB trajectory and conducted a study assessing the accuracy of the FTL motion along these clinically desirable trajectories. [Results] In the test with elementary curved paths, the maximal deviation of the robot was increased from 0.47 mm at 30 degrees turn to 1.78 mm at 180 degrees in a simple C-shaped curve. S-shaped FTL motion had lesser deviation ranging from 0.16 mm to 0.18 mm. In the phantom study, the greatest tip-deviation was 1.45 mm, and the greatest path deviation was 1.23 mm. [Conclusion] We present the application of a continuum robot with FTL motion to perform a combined ETV/ETB procedure. The validation study using human subject data indicated that the accuracy of FTL motion is relatively high. It is expected that may be useful combined ETV and ETB.

Pages