The joint effect of aging and HIV infection on microstructure of white matter bundles.

Recent evidence suggests the aging process is accelerated by HIV. Degradation of white matter (WM) has been independently associated with HIV and healthy aging. Thus, WM may be vulnerable to joint effects of HIV and aging. Diffusion-weighted imaging (DWI) was conducted with HIV-seropositive (n = 72) and HIV-seronegative (n = 34) adults. DWI data underwent tractography, which was parcellated into 18 WM tracts of interest (TOIs). Functional Analysis of Diffusion Tensor Tract Statistics (FADTTS) regression was conducted assessing the joint effect of advanced age and HIV on fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD), and radial diffusivity (RD) along TOI fibers. In addition to main effects of age and HIV on WM microstructure, the interactive effect of age and HIV was significantly related to lower FA and higher MD, AD, and RD across all TOIs. The location of findings was consistent with the clinical presentation of HIV-associated neurocognitive disorders. While older age is related to poorer WM microstructure, its detrimental effect on WM is stronger among HIV+ relative to HIV- individuals. Loss of WM integrity in the context of advancing age may place HIV+ individuals at increased risk for brain and cognitive compromise.

Superficial white matter damage in anti-NMDA receptor encephalitis.

BACKGROUND: Clinical brain MRI is normal in the majority of patients with anti-N-methyl-D-aspartate receptor (NMDAR) encephalitis. However, extensive deep white matter damage wasrecently identifiedin these patients using diffusion weighted imaging. Here, our aim was to study a particularly vulnerable brain compartment, the late myelinating superficial white matter. METHODS: Forty-six patients with anti-NMDAR encephalitis were included. Ten out of these were considered neurologically recovered (modified Rankin scale of zero), while 36 patients were non-recovered. In addition, 30 healthy controls were studied. MRI data were collected from all subjects and superficial white matter mean diffusivity derived from diffusion tensor imaging was compared between groups in whole brain, lobar and vertex-based analyses. Patients underwent comprehensive cognitive testing, and correlation analyses were performed between cognitive performance and superficial white matter integrity. RESULTS: Non-recovered patients showed widespread superficial white matter damage in comparison to recovered patients and healthy controls. Vertex-based analyses revealed that damage predominated in frontal and temporal lobes. In contrast, the superficial white matter was intact in recovered patients. Importantly, persistent cognitive impairments in working memory, verbal memory, visuospatial memory and attention significantly correlated with damage of the superficial white matter in patients. CONCLUSIONS: Anti-NMDAR encephalitis is associated with extensive superficial white matter damage in patients with incomplete recovery. The strong association with impairment in several cognitive domains highlights the clinical relevance of white matter damage in this disorder and warrants investigations of the underlying pathophysiological mechanisms.

The superficial white matter in Alzheimer's disease.

White matter abnormalities have been shown in the large deep fibers of Alzheimer's disease patients. However, the late myelinating superficial white matter comprised of intracortical myelin and short-range association fibers has not received much attention. To investigate this area, we extracted a surface corresponding to the superficial white matter beneath the cortex and then applied a cortical pattern-matching approach which allowed us to register and subsequently sample diffusivity along thousands of points at the interface between the gray matter and white matter in 44 patients with Alzheimer's disease (Age: 71.02 ± 5.84, 16M/28F) and 47 healthy controls (Age 69.23 ± 4.45, 19M/28F). In patients we found an overall increase in the axial and radial diffusivity across most of the superficial white matter (P < 0.001) with increases in diffusivity of more than 20% in the bilateral parahippocampal regions and the temporal and frontal lobes. Furthermore, diffusivity correlated with the cognitive deficits measured by the Mini-Mental State Examination scores (P < 0.001). The superficial white matter has a unique microstructure and is critical for the integration of multimodal information during brain maturation and aging. Here we show that there are major abnormalities in patients and the deterioration of these fibers relates to clinical symptoms in Alzheimer's disease.

Interactive effects of BDNF Val66Met genotype and trauma on limbic brain anatomy in childhood.

