Search

Emerging evidence suggests that repeated blast exposure (RBE) is associated with brain injury in military personnel. United States (U.S.) Special Operations Forces (SOF) personnel experience high rates of blast exposure during training and combat, but the effects of low-level RBE on brain structure and function in SOF have not been comprehensively characterized. Further, the pathophysiological link between RBE-related brain injuries and cognitive, behavioral, and physical symptoms has not been fully elucidated. We present a protocol for an observational pilot study, Long-Term Effects of Repeated Blast Exposure in U.S. SOF Personnel (ReBlast). In this exploratory study, 30 active-duty SOF personnel with RBE will participate in a comprehensive evaluation of: 1) brain network structure and function using Connectome magnetic resonance imaging (MRI) and 7 Tesla MRI; 2) neuroinflammation and tau deposition using positron emission tomography; 3) blood proteomics and metabolomics; 4) behavioral and physical symptoms using self-report measures; and 5) cognition using a battery of conventional and digitized assessments designed to detect subtle deficits in otherwise high-performing individuals. We will identify clinical, neuroimaging, and blood-based phenotypes that are associated with level of RBE, as measured by the Generalized Blast Exposure Value. Candidate biomarkers of RBE-related brain injury will inform the design of a subsequent study that will test a diagnostic assessment battery for detecting RBE-related brain injury. Ultimately, we anticipate that the ReBlast study will facilitate the development of interventions to optimize the brain health, quality of life, and battle readiness of U.S. SOF personnel.

Serial sectioning optical coherence tomography (OCT) enables accurate volumetric reconstruction of several cubic centimeters of human brain samples. We aimed to identify anatomical features of the ex vivo human brain, such as intraparenchymal blood vessels and axonal fiber bundles, from the OCT data in 3D, using intrinsic optical contrast.We developed an automatic processing pipeline to enable characterization of the intraparenchymal microvascular network in human brain samples.We demonstrated the automatic extraction of the vessels down to a 20 μm in diameter using a filtering strategy followed by a graphing representation and characterization of the geometrical properties of microvascular network in 3D. We also showed the ability to extend this processing strategy to extract axonal fiber bundles from the volumetric OCT image.This method provides a viable tool for quantitative characterization of volumetric microvascular network as well as the axonal bundle properties in normal and pathological tissues of the ex vivo human brain.

