This work presents an anatomically curated white matter atlas to enable consistent white matter tract parcellation across different populations. Leveraging a well-established computational pipeline for fiber clustering, we create a tract-based white matter atlas including information from 100 subjects. A novel anatomical annotation method is proposed that leverages population-based brain anatomical information and expert neuroanatomical knowledge to annotate and categorize the fiber clusters. A total of 256 white matter structures are annotated in the proposed atlas, which provides one of the most comprehensive tract-based white matter atlases covering the entire brain to date. These structures are composed of 58 deep white matter tracts including major long range association and projection tracts, commissural tracts, and tracts related to the brainstem and cerebellar connections, plus 198 short and medium range superficial fiber clusters organized into 16 categories according to the brain lobes they connect. Potential false positive connections are annotated in the atlas to enable their exclusion from analysis or visualization. In addition, the proposed atlas allows for a whole brain white matter parcellation into 800 fiber clusters to enable whole brain connectivity analyses. The atlas and related computational tools are open-source and publicly available. We evaluate the proposed atlas using a testing dataset of 584 diffusion MRI scans from multiple independently acquired populations, across genders, the lifespan (1 day-82 years), and different health conditions (healthy control, neuropsychiatric disorders, and brain tumor patients). Experimental results show successful white matter parcellation across subjects from different populations acquired on multiple scanners, irrespective of age, gender or disease indications. Over 99% of the fiber tracts annotated in the atlas were detected in all subjects on average. One advantage in terms of robustness is that the tract-based pipeline does not require any cortical or subcortical segmentations, which can have limited success in young children and patients with brain tumors or other structural lesions. We believe this is the first demonstration of consistent automated white matter tract parcellation across the full lifespan from birth to advanced age.

We present an algorithm for creating high resolution anatomically plausible images consistent with acquired clinical brain MRI scans with large inter-slice spacing. Although large data sets of clinical images contain a wealth of information, time constraints during acquisition result in sparse scans that fail to capture much of the anatomy. These characteristics often render computational analysis impractical as many image analysis algorithms tend to fail when applied to such images. Highly specialized algorithms that explicitly handle sparse slice spacing do not generalize well across problem domains. In contrast, we aim to enable application of existing algorithms that were originally developed for high resolution research scans to significantly undersampled scans. We introduce a generative model that captures fine-scale anatomical structure across subjects in clinical image collections and derive an algorithm for filling in the missing data in scans with large inter-slice spacing. Our experimental results demonstrate that the resulting method outperforms state-of-the-art upsampling super-resolution techniques, and promises to facilitate subsequent analysis not previously possible with scans of this quality. Our implementation is freely available at https://github.com/adalca/papago.

This work presents an automatically annotated fiber cluster (AAFC) method to enable identification of anatomically meaningful white matter structures from the whole brain tractography. The proposed method consists of 1) a study-specific whole brain white matter parcellation using a well-established data-driven groupwise fiber clustering pipeline to segment tractography into multiple fiber clusters, and 2) a novel cluster annotation method to automatically assign an anatomical tract annotation to each fiber cluster by employing cortical parcellation information across multiple subjects. The novelty of the AAFC method is that it leverages group-wise information about the fiber clusters, including their fiber geometry and cortical terminations, to compute a tract anatomical label for each cluster in an automated fashion. We demonstrate the proposed AAFC method in an application of investigating white matter abnormality in emotional processing and sensorimotor areas in major depressive disorder (MDD). Seven tracts of interest related to emotional processing and sensorimotor functions are automatically identified using the proposed AAFC method as well as a comparable method that uses a cortical parcellation alone. Experimental results indicate that our proposed method is more consistent in identifying the tracts across subjects and across hemispheres in terms of the number of fibers. In addition, we perform a between-group statistical analysis in 31 MDD patients and 62 healthy subjects on the identified tracts using our AAFC method. We find statistical differences in diffusion measures in local regions within a fiber tract (e.g. 4 fiber clusters within the identified left hemisphere cingulum bundle (consisting of 14 clusters) are significantly different between the two groups), suggesting the ability of our method in identifying potential abnormality specific to subdivisions of a white matter structure.