Childhood trauma is a major precipitating factor in psychiatric disease. Emerging data suggest that stress susceptibility is genetically determined, and that risk is mediated by changes in limbic brain circuitry. There is a need to identify markers of disease vulnerability, and it is critical that these markers be investigated in childhood and adolescence, a time when neural networks are particularly malleable and when psychiatric disorders frequently emerge. In this preliminary study, we evaluated whether a common variant in the brain-derived neurotrophic factor (BDNF) gene (Val66Met; rs6265) interacts with childhood trauma to predict limbic gray matter volume in a sample of 55 youth high in sociodemographic risk. We found trauma-by-BDNF interactions in the right subcallosal area and right hippocampus, wherein BDNF-related gray matter changes were evident in youth without histories of trauma. In youth without trauma exposure, lower hippocampal volume was related to higher symptoms of anxiety. These data provide preliminary evidence for a contribution of a common BDNF gene variant to the neural correlates of childhood trauma among high-risk urban youth. Altered limbic structure in early life may lay the foundation for longer term patterns of neural dysfunction, and hold implications for understanding the psychiatric and psychobiological consequences of traumatic stress on the developing brain.

Co-registration and distortion correction of diffusion and anatomical images based on inverse contrast normalization.

Diffusion MRI provides quantitative information about microstructural properties which can be useful in neuroimaging studies of the human brain. Echo planar imaging (EPI) sequences, which are frequently used for acquisition of diffusion images, are sensitive to inhomogeneities in the primary magnetic (B0) field that cause localized distortions in the reconstructed images. We describe and evaluate a new method for correction of susceptibility-induced distortion in diffusion images in the absence of an accurate B0 fieldmap. In our method, the distortion field is estimated using a constrained non-rigid registration between an undistorted T1-weighted anatomical image and one of the distorted EPI images from diffusion acquisition. Our registration framework is based on a new approach, INVERSION (Inverse contrast Normalization for VERy Simple registratION), which exploits the inverted contrast relationship between T1- and T2-weighted brain images to define a simple and robust similarity measure. We also describe how INVERSION can be used for rigid alignment of diffusion images and T1-weighted anatomical images. Our approach is evaluated with multiple in vivo datasets acquired with different acquisition parameters. Compared to other methods, INVERSION shows robust and consistent performance in rigid registration and shows improved alignment of diffusion and anatomical images relative to normalized mutual information for non-rigid distortion correction.

How age of acquisition influences brain architecture in bilinguals.

In the present study, we explored how Age of Acquisition (AoA) of L2 affected brain structures in bilingual individuals. Thirty-six native English speakers who were bilingual were scanned with high resolution MRI. After MRI signal intensity inhomogeneity correction, we applied both voxel-based morphometry (VBM) and surface-based morphometry (SBM) approaches to the data. VBM analysis was performed using FSL's standard VBM processing pipeline. For the SBM analysis, we utilized a semi-automated sulci delineation procedure, registered the brains to an atlas, and extracted measures of twenty four pre-selected regions of interest. We addressed three questions: (1) Which areas are more susceptible to differences in AoA? (2) How do AoA, proficiency and current level of exposure work together in predicting structural differences in the brain? And (3) What is the direction of the effect of AoA on regional volumetric and surface measures? Both VBM and SBM results suggested that earlier second language exposure was associated with larger volumes in the right parietal cortex. Consistently, SBM showed that the cortical area of the right superior parietal lobule increased as AoA decreased. In contrast, in the right pars orbitalis of the inferior frontal gyrus, AoA, proficiency, and current level of exposure are equally important in accounting for the structural differences. We interpret our results in terms of current theory and research on the effects of L2 learning on brain structures and functions.

Visualizing biological data-now and in the future.

Methods and tools for visualizing biological data have improved considerably over the last decades, but they are still inadequate for some high-throughput data sets. For most users, a key challenge is to benefit from the deluge of data without being overwhelmed by it. This challenge is still largely unfulfilled and will require the development of truly integrated and highly useable tools.

Visualization of image data from cells to organisms.

Advances in imaging techniques and high-throughput technologies are providing scientists with unprecedented possibilities to visualize internal structures of cells, organs and organisms and to collect systematic image data characterizing genes and proteins on a large scale. To make the best use of these increasingly complex and large image data resources, the scientific community must be provided with methods to query, analyze and crosslink these resources to give an intuitive visual representation of the data. This review gives an overview of existing methods and tools for this purpose and highlights some of their limitations and challenges.

