Publications

Self-diffusion in polymer solutions studied with pulsed-field gradient nuclear magnetic resonance (PFG NMR) is typically based either on a single self-diffusion coefficient, or a log-normal distribution of self-diffusion coefficients, or in some cases mixtures of these. Experimental data on polyethylene glycol (PEG) solutions and simulations were used to compare a model based on a gamma distribution of self-diffusion coefficients to more established models such as the single exponential, the stretched exponential, and the log-normal distribution model with regard to performance and consistency. Even though the gamma distribution is very similar to the log-normal distribution, its NMR signal attenuation can be written in a closed form and therefore opens up for increased computational speed. Estimates of the mean self-diffusion coefficient, the spread, and the polydispersity index that were obtained using the gamma model were in excellent agreement with estimates obtained using the log-normal model. Furthermore, we demonstrate that the gamma distribution is by far superior to the log-normal, and comparable to the two other models, in terms of computational speed. This effect is particularly striking for multi-component signal attenuation. Additionally, the gamma distribution as well as the log-normal distribution incorporates explicitly a physically plausible model for polydispersity and spread, in contrast to the single exponential and the stretched exponential. Therefore, the gamma distribution model should be preferred in many experimental situations.

In this study, we present pituitary adenoma volumetry using the free and open source medical image computing platform for biomedical research: (3D) Slicer. Volumetric changes in cerebral pathologies like pituitary adenomas are a critical factor in treatment decisions by physicians and in general the volume is acquired manually. Therefore, manual slice-by-slice segmentations in magnetic resonance imaging (MRI) data, which have been obtained at regular intervals, are performed. In contrast to this manual time consuming slice-by-slice segmentation process Slicer is an alternative which can be significantly faster and less user intensive. In this contribution, we compare pure manual segmentations of ten pituitary adenomas with semi-automatic segmentations under Slicer. Thus, physicians drew the boundaries completely manually on a slice-by-slice basis and performed a Slicer-enhanced segmentation using the competitive region-growing based module of Slicer named GrowCut. Results showed that the time and user effort required for GrowCut-based segmentations were on average about thirty percent less than the pure manual segmentations. Furthermore, we calculated the Dice Similarity Coefficient (DSC) between the manual and the Slicer-based segmentations to proof that the two are comparable yielding an average DSC of 81.97±3.39%.

BACKGROUND: Natural orifice transluminal endoscopic surgery (NOTES) is technically challenging owing to endoscopic short-sighted visualization, excessive scope flexibility and lack of adequate instrumentation. Augmented reality may overcome these difficulties. This study tested whether an image registration system for NOTES procedures (IR-NOTES) can facilitate navigation. METHODS: In three human cadavers 15 intra-abdominal organs were targeted endoscopically with and without IR-NOTES via both transgastric and transcolonic routes, by three endoscopists with different levels of expertise. Ease of navigation was evaluated objectively by kinematic analysis, and navigation complexity was determined by creating an organ access complexity score based on the same data.

Available approaches to the investigation of traumatic brain injury (TBI) are frequently hampered, to some extent, by the unsatisfactory abilities of existing methodologies to efficiently define and represent affected structural connectivity and functional mechanisms underlying TBI-related pathology. In this paper, we describe a patient-tailored framework which allows mapping and characterization of TBI-related structural damage to the brain via multimodal neuroimaging and personalized connectomics. Specifically, we introduce a graphically driven approach for the assessment of trauma-related atrophy of white matter connections between cortical structures, with relevance to the quantification of TBI chronic case evolution. This approach allows one to inform the formulation of graphical neurophysiological and neuropsychological TBI profiles based on the particular structural deficits of the affected patient. In addition, it allows one to relate the findings supplied by our workflow to the existing body of research that focuses on the functional roles of the cortical structures being targeted. A graphical means for representing patient TBI status is relevant to the emerging field of personalized medicine and to the investigation of neural atrophy.

Type 2 diabetes mellitus (T2DM) is a risk factor for Alzheimer disease (AD). Populations at risk for AD show altered brain activity in the default mode network (DMN) before cognitive dysfunction. We evaluated this brain pattern in T2DM patients. We compared T2DM patients (n = 10, age = 56 ± 2.2 years, fasting plasma glucose [FPG] = 8.4 ± 1.3 mmol/L, HbA(1c) = 7.5 ± 0.54%) with nondiabetic age-matched control subjects (n = 11, age = 54 ± 1.8 years, FPG = 4.8 ± 0.2 mmol/L) using resting-state functional magnetic resonance imaging to evaluate functional connectivity strength among DMN regions. We also evaluated hippocampal volume, cognition, and insulin sensitivity by homeostasis model assessment of insulin resistance (HOMA-IR). Control subjects showed stronger correlations versus T2DM patients in the DMN between the seed (posterior cingulate) and bilateral middle temporal gyrus (β = 0.67 vs. 0.43), the right inferior and left medial frontal gyri (β = 0.75 vs. 0.54), and the left thalamus (β = 0.59 vs. 0.37), respectively, with no group differences in cognition or hippocampal size. In T2DM patients, HOMA-IR was inversely correlated with functional connectivity in the right inferior frontal gyrus and precuneus. T2DM patients showed reduced functional connectivity in the DMN compared with control subjects, which was associated with insulin resistance in selected brain regions, but there were no group effects of brain structure or cognition.

