Publications

2018

Koppelmans V, Scott JM, Downs ME, Cassady KE, Yuan P, Pasternak O, Wood SJ, De Dios YE, Gadd NE, Kofman I, Riascos R, Reuter-Lorenz PA, Bloomberg JJ, Mulavara AP, Ploutz-Snyder LL, Seidler RD. Exercise Effects on Bed Rest-induced Brain Changes. PLoS One. 2018;13(10):e0205515.
PURPOSE: Spaceflight negatively affects sensorimotor behavior; exercise mitigates some of these effects. Head down tilt bed rest (HDBR) induces body unloading and fluid shifts, and is often used to investigate spaceflight effects. Here, we examined whether exercise mitigates effects of 70 days HDBR on the brain and if fitness and brain changes with HDBR are related. METHODS: HDBR subjects were randomized to no-exercise (n = 5) or traditional aerobic and resistance exercise (n = 5). Additionally, a flywheel exercise group was included (n = 8). Exercise protocols for exercise groups were similar in intensity, therefore these groups were pooled in statistical analyses. Pre and post-HDBR MRI (structure and structural/functional connectivity) and physical fitness measures (lower body strength, muscle cross sectional area, VO2 max, body composition) were collected. Voxel-wise permutation analyses were used to test group differences in brain changes, and their associations with fitness changes. RESULTS: Comparisons of exercisers to controls revealed that exercise led to smaller fitness deterioration with HDBR but did not affect brain volume or connectivity. Group comparisons showed that exercise modulated post-HDBR recovery of brain connectivity in somatosensory regions. Posthoc analysis showed that this was related to functional connectivity decrease with HDBR in non-exercisers but not in exercisers. Correlational analyses between fitness and brain changes showed that fitness decreases were associated with functional connectivity and volumetric increases (all r >.74), potentially reflecting compensation. Modest brain changes or even decreases in connectivity and volume were observed in subjects who maintained or showed small fitness gains. These results did not survive Bonferroni correction, but can be considered meaningful because of the large effect sizes. CONCLUSION: Exercise performed during HDBR mitigates declines in fitness and strength. Associations between fitness and brain connectivity and volume changes, although unadjusted for multiple comparisons in this small sample, suggest that supine exercise reduces compensatory HDBR-induced brain changes.
Coherence and causality measures are often used to analyze the influence of one region on another during analysis of functional brain networks. The analysis methods usually involve a regression problem, where the signal of interest is decomposed into a mixture of regressor and a residual signal. In this paper, we revisit this basic problem and present solutions that provide the minimal-entropy residuals for different types of regression filters, such as causal, instantaneously causal, and noncausal filters. Using optimal prediction theory, we derive several novel frequency-domain expressions for partial coherence, causality, and conditional causality analysis. In particular, our solution provides a more accurate estimation of the frequency-domain causality compared with the classical Geweke causality measure. Using synthetic examples and in vivo resting-state functional magnetic resonance imaging data from the human connectome project, we show that the proposed solution is more accurate at revealing frequency-domain linear dependence among high-dimensional signals.
Herrmann MD, Clunie DA, Fedorov A, Doyle SW, Pieper S, Klepeis V, Le LP, Mutter GL, Milstone DS, Schultz TJ, Kikinis R, Kotecha GK, Hwang DH, Andriole KP, Iafrate J, Brink JA, Boland GW, Dreyer KJ, Michalski M, Golden JA, Louis DN, Lennerz JK. Implementing the DICOM Standard for Digital Pathology. J Pathol Inform. 2018;9:37.
Background: Digital Imaging and Communications in Medicine (DICOM) is the standard for the representation, storage, and communication of medical images and related information. A DICOM file format and communication protocol for pathology have been defined; however, adoption by vendors and in the field is pending. Here, we implemented the essential aspects of the standard and assessed its capabilities and limitations in a multisite, multivendor healthcare network. Methods: We selected relevant DICOM attributes, developed a program that extracts pixel data and pixel-related metadata, integrated patient and specimen-related metadata, populated and encoded DICOM attributes, and stored DICOM files. We generated the files using image data from four vendor-specific image file formats and clinical metadata from two departments with different laboratory information systems. We validated the generated DICOM files using recognized DICOM validation tools and measured encoding, storage, and access efficiency for three image compression methods. Finally, we evaluated storing, querying, and retrieving data over the web using existing DICOM archive software. Results: Whole slide image data can be encoded together with relevant patient and specimen-related metadata as DICOM objects. These objects can be accessed efficiently from files or through RESTful web services using existing software implementations. Performance measurements show that the choice of image compression method has a major impact on data access efficiency. For lossy compression, JPEG achieves the fastest compression/decompression rates. For lossless compression, JPEG-LS significantly outperforms JPEG 2000 with respect to data encoding and decoding speed. Conclusion: Implementation of DICOM allows efficient access to image data as well as associated metadata. By leveraging a wealth of existing infrastructure solutions, the use of DICOM facilitates enterprise integration and data exchange for digital pathology.
Schultz V, Stern RA, Tripodis Y, Stamm J, Wrobel P, Lepage C, Weir I, Guenette JP, Chua A, Alosco ML, Baugh CM, Fritts NG, Martin BM, Chaisson CE, Coleman MJ, Lin AP, Pasternak O, Shenton ME, Koerte IK. Age at First Exposure to Repetitive Head Impacts Is Associated with Smaller Thalamic Volumes in Former Professional American Football Players. J Neurotrauma. 2018;35(2):278–85.
Thalamic atrophy has been associated with exposure to repetitive head impacts (RHI) in professional fighters. The aim of this study is to investigate whether or not age at first exposure (AFE) to RHI is associated with thalamic volume in symptomatic former National Football League (NFL) players at risk for chronic traumatic encephalopathy (CTE). Eighty-six symptomatic former NFL players (mean age = 54.9 ± 7.9 years) were included. T1-weighted data were acquired on a 3T magnetic resonance imager, and thalamic volumes were derived using FreeSurfer. Mood and behavior, psychomotor speed, and visual and verbal memory were assessed. The association between thalamic volume and AFE to playing football and to number of years playing was calculated. Decreased thalamic volume was associated with more years of play (left: p = 0.03; right: p = 0.03). Younger AFE was associated with decreased right thalamic volume (p = 0.014). This association remained significant after adjusting for total years of play. Decreased left thalamic volume was associated with worse visual memory (p = 0.014), whereas increased right thalamic volume was associated with fewer mood and behavior symptoms (p = 0.003). In our sample of symptomatic former NFL players at risk for CTE, total years of play and AFE were associated with decreased thalamic volume. The effect of AFE on right thalamic volume was almost twice as strong as the effect of total years of play. Our findings confirm previous reports of an association between thalamic volume and exposure to RHI. They suggest further that younger AFE may result in smaller thalamic volume later in life.
Pasternak O, Kelly S, Sydnor VJ, Shenton ME. Advances in Microstructural Diffusion Neuroimaging for Psychiatric Disorders. Neuroimage. 2018;182:259–82.
Understanding the neuropathological underpinnings of mental disorders such as schizophrenia, major depression, and bipolar disorder is an essential step towards the development of targeted treatments. Diffusion MRI studies utilizing the diffusion tensor imaging (DTI) model have been extremely successful to date in identifying microstructural brain abnormalities in individuals suffering from mental illness, especially in regions of white matter, although identified abnormalities have been biologically non-specific. Building on DTI’s success, in recent years more advanced diffusion MRI methods have been developed and applied to the study of psychiatric populations, with the aim of offering increased sensitivity to subtle neurological abnormalities, as well as improved specificity to candidate pathologies such as demyelination and neuroinflammation. These advanced methods, however, usually come at the cost of prolonged imaging sequences or reduced signal to noise, and they are more difficult to evaluate compared with the more simplified approach taken by the now common DTI model. To date, a limited number of advanced diffusion MRI methods have been employed to study schizophrenia, major depression and bipolar disorder populations. In this review we survey these studies, compare findings across diverse methods, discuss the main benefits and limitations of the different methods, and assess the extent to which the application of more advanced diffusion imaging approaches has led to novel and transformative information with regards to our ability to better understand the etiology and pathology of mental disorders.
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.
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.