Diffusion magnetic resonance imaging (dMRI) is an important method for studying white matter connectivity in the brain in vivo in both healthy and clinical populations. Improvements in dMRI tractography algorithms, which reconstruct macroscopic three-dimensional white matter fiber pathways, have allowed for methodological advances in the study of white matter; however, insufficient attention has been paid to comparing post-tractography methods that extract white matter fiber tracts of interest from whole-brain tractography. Here we conduct a comparison of three representative and conceptually distinct approaches to fiber tract delineation: 1) a manual multiple region of interest-based approach, 2) an atlas-based approach, and 3) a groupwise fiber clustering approach, by employing methods that exemplify these approaches to delineate the arcuate fasciculus, the middle longitudinal fasciculus, and the uncinate fasciculus in 10 healthy male subjects. We enable qualitative comparisons across methods, conduct quantitative evaluations of tract volume, tract length, mean fractional anisotropy, and true positive and true negative rates, and report measures of intra-method and inter-method agreement. We discuss methodological similarities and differences between the three approaches and the major advantages and drawbacks of each, and review research and clinical contexts for which each method may be most apposite. Emphasis is given to the means by which different white matter fiber tract delineation approaches may systematically produce variable results, despite utilizing the same input tractography and reliance on similar anatomical knowledge.

We propose a new iterative segmentation model which can be accurately learned from a small dataset. A common approach is to train a model to directly segment an image, requiring a large collection of manually annotated images to capture the anatomical variability in a cohort. In contrast, we develop a segmentation model that recursively evolves a segmentation in several steps, and implement it as a recurrent neural network. We learn model parameters by optimizing the intermediate steps of the evolution in addition to the final segmentation. To this end, we train our segmentation propagation model by presenting incomplete and/or inaccurate input segmentations paired with a recommended next step. Our work aims to alleviate challenges in segmenting heart structures from cardiac MRI for patients with congenital heart disease (CHD), which encompasses a range of morphological deformations and topological changes. We demonstrate the advantages of this approach on a dataset of 20 images from CHD patients, learning a model that accurately segments individual heart chambers and great vessels. Compared to direct segmentation, the iterative method yields more accurate segmentation for patients with the most severe CHD malformations.

This paper presents a novel approach to modeling the pos terior distribution in image registration that is computationally efficient for large deformation diffeomorphic metric mapping (LDDMM). We develop a Laplace approximation of Bayesian registration models entirely in a bandlimited space that fully describes the properties of diffeomorphic transformations. In contrast to current methods, we compute the inverse Hessian at the mode of the posterior distribution of diffeomorphisms directly in the low dimensional frequency domain. This dramatically reduces the computational complexity of approximating posterior marginals in the high dimensional imaging space. Experimental results show that our method is significantly faster than the state-of-the-art diffeomorphic image registration uncertainty quantification algorithms, while producing comparable results. The efficiency of our method strengthens the feasibility in prospective clinical applications, e.g., real- time image-guided navigation for brain surgery.

A reliable Ultrasound (US)-to-US registration method to compensate for brain shift would substantially improve Image-Guided Neurological Surgery. Developing such a registration method is very challenging, due to factors such as the tumor resection, the complexity of brain pathology and the demand for fast computation. We propose a novel feature-driven active registration framework. Here, landmarks and their displacement are first estimated from a pair of US images using corresponding local image features. Subsequently, a Gaussian Process (GP) model is used to interpolate a dense deformation field from the sparse landmarks. Kernels of the GP are estimated by using variograms and a discrete grid search method. If necessary, the user can actively add new landmarks based on the image context and visualization of the uncertainty measure provided by the GP to further improve the result. We retrospectively demonstrate our registration framework as a robust and accurate brain shift compensation solution on clinical data.