BrainSuite BIDS App: Containerized Workflows for MRI Analysis.

There has been a concerted effort by the neuroimaging community to establish standards for computational methods for data analysis that promote reproducibility and portability. In particular, the Brain Imaging Data Structure (BIDS) specifies a standard for storing imaging data, and the related BIDS App methodology provides a standard for implementing containerized processing environments that include all necessary dependencies to process BIDS datasets using image processing workflows. We present the BrainSuite BIDS App, which encapsulates the core MRI processing functionality of BrainSuite within the BIDS App framework. Specifically, the BrainSuite BIDS App implements a participant-level workflow comprising three pipelines and a corresponding set of group-level analysis workflows for processing the participant-level outputs. The BrainSuite Anatomical Pipeline (BAP) extracts cortical surface models from a T1-weighted (T1w) MRI. It then performs surface-constrained volumetric registration to align the T1w MRI to a labeled anatomical atlas, which is used to delineate anatomical regions of interest in the MRI brain volume and on the cortical surface models. The BrainSuite Diffusion Pipeline (BDP) processes diffusion-weighted imaging (DWI) data, with steps that include coregistering the DWI data to the T1w scan, correcting for geometric image distortion, and fitting diffusion models to the DWI data. The BrainSuite Functional Pipeline (BFP) performs fMRI processing using a combination of FSL, AFNI, and BrainSuite tools. BFP coregisters the fMRI data to the T1w image, then transforms the data to the anatomical atlas space and to the Human Connectome Project's grayordinate space. Each of these outputs can then be processed during group-level analysis. The outputs of BAP and BDP are analyzed using the BrainSuite Statistics in R (bssr) toolbox, which provides functionality for hypothesis testing and statistical modeling. The outputs of BFP can be analyzed using atlas-based or atlas-free statistical methods during group-level processing. These analyses include the application of BrainSync, which synchronizes the time-series data temporally and enables comparison of resting-state or task-based fMRI data across scans. We also present the BrainSuite Dashboard quality control system, which provides a browser-based interface for reviewing the outputs of individual modules of the participant-level pipelines across a study in real-time as they are generated. BrainSuite Dashboard facilitates rapid review of intermediate results, enabling users to identify processing errors and make adjustments to processing parameters if necessary. The comprehensive functionality included in the BrainSuite BIDS App provides a mechanism for rapidly deploying the BrainSuite workflows into new environments to perform large-scale studies. We demonstrate the capabilities of the BrainSuite BIDS App using structural, diffusion, and functional MRI data from the Amsterdam Open MRI Collection's Population Imaging of Psychology dataset.

A hybrid high-resolution anatomical MRI atlas with sub-parcellation of cortical gyri using resting fMRI.

We present a new high-quality, single-subject atlas with sub-millimeter voxel resolution, high SNR, and excellent gray-white tissue contrast to resolve fine anatomical details. The atlas is labeled into two parcellation schemes: 1) the anatomical BCI-DNI atlas, which is manually labeled based on known morphological and anatomical features, and 2) the hybrid USCBrain atlas, which incorporates functional information to guide the sub-parcellation of cerebral cortex. In both cases, we provide consistent volumetric and cortical surface-based parcellation and labeling. The intended use of the atlas is as a reference template for structural coregistration and labeling of individual brains. A single-subject T1-weighted image was acquired five times at a resolution of 0.547 mm × 0.547 mm × 0.800 mm and averaged. Images were processed by an expert neuroanatomist using semi-automated methods in BrainSuite to extract the brain, classify tissue-types, and render anatomical surfaces. Sixty-six cortical and 29 noncortical regions were manually labeled to generate the BCI-DNI atlas. The cortical regions were further sub-parcellated into 130 cortical regions based on multi-subject connectivity analysis using resting fMRI (rfMRI) data from the Human Connectome Project (HCP) database to produce the USCBrain atlas. In addition, we provide a delineation between sulcal valleys and gyral crowns, which offer an additional set of 26 sulcal subregions per hemisphere. Lastly, a probabilistic map is provided to give users a quantitative measure of reliability for each gyral subdivision. Utility of the atlas was assessed by computing Adjusted Rand Indices (ARIs) between individual sub-parcellations obtained through structural-only coregistration to the USCBrain atlas and sub-parcellations obtained directly from each subject's resting fMRI data. Both atlas parcellations can be used with the BrainSuite, FreeSurfer, and FSL software packages.