Background: Neurofibrillary tangle (NFT) accumulation in the entorhinal cortex (EC) precedes the transformation from cognitive controls to mild cognitive impairment and Alzheimer’s disease (AD). While tauopathy has been described in the EC before, the order and degree to which the individual subfields within the EC are engulfed by NFTs in aging and the preclinical AD stage is unknown. Objective: We aimed to investigate substructures within the EC to map the populations of cortical neurons most vulnerable to tau pathology in aging and the preclinical AD stage. Methods: We characterized phosphorylated tau (CP13) in 10 cases at eight well-defined anterior-posterior levels and assessed NFT density within the eight entorhinal subfields (described by Insausti and colleagues) at the preclinical stages of AD. We validated with immunohistochemistry and labeled the NFT density ratings on ex vivo MRIs. We measured subfield cortical thickness and reconstructed the labels as three-dimensional isosurfaces, resulting in anatomically comprehensive, histopathologically validated tau “heat maps.” Results: We found the lateral EC subfields ELc, ECL, and ECs (lateral portion) to have the highest tau density in semi-quantitative scores and quantitative measurements. We observed significant stepwise higher tau from anterior to posterior levels (p

Optical coherence tomography (OCT) is an emerging 3D imaging technique that allows quantification of intrinsic optical properties such as scattering coefficient and back-scattering coefficient, and has proved useful in distinguishing delicate microstructures in the human brain. The origins of scattering in brain tissues are contributed by the myelin content, neuron size and density primarily; however, no quantitative relationships between them have been reported, which hampers the use of OCT in fundamental studies of architectonic areas in the human brain and the pathological evaluations of diseases. Here, we built a generalized linear model based on Mie scattering theory that quantitatively links tissue scattering to myelin content and neuron density in the human brain. We report a strong linear relationship between scattering coefficient and the myelin content that is retained across different regions of the brain. Neuronal cell body turns out to be a secondary contribution to the overall scattering. The optical property of OCT provides a label-free solution for quantifying volumetric myelin content and neuron cells in the human brain.

A tool was developed to automatically segment several subcortical limbic structures (nucleus accumbens, basal forebrain, septal nuclei, hypothalamus without mammillary bodies, the mammillary bodies, and fornix) using only a T1-weighted MRI as input. This tool fills an unmet need as there are few, if any, publicly available tools to segment these clinically relevant structures. A U-Net with spatial, intensity, contrast, and noise augmentation was trained using 39 manually labeled MRI data sets. In general, the Dice scores, true positive rates, false discovery rates, and manual-automatic volume correlation were very good relative to comparable tools for other structures. A diverse data set of 698 subjects were segmented using the tool; evaluation of the resulting labelings showed that the tool failed in less than 1% of cases. Test-retest reliability of the tool was excellent. The automatically segmented volume of all structures except mammillary bodies showed effectiveness at detecting either clinical AD effects, age effects, or both. This tool will be publicly released with FreeSurfer (surfer.nmr.mgh.harvard.edu/fswiki/ScLimbic). Together with the other cortical and subcortical limbic segmentations, this tool will allow FreeSurfer to provide a comprehensive view of the limbic system in an automated way., Graphical Abstract

Joint registration of a stack of 2D histological sections to recover 3D structure (``3D histology reconstruction'') finds application in areas such as atlas building and validation of \emph{in vivo} imaging. Straightforward pairwise registration of neighbouring sections yields smooth reconstructions but has well-known problems such as ``banana effect'' (straightening of curved structures) and ``z-shift'' (drift). While these problems can be alleviated with an external, linearly aligned reference (e.g., Magnetic Resonance (MR) images), registration is often inaccurate due to contrast differences and the strong nonlinear distortion of the tissue, including artefacts such as folds and tears. In this paper, we present a probabilistic model of spatial deformation that yields reconstructions for multiple histological stains that that are jointly smooth, robust to outliers, and follow the reference shape. The model relies on a spanning tree of latent transforms connecting all the sections and slices of the reference volume, and assumes that the registration between any pair of images can be see as a noisy version of the composition of (possibly inverted) latent transforms connecting the two images. Bayesian inference is used to compute the most likely latent transforms given a set of pairwise registrations between image pairs within and across modalities. The framework is used for accurate 3D reconstruction of two stains (Nissl and parvalbumin) from the Allen human brain atlas, showing its benefits on real data with severe distortions. Moreover, we also provide the registration of the reconstructed volume to MNI space, bridging the gaps between two of the most widely used atlases in histology and MRI. The 3D reconstructed volumes and atlas registration can be downloaded from https://openneuro.org/datasets/ds003590. The code is freely available at https://github.com/acasamitjana/3dhirest., Comment: Medical Image Analysis (Accepted on 07/10/2021). 22 pages, 11 figures

The removal of non-brain signal from magnetic resonance imaging (MRI) data, known as skull-stripping, is an integral component of many neuroimage analysis streams. Despite their abundance, popular classical skull-stripping methods are usually tailored to images with specific acquisition properties, namely near-isotropic resolution and T1-weighted (T1w) MRI contrast, which are prevalent in research settings. As a result, existing tools tend to adapt poorly to other image types, such as stacks of thick slices acquired with fast spin-echo (FSE) MRI that are common in the clinic. While learning-based approaches for brain extraction have gained traction in recent years, these methods face a similar burden, as they are only effective for image types seen during the training procedure. To achieve robust skull-stripping across a landscape of imaging protocols, we introduce SynthStrip, a rapid, learning-based brain-extraction tool. By leveraging anatomical segmentations to generate an entirely synthetic training dataset with anatomies, intensity distributions, and artifacts that far exceed the realistic range of medical images, SynthStrip learns to successfully generalize to a variety of real acquired brain images, removing the need for training data with target contrasts. We demonstrate the efficacy of SynthStrip for a diverse set of image acquisitions and resolutions across subject populations, ranging from newborn to adult. We show substantial improvements in accuracy over popular skull-stripping baselines -- all with a single trained model. Our method and labeled evaluation data are available at https://w3id.org/synthstrip., 19 pages, 9 figures, 7 tables, skull stripping, brain extraction, image synthesis, MRI-contrast agnosticism, deep learning, final published version

BackgroundStructural brain connectivity has been shown to be sensitive to the changes that the brain undergoes during Alzheimer’s disease (AD) progression.MethodsIn this work, we used our recently proposed structural connectivity quantification measure derived from diffusion MRI, which accounts for both direct and indirect pathways, to quantify brain connectivity in dementia. We analyzed data from the ADNI-2 and OASIS-3 datasets to derive relevant information for the study of the changes that the brain undergoes in AD. We also compared these datasets to the HCP dataset, as a reference, and eventually validated externally on two cohorts of the EDSD database.ResultsOur analysis shows expected trends of mean conductance with respect to age and cognitive scores, significant age prediction values in aging data, and regional effects centered among sub-cortical regions, and cingulate and temporal cortices.DiscussionResults indicate that the conductance measure has prediction potential, especially for age, that age and cognitive scores largely overlap, and that this measure could be used to study effects such as anti-correlation in structural connections.Impact statementThis work presents a methodology and a set of analyses that open new possibilities in the study of healthy and pathological aging. The methodology used here is sensitive to direct and indirect pathways in deriving brain connectivity measures from dMRI, and therefore provides information that many state-of-the-art methods do not account for. As a result, this technique may provide the research community with ways to detect subtle effects of healthy aging and AD.

The default mode network (DMN) mediates self-awareness and introspection, core components of human consciousness. Therapies to restore consciousness in patients with severe brain injuries have historically targeted subcortical sites in the brainstem, thalamus, hypothalamus, basal forebrain, and basal ganglia, with the goal of reactivating cortical DMN nodes. However, the subcortical connectivity of the DMN has not been fully mapped and optimal subcortical targets for therapeutic neuromodulation of consciousness have not been identified. In this work, we created a comprehensive map of DMN subcortical connectivity by combining high-resolution functional and structural datasets with advanced signal processing methods. We analyzed 7 Tesla resting-state functional MRI (rs-fMRI) data from 168 healthy volunteers acquired in the Human Connectome Project. The rs-fMRI blood-oxygen-level-dependent (BOLD) data were temporally synchronized across subjects using the BrainSync algorithm. Cortical and subcortical DMN nodes were jointly analyzed and identified at the group level by applying a novel Nadam-Accelerated SCAlable and Robust (NASCAR) tensor decomposition method to the synchronized dataset. The subcortical connectivity map was then overlaid on a 7 Tesla 100 micron ex vivo MRI dataset for neuroanatomic analysis using automated segmentation of nuclei within the brainstem, thalamus, hypothalamus, basal forebrain, and basal ganglia. We further compared the NASCAR subcortical connectivity map with its counterpart generated from canonical seed-based correlation analyses. The NASCAR method revealed that BOLD signal in the central lateral nucleus of the thalamus and ventral tegmental area of the midbrain is strongly correlated with that of the DMN. In an exploratory analysis, additional subcortical sites in the median and dorsal raphe, lateral hypothalamus, and caudate nuclei were correlated with the cortical DMN. We also found that the putamen and globus pallidus are negatively correlated (i.e., anti-correlated) with the DMN, providing rs-fMRI evidence for the mesocircuit hypothesis of human consciousness, whereby a striatopallidal feedback system modulates anterior forebrain function via disinhibition of the central thalamus. Seed-based analyses yielded similar subcortical DMN connectivity, but the NASCAR result showed stronger contrast and better spatial alignment with dopamine immunostaining data. The DMN subcortical connectivity map identified here advances understanding of the subcortical regions that contribute to human consciousness and can be used to inform the selection of therapeutic targets in clinical trials for patients with disorders of consciousness.

The neuropathology associated with cognitive decline in military personnel exposed to traumatic brain injury (TBI) and chronic stress is incompletely understood. Few studies have examined clinicopathologic correlations between phosphorylated-tau neurofibrillary tangles, β-amyloid neuritic plaques, neuroinflammation, or white matter (WM) lesions, and neuropsychiatric disorders in veterans. We describe clinicopathologic findings in 4 military veterans with early-onset dementia (EOD) who had varying histories of blunt- and blast-TBI, cognitive decline, behavioral abnormalities, post-traumatic stress disorder, suicidal ideation, and suicide. We found that pathologic lesions in these military-EOD cases could not be categorized as classic Alzheimer’s disease (AD), chronic traumatic encephalopathy, traumatic axonal injury, or other well-characterized clinicopathologic entities. Rather, we observed a mixture of polypathology with unusual patterns compared with pathologies found in AD or other dementias. Also, ultrahigh resolution ex vivo MRI in 2 of these 4 brains revealed unusual patterns of periventricular WM injury. These findings suggest that military-EOD cases are associated with atypical combinations of brain lesions and distribution rarely seen in nonmilitary populations. Future prospective studies that acquire neuropsychiatric data before and after deployments, as well as genetic and environmental exposure data, are needed to further elucidate clinicopathologic correlations in military-EOD.

Continued improvement in MRI acquisition technology has made functional MRI (fMRI) with small isotropic voxel sizes down to 1 mm and below more commonly available. Although many conventional fMRI studies seek to investigate regional patterns of cortical activation for which conventional voxel sizes of 3 mm and larger provide sufficient spatial resolution, smaller voxels can help avoid contamination from adjacent white matter (WM) and cerebrospinal fluid (CSF), and thereby increase the specificity of fMRI to signal changes within the gray matter. Unfortunately, temporal signal-to-noise ratio (tSNR), a metric of fMRI sensitivity, is reduced in high-resolution acquisitions, which offsets the benefits of small voxels. Here we introduce a framework that combines small, isotropic fMRI voxels acquired at 7 T field strength with a novel anatomically-informed, surface mesh-navigated spatial smoothing that can provide both higher detection power and higher resolution than conventional voxel sizes. Our smoothing approach uses a family of intracortical surface meshes and allows for kernels of various shapes and sizes, including curved 3D kernels that adapt to and track the cortical folding pattern. Our goal is to restrict smoothing to the cortical gray matter ribbon and avoid noise contamination from CSF and signal dilution from WM via partial volume effects. We found that the intracortical kernel that maximizes tSNR does not maximize percent signal change (ΔS/S), and therefore the kernel configuration that optimizes detection power cannot be determined from tSNR considerations alone. However, several kernel configurations provided a favorable balance between boosting tSNR and ΔS/S, and allowed a 1.1-mm isotropic fMRI acquisition to have higher performance after smoothing (in terms of both detection power and spatial resolution) compared to an unsmoothed 3.0-mm isotropic fMRI acquisition. Overall, the results of this study support the strategy of acquiring voxels smaller than the cortical thickness, even for studies not requiring high spatial resolution, and smoothing them down within the cortical ribbon with a kernel of an appropriate shape to achieve the best performance—thus decoupling the choice of fMRI voxel size from the spatial resolution requirements of the particular study. The improvement of this new intracortical smoothing approach over conventional surface-based smoothing is expected to be modest for conventional resolutions, however the improvement is expected to increase with higher resolutions. This framework can also be applied to anatomically-informed intracortical smoothing of higher-resolution data (e.g. along columns and layers) in studies with prior information about the spatial structure of activation.

Connectomics has proved promising in quantifying and understanding the effects of development, aging and an array of diseases on the brain. In this work, we propose a new structural connectivity measure from diffusion MRI that allows us to incorporate direct brain connections, as well as indirect ones that would not be otherwise accounted for by standard techniques and that may be key for the better understanding of function from structure. From our experiments on the Human Connectome Project dataset, we find that our measure of structural connectivity better correlates with functional connectivity than streamline tractography does, meaning that it provides new structural information related to function. Through additional experiments on the ADNI-2 dataset, we demonstrate the ability of this new measure to better discriminate different stages of Alzheimer’s disease. Our findings suggest that this measure is useful in the study of the normal brain structure, and for quantifying the effects of disease on the brain structure.

Optical coherence tomography (OCT) images of ex vivo human brain tissue are corrupted by multiplicative speckle noise that degrades the contrast to noise ratio (CNR) of microstructural compartments. This work proposes a novel algorithm to reduce noise corruption in OCT images that minimizes the penalized negative log likelihood of gamma distributed speckle noise. The proposed method is formulated as a majorize-minimize problem that reduces to solving an iterative regularized least squares optimization. We demonstrate the usefulness of the proposed method by removing speckle in simulated data, phantom data and real OCT images of human brain tissue. We compare the proposed method with state of the art filtering and non-local means based denoising methods. We demonstrate that our approach removes speckle accurately, improves CNR between different tissue types and better preserves small features and edges in human brain tissue.

Diffusion tractography MRI can infer changes in network connectivity in patients with traumatic brain injury (TBI), but pathological substrates of disconnected tracts have not been well-defined due to a lack of high-resolution imaging with histopathological validation. We developed an ex vivo MRI protocol to analyze tract terminations at 750 μm resolution, followed by histopathologic evaluation of white matter pathology, and applied these methods to a 60-year-old man who died 26 days after TBI. Analysis of 74 cerebral hemispheric white matter regions revealed a heterogeneous distribution of tract disruptions. Associated histopathology identified variable white matter injury with patchy deposition of amyloid precursor protein and loss of neurofilament-positive axonal processes, myelin dissolution, astrogliosis, microgliosis, and perivascular hemosiderin-laden macrophages. Multiple linear regression revealed that tract disruption strongly correlated with neurofilament loss. Ex vivo diffusion MRI can detect tract disruptions in the human brain that reflect axonal injury.

The MarkVCID consortium was formed under cooperative agreements with the National Institute of Neurologic Diseases and Stroke (NINDS) and National Institute on Aging (NIA) in 2016 with the goals of developing and validating biomarkers for the cerebral small vessel diseases associated with the vascular contributions to cognitive impairment and dementia (VCID). Rigorously validated biomarkers have consistently been identified as crucial for multicenter studies to identify effective strategies to prevent and treat VCID, specifically to detect increased VCID risk, diagnose the presence of small vessel disease and its subtypes, assess prognosis for disease progression or response to treatment, demonstrate target engagement or mechanism of action for candidate interventions, and monitor disease progression during treatment. The seven project sites and central coordinating center comprising MarkVCID, working with NINDS and NIA, identified a panel of 11 candidate fluid- and neuroimaging-based biomarker kits and established harmonized multicenter study protocols (see companion paper "MarkVCID cerebral small vessel consortium: I. Enrollment, clinical, fluid protocols" for full details). Here we describe the MarkVCID neuroimaging protocols with specific focus on validating their application to future multicenter trials. MarkVCID procedures for participant enrollment; clinical and cognitive evaluation; and collection, handling, and instrumental validation of fluid samples are described in detail in a companion paper. Magnetic resonance imaging (MRI) has long served as the neuroimaging modality of choice for cerebral small vessel disease and VCID because of its sensitivity to a wide range of brain properties, including small structural lesions, connectivity, and cerebrovascular physiology. Despite MRI's widespread use in the VCID field, there have been relatively scant data validating the repeatability and reproducibility of MRI-based biomarkers across raters, scanner types, and time intervals (collectively defined as instrumental validity). The MRI protocols described here address the core MRI sequences for assessing cerebral small vessel disease in future research studies, specific sequence parameters for use across various research scanner types, and rigorous procedures for determining instrumental validity. Another candidate neuroimaging modality considered by MarkVCID is optical coherence tomography angiography (OCTA), a non-invasive technique for directly visualizing retinal capillaries as a marker of the cerebral capillaries. OCTA has theoretical promise as a unique opportunity to visualize small vessels derived from the cerebral circulation, but at a considerably earlier stage of development than MRI. The additional OCTA protocols described here address procedures for determining OCTA instrumental validity, evaluating sources of variability such as pupil dilation, and handling data to maintain participant privacy. MRI protocol and instrumental validation The core sequences selected for the MarkVCID MRI protocol are three-dimensional T1-weighted multi-echo magnetization-prepared rapid-acquisition-of-gradient-echo (ME-MPRAGE), three-dimensional T2-weighted fast spin echo fluid-attenuated-inversion-recovery (FLAIR), two-dimensional diffusion-weighted spin-echo echo-planar imaging (DWI), three-dimensional T2*-weighted multi-echo gradient echo (3D-GRE), three-dimensional T2 -weighted fast spin-echo imaging (T2w), and two-dimensional T2*-weighted gradient echo echo-planar blood-oxygenation-level-dependent imaging with brief periods of CO2 inhalation (BOLD-CVR). Harmonized parameters for each of these core sequences were developed for four 3 Tesla MRI scanner models in widespread use at academic medical centers. MarkVCID project sites are trained and certified for their instantiation of the consortium MRI protocols. Sites are required to perform image quality checks every 2 months using the Alzheimer's Disease Neuroimaging Initiative phantom. Instrumental validation for MarkVCID MRI-based biomarkers is operationally defined as inter-rater reliability, test-retest repeatability, and inter-scanner reproducibility. Assessments of these instrumental properties are performed on individuals representing a range of cerebral small vessel disease from mild to severe. Inter-rater reliability is determined by distribution of an independent dataset of MRI scans to each analysis site. Test-retest repeatability is determined by repeat MRI scans performed on individual participants on a single MRI scanner after a short (1- to 14-day) interval. Inter-scanner reproducibility is determined by repeat MRI scans performed on individuals performed across four MRI scanner models. OCTA protocol and instrumental validation The MarkVCID OCTA protocol uses a commercially available, Food and Drug Administration-approved OCTA apparatus. Imaging is performed on one dilated and one undilated eye to assess the need for dilation. Scans are performed in quadruplicate. MarkVCID project sites participating in OCTA validation are trained and certified by this biomarker's lead investigator. Inter-rater reliability for OCTA is assessed by distribution of OCTA datasets to each analysis site. Test-retest repeatability is assessed by repeat OCTA imaging on individuals on the same day as their baseline OCTA and a different-day repeat session after a short (1- to 14-day) interval. Methods were developed to allow the OCTA data to be de-identified by the sites before transmission to the central data management system. The MarkVCID neuroimaging protocols, like the other MarkVCID procedures, are designed to allow translation to multicenter trials and as a template for outside groups to generate directly comparable neuroimaging data. The MarkVCID neuroimaging protocols are available to the biomedical community and intended to be shared. In addition to the instrumental validation procedures described here, each of the neuroimaging MarkVCID kits will undergo biological validation to determine its ability to measure important aspects of VCID such as cognitive function. The analytic methods for the neuroimaging-based kits and the results of these validation studies will be published separately. The results will ultimately determine the neuroimaging kits' potential usefulness for multicenter interventional trials in small vessel disease-related VCID.

Accurate and reliable whole-brain segmentation is critical to longitudinal neuroimaging studies. We undertake a comparative analysis of two subcortical segmentation methods, Automatic Segmentation (ASEG) and Sequence Adaptive Multimodal Segmentation (SAMSEG), recently provided in the open-source neuroimaging package FreeSurfer 7.1, with regard to reliability, bias, sensitivity to detect longitudinal change, and diagnostic sensitivity to Alzheimer’s disease. First, we assess intra- and inter-scanner reliability for eight bilateral subcortical structures: amygdala, caudate, hippocampus, lateral ventricles, nucleus accumbens, pallidum, putamen and thalamus. For intra-scanner analysis we use a large sample of participants (n = 1629) distributed across the lifespan (age range = 4–93 years) and acquired on a 1.5T Siemens Avanto (n = 774) and a 3T Siemens Skyra (n = 855) scanners. For inter-scanner analysis we use a sample of 24 participants scanned on the day with three models of Siemens scanners: 1.5T Avanto, 3T Skyra and 3T Prisma. Second, we test how each method detects volumetric age change using longitudinal follow up scans (n = 491 for Avanto and n = 245 for Skyra; interscan interval = 1–10 years). Finally, we test sensitivity to clinically relevant change. We compare annual rate of hippocampal atrophy in cognitively normal older adults (n = 20), patients with mild cognitive impairment (n = 20) and Alzheimer’s disease (n = 20). We find that both ASEG and SAMSEG are reliable and lead to the detection of within-person longitudinal change, although with notable differences between age-trajectories for most structures, including hippocampus and amygdala. In summary, SAMSEG yields significantly lower differences between repeated measures for intra- and inter-scanner analysis without compromising sensitivity to changes and demonstrating ability to detect clinically relevant longitudinal changes.

Traditional neuroimage analysis pipelines involve computationally intensive, time-consuming optimization steps, and thus, do not scale well to large cohort studies with thousands or tens of thousands of individuals. In this work we propose a fast and accurate deep learning based neuroimaging pipeline for the automated processing of structural human brain MRI scans, replicating FreeSurfer's anatomical segmentation including surface reconstruction and cortical parcellation. To this end, we introduce an advanced deep learning architecture capable of whole brain segmentation into 95 classes. The network architecture incorporates local and global competition via competitive dense blocks and competitive skip pathways, as well as multi-slice information aggregation that specifically tailor network performance towards accurate segmentation of both cortical and sub-cortical structures. Further, we perform fast cortical surface reconstruction and thickness analysis by introducing a spectral spherical embedding and by directly mapping the cortical labels from the image to the surface. This approach provides a full FreeSurfer alternative for volumetric analysis (in under 1 minute) and surface-based thickness analysis (within only around 1h runtime). For sustainability of this approach we perform extensive validation: we assert high segmentation accuracy on several unseen datasets, measure generalizability and demonstrate increased test-retest reliability, and high sensitivity to group differences in dementia., Comment: Submitted to NeuroImage

The development of automated tools for brain morphometric analysis in infants has lagged significantly behind analogous tools for adults. This gap reflects the greater challenges in this domain due to: 1) a smaller-scaled region of interest, 2) increased motion corruption, 3) regional changes in geometry due to heterochronous growth, and 4) regional variations in contrast properties corresponding to ongoing myelination and other maturation processes. Nevertheless, there is a great need for automated image-processing tools to quantify differences between infant groups and other individuals, because aberrant cortical morphologic measurements (including volume, thickness, surface area, and curvature) have been associated with neuropsychiatric, neurologic, and developmental disorders in children. In this paper we present an automated segmentation and surface extraction pipeline designed to accommodate clinical MRI studies of infant brains in a population 0-2 year-olds. The algorithm relies on a single channel of T1-weighted MR images to achieve automated segmentation of cortical and subcortical brain areas, producing volumes of subcortical structures and surface models of the cerebral cortex. We evaluated the algorithm both qualitatively and quantitatively using manually labeled datasets, relevant comparator software solutions cited in the literature, and expert evaluations. The computational tools and atlases described in this paper will be distributed to the research community as part of the FreeSurfer image analysis package., Comment: 49 pages, 25 figures, submitted to NeuroImage

