MaPPeRTrac: A Massively Parallel, Portable, and Reproducible Tractography Pipeline

Large-scale diffusion MRI tractography remains a significant challenge. Users must orchestrate a complex sequence of instructions that requires many software packages with complex dependencies and high computational costs. We developed MaPPeRTrac, an edge-centric tractography pipeline that simplifies and accelerates this process in a wide range of high-performance computing (HPC) environments. It fully automates either probabilistic or deterministic tractography, starting from a subject's magnetic resonance imaging (MRI) data, including structural and diffusion MRI images, to the edge density image (EDI) of their structural connectomes. Dependencies are containerized with Singularity (now called Apptainer) and decoupled from code to enable rapid prototyping and modification. Data derivatives are organized with the Brain Imaging Data Structure (BIDS) to ensure that they are findable, accessible, interoperable, and reusable following FAIR principles. The pipeline takes full advantage of HPC resources using the Parsl parallel programming framework, resulting in the creation of connectome datasets of unprecedented size. MaPPeRTrac is publicly available and tested on commercial and scientific hardware, so it can accelerate brain connectome research for a broader user community. MaPPeRTrac is available at: https://github.com/LLNL/mappertrac .

Connectome mapping with edge density imaging differentiates pediatric mild traumatic brain injury from typically developing controls: proof of concept

BACKGROUND

Although acute neurologic impairment might be transient, other long-term effects can be observed with mild traumatic brain injury. However, when pediatric patients with mild traumatic brain injury present for medical care, conventional imaging with CT and MR imaging often does not reveal abnormalities.

OBJECTIVE

To determine whether edge density imaging can separate pediatric mild traumatic brain injury from typically developing controls.

MATERIALS AND METHODS

Subjects were recruited as part of the "Therapeutic Resources for Attention Improvement using Neuroimaging in Traumatic Brain Injury" (TRAIN-TBI) study. We included 24 adolescents (χ=14.1 years of age, σ=1.6 years, range 10-16 years), 14 with mild traumatic brain injury (TBI) and 10 typically developing controls. Neurocognitive assessments included the pediatric version of the California Verbal Learning Test (CVLT) and the Attention Network Task (ANT). Diffusion MR imaging was acquired on a 3-tesla (T) scanner. Edge density images were computed utilizing fiber tractography. Principal component analysis (PCA) and support vector machines (SVM) were used in an exploratory analysis to separate mild TBI and control groups. The diagnostic accuracy of edge density imaging, neurocognitive tests, and fractional anisotropy (FA) from diffusion tensor imaging (DTI) was computed with two-sample t-tests and receiver operating characteristic (ROC) metrics.

RESULTS

Support vector machine-principal component analysis of edge density imaging maps identified three white matter regions distinguishing pediatric mild TBI from controls. The bilateral tapetum, sagittal stratum, and callosal splenium identified mild TBI subjects with sensitivity of 79% and specificity of 100%. Accuracy from the area under the ROC curve (AUC) was 94%. Neurocognitive testing provided an AUC of 61% (CVLT) and 71% (ANT). Fractional anisotropy yielded an AUC of 48%.

CONCLUSION

In this proof-of-concept study, we show that edge density imaging is a new form of connectome mapping that provides better diagnostic delineation between pediatric mild TBI and healthy controls than DTI or neurocognitive assessments of memory or attention.

Diffusion tensor tractography in children with sensory processing disorder: Potentials for devising machine learning classifiers

The "sensory processing disorder" (SPD) refers to brain's inability to organize sensory input for appropriate use. In this study, we determined the diffusion tensor imaging (DTI) microstructural and connectivity correlates of SPD, and apply machine learning algorithms for identification of children with SPD based on DTI/tractography metrics. A total of 44 children with SPD and 41 typically developing children (TDC) were prospectively recruited and scanned. In addition to fractional anisotropy (FA), mean diffusivity (MD), and radial diffusivity (RD), we applied probabilistic tractography to generate edge density (ED) and track density (TD) from DTI maps. For identification of children with SPD, accurate classification rates from a combination of DTI microstructural (FA, MD, AD, and RD), connectivity (TD) and connectomic (ED) metrics with different machine learning algorithms - including naïve Bayes, random forest, support vector machine, and neural networks - were determined. In voxel-wise analysis, children with SPD had lower FA, ED, and TD but higher MD and RD compared to TDC - predominantly in posterior white matter tracts including posterior corona radiata, posterior thalamic radiation, and posterior body and splenium of corpus callosum. In stepwise penalized logistic regression, the only independent variable distinguishing children with SPD from TDC was the average TD in the splenium (p < 0.001). Among different combinations of machine learning algorithms and DTI/connectivity metrics, random forest models using tract-based TD yielded the highest accuracy in classification of SPD - 77.5% accuracy, 73.8% sensitivity, and 81.6% specificity. Our findings demonstrate impaired microstructural and connectivity/connectomic integrity in children with SPD, predominantly in posterior white matter tracts, and with reduced TD of the splenium of corpus callosum as the most distinctive pattern. Applying machine learning algorithms, these connectivity metrics can be used to devise novel imaging biomarkers for neurodevelopmental disorders.