Real time and delayed effects of subcortical low intensity focused ultrasound.

Deep brain nuclei are integral components of large-scale circuits mediating important cognitive and sensorimotor functions. However, because they fall outside the domain of conventional non-invasive neuromodulatory techniques, their study has been primarily based on neuropsychological models, limiting the ability to fully characterize their role and to develop interventions in cases where they are damaged. To address this gap, we used the emerging technology of non-invasive low-intensity focused ultrasound (LIFU) to directly modulate left lateralized basal ganglia structures in healthy volunteers. During sonication, we observed local and distal decreases in blood oxygenation level dependent (BOLD) signal in the targeted left globus pallidus (GP) and in large-scale cortical networks. We also observed a generalized decrease in relative perfusion throughout the cerebrum following sonication. These results show, for the first time using functional MRI data, the ability to modulate deep-brain nuclei using LIFU while measuring its local and global consequences, opening the door for future applications of subcortical LIFU.

Genetics of brain fiber architecture and intellectual performance.

The study is the first to analyze genetic and environmental factors that affect brain fiber architecture and its genetic linkage with cognitive function. We assessed white matter integrity voxelwise using diffusion tensor imaging at high magnetic field (4 Tesla), in 92 identical and fraternal twins. White matter integrity, quantified using fractional anisotropy (FA), was used to fit structural equation models (SEM) at each point in the brain, generating three-dimensional maps of heritability. We visualized the anatomical profile of correlations between white matter integrity and full-scale, verbal, and performance intelligence quotients (FIQ, VIQ, and PIQ). White matter integrity (FA) was under strong genetic control and was highly heritable in bilateral frontal (a(2)=0.55, p=0.04, left; a(2)=0.74, p=0.006, right), bilateral parietal (a(2)=0.85, p<0.001, left; a(2)=0.84, p<0.001, right), and left occipital (a(2)=0.76, p=0.003) lobes, and was correlated with FIQ and PIQ in the cingulum, optic radiations, superior fronto-occipital fasciculus, internal capsule, callosal isthmus, and the corona radiata (p=0.04 for FIQ and p=0.01 for PIQ, corrected for multiple comparisons). In a cross-trait mapping approach, common genetic factors mediated the correlation between IQ and white matter integrity, suggesting a common physiological mechanism for both, and common genetic determination. These genetic brain maps reveal heritable aspects of white matter integrity and should expedite the discovery of single-nucleotide polymorphisms affecting fiber connectivity and cognition.

Comparison of landmark-based and automatic methods for cortical surface registration.

Group analysis of structure or function in cerebral cortex typically involves, as a first step, the alignment of cortices. A surface-based approach to this problem treats the cortex as a convoluted surface and coregisters across subjects so that cortical landmarks or features are aligned. This registration can be performed using curves representing sulcal fundi and gyral crowns to constrain the mapping. Alternatively, registration can be based on the alignment of curvature metrics computed over the entire cortical surface. The former approach typically involves some degree of user interaction in defining the sulcal and gyral landmarks while the latter methods can be completely automated. Here we introduce a cortical delineation protocol consisting of 26 consistent landmarks spanning the entire cortical surface. We then compare the performance of a landmark-based registration method that uses this protocol with that of two automatic methods implemented in the software packages FreeSurfer and BrainVoyager. We compare performance in terms of discrepancy maps between the different methods, the accuracy with which regions of interest are aligned, and the ability of the automated methods to correctly align standard cortical landmarks. Our results show similar performance for ROIs in the perisylvian region for the landmark-based method and FreeSurfer. However, the discrepancy maps showed larger variability between methods in occipital and frontal cortex and automated methods often produce misalignment of standard cortical landmarks. Consequently, selection of the registration approach should consider the importance of accurate sulcal alignment for the specific task for which coregistration is being performed. When automatic methods are used, the users should ensure that sulci in regions of interest in their studies are adequately aligned before proceeding with subsequent analysis.