BackgroundWe sought to develop an automatable score to predict hospitalization, critical illness, or death in patients at risk for COVID-19 presenting for urgent care during the Massachusetts outbreak.MethodsSingle-center study of adult outpatients seen in respiratory illness clinics (RICs) or the emergency department (ED), including development (n = 9381, March 7-May 2) and prospective (n = 2205, May 3-14) cohorts. Data was queried from Partners Enterprise Data Warehouse. Outcomes were hospitalization, critical illness or death within 7 days. We developed the COVID-19 Acuity Score (CoVA) using automatically extracted data from the electronic medical record and learning-to-rank ordinal logistic regression modeling. Calibration was assessed using predicted-to-observed ratio (E/O). Discrimination was assessed by C-statistics (AUC).ResultsIn the development cohort, 27.3%, 7.2%, and 1.1% of patients experienced hospitalization, critical illness, or death, respectively; and in the prospective cohort, 26.1%, 6.3%, and 0.5%. CoVA showed excellent performance in the development cohort (concurrent validation) for hospitalization (E/O: 1.00, AUC: 0.80); for critical illness (E/O: 1.00, AUC: 0.82); and for death (E/O: 1.00, AUC: 0.87). Performance in the prospective cohort (prospective validation) was similar for hospitalization (E/O: 1.01, AUC: 0.76); for critical illness (E/O 1.03, AUC: 0.79); and for death (E/O: 1.63, AUC=0.93). Among 30 predictors, the top five were age, diastolic blood pressure, blood oxygen saturation, COVID-19 testing status, and respiratory rate.ConclusionsCoVA is a prospectively validated automatable score to assessing risk for adverse outcomes related to COVID-19 infection in the outpatient setting.