Free Water Imaging is a novel diffusion magnetic resonance (MR) imaging method that is able to separate changes affecting the extracellular space from those that reflect changes in neuronal cells and processes. A previous Free Water Imaging study in schizophrenia identified significantly greater extracellular water volume in the early stages of the disorder; however, its clinical and functional sequelae have not yet been investigated. Here, we applied Free Water Imaging to a larger cohort of 63 first-episode patients with psychosis and 70 healthy matched controls to better understand the functional significance of greater extracellular water. We used diffusion MR imaging data and the Tract-Based Spatial Statistics analytic pipeline to first analyze fractional anisotropy (FA), the most commonly employed metric for assessing white matter. This comparison was then followed by Free Water Imaging analysis, where two parameters, the fractional volume of extracellular free-water (FW) and cellular tissue FA (FA-t), were estimated and compared across the entire white matter skeleton between groups, and correlated with cognitive measures at baseline and following 12 weeks of antipsychotic treatment. Our results indicated lower FA across the whole brain in patients compared with healthy controls that overlap with significant increases in FW, with only limited decreases in FA-t. In addition, higher FW correlated with better neurocognitive functioning following 12 weeks of antipsychotic treatment. We believe this is the first study to suggest that an extracellular water increase during the first-episode of psychosis, which may be indicative of an acute neuroinflammatory process, and/or cerebral edema may predict better functional outcome.

Mild traumatic brain injuries (mTBIs) are often associated with posttraumatic stress disorder (PTSD). In cases of chronic mTBI, accurate diagnosis can be challenging due to the overlapping symptoms this condition shares with PTSD. Furthermore, mTBIs are heterogeneous and not easily observed using conventional neuroimaging tools, despite the fact that diffuse axonal injuries are the most common injury. Diffusion tensor imaging (DTI) is sensitive to diffuse axonal injuries and is thus more likely to detect mTBIs, especially when analyses account for the inter-individual variability of these injuries. Using a subject-specific approach, we compared fractional anisotropy (FA) abnormalities between groups with a history of mTBI (n = 35), comorbid mTBI and PTSD (mTBI + PTSD; n = 22), and healthy controls (n = 37). We compared all three groups on the number of abnormal FA clusters derived from subject-specific injury profiles (i.e., individual z-score maps) along a common white matter skeleton. The mTBI + PTSD group evinced a greater number of abnormally low FA clusters relative to both the healthy controls and the mTBI group without PTSD (p < .05). Across the groups with a history of mTBI, increased numbers of abnormally low FA clusters were significantly associated with PTSD symptom severity, depression, post-concussion symptoms, and reduced information processing speed (p < .05). These findings highlight the utility of subject-specific microstructural analyses when searching for mTBI-related brain abnormalities, particularly in patients with PTSD. This study also suggests that patients with a history of mTBI and comorbid PTSD, relative to those without PTSD, are at increased risk of FA abnormalities.

Mastering the technical skills required to perform pediatric cardiac valve surgery is challenging in part due to limited opportunity for practice. Transformation of 3D echocardiographic (echo) images of congenitally abnormal heart valves to realistic physical models could allow patient-specific simulation of surgical valve repair. We compared materials, processes, and costs for 3D printing and molding of patient-specific models for visualization and surgical simulation of congenitally abnormal heart valves. Pediatric atrioventricular valves (mitral, tricuspid, and common atrioventricular valve) were modeled from transthoracic 3D echo images using semi-automated methods implemented as custom modules in 3D Slicer. Valve models were then both 3D printed in soft materials and molded in silicone using 3D printed "negative" molds. Using pre-defined assessment criteria, valve models were evaluated by congenital cardiac surgeons to determine suitability for simulation. Surgeon assessment indicated that the molded valves had superior material properties for the purposes of simulation compared to directly printed valves (p < 0.01). Patient-specific, 3D echo-derived molded valves are a step toward realistic simulation of complex valve repairs but require more time and labor to create than directly printed models. Patient-specific simulation of valve repair in children using such models may be useful for surgical training and simulation of complex congenital cases.