Sulcal set optimization for cortical surface registration.

Flat mapping based cortical surface registration constrained by manually traced sulcal curves has been widely used for inter subject comparisons of neuroanatomical data. Even for an experienced neuroanatomist, manual sulcal tracing can be quite time consuming, with the cost increasing with the number of sulcal curves used for registration. We present a method for estimation of an optimal subset of size N(C) from N possible candidate sulcal curves that minimizes a mean squared error metric over all combinations of N(C) curves. The resulting procedure allows us to estimate a subset with a reduced number of curves to be traced as part of the registration procedure leading to optimal use of manual labeling effort for registration. To minimize the error metric we analyze the correlation structure of the errors in the sulcal curves by modeling them as a multivariate Gaussian distribution. For a given subset of sulci used as constraints in surface registration, the proposed model estimates registration error based on the correlation structure of the sulcal errors. The optimal subset of constraint curves consists of the N(C) sulci that jointly minimize the estimated error variance for the subset of unconstrained curves conditioned on the N(C) constraint curves. The optimal subsets of sulci are presented and the estimated and actual registration errors for these subsets are computed.

Online resource for validation of brain segmentation methods.

One key issue that must be addressed during the development of image segmentation algorithms is the accuracy of the results they produce. Algorithm developers require this so they can see where methods need to be improved and see how new developments compare with existing ones. Users of algorithms also need to understand the characteristics of algorithms when they select and apply them to their neuroimaging analysis applications. Many metrics have been proposed to characterize error and success rates in segmentation, and several datasets have also been made public for evaluation. Still, the methodologies used in analyzing and reporting these results vary from study to study, so even when studies use the same metrics their numerical results may not necessarily be directly comparable. To address this problem, we developed a web-based resource for evaluating the performance of skull-stripping in T1-weighted MRI. The resource provides both the data to be segmented and an online application that performs a validation study on the data. Users may download the test dataset, segment it using whichever method they wish to assess, and upload their segmentation results to the server. The server computes a series of metrics, displays a detailed report of the validation results, and archives these for future browsing and analysis. We applied this framework to the evaluation of 3 popular skull-stripping algorithms--the Brain Extraction Tool [Smith, S.M., 2002. Fast robust automated brain extraction. Hum. Brain Mapp. 17 (3),143-155 (Nov)], the Hybrid Watershed Algorithm [Ségonne, F., Dale, A.M., Busa, E., Glessner, M., Salat, D., Hahn, H.K., Fischl, B., 2004. A hybrid approach to the skull stripping problem in MRI. NeuroImage 22 (3), 1060-1075 (Jul)], and the Brain Surface Extractor [Shattuck, D.W., Sandor-Leahy, S.R., Schaper, K.A., Rottenberg, D.A., Leahy, R.M., 2001. Magnetic resonance image tissue classification using a partial volume model. NeuroImage 13 (5), 856-876 (May) under several different program settings. Our results show that with proper parameter selection, all 3 algorithms can achieve satisfactory skull-stripping on the test data.

Active fibers: matching deformable tract templates to diffusion tensor images.

Reliable quantitative analysis of white matter connectivity in the brain is an open problem in neuroimaging, with common solutions requiring tools for fiber tracking, tractography segmentation and estimation of intersubject correspondence. This paper proposes a novel, template matching approach to the problem. In the proposed method, a deformable fiber-bundle model is aligned directly with the subject tensor field, skipping the fiber tracking step. Furthermore, the use of a common template eliminates the need for tractography segmentation and defines intersubject shape correspondence. The method is validated using phantom DTI data and applications are presented, including automatic fiber-bundle reconstruction and tract-based morphometry.

Mean diffusivity: a biomarker for CSF-related disease and genetic liability effects in schizophrenia.

Mean diffusivity (MD), the rotationally invariant magnitude of water diffusion that is greater in cerebrospinal fluid (CSF) and smaller in organized brain tissue, has been suggested to reflect schizophrenia-associated cortical atrophy. Regional changes, associations with CSF, and the effects of genetic predisposition towards schizophrenia, however, remain uncertain. Six-direction diffusion tensor imaging DTI and high-resolution structural images were obtained from 26 schizophrenia patients, 36 unaffected first-degree patient relatives, 20 control subjects and 32 control relatives (N=114). Registration procedures aligned diffusion tensor imaging (DTI) data across imaging modalities. MD was averaged within lobar regions and the cingulate and superior temporal gyri. CSF volume and MD were highly correlated. Significant bilateral temporal, and superior temporal MD increases were observed in schizophrenia compared with unrelated control probands. First-degree relatives of schizophrenia probands showed larger MD measures compared with controls within bilateral superior temporal regions with CSF volume correction. Superior temporal lobe brain tissue deficits and proximal CSF enlargements are widely documented in schizophrenia. Larger MD indices in patients and their relatives may thus reflect similar pathophysiological mechanisms. However, persistence of regional MD effects after controlling for CSF volume, suggests that MD is a sensitive biological marker of disease and genetic liability, characterizing at least partially distinct aspects of brain structural integrity.

A parameterization-based numerical method for isotropic and anisotropic diffusion smoothing on non-flat surfaces.

Neuroimaging data, such as 3-D maps of cortical thickness or neural activation, can often be analyzed more informatively with respect to the cortical surface rather than the entire volume of the brain. Any cortical surface-based analysis should be carried out using computations in the intrinsic geometry of the surface rather than using the metric of the ambient 3-D space. We present parameterization-based numerical methods for performing isotropic and anisotropic filtering on triangulated surface geometries. In contrast to existing FEM-based methods for triangulated geometries, our approach accounts for the metric of the surface. In order to discretize and numerically compute the isotropic and anisotropic geometric operators, we first parameterize the surface using a p-harmonic mapping. We then use this parameterization as our computational domain and account for the surface metric while carrying out isotropic and anisotropic filtering. To validate our method, we compare our numerical results to the analytical expression for isotropic diffusion on a spherical surface. We apply these methods to smoothing of mean curvature maps on the cortical surface, a step commonly required for analysis of gyrification or for registering surface-based maps across subjects.

Multiuser virtual reality environment for visualising neuroimaging data.

The recent advent of high-performance consumer virtual reality (VR) systems has opened new possibilities for immersive visualisation of numerous types of data. Medical imaging has long made use of advanced visualisation techniques, and VR offers exciting new opportunities for data exploration. The author presents a new framework for interacting with neuroimaging data, including MRI volumes, neuroanatomical surface models, diffusion tensors, and streamline tractography, as well as text-based annotations. The system was developed for the HTC Vive using C++, OpenGL, and the OpenVR software development kit. The author developed custom GLSL shaders for each type of data to provide high-performance real-time rendering suitable for use in a VR environment. These are integrated with an interface that enables the user to manipulate the scene through the Vive controllers and perform operations such as volume slicing, fibre track selection, and structural queries. The software can read data generated by existing automated brain MRI analysis packages, enabling the rapid development of subject-specific visualisations of multimodal data or annotated atlases. The system can also support multiple simultaneous users, placing them in the same virtual space to interact with each other while visualising the same datasets, opening new possibilities for teaching and for collaborative exploration of neuroimaging data.

VOTING-BASED SEGMENTATION OF OVERLAPPING NUCLEI IN CLARITY IMAGES.

New tissue-clearing techniques and improvements in optical microscopy have rapidly advanced capabilities to acquire volumetric imagery of neural tissue at resolutions of one micron or better. As sizes for data collections increase, accurate automatic segmentation of cell nuclei becomes increasingly important for quantitative analysis of imaged tissue. We present a cell nucleus segmentation method that is formulated as a parameter estimation problem with the goal of determining the count, shapes, and locations of nuclei that most accurately describe an image. We applied our new voting-based approach to fluorescence confocal microscopy images of neural tissue stained with DAPI, which highlights nuclei. Compared to manual counting of cells in three DAPI images, our method outperformed three existing approaches. On a manually labeled high-resolution DAPI image, our method also outperformed those methods and achieved a cell count accuracy of 98.99% and mean Dice coefficient of 0.6498.