The human hippocampus is vulnerable to a range of degenerative conditions and as such, accurate in vivo measurement of the hippocampus and hippocampal substructures via neuroimaging is of great interest for understanding mechanisms of disease as well as for use as a biomarker in clinical trials of novel therapeutics. Although total hippocampal volume can be measured relatively reliably, it is critical to understand how this reliability is affected by acquisition on different scanners, as multiple scanning platforms would likely be utilized in large-scale clinical trials. This is particularly true for hippocampal subregional measurements, which have only relatively recently been measurable through common image processing platforms such as FreeSurfer. Accurate segmentation of these subregions is challenging due to their small size, magnetic resonance imaging (MRI) signal loss in medial temporal regions of the brain, and lack of contrast for delineation from standard neuroimaging procedures.Here, we assess the test-retest reliability of the FreeSurfer automated hippocampal subfield segmentation procedure using two Siemens model scanners (a Siemens Trio and Prismafit Trio upgrade). T1-weighted images were acquired for 11 generally healthy younger participants (two scans on the Trio and one scan on the Prismafit). Each scan was processed through the standard cross-sectional stream and the recently released longitudinal pipeline in FreeSurfer v6.0 for hippocampal segmentation. Test-retest reliability of the volumetric measures was examined for individual subfields as well as percent volume difference and Dice overlap among scans and intra-class correlation coefficients (ICC). Reliability was high in the molecular layer, dentate gyrus, and whole hippocampus with the inclusion of three time points with mean volume differences among scans less than 3%, overlap greater than 80%, and ICC >0.95. The parasubiculum and hippocampal fissure showed the least improvement in reliability with mean volume difference greater than 5%, overlap less than 70%, and ICC scores ranging from 0.78 to 0.89. Other subregions, including the CA regions, were stable in their mean volume difference and overlap (75% respectively) and showed improvement in reliability with the inclusion of three scans (ICC ​> ​0.9). Reliability was generally higher within scanner (Trio-Trio), however, Trio-Prismafit reliability was also high and did not exhibit an obvious bias. These results suggest that the FreeSurfer automated segmentation procedure is a reliable method to measure total as well as hippocampal subregional volumes and may be useful in clinical applications including as an endpoint for future clinical trials of conditions affecting the hippocampus.

Non-rigid cortical registration is an important and challenging task due to the geometric complexity of the human cortex and the high degree of inter-subject variability. A conventional solution is to use a spherical representation of surface properties and perform registration by aligning cortical folding patterns in that space. This strategy produces accurate spatial alignment but often requires a high computational cost. Recently, convolutional neural networks (CNNs) have demonstrated the potential to dramatically speed up volumetric registration. However, due to distortions introduced by projecting a sphere to a 2D plane, a direct application of recent learning-based methods to surfaces yields poor results. In this study, we present SphereMorph, a diffeomorphic registration framework for cortical surfaces using deep networks that addresses these issues. SphereMorph uses a UNet-style network associated with a spherical kernel to learn the displacement field and warps the sphere using a modified spatial transformer layer. We propose a resampling weight in computing the data fitting loss to account for distortions introduced by polar projection, and demonstrate the performance of our proposed method on two tasks, including cortical parcellation and group-wise functional area alignment. The experiments show that the proposed SphereMorph is capable of modeling the geometric registration problem in a CNN framework and demonstrate superior registration accuracy and computational efficiency. The source code of SphereMorph will be released to the public upon acceptance of this manuscript at https://github.com/voxelmorph/spheremorph., Comment: cortical surface registration, deep network, unsupervised learning, registration, deep learning, cortical, spherical, invertible

Tuberous sclerosis complex (TSC), a heritable neurodevelopmental disorder, is caused by mutations in the TSC1 or TSC2 genes. To date, there has been little work to elucidate regional TSC1 and TSC2 gene expression within the human brain, how it changes with age, and how it may influence disease. Using a publicly available microarray dataset, we found that TSC1 and TSC2 gene expression was highest within the adult neo-cerebellum and that this pattern of increased cerebellar expression was maintained throughout postnatal development. During mid-gestational fetal development, however, TSC1 and TSC2 expression was highest in the cortical plate. Using a bioinformatics approach to explore protein and genetic interactions, we confirmed extensive connections between TSC1/TSC2 and the other genes that comprise the mammalian target of rapamycin (mTOR) pathway, and show that the mTOR pathway genes with the highest connectivity are also selectively expressed within the cerebellum. Finally, compared to age-matched controls, we found increased cerebellar volumes in pediatric TSC patients without current exposure to antiepileptic drugs. Considered together, these findings suggest that the cerebellum may play a central role in TSC pathogenesis and may contribute to the cognitive impairment, including the high incidence of autism spectrum disorder, observed in the TSC population.

The false positive rates (FPR) for surface-based group analysis of cortical thickness, surface area, and volume were evaluated for parametric and non-parametric clusterwise correction for multiple comparisons for a range of smoothing levels and cluster-forming thresholds (CFT) using real data under group assignments that should not yield significant results. For whole cortical surface analysis, thickness showed modest inflation in parametric FPRs above the nominal level (10% versus 5%). Surface area and volume FPRs were much higher (20–30%). In the analysis of interhemispheric thickness asymmetries, FPRs were well controlled by parametric correction, but FPRs for surface area and volume asymmetries were still inflated. In all cases, non-parametric permutation adequately controlled the FPRs. It was found that inflated parametric FPRs were caused by violations in the parametric assumptions, namely a heavier-than-Gaussian spatial correlation. The non-Gaussian spatial correlation originates from anatomical features unique to individuals (e.g., a patch of cortex slightly thicker or thinner than average) and is not a by-product of scanning or processing. Thickness performed better than surface area and volume because thickness does not require a Jacobian correction.

Mesopontine tegmental nuclei such as the cuneiform, pedunculotegmental, oral pontine reticular, paramedian raphe and caudal linear raphe nuclei, are deep brain structures involved in arousal and motor function. Dysfunction of these nuclei is implicated in the pathogenesis of disorders of consciousness and sleep, as well as in neurodegenerative diseases. However, their localization in conventional neuroimages of living humans is difficult due to limited image sensitivity and contrast, and a stereotaxic probabilistic neuroimaging template of these nuclei in humans does not exist. We used semi-automatic segmentation of single-subject 1.1 mm-isotropic 7 Tesla diffusion-fractional-anisotropy and T(2)-weighted images in healthy adults to generate an in vivo probabilistic neuroimaging structural template of these nuclei in standard stereotaxic (Montreal Neurological Institute, MNI) space. The template was validated through independent manual delineation, as well as leave-one-out validation and evaluation of nuclei volumes. This template can enable localization of five mesopontine tegmental nuclei in conventional images (e.g. 1.5 Tesla, 3 Tesla) in future studies of arousal and motor physiology (e.g. sleep, anesthesia, locomotion) and pathology (e.g. disorders of consciousness, sleep disorders, Parkinson’s disease). The 7 Tesla magnetic resonance imaging procedure for single-subject delineation of these nuclei may also prove useful for future 7 Tesla studies of arousal and motor mechanisms.

Functional MRI (fMRI) benefits from both increased sensitivity and specificity with increasing magnetic field strength, making it a key application for Ultra-High Field (UHF) MRI scanners. Most UHF-fMRI studies utilize the dramatic increases in sensitivity and specificity to acquire high-resolution data reaching sub-millimeter scales, which enable new classes of experiments to probe the functional organization of the human brain. This review article surveys advanced data analysis strategies developed for high-resolution fMRI at UHF. These include strategies designed to mitigate distortion and artifacts associated with higher fields in ways that attempt to preserve spatial resolution of the fMRI data, as well as recently introduced analysis techniques that are enabled by these extremely high-resolution data. Particular focus is placed on anatomically-informed analyses, including cortical surface-based analysis, which are powerful techniques that can guide each step of the analysis from preprocessing to statistical analysis to interpretation and visualization. New intracortical analysis techniques for laminar and columnar fMRI are also reviewed and discussed. Prospects for single-subject individualized analyses are also presented and discussed. Altogether, there are both specific challenges and opportunities presented by UHF-fMRI, and the use of proper analysis strategies can help these valuable data reach their full potential.