2017

Sastry R, Bi WL, Pieper S, Frisken S, Kapur T, Wells WM III, Golby AJ. Applications of Ultrasound in the Resection of Brain Tumors. J Neuroimaging. 2017;27(1):5–15.

Neurosurgery makes use of preoperative imaging to visualize pathology, inform surgical planning, and evaluate the safety of selected approaches. The utility of preoperative imaging for neuronavigation, however, is diminished by the well-characterized phenomenon of brain shift, in which the brain deforms intraoperatively as a result of craniotomy, swelling, gravity, tumor resection, cerebrospinal fluid (CSF) drainage, and many other factors. As such, there is a need for updated intraoperative information that accurately reflects intraoperative conditions. Since 1982, intraoperative ultrasound has allowed neurosurgeons to craft and update operative plans without ionizing radiation exposure or major workflow interruption. Continued evolution of ultrasound technology since its introduction has resulted in superior imaging quality, smaller probes, and more seamless integration with neuronavigation systems. Furthermore, the introduction of related imaging modalities, such as 3-dimensional ultrasound, contrast-enhanced ultrasound, high-frequency ultrasound, and ultrasound elastography, has dramatically expanded the options available to the neurosurgeon intraoperatively. In the context of these advances, we review the current state, potential, and challenges of intraoperative ultrasound for brain tumor resection. We begin by evaluating these ultrasound technologies and their relative advantages and disadvantages. We then review three specific applications of these ultrasound technologies to brain tumor resection: (1) intraoperative navigation, (2) assessment of extent of resection, and (3) brain shift monitoring and compensation. We conclude by identifying opportunities for future directions in the development of ultrasound technologies.

