Neurosurgery Research Publications

Yao S, Lin P, Vera M, Akter F, Zhang R-Y, Zeng A, Golby AJ, Xu G, Tie Y, Song J. Hormone Levels are Related to Functional Compensation in Prolactinomas: A Resting-state fMRI Study. J Neurol Sci. 2020;411 :116720.Abstract
Prolactinomas are tumors of the pituitary gland, which overproduces prolactin leading to dramatic fluctuations of endogenous hormone levels throughout the body. While it is not fully understood how endogenous hormone disorders affect a patient's brain, it is well known that fluctuating hormone levels can have negative neuropsychological effects. Using resting-state functional magnetic resonance imaging (rs-fMRI), we investigated whole-brain functional connectivity (FC) and its relationship with hormone levels in prolactinomas. By performing seed-based FC analyses, we compared FC metrics between 33 prolactinoma patients and 31 healthy controls matched for age, sex, and hand dominance. We then carried out a partial correlation analysis to examine the relationship between FC metrics and hormone levels. Compared to healthy controls, prolactinoma patients showed significantly increased thalamocortical and cerebellar-cerebral FC. Endogenous hormone levels were also positively correlated with increased FC metrics, and these hormone-FC relationships exhibited sex differences in prolactinoma patients. Our study is the first to reveal altered FC patterns in prolactinomas and to quantify the hormone-FC relationships. These results indicate the importance of endogenous hormones on functional compensation of the brain in patients with prolactinomas.
Randall EC, Lopez BGC, Peng S, Regan MS, Abdelmoula WM, Basu SS, Santagata S, Yoon H, Haigis MC, Agar JN, et al. Localized Metabolomic Gradients in Patient-Derived Xenograft Models of Glioblastoma. Cancer Res. 2020;80 (6) :1258-67.Abstract
Glioblastoma (GBM) is increasingly recognized as a disease involving dysfunctional cellular metabolism. GBMs are known to be complex heterogeneous systems containing multiple distinct cell populations and are supported by an aberrant network of blood vessels. A better understanding of GBM metabolism, its variation with respect to the tumor microenvironment, and resulting regional changes in chemical composition is required. This may shed light on the observed heterogeneous drug distribution, which cannot be fully described by limited or uneven disruption of the blood-brain barrier. In this work, we used mass spectrometry imaging (MSI) to map metabolites and lipids in patient-derived xenograft models of GBM. A data analysis workflow revealed that distinctive spectral signatures were detected from different regions of the intracranial tumor model. A series of long-chain acylcarnitines were identified and detected with increased intensity at the tumor edge. A 3D MSI dataset demonstrated that these molecules were observed throughout the entire tumor/normal interface and were not confined to a single plane. mRNA sequencing demonstrated that hallmark genes related to fatty acid metabolism were highly expressed in samples with higher acylcarnitine content. These data suggest that cells in the core and the edge of the tumor undergo different fatty acid metabolism, resulting in different chemical environments within the tumor. This may influence drug distribution through changes in tissue drug affinity or transport and constitute an important consideration for therapeutic strategies in the treatment of GBM. SIGNIFICANCE: GBM tumors exhibit a metabolic gradient that should be taken into consideration when designing therapeutic strategies for treatment..
Bunevicius A, McDannold NJ, Golby AJ. Focused Ultrasound Strategies for Brain Tumor Therapy. Oper Neurosurg (Hagerstown). 2020;19 (1) :9-18.Abstract
BACKGROUND: A key challenge in the medical treatment of brain tumors is the limited penetration of most chemotherapeutic agents across the blood-brain barrier (BBB) into the tumor and the infiltrative margin around the tumor. Magnetic resonance-guided focused ultrasound (MRgFUS) is a promising tool to enhance the delivery of chemotherapeutic agents into brain tumors. OBJECTIVE: To review the mechanism of FUS, preclinical evidence, and clinical studies that used low-frequency FUS for a BBB opening in gliomas. METHODS: Literature review. RESULTS: The potential of externally delivered low-intensity ultrasound for a temporally and spatially precise and predictable disruption of the BBB has been investigated for over a decade, yielding extensive preclinical literature demonstrating that FUS can disrupt the BBB in a spatially targeted and temporally reversible manner. Studies in animal models documented that FUS enhanced the delivery of numerous chemotherapeutic and investigational agents across the BBB and into brain tumors, including temozolomide, bevacizumab, 1,3-bis (2-chloroethyl)-1-nitrosourea, doxorubicin, viral vectors, and cells. Chemotherapeutic interventions combined with FUS slowed tumor progression and improved animal survival. Recent advances of MRgFUS systems allow precise, temporally and spatially controllable, and safe transcranial delivery of ultrasound energy. Initial clinical evidence in glioma patients has shown the efficacy of MRgFUS in disrupting the BBB, as demonstrated by an enhanced gadolinium penetration. CONCLUSION: Thus far, a temporary disruption of the BBB followed by the administration of chemotherapy has been both feasible and safe. Further studies are needed to determine the actual drug delivery, including the drug distribution at a tissue-level scale, as well as effects on tumor growth and patient prognosis.
Freedman RA, Gelman RS, Agar NYR, Santagata S, Randall EC, Gimenez-Cassina Lopez B, Connolly RM, Dunn IF, Van Poznak CH, Anders CK, et al. Pre- and Postoperative Neratinib for HER2-Positive Breast Cancer Brain Metastases: Translational Breast Cancer Research Consortium 022. Clin Breast Cancer. 2020;20 (2) :145-51.Abstract
PURPOSE: This pilot study was performed to test our ability to administer neratinib monotherapy before clinically recommended craniotomy in patients with HER2-positive metastatic breast cancer to the central nervous system, to examine neratinib's central nervous system penetration at craniotomy, and to examine postoperative neratinib maintenance. PATIENTS AND METHODS: Patients with HER2-positive brain metastases undergoing clinically indicated cranial resection of a parenchymal tumor received neratinib 240 mg orally once a day for 7 to 21 days preoperatively, and resumed therapy postoperatively in 28-day cycles. Exploratory evaluations of time to disease progression, survival, and correlative tissue, cerebrospinal fluid (CSF), and blood-based analyses examining neratinib concentrations were planned. The study was registered at under number NCT01494662. RESULTS: We enrolled 5 patients between May 22, 2013, and October 18, 2016. As of March 1, 2019, patients had remained on the study protocol for 1 to 75+ postoperative cycles pf therapy. Two patients had grade 3 diarrhea. Evaluation of the CSF showed low concentrations of neratinib; nonetheless, 2 patients continued to receive therapy without disease progression for at least 13 cycles, with one on-study treatment lasting for nearly 6 years. Neratinib distribution in surgical tissue was variable for 1 patient, while specimens from 2 others did not produce conclusive results as a result of limited available samples. CONCLUSION: Neratinib resulted in expected rates of diarrhea in this small cohort, with 2 of 5 patients receiving the study treatment for durable periods. Although logistically challenging, we were able to test a limited number of CSF- and parenchymal-based neratinib concentrations. Our findings from resected tumor tissue in one patient revealed heterogeneity in drug distribution and tumor histopathology.
Xie G, Zhang F, Leung L, Mooney MA, Epprecht L, Norton I, Rathi Y, Kikinis R, Al-Mefty O, Makris N, et al. Anatomical Assessment of Trigeminal Nerve Tractography Using Diffusion MRI: A Comparison of Acquisition B-Values and Single- and Multi-Fiber Tracking Strategies. Neuroimage Clin. 2020;25 :102160.Abstract
BACKGROUND: The trigeminal nerve (TGN) is the largest cranial nerve and can be involved in multiple inflammatory, compressive, ischemic or other pathologies. Currently, imaging-based approaches to identify the TGN mostly rely on T2-weighted magnetic resonance imaging (MRI), which provides localization of the cisternal portion of the TGN where the contrast between nerve and cerebrospinal fluid (CSF) is high enough to allow differentiation. The course of the TGN within the brainstem as well as anterior to the cisternal portion, however, is more difficult to display on traditional imaging sequences. An advanced imaging technique, diffusion MRI (dMRI), enables tracking of the trajectory of TGN fibers and has the potential to visualize anatomical regions of the TGN not seen on T2-weighted imaging. This may allow a more comprehensive assessment of the nerve in the context of pathology. To date, most work in TGN tracking has used clinical dMRI acquisitions with a b-value of 1000 s/mm and conventional diffusion tensor MRI (DTI) tractography methods. Though higher b-value acquisitions and multi-tensor tractography methods are known to be beneficial for tracking brain white matter fiber tracts, there have been no studies conducted to evaluate the performance of these advanced approaches on nerve tracking of the TGN, in particular on tracking different anatomical regions of the TGN. OBJECTIVE: We compare TGN tracking performance using dMRI data with different b-values, in combination with both single- and multi-tensor tractography methods. Our goal is to assess the advantages and limitations of these different strategies for identifying the anatomical regions of the TGN. METHODS: We proposed seven anatomical rating criteria including true and false positive structures, and we performed an expert rating study of over 1000 TGN visualizations, as follows. We tracked the TGN using high-quality dMRI data from 100 healthy adult subjects from the Human Connectome Project (HCP). TGN tracking performance was compared across dMRI acquisitions with b = 1000 s/mm, b = 2000 s/mm and b = 3000 s/mm, using single-tensor (1T) and two-tensor (2T) unscented Kalman filter (UKF) tractography. This resulted in a total of six tracking strategies. The TGN was identified using an anatomical region-of-interest (ROI) selection approach. First, in a subset of the dataset we identified ROIs that provided good TGN tracking performance across all tracking strategies. Using these ROIs, the TGN was then tracked in all subjects using the six tracking strategies. An expert rater (GX) visually assessed and scored each TGN based on seven anatomical judgment criteria. These criteria included the presence of multiple expected anatomical segments of the TGN (true positive structures), specifically branch-like structures, cisternal portion, mesencephalic trigeminal tract, and spinal cord tract of the TGN. False positive criteria included the presence of any fibers entering the temporal lobe, the inferior cerebellar peduncle, or the middle cerebellar peduncle. Expert rating scores were analyzed to compare TGN tracking performance across the six tracking strategies. Intra- and inter-rater validation was performed to assess the reliability of the expert TGN rating result. RESULTS: The TGN was selected using two anatomical ROIs (Meckel's Cave and cisternal portion of the TGN). The two-tensor tractography method had significantly better performance on identifying true positive structures, while generating more false positive streamlines in comparison to the single-tensor tractography method. TGN tracking performance was significantly different across the three b-values for almost all structures studied. Tracking performance was reported in terms of the percentage of subjects achieving each anatomical rating criterion. Tracking of the cisternal portion and branching structure of the TGN was generally successful, with the highest performance of over 98% using two-tensor tractography and b = 1000 or b = 2000. However, tracking the smaller mesencephalic and spinal cord tracts of the TGN was quite challenging (highest performance of 37.5% and 57.07%, using two-tensor tractography with b = 1000 and b = 2000, respectively). False positive connections to the temporal lobe (over 38% of subjects for all strategies) and cerebellar peduncles (100% of subjects for all strategies) were prevalent. High joint probability of agreement was obtained in the inter-rater (on average 83%) and intra-rater validation (on average 90%), showing a highly reliable expert rating result. CONCLUSIONS: Overall, the results of the study suggest that researchers and clinicians may benefit from tailoring their acquisition and tracking methodology to the specific anatomical portion of the TGN that is of the greatest interest. For example, tracking of branching structures and TGN-T2 overlap can be best achieved with a two-tensor model and an acquisition using b = 1000 or b = 2000. In general, b = 1000 and b = 2000 acquisitions provided the best-rated tracking results. Further research is needed to improve both sensitivity and specificity of the depiction of the TGN anatomy using dMRI.
Rigolo L, Essayed WI, Tie Y, Norton I, Mukundan S, Golby A. Intraoperative Use of Functional MRI for Surgical Decision Making after Limited or Infeasible Electrocortical Stimulation Mapping. J Neuroimaging. 2020;30 (2) :184-91.Abstract
BACKGROUND AND PURPOSE: Functional magnetic resonance imaging (fMRI) is becoming widely recognized as a key component of preoperative neurosurgical planning, although intraoperative electrocortical stimulation (ECS) is considered the gold standard surgical brain mapping method. However, acquiring and interpreting ECS results can sometimes be challenging. This retrospective study assesses whether intraoperative availability of fMRI impacted surgical decision-making when ECS was problematic or unobtainable. METHODS: Records were reviewed for 191 patients who underwent presurgical fMRI with fMRI loaded into the neuronavigation system. Four patients were excluded as a bur-hole biopsy was performed. Imaging was acquired at 3 Tesla and analyzed using the general linear model with significantly activated pixels determined via individually determined thresholds. fMRI maps were displayed intraoperatively via commercial neuronavigation systems. RESULTS: Seventy-one cases were planned ECS; however, 18 (25.35%) of these procedures were either not attempted or aborted/limited due to: seizure (10), patient difficulty cooperating with the ECS mapping (4), scarring/limited dural opening (3), or dural bleeding (1). In all aborted/limited ECS cases, the surgeon continued surgery using fMRI to guide surgical decision-making. There was no significant difference in the incidence of postoperative deficits between cases with completed ECS and those with limited/aborted ECS. CONCLUSIONS: Preoperative fMRI allowed for continuation of surgery in over one-fourth of patients in which planned ECS was incomplete or impossible, without a significantly different incidence of postoperative deficits compared to the patients with completed ECS. This demonstrates additional value of fMRI beyond presurgical planning, as fMRI data served as a backup method to ECS.
Bunevicius A, Schregel K, Sinkus R, Golby A, Patz S. REVIEW: MR Elastography of Brain Tumors. Neuroimage Clin. 2020;25 :102109.Abstract
MR elastography allows non-invasive quantification of the shear modulus of tissue, i.e. tissue stiffness and viscosity, information that offers the potential to guide presurgical planning for brain tumor resection. Here, we review brain tumor MRE studies with particular attention to clinical applications. Studies that investigated MRE in patients with intracranial tumors, both malignant and benign as well as primary and metastatic, were queried from the Pubmed/Medline database in August 2018. Reported tumor and normal appearing white matter stiffness values were extracted and compared as a function of tumor histopathological diagnosis and MRE vibration frequencies. Because different studies used different elastography hardware, pulse sequences, reconstruction inversion algorithms, and different symmetry assumptions about the mechanical properties of tissue, effort was directed to ensure that similar quantities were used when making inter-study comparisons. In addition, because different methodologies and processing pipelines will necessarily bias the results, when pooling data from different studies, whenever possible, tumor values were compared with the same subject's contralateral normal appearing white matter to minimize any study-dependent bias. The literature search yielded 10 studies with a total of 184 primary and metastatic brain tumor patients. The group mean tumor stiffness, as measured with MRE, correlated with intra-operatively assessed stiffness of meningiomas and pituitary adenomas. Pooled data analysis showed significant overlap between shear modulus values across brain tumor types. When adjusting for the same patient normal appearing white matter shear modulus values, meningiomas were the stiffest tumor-type. MRE is increasingly being examined for potential in brain tumor imaging and might have value for surgical planning. However, significant overlap of shear modulus values between a number of different tumor types limits applicability of MRE for diagnostic purposes. Thus, further rigorous studies are needed to determine specific clinical applications of MRE for surgical planning, disease monitoring and molecular stratification of brain tumors.
Lee TC, Guenette JP, Moses ZB, Chi JH. MRI-Guided Cryoablation of Epidural Malignancies in the Spinal Canal Resulting in Neural Decompression and Regrowth of Bone. AJR Am J Roentgenol. 2019;212 (1) :205-8.Abstract
OBJECTIVE: The purpose of this article is to describe the use of MRI to safely monitor cryoablation for the treatment of spinal epidural malignancies. CONCLUSION: Use of MRI guidance to monitor percutaneous cryoablation allows ablation margins more distinct than those allowed by heat-based ablation modalities. MRI-guided cryoablation is a feasible option for treating epidural tumors involving the spinal canal, resulting in successful decompression of the tumor away from the spinal cord with regrowth of previously eroded bone around the spinal canal.
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-7.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.
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.
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.
Frisken S, Luo M, Machado I, Unadkat P, Juvekar P, Bunevicius A, Toews M, Wells WM, Miga MI, Golby AJ. Preliminary Results Comparing Thin Plate Splines with Finite Element Methods for Modeling Brain Deformation during Neurosurgery using Intraoperative Ultrasound. Proc SPIE Int Soc Opt Eng. 2019;10951.Abstract
Brain shift compensation attempts to model the deformation of the brain which occurs during the surgical removal of brain tumors to enable mapping of presurgical image data into patient coordinates during surgery and thus improve the accuracy and utility of neuro-navigation. We present preliminary results from clinical tumor resections that compare two methods for modeling brain deformation, a simple thin plate spline method that interpolates displacements and a more complex finite element method (FEM) that models physical and geometric constraints of the brain and its material properties. Both methods are driven by the same set of displacements at locations surrounding the tumor. These displacements were derived from sets of corresponding matched features that were automatically detected using the SIFT-Rank algorithm. The deformation accuracy was tested using a set of manually identified landmarks. The FEM method requires significantly more preprocessing than the spline method but both methods can be used to model deformations in the operating room in reasonable time frames. Our preliminary results indicate that the FEM deformation model significantly out-performs the spline-based approach for predicting the deformation of manual landmarks. While both methods compensate for brain shift, this work suggests that models that incorporate biophysics and geometric constraints may be more accurate.
Lampinen B, Szczepankiewicz F, Novén M, van Westen D, Hansson O, Englund E, Mårtensson J, Westin C-F, Nilsson M. Searching for the Neurite Density with Diffusion MRI: Challenges for Biophysical Modeling. Hum Brain Mapp. 2019;40 (8) :2529-45.Abstract
In vivo mapping of the neurite density with diffusion MRI (dMRI) is a high but challenging aim. First, it is unknown whether all neurites exhibit completely anisotropic ("stick-like") diffusion. Second, the "density" of tissue components may be confounded by non-diffusion properties such as T2 relaxation. Third, the domain of validity for the estimated parameters to serve as indices of neurite density is incompletely explored. We investigated these challenges by acquiring data with "b-tensor encoding" and multiple echo times in brain regions with low orientation coherence and in white matter lesions. Results showed that microscopic anisotropy from b-tensor data is associated with myelinated axons but not with dendrites. Furthermore, b-tensor data together with data acquired for multiple echo times showed that unbiased density estimates in white matter lesions require data-driven estimates of compartment-specific T2 values. Finally, the "stick" fractions of different biophysical models could generally not serve as neurite density indices across the healthy brain and white matter lesions, where outcomes of comparisons depended on the choice of constraints. In particular, constraining compartment-specific T2 values was ambiguous in the healthy brain and had a large impact on estimated values. In summary, estimating neurite density generally requires accounting for different diffusion and/or T2 properties between axons and dendrites. Constrained "index" parameters could be valid within limited domains that should be delineated by future studies.
Zhang F, Ning L, O'Donnell LJ, Pasternak O. MK-curve - Characterizing the Relation between Mean Kurtosis and Alterations in the Diffusion MRI Signal. Neuroimage. 2019;196 :68-80.Abstract
Diffusion kurtosis imaging (DKI) is a diffusion MRI (dMRI) technique to quantify brain microstructural properties. While DKI measures are sensitive to tissue alterations, they are also affected by signal alterations caused by imaging artifacts such as noise, motion and Gibbs ringing. Consequently, DKI often yields output parameter values (e.g. mean kurtosis; MK) that are implausible. These include implausible values that are outside of the range dictated by physics/biology, and visually apparent implausible values that form unexpected discontinuities, being too high or too low comparing with their neighborhood. These implausible values will introduce bias into any following data analyses (e.g. between-population statistical computation). Existing studies have attempted to correct implausible DKI parameter values in multiple ways; however, these approaches are not always effective. In this study, we propose a novel method for detecting and correcting voxels with implausible values to enable improved DKI parameter estimation. In particular, we focus on MK parameter estimation. We first characterize the relation between MK and alterations in the dMRI signal including diffusion weighted images (DWIs) and the baseline (b0) images. This is done by calculating MK for a range of synthetic DWI or b0 for each voxel, and generating curves (MK-curve) representing how alterations to the input dMRI signals affect the resulting output MK. We find that voxels with implausible MK values are more likely caused by artifacts in the b0 images than artifacts in DWIs with higher b-values. Accordingly, two characteristic b0 values, which define a range of synthetic b0 values that generate implausible MK values, are identified on the MK-curve. Based on this characterization, we propose an automatic approach for detection of voxels with implausible MK values by comparing a voxel's original b0 signal to the identified two characteristic b0 values, along with a correction strategy to replace the original b0 in each detected implausible voxel with a synthetic b0 value computed from the MK-curve. We evaluate the method on a DKI phantom dataset and dMRI datasets from the Human Connectome Project (HCP), and we compare the proposed correction method with other previously proposed correction methods. Results show that our proposed method is able to identify and correct most voxels with implausible DKI parameter values as well as voxels with implausible diffusion tensor parameter values.
Nitsch J, Klein J, Dammann P, Wrede K, Gembruch O, Moltz JH, Meine H, Sure U, Kikinis R, Miller D. Automatic and Efficient MRI-US Segmentations for Improving Intraoperative Image Fusion in Image-guided Neurosurgery. Neuroimage Clin. 2019;22 :101766.Abstract
Knowledge of the exact tumor location and structures at risk in its vicinity are crucial for neurosurgical interventions. Neuronavigation systems support navigation within the patient's brain, based on preoperative MRI (preMRI). However, increasing tissue deformation during the course of tumor resection reduces navigation accuracy based on preMRI. Intraoperative ultrasound (iUS) is therefore used as real-time intraoperative imaging. Registration of preMRI and iUS remains a challenge due to different or varying contrasts in iUS and preMRI. Here, we present an automatic and efficient segmentation of B-mode US images to support the registration process. The falx cerebri and the tentorium cerebelli were identified as examples for central cerebral structures and their segmentations can serve as guiding frame for multi-modal image registration. Segmentations of the falx and tentorium were performed with an average Dice coefficient of 0.74 and an average Hausdorff distance of 12.2 mm. The subsequent registration incorporates these segmentations and increases accuracy, robustness and speed of the overall registration process compared to purely intensity-based registration. For validation an expert manually located corresponding landmarks. Our approach reduces the initial mean Target Registration Error from 16.9 mm to 3.8 mm using our intensity-based registration and to 2.2 mm with our combined segmentation and registration approach. The intensity-based registration reduced the maximum initial TRE from 19.4 mm to 5.6 mm, with the approach incorporating segmentations this is reduced to 3.0 mm. Mean volumetric intensity-based registration of preMRI and iUS took 40.5 s, including segmentations 12.0 s.
Zhang F, Wu Y, Norton I, Rathi Y, Golby AJ, O'Donnell LJ. Test-retest Reproducibility of White Matter Parcellation using Diffusion MRI Tractography Fiber Clustering. Hum Brain Mapp. 2019;40 (10) :3041-57.Abstract
There are two popular approaches for automated white matter parcellation using diffusion MRI tractography, including fiber clustering strategies that group white matter fibers according to their geometric trajectories and cortical-parcellation-based strategies that focus on the structural connectivity among different brain regions of interest. While multiple studies have assessed test-retest reproducibility of automated white matter parcellations using cortical-parcellation-based strategies, there are no existing studies of test-retest reproducibility of fiber clustering parcellation. In this work, we perform what we believe is the first study of fiber clustering white matter parcellation test-retest reproducibility. The assessment is performed on three test-retest diffusion MRI datasets including a total of 255 subjects across genders, a broad age range (5-82 years), health conditions (autism, Parkinson's disease and healthy subjects), and imaging acquisition protocols (three different sites). A comprehensive evaluation is conducted for a fiber clustering method that leverages an anatomically curated fiber clustering white matter atlas, with comparison to a popular cortical-parcellation-based method. The two methods are compared for the two main white matter parcellation applications of dividing the entire white matter into parcels (i.e., whole brain white matter parcellation) and identifying particular anatomical fiber tracts (i.e., anatomical fiber tract parcellation). Test-retest reproducibility is measured using both geometric and diffusion features, including volumetric overlap (wDice) and relative difference of fractional anisotropy. Our experimental results in general indicate that the fiber clustering method produced more reproducible white matter parcellations than the cortical-parcellation-based method.
Randall EC, Zadra G, Chetta P, Lopez BGC, Syamala S, Basu SS, Agar JN, Loda M, Tempany CM, Fennessy FM, et al. Molecular Characterization of Prostate Cancer with Associated Gleason Score using Mass Spectrometry Imaging. Mol Cancer Res. 2019;17 (5) :1155-65.Abstract
Diagnosis of prostate cancer is based on histological evaluation of tumor architecture using a system known as the 'Gleason score'. This diagnostic paradigm, while the standard of care, is time-consuming, shows intra-observer variability and provides no information about the altered metabolic pathways, which result in altered tissue architecture. Characterization of the molecular composition of prostate cancer and how it changes with respect to the Gleason score (GS) could enable a more objective and faster diagnosis. It may also aid in our understanding of disease onset and progression. In this work, we present mass spectrometry imaging for identification and mapping of lipids and metabolites in prostate tissue from patients with known prostate cancer with GS from 6 to 9. A gradient of changes in the intensity of various lipids was observed, which correlated with increasing GS. Interestingly, these changes were identified in both regions of high tumor cell density, and in regions of tissue that appeared histologically benign, possibly suggestive of pre-cancerous metabolomic changes. A total of 31 lipids, including several phosphatidylcholines, phosphatidic acids, phosphatidylserines, phosphatidylinositols and cardiolipins were detected with higher intensity in GS (4+3) compared with GS (3+4), suggesting they may be markers of prostate cancer aggression. Results obtained through mass spectrometry imaging studies were subsequently correlated with a fast, ambient mass spectrometry method for potential use as a clinical tool to support image-guided prostate biopsy. Implications: In this study we suggest that metabolomic differences between prostate cancers with different Gleason scores can be detected by mass spectrometry imaging.
Stojanovski S, Felsky D, Viviano JD, Shahab S, Bangali R, Burton CL, Devenyi GA, O'Donnell LJ, Szatmari P, Chakravarty MM, et al. Polygenic Risk and Neural Substrates of Attention-Deficit/Hyperactivity Disorder Symptoms in Youths With a History of Mild Traumatic Brain Injury. Biol Psychiatry. 2019;85 (5) :408-16.Abstract
BACKGROUND: Attention-deficit/hyperactivity disorder (ADHD) is a major sequela of traumatic brain injury (TBI) in youths. The objective of this study was to examine whether ADHD symptoms are differentially associated with genetic risk and brain structure in youths with and without a history of TBI. METHODS: Medical history, ADHD symptoms, genetic data, and neuroimaging data were obtained from a community sample of youths. ADHD symptom severity was compared between those with and without TBI (TBI n = 418, no TBI n = 3193). The relationship of TBI history, genetic vulnerability, brain structure, and ADHD symptoms was examined by assessing 1) ADHD polygenic score (discovery sample ADHD n = 19,099, control sample n = 34,194), 2) basal ganglia volumes, and 3) fractional anisotropy in the corpus callosum and corona radiata. RESULTS: Youths with TBI reported greater ADHD symptom severity compared with those without TBI. Polygenic score was positively associated with ADHD symptoms in youths without TBI but not in youths with TBI. The negative association between the caudate volume and ADHD symptoms was not moderated by a history of TBI. However, the relationship between ADHD symptoms and structure of the genu of the corpus callosum was negative in youths with TBI and positive in youths without TBI. CONCLUSIONS: The identification of distinct ADHD etiology in youths with TBI provides neurobiological insight into the clinical heterogeneity in the disorder. Results indicate that genetic predisposition to ADHD does not increase the risk for ADHD symptoms associated with TBI. ADHD symptoms associated with TBI may be a result of a mechanical insult rather than neurodevelopmental factors.
Unadkat P, Fumagalli L, Rigolo L, Vangel MG, Young GS, Huang R, Mukundan S, Golby A, Tie Y. Functional MRI Task Comparison for Language Mapping in Neurosurgical Patients. J Neuroimaging. 2019;29 (3) :348-56.Abstract
BACKGROUND AND PURPOSE: Language task-based functional MRI (fMRI) is increasingly used for presurgical planning in patients with brain lesions. Different paradigms elicit activations of different components of the language network. The aim of this study is to optimize fMRI clinical usage by comparing the effectiveness of three language tasks for language lateralization and localization in a large group of patients with brain lesions. METHODS: We analyzed fMRI data from a sequential retrospective cohort of 51 patients with brain lesions who underwent presurgical fMRI language mapping. We compared the effectiveness of three language tasks (Antonym Generation, Sentence Completion (SC), and Auditory Naming) for lateralizing language function and for activating cortex within patient-specific regions-of-interest representing eloquent language areas, and assessed the degree of spatial overlap of the areas of activation elicited by each task. RESULTS: The tasks were similarly effective for lateralizing language within the anterior language areas. The SC task produced higher laterality indices within the posterior language areas and had a significantly higher agreement with the clinical report. Dice coefficients between the task pairs were in the range of .351-.458, confirming substantial variation in the components of the language network activated by each task. CONCLUSIONS: SC task consistently produced large activations within the dominant hemisphere and was more effective for lateralizing language within the posterior language areas. The low degree of spatial overlap among the tasks strongly supports the practice of using a battery of tasks to help the surgeon to avoid eloquent language areas.
Sjölund J, Eklund A, Özarslan E, Herberthson M, Bånkestad M, Knutsson H. Bayesian Uncertainty Quantification in Linear Models for Diffusion MRI. Neuroimage. 2018;175 :272-85.Abstract
Diffusion MRI (dMRI) is a valuable tool in the assessment of tissue microstructure. By fitting a model to the dMRI signal it is possible to derive various quantitative features. Several of the most popular dMRI signal models are expansions in an appropriately chosen basis, where the coefficients are determined using some variation of least-squares. However, such approaches lack any notion of uncertainty, which could be valuable in e.g. group analyses. In this work, we use a probabilistic interpretation of linear least-squares methods to recast popular dMRI models as Bayesian ones. This makes it possible to quantify the uncertainty of any derived quantity. In particular, for quantities that are affine functions of the coefficients, the posterior distribution can be expressed in closed-form. We simulated measurements from single- and double-tensor models where the correct values of several quantities are known, to validate that the theoretically derived quantiles agree with those observed empirically. We included results from residual bootstrap for comparison and found good agreement. The validation employed several different models: Diffusion Tensor Imaging (DTI), Mean Apparent Propagator MRI (MAP-MRI) and Constrained Spherical Deconvolution (CSD). We also used in vivo data to visualize maps of quantitative features and corresponding uncertainties, and to show how our approach can be used in a group analysis to downweight subjects with high uncertainty. In summary, we convert successful linear models for dMRI signal estimation to probabilistic models, capable of accurate uncertainty quantification.