The surface of the human cerebellar cortex is much more tightly folded than the cerebral cortex. Volumetric analysis of cerebellar morphometry in magnetic resonance imaging studies suffers from insufficient resolution, and therefore has had limited impact on disease assessment. Automatic serial polarization-sensitive optical coherence tomography (as-PSOCT) is an emerging technique that offers the advantages of microscopic resolution and volumetric reconstruction of large-scale samples. In this study, we reconstructed multiple cubic centimeters of ex vivo human cerebellum tissue using as-PSOCT. The morphometric and optical properties of the cerebellar cortex across five subjects were quantified. While the molecular and granular layers exhibited similar mean thickness in the five subjects, the thickness varied greatly in the granular layer within subjects. Layer-specific optical property remained homogenous within individual subjects but showed higher cross-subject variability than layer thickness. High-resolution volumetric morphometry and optical property maps of human cerebellar cortex revealed by as-PSOCT have great potential to advance our understanding of cerebellar function and diseases.

While many useful microstructural indices, as well as orientation distribution functions, can be obtained from multi-shell dMRI data, there is growing interest in exploring the richer set of microstructural features that can be extracted from the full ensemble average propagator (EAP). The EAP can be readily computed from diffusion spectrum imaging (DSI) data, at the cost of a very lengthy acquisition. Compressed sensing (CS) has been used to make DSI more practical by reducing its acquisition time. CS applied to DSI (CS-DSI) attempts to reconstruct the EAP from significantly undersampled q-space data. We present a post mortem validation study where we evaluate the ability of CS-DSI to approximate not only fully sampled DSI but also multi-shell acquisitions with high fidelity. Human brain samples are imaged with high-resolution DSI at 9.4T and with polarization-sensitive optical coherence tomography (PSOCT). The latter provides direct measurements of axonal orientations at microscopic resolutions, allowing us to evaluate the mesoscopic orientation estimates obtained from diffusion MRI, in terms of their angular error and the presence of spurious peaks. We test two fast, dictionary-based, L2-regularized algorithms for CS-DSI reconstruction. We find that, for a CS acceleration factor of R=3, i.e., an acquisition with 171 gradient directions, one of these methods is able to achieve both low angular error and low number of spurious peaks. With a scan length similar to that of high angular resolution multi-shell acquisition schemes, this CS-DSI approach is able to approximate both fully sampled DSI and multi-shell data with high accuracy. Thus it is suitable for orientation reconstruction and microstructural modeling techniques that require either grid- or shell-based acquisitions. We find that the signal-to-noise ratio (SNR) of the training data used to construct the dictionary can have an impact on the accuracy of CS-DSI, but that there is substantial robustness to loss of SNR in the test data. Finally, we show that, as the CS acceleration factor increases beyond R=3, the accuracy of these reconstruction methods degrade, either in terms of the angular error, or in terms of the number of spurious peaks. Our results provide useful benchmarks for the future development of even more efficient q-space acceleration techniques.

Most existing algorithms for automatic 3D morphometry of human brain MRI scans are designed for data with near-isotropic voxels at approximately 1 mm resolution, and frequently have contrast constraints as well-typically requiring T1-weighted images (e.g., MP-RAGE scans). This limitation prevents the analysis of millions of MRI scans acquired with large inter-slice spacing in clinical settings every year. In turn, the inability to quantitatively analyze these scans hinders the adoption of quantitative neuro imaging in healthcare, and also precludes research studies that could attain huge sample sizes and hence greatly improve our understanding of the human brain. Recent advances in convolutional neural networks (CNNs) are producing outstanding results in super-resolution and contrast synthesis of MRI. However, these approaches are very sensitive to the specific combination of contrast, resolution and orientation of the input images, and thus do not generalize to diverse clinical acquisition protocols – even within sites. In this article, we present SynthSR, a method to train a CNN that receives one or more scans with spaced slices, acquired with different contrast, resolution and orientation, and produces an isotropic scan of canonical contrast (typically a 1 mm MP-RAGE). The presented method does not require any preprocessing, beyond rigid coregistration of the input scans. Crucially, SynthSR trains on synthetic input images generated from 3D segmentations, and can thus be used to train CNNs for any combination of contrasts, resolutions and orientations without high-resolution real images of the input contrasts. We test the images generated with SynthSR in an array of common downstream analyses, and show that they can be reliably used for subcortical segmentation and volumetry, image registration (e.g., for tensor-based morphometry), and, if some image quality requirements are met, even cortical thickness morphometry. The source code is publicly available at https://github.com/BBillot/SynthSR .

There is great value to use of structural neuroimaging in the assessment of Alzheimer's disease (AD). However, to date, predictive value of structural imaging tend to range between 80% and 90% in accuracy and it is unclear why this is the case given that structural imaging should parallel the pathologic processes of AD. There is a possibility that clinical misdiagnosis relative to the gold standard pathologic diagnosis and/or additional brain pathologies are confounding factors contributing to reduced structural imaging classification accuracy. We examined potential factors contributing to misclassification of individuals with clinically diagnosed AD purely from cortical thickness measures. Correctly classified and incorrectly classified groups were compared across a range of demographic, biological, and neuropsychological data including cerebrospinal fluid biomarkers, amyloid imaging, white matter hyperintensity (WMH) volume, cognitive, and genetic factors. Individual subject analyses suggested that at least a portion of the control individuals misclassified as AD from structural imaging additionally harbor substantial AD biomarker pathology and risk, yet are relatively resistant to cognitive symptoms, likely due to "cognitive reserve," and therefore clinically unimpaired. In contrast, certain clinical control individuals misclassified as AD from cortical thickness had increased WMH volume relative to other controls in the sample, suggesting that vascular conditions may contribute to classification accuracy from cortical thickness measures. These results provide examples of factors that contribute to the accuracy of structural imaging in predicting a clinical diagnosis of AD, and provide important information about considerations for future work aimed at optimizing structural based diagnostic classifiers for AD.

Background White matter signal abnormalities (WMSA) (also known as 'hyperintensities') on MRI are commonly seen in normal aging and increases have been noted in Alzheimer's disease (AD), but whether there is a spatial specificity to these increases is unknown. Objective To discern whether or not there is a spatial pattern of WMSA in the brains of individuals with AD that differs from those who exhibit cognitively healthy aging. Method Structural MRI data from the Alzheimer's Disease Neuroimaging Initiative public database were used to quantify WMSA in 35 regions of interest (ROIs). Regional measures were compared between cognitively healthy older controls (OC; n = 107) and individuals with a clinical diagnosis of AD (n = 127). Regional WMSA volume was also assessed in individuals with mild cognitive impairment (MCI; n = 74) who were 6, 12, and 24 months away from AD conversion. Results WMSA volume was significantly greater in AD compared to OC in 24 out of 35 ROIs after controlling for age, and nine were significantly higher after normalizing for total WMSA. Regions with greater WMSA volume in AD included rostral frontal, inferior temporal, and inferior parietal WM. In MCI, frontal and temporal regions demonstrated significantly greater WMSA volume with decreasing time-to-AD-conversion. Discussion Individuals with AD have greater regional volume of WMSA compared to OC regardless of age or total WMSA volume. Accumulation of regional WMSA is linked to time to AD conversion in individuals with MCI. These findings indicate WMSA is an important pathological component of AD development.

White matter lesions, quantified as ‘white matter signal abnormalities’ (WMSA) on neuroimaging, are common incidental findings on brain images of older adults. This tissue damage is linked to cerebrovascular dysfunction and is associated with cognitive decline. The regional distribution of WMSA throughout the cerebral white matter has been described at a gross scale; however, to date no prior study has described regional patterns relative to cortical gyral landmarks which may be important for understanding functional impact. Additionally, no prior study has described how regional WMSA volume scales with total global WMSA. Such information could be used in the creation of a pathologic ‘staging’ of WMSA through a detailed regional characterization at the individual level. Magnetic resonance imaging data from 97 cognitively-healthy older individuals (OC) aged 52–90 from the Alzheimer's Disease Neuroimaging Initiative (ADNI) study were processed using a novel WMSA labeling procedure described in our prior work. WMSA were quantified regionally using a procedure that segments the cerebral white matter into 35 bilateral units based on proximity to landmarks in the cerebral cortex. An initial staging was performed by quantifying the regional WMSA volume in four groups based on quartiles of total WMSA volume (quartiles I–IV). A consistent spatial pattern of WMSA accumulation was observed with increasing quartile. A clustering procedure was then used to distinguish regions based on patterns of scaling of regional WMSA to global WMSA. Three patterns were extracted that showed high, medium, and non-scaling with global WMSA. Regions in the high-scaling cluster included periventricular, caudal and rostral middle frontal, inferior and superior parietal, supramarginal, and precuneus white matter. A data-driven staging procedure was then created based on patterns of WMSA scaling and specific regional cut-off values from the quartile analyses. Individuals with Alzheimer's disease (AD) and mild cognitive impairment (MCI) were then additionally staged, and significant differences in the percent of each diagnostic group in Stages I and IV were observed, with more AD individuals residing in Stage IV and more OC and MCI individuals residing in Stage I. These data demonstrate a consistent regional scaling relationship between global and regional WMSA that can be used to classify individuals into one of four stages of white matter disease. White matter staging could play an important role in a better understanding and the treatment of cerebrovascular contributions to brain aging and dementia., Highlights • White signal abnormalities (WMSA) manifest in a stereotyped juxtacortical spatial pattern in cognitively healthy aging • We present a procedure for staging WMSA using cut-off values derived from statistical analyses • We present a method to assess how regional WMSA scales with global WMSA to identify vulnerable white matter regions • Regions at cerebral perfusion boundary zones scale higher with global WMSA load than regions with a dual blood supply • There are more individuals with Alzheimer's disease than cognitively heathy older controls in higher stages of WMSA