INTRODUCTION: Neuroanatomical and cognitive alterations typical of schizophrenia (SZ) patients are observed to a lesser extent in their adolescent and adult first-degree relatives, likely reflecting neurodevelopmental abnormalities associated with genetic risk for the illness. The anatomical pathways for language are hypothesized to be abnormal and to underlie the positive symptoms of schizophrenia. Examining non-psychotic relatives at high familial risk (FHR) for schizophrenia may clarify if these deficits represent trait markers associated with genetic vulnerability, rather than specific markers resulting from the pathological process underlying schizophrenia.

White matter (WM) mapping of the human brain using neuroimaging techniques has gained considerable interest in the neuroscience community. Using diffusion weighted (DWI) and magnetic resonance imaging (MRI), WM fiber pathways between brain regions may be systematically assessed to make inferences concerning their role in normal brain function, influence on behavior, as well as concerning the consequences of network-level brain damage. In this paper, we investigate the detailed connectomics in a noted example of severe traumatic brain injury (TBI) which has proved important to and controversial in the history of neuroscience. We model the WM damage in the notable case of Phineas P. Gage, in whom a "tamping iron" was accidentally shot through his skull and brain, resulting in profound behavioral changes. The specific effects of this injury on Mr. Gage’s WM connectivity have not previously been considered in detail. Using computed tomography (CT) image data of the Gage skull in conjunction with modern anatomical MRI and diffusion imaging data obtained in contemporary right handed male subjects (aged 25-36), we computationally simulate the passage of the iron through the skull on the basis of reported and observed skull fiducial landmarks and assess the extent of cortical gray matter (GM) and WM damage. Specifically, we find that while considerable damage was, indeed, localized to the left frontal cortex, the impact on measures of network connectedness between directly affected and other brain areas was profound, widespread, and a probable contributor to both the reported acute as well as long-term behavioral changes. Yet, while significantly affecting several likely network hubs, damage to Mr. Gage’s WM network may not have been more severe than expected from that of a similarly sized "average" brain lesion. These results provide new insight into the remarkable brain injury experienced by this noteworthy patient.

A method for automated location of shape differences in diseased anatomical structures via high resolution biomedical atlases annotated with labels from formal ontologies is described. In particular, a high resolution magnetic resonance image of the myocardium of the human left ventricle was segmented and annotated with structural terms from an extracted subset of the Foundational Model of Anatomy ontology. The atlas was registered to the end systole template of a previous study of left ventricular remodeling in cardiomyopathy using a diffeomorphic registration algorithm. The previous study used thresholding and visual inspection to locate a region of statistical significance which distinguished patients with ischemic cardiomyopathy from those with nonischemic cardiomyopathy. Using semantic technologies and the deformed annotated atlas, this location was more precisely found. Although this study used only a cardiac atlas, it provides a proof-of-concept that ontologically labeled biomedical atlases of any anatomical structure can be used to automate location-based inferences.

The validity of parametric functional magnetic resonance imaging (fMRI) analysis has only been reported for simulated data. Recent advances in computer science and data sharing make it possible to analyze large amounts of real fMRI data. In this study, 1484 rest datasets have been analyzed in SPM8, to estimate true familywise error rates. For a familywise significance threshold of 5%, significant activity was found in 1%-70% of the 1484 rest datasets, depending on repetition time, paradigm and parameter settings. This means that parametric significance thresholds in SPM both can be conservative or very liberal. The main reason for the high familywise error rates seems to be that the global AR(1) auto correlation correction in SPM fails to model the spectra of the residuals, especially for short repetition times. The findings that are reported in this study cannot be generalized to parametric fMRI analysis in general, other software packages may give different results. By using the computational power of the graphics processing unit (GPU), the 1484 rest datasets were also analyzed with a random permutation test. Significant activity was then found in 1%-19% of the datasets. These findings speak to the need for a better model of temporal correlations in fMRI timeseries.

Heart wall myofibers wind as helices around the ventricles, strengthening them in a manner analogous to the reinforcement of concrete cylindrical columns by spiral steel cables [Richart FE, et al. (1929) Univ of Illinois, Eng Exp Stn Bull 190]. A multitude of such fibers, arranged smoothly and regularly, contract and relax as an integrated functional unit as the heart beats. To orchestrate this motion, fiber tangling must be avoided and pumping should be efficient. Current models of myofiber orientation across the heart wall suggest groupings into sheets or bands, but the precise geometry of bundles of myofibers is unknown. Here we show that this arrangement takes the form of a special minimal surface, the generalized helicoid [Blair DE, Vanstone JR (1978) Minimal Submanifolds and Geodesics 13-16], closing the gap between individual myofibers and their collective wall structure. The model holds across species, with a smooth variation in its three curvature parameters within the myocardial wall providing tight fits to diffusion magnetic resonance images from the rat, the dog, and the human. Mathematically it explains how myofibers are bundled in the heart wall while economizing fiber length and optimizing ventricular ejection volume as they contract. The generalized helicoid provides a unique foundation for analyzing the fibrous composite of the heart wall and should therefore find applications in heart tissue engineering and in the study of heart muscle diseases.