Diffusion imaging and tractography enable mapping structural connections in the human brain, in-vivo. Linear Fascicle Evaluation (LiFE) is a state-of-the-art approach for pruning spurious connections in the estimated structural connectome, by optimizing its fit to the measured diffusion data. Yet, LiFE imposes heavy demands on computing time, precluding its use in analyses of large connectome databases. Here, we introduce a GPU-based implementation of LiFE that achieves 50-100x speedups over conventional CPU-based implementations for connectome sizes of up to several million fibers. Briefly, the algorithm accelerates generalized matrix multiplications on a compressed tensor through efficient GPU kernels, while ensuring favorable memory access patterns. Leveraging these speedups, we advance LiFE’s algorithm by imposing a regularization constraint on estimated fiber weights during connectome pruning. Our regularized, accelerated, LiFE algorithm (“ReAl-LiFE”) estimates sparser connectomes that also provide more accurate fits to the underlying diffusion signal. We demonstrate the utility of our approach by classifying pathological signatures of structural connectivity in patients with Alzheimer’s Disease (AD). We estimated million fiber whole-brain connectomes, followed by pruning with ReAl-LiFE, for 90 individuals (45 AD patients and 45 healthy controls). Linear classifiers, based on support vector machines, achieved over 80% accuracy in classifying AD patients from healthy controls based on their ReAl-LiFE pruned structural connectomes alone. Moreover, classification based on the ReAl-LiFE pruned connectome outperformed both the unpruned connectome, as well as the LiFE pruned connectome, in terms of accuracy. We propose our GPU-accelerated approach as a widely relevant tool for non-negative least-squares optimization, across many domains.
more »
« less
CHIASM, the human brain albinism and achiasma MRI dataset
Abstract We describe a collection of T1-, diffusion- and functional T2*-weighted magnetic resonance imaging data from human individuals with albinism and achiasma. This repository can be used as a test-bed to develop and validate tractography methods like diffusion-signal modeling and fiber tracking as well as to investigate the properties of the human visual system in individuals with congenital abnormalities. The MRI data is provided together with tools and files allowing for its preprocessing and analysis, along with the data derivatives such as manually curated masks and regions of interest for performing tractography.
more »
« less
- NSF-PAR ID:
- 10304441
- Date Published:
- Journal Name:
- Scientific Data
- Volume:
- 8
- Issue:
- 1
- ISSN:
- 2052-4463
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We propose a novel matrix autoencoder to map functional connectomes from resting state fMRI (rs-fMRI) to structural connectomes from Diffusion Tensor Imaging (DTI), as guided by subject-level phenotypic measures. Our specialized autoencoder infers a low dimensional manifold embedding for the rs-fMRI correlation matrices that mimics a canonical outer-product decomposition. The embedding is simultaneously used to reconstruct DTI tractography matrices via a second manifold alignment decoder and to predict inter-subject phenotypic variability via an artificial neural network. We validate our framework on a dataset of 275 healthy individuals from the Human Connectome Project database and on a second clinical dataset consisting of 57 subjects with Autism Spectrum Disorder. We demonstrate that the model reliably recovers structural connectivity patterns across individuals, while robustly extracting predictive and interpretable brain biomarkers in a cross-validated setting. Finally, our framework outperforms several baselines at predicting behavioral phenotypes in both real-world datasets.more » « less
-
The human sense of smell plays an important role in appetite and food intake, detecting environmental threats, social interactions, and memory processing. However, little is known about the neural circuity supporting its function. The olfactory tracts project from the olfactory bulb along the base of the frontal cortex, branching into several striae to meet diverse cortical regions. Historically, using diffusion magnetic resonance imaging (dMRI) to reconstruct the human olfactory tracts has been prevented by susceptibility and motion artifacts. Here, we used a dMRI method with readout segmentation of long variable echo-trains (RESOLVE) to minimize image distortions and characterize the human olfactory tracts in vivo . We collected high-resolution dMRI data from 25 healthy human participants (12 male and 13 female) and performed probabilistic tractography using constrained spherical deconvolution (CSD). At the individual subject level, we identified the lateral, medial, and intermediate striae with their respective cortical connections to the piriform cortex and amygdala (AMY), olfactory tubercle (OT), and anterior olfactory nucleus (AON). We combined individual results across subjects to create a normalized, probabilistic atlas of the olfactory tracts. We then investigated the relationship between olfactory perceptual scores and measures of white matter integrity, including mean diffusivity (MD). Importantly, we found that olfactory tract MD negatively correlated with odor discrimination performance. In summary, our results provide a detailed characterization of the connectivity of the human olfactory tracts and demonstrate an association between their structural integrity and olfactory perceptual function. SIGNIFICANCE STATEMENT This study provides the first detailed in vivo description of the cortical connectivity of the three olfactory tract striae in the human brain, using diffusion magnetic resonance imaging (dMRI). Additionally, we show that tract microstructure correlates with performance on an odor discrimination task, suggesting a link between the structural integrity of the olfactory tracts and odor perception. Lastly, we generated a normalized probabilistic atlas of the olfactory tracts that may be used in future research to study its integrity in health and disease.more » « less
-
more » « less
Background A lack of in utero imaging data hampers our understanding of the connections in the human fetal brain. Generalizing observations from postmortem subjects and premature newborns is inaccurate due to technical and biological differences.
Purpose To evaluate changes in fetal brain structural connectivity between 23 and 35 weeks postconceptional age using a spatiotemporal atlas of diffusion tensor imaging (DTI).
Study Type Retrospective.
Population Publicly available diffusion atlases, based on 60 healthy women (age 18–45 years) with normal prenatal care, from 23 and 35 weeks of gestation.
Field Strength/Sequence 3.0 Tesla/DTI acquired with diffusion‐weighted echo planar imaging (EPI).
Assessment We performed whole‐brain fiber tractography from DTI images. The cortical plate of each diffusion atlas was segmented and parcellated into 78 regions derived from the Edinburgh Neonatal Atlas (ENA33). Connectivity matrices were computed, representing normalized fiber connections between nodes. We examined the relationship between global efficiency (GE), local efficiency (LE), small‐worldness (SW), nodal efficiency (NE), and betweenness centrality (BC) with gestational age (GA) and with laterality.
Statistical Tests Linear regression was used to analyze changes in GE, LE, NE, and BC throughout gestation, and to assess changes in laterality. The
t ‐tests were used to assess SW.P ‐values were corrected using Holm‐Bonferroni method. A correctedP ‐value <0.05 was considered statistically significant.Results Network analysis revealed a significant weekly increase in GE (5.83%/week, 95% CI 4.32–7.37), LE (5.43%/week, 95% CI 3.63–7.25), and presence of SW across GA. No significant hemisphere differences were found in GE (
P = 0.971) or LE (P = 0.458). Increasing GA was significantly associated with increasing NE in 41 nodes, increasing BC in 3 nodes, and decreasing BC in 2 nodes.Data Conclusion Extensive network development and refinement occur in the second and third trimesters, marked by a rapid increase in global integration and local segregation.
Level of Evidence 3
Technical Efficacy Stage 2
-
Yap, Pew-Thian (Ed.)Diffusion weighted imaging (DWI) with multiple, high b-values is critical for extracting tissue microstructure measurements; however, high b-value DWI images contain high noise levels that can overwhelm the signal of interest and bias microstructural measurements. Here, we propose a simple denoising method that can be applied to any dataset, provided a low-noise, single-subject dataset is acquired using the same DWI sequence. The denoising method uses a one-dimensional convolutional neural network (1D-CNN) and deep learning to learn from a low-noise dataset, voxel-by-voxel. The trained model can then be applied to high-noise datasets from other subjects. We validated the 1D-CNN denoising method by first demonstrating that 1D-CNN denoising resulted in DWI images that were more similar to the noise-free ground truth than comparable denoising methods, e.g., MP-PCA, using simulated DWI data. Using the same DWI acquisition but reconstructed with two common reconstruction methods, i.e. SENSE1 and sum-of-square, to generate a pair of low-noise and high-noise datasets, we then demonstrated that 1D-CNN denoising of high-noise DWI data collected from human subjects showed promising results in three domains: DWI images, diffusion metrics, and tractography. In particular, the denoised images were very similar to a low-noise reference image of that subject, more than the similarity between repeated low-noise images (i.e. computational reproducibility). Finally, we demonstrated the use of the 1D-CNN method in two practical examples to reduce noise from parallel imaging and simultaneous multi-slice acquisition. We conclude that the 1D-CNN denoising method is a simple, effective denoising method for DWI images that overcomes some of the limitations of current state-of-the-art denoising methods, such as the need for a large number of training subjects and the need to account for the rectified noise floor.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1038/s41597-021-01080-w