We present an ultra-high resolution MRI dataset of an ex vivo human brain specimen. The brain specimen was donated by a 58-year-old woman who had no history of neurological disease and died of non-neurological causes. After fixation in 10% formalin, the specimen was imaged on a 7 Tesla MRI scanner at 100 µm isotropic resolution using a custom-built 31-channel receive array coil. Single-echo multi-flip Fast Low-Angle SHot (FLASH) data were acquired over 100 hours of scan time (25 hours per flip angle), allowing derivation of synthesized FLASH volumes. This dataset provides an unprecedented view of the three-dimensional neuroanatomy of the human brain. To optimize the utility of this resource, we warped the dataset into standard stereotactic space. We now distribute the dataset in both native space and stereotactic space to the academic community via multiple platforms. We envision that this dataset will have a broad range of investigational, educational, and clinical applications that will advance understanding of human brain anatomy in health and disease., Measurement(s)nuclear magnetic resonance assayTechnology Type(s)MRI ScannerFactor Type(s)flip angleSample Characteristic - OrganismHomo sapiens Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.9958688

With the advent of convolutional neural networks~(CNN), supervised learning methods are increasingly being used for whole brain segmentation. However, a large, manually annotated training dataset of labeled brain images required to train such supervised methods is frequently difficult to obtain or create. In addition, existing training datasets are generally acquired with a homogeneous magnetic resonance imaging~(MRI) acquisition protocol. CNNs trained on such datasets are unable to generalize on test data with different acquisition protocols. Modern neuroimaging studies and clinical trials are necessarily multi-center initiatives with a wide variety of acquisition protocols. Despite stringent protocol harmonization practices, it is very difficult to standardize the gamut of MRI imaging parameters across scanners, field strengths, receive coils etc., that affect image contrast. In this paper we propose a CNN-based segmentation algorithm that, in addition to being highly accurate and fast, is also resilient to variation in the input acquisition. Our approach relies on building approximate forward models of pulse sequences that produce a typical test image. For a given pulse sequence, we use its forward model to generate plausible, synthetic training examples that appear as if they were acquired in a scanner with that pulse sequence. Sampling over a wide variety of pulse sequences results in a wide variety of augmented training examples that help build an image contrast invariant model. Our method trains a single CNN that can segment input MRI images with acquisition parameters as disparate as $T_1$-weighted and $T_2$-weighted contrasts with only $T_1$-weighted training data. The segmentations generated are highly accurate with state-of-the-art results~(overall Dice overlap$=0.94$), with a fast run time~($\approx$ 45 seconds), and consistent across a wide range of acquisition protocols., Typo in author name corrected. Greves -> Greve

No abstract available.

Online publication 15/07/2016 The hippocampal formation is a complex, heterogeneous structure that consists of a number of distinct, interacting subregions. Atrophy of these subregions is implied in a variety of neurodegenerative diseases, most prominently in Alzheimer’s disease (AD). Thanks to the increasing resolution ofMRimages and computational atlases, automatic segmentation of hippocampal subregions is becoming feasible in MRI scans. Here we introduce a generative model for dedicated longitudinal segmentation that relies on subject-specific atlases. The segmentations of the scans at the different time points are jointly computed using Bayesian inference. All time points are treated the same to avoid processing bias. We evaluate this approach using over 4700 scans from two publicly available datasets (ADNI and MIRIAD). In test–retest reliability experiments, the proposed method yielded significantly lower volume differences and significantly higher Dice overlaps than the cross-sectional approach for nearly every subregion (average across subregions: 4.5% vs. 6.5%, Dice overlap: 81.8% vs. 75.4%). The longitudinal algorithm also demonstrated increased sensitivity to group differences: in MIRIAD (69 subjects: 46 with AD and 23 controls), it found differences in atrophy rates between AD and controls that the cross sectional method could not detect in a number of subregions: right parasubiculum, left and right presubiculum, right subiculum, left dentate gyrus, left CA4, left HATA and right tail. In ADNI (836 subjects: 369 with AD, 215 with early cognitive impairment — eMCI — and 252 controls), all methods found significant differences between AD and controls, but the proposed longitudinal algorithm detected differences between controls and eMCI and differences between eMCI and AD that the cross sectional method could not find: left presubiculum, right subiculum, left and right parasubiculum, left and right HATA. Moreover, many of the differences that the cross-sectional method already found were detected with higher significance. The presented algorithm will be made available as part of the open-source neuroimaging package FreeSurfer. This project has received funding from the European Union's Horizon 2020 research and innovation program under the Marie Sklodowska-Curie grant agreement No.654911 (project “THALAMODEL”), and also from the Spanish Ministry of Economy and Competitiveness (MINECO, reference TEC2014-51882-P). Support for this research was also provided in part by the National Cancer Institute (1K25-CA181632-01), the Genentech Foundation (G-40819) and the Nvidia corporation, which donated a Titan X GPU. Further support was provided by the A.A. Martinos Center for Biomedical Imaging (P41RR014075, P41EB015896, U24RR021382), and was made possible by the resources provided by Shared Instrumentation Grants1S10RR023401, 1S10RR019307, and 1S10RR023043. Support was also provided by the National Institute for Biomedical Imaging and Bioengineering (R01EB006758, R21EB018907, R01EB019956), the National Institute on Aging (5R01AG008122, R01AG016495), the National Institute for Neurological Disorders and Stroke (R01NS0525851, R21NS072652, R01NS070963, R01NS083534, 5U01NS086625) and the Lundbeck Foundation (R141-2013-13117), Additional support was provided by the NIH Blueprint for Neuroscience Research (5U01-MH093765), as part of the multi-institutional Human Connectome Project. In addition, BF has a financial interest in CorticoMetrics, a company whose medical pursuits focus on brain imaging and measurement technologies. BF's interests were reviewed and are managed by Massachusetts General Hospital and Partners HealthCare in accordance with their conflict of interest policies. The collection and sharing of the MRI data used in the group study based on ADNI was funded by the Alzheimer's Disease Neuroimaging Initiative (National Institutes of Health Grant U01 AG024904) and DOD ADNI (U.S. Department of Defense award number W81XWH-12-2-0012). ADNI is funded by the National Institute on Aging, the National Institute of Biomedical Imaging and Bioengineering, and through generous contributions from the following: Alzheimer's Association; Alzheimer's Drug Discovery Foundation; BioClinica, Inc.; Biogen Idec Inc.; Bristol-Myers Squibb Company; Eisai Inc.; Elan Pharmaceuticals, Inc.; Eli Lilly and Company; F. Hoffmann-La Roche Ltd and its affiliated company Genentech, Inc.; GE Healthcare; Innogenetics, N.V.; IXICO Ltd.; Janssen Alzheimer Immunotherapy Research & Development, LLC.; Johnson & Johnson Pharmaceutical Research & Development LLC.; Medpace, Inc.; Merck & Co., Inc.; Meso Scale Diagnostics, LLC.; NeuroRx Research; Novartis Pharmaceuticals Corporation; Pfizer Inc.; Piramal Imaging; Servier; Synarc Inc.; and Takeda Pharmaceutical Company. The Canadian Institutes of Health Research is providing funds to support ADNI clinical sites in Canada. Private sector contributions are facilitated by the Foundation for the National Institutes of Health (www.fnih.org). The grantee organization is the Northern California Institute for Research and Education, and the study is coordinated by the Alzheimer's Disease Cooperative Study at the University of California, San Diego. ADNI data are disseminated by the Laboratory for Neuro Imaging at the University of Southern California.

Diffusion magnetic resonance imaging (dMRI) is a unique technology that allows the noninvasive quantification of microstructural tissue properties of the human brain in healthy subjects as well as the probing of disease-induced variations. Population studies of dMRI data have been essential in identifying pathological structural changes in various conditions, such as Alzheimer's and Huntington's diseases (Salat et al., 2010; Rosas et al., 2006). The most common form of dMRI involves fitting a tensor to the underlying imaging data (known as diffusion tensor imaging, or DTI), then deriving parametric maps, each quantifying a different aspect of the underlying microstructure, e.g. fractional anisotropy and mean diffusivity. To date, the statistical methods utilized in most DTI population studies either analyzed only one such map or analyzed several of them, each in isolation. However, it is most likely that variations in the microstructure due to pathology or normal variability would affect several parameters simultaneously, with differing variations modulating the various parameters to differing degrees. Therefore, joint analysis of the available diffusion maps can be more powerful in characterizing histopathology and distinguishing between conditions than the widely used univariate analysis. In this article, we propose a multivariate approach for statistical analysis of diffusion parameters that uses partial least squares correlation (PLSC) analysis and permutation testing as building blocks in a voxel-wise fashion. Stemming from the common formulation, we present three different multivariate procedures for group analysis, regressing-out nuisance parameters and comparing effects of different conditions. We used the proposed procedures to study the effects of non-demented aging, Alzheimer's disease and mild cognitive impairment on the white matter. Here, we present results demonstrating that the proposed PLSC-based approach can differentiate between effects of different conditions in the same region as well as uncover spatial variations of effects across the white matter. The proposed procedures were able to answer questions on structural variations such as: "are there regions in the white matter where Alzheimer's disease has a different effect than aging or similar effect as aging?" and "are there regions in the white matter that are affected by both mild cognitive impairment and Alzheimer's disease but with differing multivariate effects?"