Ning L, Setsompop K, Westin CF, Rathi Y. New Insights about Time-varying Diffusivity and its Estimation from Diffusion MRI. Magn Reson Med. 2017;78(2):763–74.

PURPOSE: Characterizing the relation between the applied gradient sequences and the measured diffusion MRI signal is important for estimating the time-dependent diffusivity, which provides important information about the microscopic tissue structure. THEORY AND METHODS: In this article, we extend the classical theory of Stepi\v snik for measuring time-dependent diffusivity under the Gaussian phase approximation. In particular, we derive three novel expressions which represent the diffusion MRI signal in terms of the mean-squared displacement, the instantaneous diffusivity, and the velocity autocorrelation function. We present the explicit signal expressions for the case of single diffusion encoding and oscillating gradient spin-echo sequences. Additionally, we also propose three different models to represent time-varying diffusivity and test them using Monte-Carlo simulations and in vivo human brain data. RESULTS: The time-varying diffusivities are able to distinguish the synthetic structures in the Monte-Carlo simulations. There is also strong statistical evidence about time-varying diffusivity from the in vivo human data set. CONCLUSION: The proposed theory provides new insights into our understanding of the time-varying diffusivity using different gradient sequences. The proposed models for representing time-varying diffusivity can be utilized to study time-varying diffusivity using in vivo human brain diffusion MRI data. 

Ning L, Özarslan E, Westin CF, Rathi Y. Precise Inference and Characterization of Structural Organization (PICASO) of Tissue from Molecular Diffusion. Neuroimage. 2017;146:452–73.

Inferring the microstructure of complex media from the diffusive motion of molecules is a challenging problem in diffusion physics. In this paper, we introduce a novel representation of diffusion MRI (dMRI) signal from tissue with spatially-varying diffusivity using a diffusion disturbance function. This disturbance function contains information about the (intra-voxel) spatial fluctuations in diffusivity due to restrictions, hindrances and tissue heterogeneity of the underlying tissue substrate. We derive the short- and long-range disturbance coefficients from this disturbance function to characterize the tissue structure and organization. Moreover, we provide an exact relation between the disturbance coefficients and the time-varying moments of the diffusion propagator, as well as their relation to specific tissue microstructural information such as the intra-axonal volume fraction and the apparent axon radius. The proposed approach is quite general and can model dMRI signal for any type of gradient sequence (rectangular, oscillating, etc.) without using the Gaussian phase approximation. The relevance of the proposed PICASO model is explored using Monte-Carlo simulations and in-vivo dMRI data. The results show that the estimated disturbance coefficients can distinguish different types of microstructural organization of axons.