Recent work has demonstrated that subject motion produces systematic biases in the metrics computed by widely used morphometry software packages, even when the motion is too small to produce noticeable image artefacts. In the common situation where the control population exhibits different behaviours in the scanner when compared to the experimental population, these systematic measurement biases may produce significant confounds for between-group analyses, leading to erroneous conclusions about group differences. While previous work has shown that prospective motion correction can improve perceived image quality, here we demonstrate that, in healthy subjects performing a variety of directed motions, the use of the volumetric navigator (vNav) prospective motion correction system significantly reduces the motion-induced bias and variance in morphometry.

In the first study comparing high angular resolution diffusion MRI (dMRI) in the human brain to axonal orientation measurements from polarization-sensitive optical coherence tomography (PSOCT), we compare the accuracy of orientation estimates from various dMRI sampling schemes and reconstruction methods. We find that, if the reconstruction approach is chosen carefully, single-shell dMRI data can yield the same accuracy as multi-shell data, and only moderately lower accuracy than a full Cartesian-grid sampling scheme. Our results suggest that current dMRI reconstruction approaches do not benefit substantially from ultra-high b-values or from very large numbers of diffusion-encoding directions. We also show that accuracy remains stable across dMRI voxel sizes of 1 ​mm or smaller but degrades at 2 ​mm, particularly in areas of complex white-matter architecture. We also show that, as the spatial resolution is reduced, axonal configurations in a dMRI voxel can no longer be modeled as a small set of distinct axon populations, violating an assumption that is sometimes made by dMRI reconstruction techniques. Our findings have implications for in vivo studies and illustrate the value of PSOCT as a source of ground-truth measurements of white-matter organization that does not suffer from the distortions typical of histological techniques.

The Human Connectome Projects in Development (HCP-D) and Aging (HCP-A) are two large-scale brain imaging studies that will extend the recently completed HCP Young-Adult (HCP-YA) project to nearly the full lifespan, collecting structural, resting-state fMRI, task-fMRI, diffusion, and perfusion MRI in participants from 5 to 100+years of age. HCP-D is enrolling 1300+healthy children, adolescents, and young adults (ages 5-21), and HCP-A is enrolling 1200+healthy adults (ages 36-100+), with each study collecting longitudinal data in a subset of individuals at particular age ranges. The imaging protocols of the HCP-D and HCP-A studies are very similar, differing primarily in the selection of different task-fMRI paradigms. We strove to harmonize the imaging protocol to the greatest extent feasible with the completed HCP-YA (1200+participants, aged 22-35), but some imaging-related changes were motivated or necessitated by hardware changes, the need to reduce the total amount of scanning per participant, and/or the additional challenges of working with young and elderly populations. Here, we provide an overview of the common HCP-D/A imaging protocol including data and rationales for protocol decisions and changes relative to HCP-YA. The result will be a large, rich, multi-modal, and freely available set of consistently acquired data for use by the scientific community to investigate and define normative developmental and aging related changes in the healthy human brain.

The original Human Connectome Project yielded a rich data set on structural and functional connectivity in a large sample of healthy young adults using improved methods of data acquisition, analysis, and sharing. More recent efforts are extending this approach to include infants, children, older adults, and brain disorders. This paper introduces and describes the Human Connectome Project in Aging (HCP-A), which is currently recruiting 1200+healthy adults aged 36 to 100+, with a subset of 600+participants returning for longitudinal assessment. Four acquisition sites using matched Siemens Prisma 3T MRI scanners with centralized quality control and data analysis are enrolling participants. Data are acquired across multimodal imaging and behavioral domains with a focus on factors known to be altered in advanced aging. MRI acquisitions include structural (whole brain and high resolution hippocampal) plus multiband resting state functional (rfMRI), task fMRI (tfMRI), diffusion MRI (dMRI), and arterial spin labeling (ASL). Behavioral characterization includes cognitive (such as processing speed and episodic memory), psychiatric, metabolic, and socioeconomic measures as well as assessment of systemic health (with a focus on menopause via hormonal assays). This dataset will provide a unique resource for examining how brain organization and connectivity changes across typical aging, and how these differences relate to key characteristics of aging including alterations in hormonal status and declining memory and general cognition. A primary goal of the HCP-A is to make these data freely available to the scientific community, supported by the Connectome Coordination Facility (CCF) platform for data quality assurance, preprocessing and basic analysis, and shared via the NIMH Data Archive (NDA). Here we provide the rationale for our study design and sufficient details of the resource for scientists to plan future analyses of these data. A companion paper describes the related Human Connectome Project in Development (HCP-D, Somerville etal., 2018), and the image acquisition protocol common to both studies (Harms etal., 2018).

Developmental neuroimaging studies report the emergence of increasingly diverse cognitive functions as closely entangled with a rise-fall modulation of cortical thickness (CTh), structural cortical and white-matter connectivity, and a time-course for the experience-dependent selective elimination of the overproduced synapses. We examine which of two visual processing networks, the dorsal (DVN; prefrontal, parietal nodes) or ventral (VVN; frontal-temporal, fusiform nodes) matures first, thus leading the neuro-cognitive developmental trajectory. Three age-dependent measures are reported: (i) the CTh at network nodes; (ii) the matrix of intra-network structural connectivity (edges); and (iii) the proficiency in network-related neuropsychological tests. Typically developing children (age ~6 years), adolescents (~11 years), and adults (~21 years) were tested using multiple-acquisition structural T1-weighted magnetic resonance imaging (MRI) and neuropsychology. MRI images reconstructed into a gray/white/pial matter boundary model were used for CTh evaluation. No significant group differences in CTh and in the matrix of edges were found for DVN (except for the left prefrontal), but a significantly thicker cortex in children for VVN with reduced prefrontal ventral-fusiform connectivity and with an abundance of connections in adolescents. The higher performance in children on tests related to DVN corroborates the age-dependent MRI structural connectivity findings. The current findings are consistent with an earlier maturational course of DVN.

Image segmentation is a critical step in numerous medical imaging studies, which can be facilitated by automatic computational techniques. Supervised methods, although highly effective, require large training datasets of manually labeled images that are labor-intensive to produce. Unsupervised methods, on the contrary, can be used in the absence of training data to segment new images. We introduce a new approach to unsupervised image segmentation that is based on the computation of the local center of mass. We propose an efficient method to group the pixels of a one-dimensional signal, which we then use in an iterative algorithm for two- and three-dimensional image segmentation. We validate our method on a 2D X-ray image, a 3D abdominal magnetic resonance (MR) image and a dataset of 3D cardiovascular MR images.

PURPOSE: To integrate markerless head motion tracking with prospectively-corrected neuroanatomical MRI sequences and to investigate high-frequency motion correction during imaging echo-trains. METHODS: A commercial 3D surface tracking system, which estimates head motion by registering point cloud reconstructions of the face, was used to adapt the imaging field-of-view (FOV) based on head movement during MPRAGE and T2-SPACE (3D variable-flip-angle TSE) sequences. The FOV position and orientation were updated every 6 lines of k-space (

Epidemiological studies suggest that a single moderate-to-severe traumatic brain injury (TBI) is associated with an increased risk of neurodegenerative disease, including Alzheimer's disease (AD) and Parkinson's disease (PD). Histopathological studies describe complex neurodegenerative pathologies in individuals exposed to single moderate-to-severe TBI or repetitive mild TBI, including chronic traumatic encephalopathy (CTE). However, the clinicopathological links between TBI and post-traumatic neurodegenerative diseases such as AD, PD, and CTE remain poorly understood. Here, we describe the methodology of the Late Effects of TBI (LETBI) study, whose goals are to characterize chronic post-traumatic neuropathology and to identify in vivo biomarkers of post-traumatic neurodegeneration. LETBI participants undergo extensive clinical evaluation using National Institutes of Health TBI Common Data Elements, proteomic and genomic analysis, structural and functional magnetic resonance imaging (MRI), and prospective consent for brain donation. Selected brain specimens undergo ultra-high resolution ex vivo MRI and histopathological evaluation including whole–mount analysis. Co-registration of ex vivo and in vivo MRI data enables identification of ex vivo lesions that were present during life. In vivo signatures of postmortem pathology are then correlated with cognitive and behavioral data to characterize the clinical phenotype(s) associated with pathological brain lesions. We illustrate the study methods and demonstrate proof of concept for this approach by reporting results from the first LETBI participant, who despite the presence of multiple in vivo and ex vivo pathoanatomic lesions had normal cognition and was functionally independent until her mid-80s. The LETBI project represents a multidisciplinary effort to characterize post-traumatic neuropathology and identify in vivo signatures of postmortem pathology in a prospective study.

The results of most neuroimaging studies are reported in volumetric (e.g., MNI152) or surface (e.g., fsaverage) coordinate systems. Accurate mappings between volumetric and surface coordinate systems can facilitate many applications, such as projecting fMRI group analyses from MNI152/Colin27 to fsaverage for visualization, or projecting resting-state fMRI parcellations from fsaverage to MNI152/Colin27 for volumetric analysis of new data. However, there has been surprisingly little research on this topic. Here, we evaluated three approaches for mapping data between MNI152/Colin27 and fsaverage coordinate systems by simulating the above applications: projection of group-average data from MNI152/Colin27 to fsaverage and projection of fsaverage parcellations to MNI152/Colin27. Two of the approaches are currently widely used. A third approach (registration fusion) was previously proposed, but not widely adopted. Two implementations of the registration fusion (RF) approach were considered, with one implementation utilizing the Advanced Normalization Tools (ANTs). We found that RF-ANTs performed the best for mapping between fsaverage and MNI152/Colin27, even for new subjects registered to MNI152/Colin27 using a different software tool (FSL FNIRT). This suggests that RF-ANTs would be useful even for researchers not using ANTs. Finally, it is worth emphasizing that the most optimal approach for mapping data to a coordinate system (e.g., fsaverage) is to register individual subjects directly to the coordinate system, rather than via another coordinate system. Only in scenarios where the optimal approach is not possible (e.g., mapping previously published results from MNI152 to fsaverage), should the approaches evaluated in this manuscript be considered. In these scenarios, we recommend RF-ANTs (https://github.com/ThomasYeoLab/CBIG/tree/master/stable_projects/registration/Wu2017_RegistrationFusion).

Nonlinear registration of 2D histological sections with corresponding slices of MRI data is a critical step of 3D histology reconstruction algorithms. This registration is difficult due to the large differences in image contrast and resolution, as well as the complex nonrigid deformations and artefacts produced when sectioning the sample and mounting it on the glass slide. It has been shown in brain MRI registration that better spatial alignment across modalities can be obtained by synthesising one modality from the other and then using intra-modality registration metrics, rather than by using information theory based metrics to solve the problem directly. However, such an approach typically requires a database of aligned images from the two modalities, which is very difficult to obtain for histology and MRI. Here, we overcome this limitation with a probabilistic method that simultaneously solves for deformable registration and synthesis directly on the target images, without requiring any training data. The method is based on a probabilistic model in which the MRI slice is assumed to be a contrast-warped, spatially deformed version of the histological section. We use approximate Bayesian inference to iteratively refine the probabilistic estimate of the synthesis and the registration, while accounting for each other's uncertainty. Moreover, manually placed landmarks can be seamlessly integrated in the framework for increased performance and robustness. Experiments on a synthetic dataset of MRI slices show that, compared with mutual information based registration, the proposed method makes it possible to use a much more flexible deformation model in the registration to improve its accuracy, without compromising robustness. Moreover, our framework also exploits information in manually placed landmarks more efficiently than mutual information: landmarks constrain the deformation field in both methods, but in our algorithm, it also has a positive effect on the synthesis - which further improves the registration. We also show results on two real, publicly available datasets: the Allen and BigBrain atlases. In both of them, the proposed method provides a clear improvement over mutual information based registration, both qualitatively (visual inspection) and quantitatively (registration error measured with pairs of manually annotated landmarks). ispartof: MEDICAL IMAGE ANALYSIS vol:50 pages:127-144 ispartof: location:Netherlands status: published

Recent advances in MR technology have enabled increased spatial resolution for routine functional and anatomical imaging, which has created demand for software tools that are able to process these data. The availability of high-resolution data also raises the question of whether higher resolution leads to substantial gains in accuracy of quantitative morphometric neuroimaging procedures, in particular the cortical surface reconstruction and cortical thickness estimation. In this study we adapted the FreeSurfer cortical surface reconstruction pipeline to process structural data at native submillimeter resolution. We then quantified the differences in surface placement between meshes generated from (0.75 mm)3 isotropic resolution data acquired in 39 volunteers and the same data downsampled to the conventional 1 mm3 voxel size. We find that when processed at native resolution, cortex is estimated to be thinner in most areas, but thicker around the Cingulate and the Calcarine sulci as well as in the posterior bank of the Central sulcus. Thickness differences are driven by two kinds of effects. First, the gray–white surface is found closer to the white matter, especially in cortical areas with high myelin content, and thus low contrast, such as the Calcarine and the Central sulci, causing local increases in thickness estimates. Second, the gray–CSF surface is placed more interiorly, especially in the deep sulci, contributing to local decreases in thickness estimates. We suggest that both effects are due to reduced partial volume effects at higher spatial resolution. Submillimeter voxel sizes can therefore provide improved accuracy for measuring cortical thickness.

We introduce BrainPrint, a compact and discriminative representation of brain morphology. BrainPrint captures shape information of an ensemble of cortical and subcortical structures by solving the eigenvalue problem of the 2D and 3D Laplace–Beltrami operator on triangular (boundary) and tetrahedral (volumetric) meshes. This discriminative characterization enables new ways to study the similarity between brains; the focus can either be on a specific brain structure of interest or on the overall brain similarity. We highlight four applications for BrainPrint in this article: (i) subject identification, (ii) age and sex prediction, (iii) brain asymmetry analysis, and (iv) potential genetic influences on brain morphology. The properties of BrainPrint require the derivation of new algorithms to account for the heterogeneous mix of brain structures with varying discriminative power. We conduct experiments on three datasets, including over 3000 MRI scans from the ADNI database, 436 MRI scans from the OASIS dataset, and 236 MRI scans from the VETSA twin study. All processing steps for obtaining the compact representation are fully automated, making this processing framework particularly attractive for handling large datasets., National Cancer Institute (U.S.) (1K25-CA181632-01), Athinoula A. Martinos Center for Biomedical Imaging (P41-RR014075), Athinoula A. Martinos Center for Biomedical Imaging (P41-EB015896), National Alliance for Medical Image Computing (U.S.) (U54-EB005149), Neuroimaging Analysis Center (U.S.) (P41-EB015902), National Center for Research Resources (U.S.) (U24 RR021382), National Institute of Biomedical Imaging and Bioengineering (U.S.) (5P41EB015896-15), National Institute of Biomedical Imaging and Bioengineering (U.S.) (R01EB006758), National Institute on Aging (AG022381), National Institute on Aging (5R01AG008122-22), National Institute on Aging (AG018344), National Institute on Aging (AG018386), National Center for Complementary and Alternative Medicine (U.S.) (RC1 AT005728-01), National Institute of Neurological Diseases and Stroke (U.S.) (R01 NS052585-01), National Institute of Neurological Diseases and Stroke (U.S.) (1R21NS072652-01), National Institute of Neurological Diseases and Stroke (U.S.) (1R01NS070963), National Institute of Neurological Diseases and Stroke (U.S.) (R01NS083534), National Institutes of Health (U.S.) ((5U01-MH093765)

In this paper we present a novel label fusion algorithm suited for scenarios in which different manual delineation protocols with potentially disparate structures have been used to annotate the training scans (hereafter referred to as "atlases"). Such scenarios arise when atlases have missing structures, when they have been labeled with different levels of detail, or when they have been taken from different heterogeneous databases. The proposed algorithm can be used to automatically label a novel scan with any of the protocols from the training data. Further, it enables us to generate new labels that are not present in any delineation protocol by defining intersections on the underling labels. We first use probabilistic models of label fusion to generalize three popular label fusion techniques to the multi-protocol setting: majority voting, semi-locally weighted voting and STAPLE. Then, we identify some shortcomings of the generalized methods, namely the inability to produce meaningful posterior probabilities for the different labels (majority voting, semi-locally weighted voting) and to exploit the similarities between the atlases (all three methods). Finally, we propose a novel generative label fusion model that can overcome these drawbacks. We use the proposed method to combine four brain MRI datasets labeled with different protocols (with a total of 102 unique labeled structures) to produce segmentations of 148 brain regions. Using cross-validation, we show that the proposed algorithm outperforms the generalizations of majority voting, semi-locally weighted voting and STAPLE (mean Dice score 83%, vs. 77%, 80% and 79%, respectively). We also evaluated the proposed algorithm in an aging study, successfully reproducing some well-known results in cortical and subcortical structures. (C) 2014 The Authors. Published by Elsevier Inc.

The choice of a reference image typically influences the results of deformable image registration, thereby making it asymmetric. This is a consequence of a spatially non-uniform weighting in the cost function integral that leads to general registration inaccuracy. The inhomogeneous integral measure – which is the local volume change in the transformation, thus varying through the course of the registration – causes image regions to contribute differently to the objective function. More importantly, the optimization algorithm is allowed to minimize the cost function by manipulating the volume change, instead of aligning the images. The approaches that restore symmetry to deformable registration successfully achieve inverse-consistency, but do not eliminate the regional bias that is the source of the error. In this work, we address the root of the problem: the non-uniformity of the cost function integral. We introduce a new quasi-volume-preserving constraint that allows for volume change only in areas with well-matching image intensities, and show that such a constraint puts a bound on the error arising from spatial non-uniformity. We demonstrate the advantages of adding the proposed constraint to standard (asymmetric and symmetrized) demons and diffeomorphic demons algorithms through experiments on synthetic images, and real X-ray and 2D/3D brain MRI data. Specifically, the results show that our approach leads to image alignment with more accurate matching of manually defined neuroanatomical structures, better tradeoff between image intensity matching and registration-induced distortion, improved native symmetry, and lower susceptibility to local optima. In summary, the inclusion of this space- and time-varying constraint leads to better image registration along every dimension that we have measured it.