Yogesh Rathi

YNICL Journal 2026 Journal Article

Study of sex differences in the whole brain white matter using diffusion MRI tractography and suprathreshold fiber cluster statistics

• Fan Zhang
• Jarrett Rushmore
• Yijie Li
• Suheyla Cetin-Karayumak
• Yang Song
• Weidong Cai
• Carl-Fredrik Westin
• James J. Levitt

Sex-specific characteristics demonstrate a substantial influence on the brain white matter, suggesting distinct structural connectivity patterns between females and males. Diffusion MRI tractography is an important tool in assessing white matter connectivity and brain tissue microstructure across different populations. Whole brain tractography analysis for group statistical comparison is a challenging task due to the large number of white matter connections. This work studies whole-brain white matter connectivity differences between females and males using dMRI tractography. We study a large cohort of 707 subjects from the Human Connectome Project Young Adult dataset. By applying a well-established fiber clustering pipeline and a suprathreshold fiber cluster statistical method, we analyze tracts in the cerebral cortex and understudied pathways like those connecting to the cerebellum. We identify several tracts with significant sex differences in terms of their fractional anisotropy and/or mean diffusivity. These include deep tracts like the arcuate fasciculus, corticospinal tract, and corpus callosum, superficial tracts in the frontal lobe, and cerebellar tracts. Finally, a canonical correlation analysis (CCA) that identifies covariance patterns between white matter and behavior measures reveals that these white matter differences are associated with a range of neurobehavioral measures, with the strongest and most consistent associations observed for motor function, suggesting motor circuits as a potential key focus for future research.

YNIMG Journal 2025 Journal Article

DeepNuParc: A novel deep clustering framework for fine-scale parcellation of brain nuclei using diffusion MRI tractography

• Haolin He
• Ce Zhu
• Le Zhang
• Yipeng Liu
• Xiao Xu
• Yuqian Chen
• Leo Zekelman
• Jarrett Rushmore

Brain nuclei are clusters of anatomically distinct neurons that serve as important hubs for processing and relaying information in various neural circuits. Fine-scale parcellation of the brain nuclei is vital for a comprehensive understanding of their anatomico-functional correlations. Diffusion MRI tractography is an advanced imaging technique that can estimate the brain's white matter structural connectivity to potentially reveal the topography of the nuclei of interest for studying their subdivisions. In this work, we present a deep clustering pipeline, namely DeepNuParc, to perform automated, fine-scale parcellation of brain nuclei using diffusion MRI tractography. First, we incorporate a newly proposed deep learning approach to enable accurate segmentation of the nuclei of interest directly on the dMRI data. Next, we design a novel streamline clustering-based structural connectivity feature for a robust representation of voxels within the nuclei. Finally, we improve the popular joint dimensionality reduction and k-means clustering approach to enable nuclei parcellation at a finer scale. We demonstrate DeepNuParc on two important brain structures, i.e. the amygdala and the thalamus, that are known to have multiple anatomically and functionally distinct nucleus subdivisions. Experimental results show that DeepNuParc enables consistent parcellation of the nuclei into multiple parcels across multiple subjects and achieves good correspondence with the widely used coarse-scale atlases. Our code is available at https://github.com/HarlandZZC/deep_nuclei_parcellation.

YNIMG Journal 2023 Journal Article

Deep fiber clustering: Anatomically informed fiber clustering with self-supervised deep learning for fast and effective tractography parcellation

• Yuqian Chen
• Chaoyi Zhang
• Tengfei Xue
• Yang Song
• Nikos Makris
• Yogesh Rathi
• Weidong Cai
• Fan Zhang

White matter fiber clustering is an important strategy for white matter parcellation, which enables quantitative analysis of brain connections in health and disease. In combination with expert neuroanatomical labeling, data-driven white matter fiber clustering is a powerful tool for creating atlases that can model white matter anatomy across individuals. While widely used fiber clustering approaches have shown good performance using classical unsupervised machine learning techniques, recent advances in deep learning reveal a promising direction toward fast and effective fiber clustering. In this work, we propose a novel deep learning framework for white matter fiber clustering, Deep Fiber Clustering (DFC), which solves the unsupervised clustering problem as a self-supervised learning task with a domain-specific pretext task to predict pairwise fiber distances. This process learns a high-dimensional embedding feature representation for each fiber, regardless of the order of fiber points reconstructed during tractography. We design a novel network architecture that represents input fibers as point clouds and allows the incorporation of additional sources of input information from gray matter parcellation. Thus, DFC makes use of combined information about white matter fiber geometry and gray matter anatomy to improve the anatomical coherence of fiber clusters. In addition, DFC conducts outlier removal naturally by rejecting fibers with low cluster assignment probability. We evaluate DFC on three independently acquired cohorts, including data from 220 individuals across genders, ages (young and elderly adults), and different health conditions (healthy control and multiple neuropsychiatric disorders). We compare DFC to several state-of-the-art white matter fiber clustering algorithms. Experimental results demonstrate superior performance of DFC in terms of cluster compactness, generalization ability, anatomical coherence, and computational efficiency.

YNIMG Journal 2022 Journal Article

Cross-site harmonization of multi-shell diffusion MRI measures based on rotational invariant spherical harmonics (RISH)

• Alberto de Luca
• Suheyla Cetin Karayumak
• Alexander Leemans
• Yogesh Rathi
• Stephan Swinnen
• Jolien Gooijers
• Amanda Clauwaert
• Roald Bahr

Quantification methods based on the acquisition of diffusion magnetic resonance imaging (dMRI) with multiple diffusion weightings (e.g., multi-shell) are becoming increasingly applied to study the in-vivo brain. Compared to single-shell data for diffusion tensor imaging (DTI), multi-shell data allows to apply more complex models such as diffusion kurtosis imaging (DKI), which attempts to capture both diffusion hindrance and restriction effects, or biophysical models such as NODDI, which attempt to increase specificity by separating biophysical components. Because of the strong dependence of the dMRI signal on the measurement hardware, DKI and NODDI metrics show scanner and site differences, much like other dMRI metrics. These effects limit the implementation of multi-shell approaches in multicenter studies, which are needed to collect large sample sizes for robust analyses. Recently, a post-processing technique based on rotation invariant spherical harmonics (RISH) features was introduced to mitigate cross-scanner differences in DTI metrics. Unlike statistical harmonization methods, which require repeated application to every dMRI metric of choice, RISH harmonization is applied once on the raw data, and can be followed by any analysis. RISH features harmonization has been tested on DTI features but not its generalizability to harmonize multi-shell dMRI. In this work, we investigated whether performing the RISH features harmonization of multi-shell dMRI data removes cross-site differences in DKI and NODDI metrics while retaining longitudinal effects. To this end, 46 subjects underwent a longitudinal (up to 3 time points) two-shell dMRI protocol at 3 imaging sites. DKI and NODDI metrics were derived before and after harmonization and compared both at the whole brain level and at the voxel level. Then, the harmonization effects on cross-sectional and on longitudinal group differences were evaluated. RISH features averaged for each of the 3 sites exhibited prominent between-site differences in the frontal and posterior part of the brain. Statistically significant differences in fractional anisotropy, mean diffusivity and mean kurtosis were observed both at the whole brain and voxel level between all the acquisition sites before harmonization, but not after. The RISH method also proved effective to harmonize NODDI metrics, particularly in white matter. The RISH based harmonization maintained the magnitude and variance of longitudinal changes as compared to the non-harmonized data of all considered metrics. In conclusion, the application of RISH feature based harmonization to multi-shell dMRI data can be used to remove cross-site differences in DKI metrics and NODDI analyses, while retaining inherent relations between longitudinal acquisitions.

YNIMG Journal 2022 Journal Article

White matter association tracts underlying language and theory of mind: An investigation of 809 brains from the Human Connectome Project

• Leo R. Zekelman
• Fan Zhang
• Nikos Makris
• Jianzhong He
• Yuqian Chen
• Tengfei Xue
• Daniela Liera
• Daniel L. Drane

Language and theory of mind (ToM) are the cognitive capacities that allow for the successful interpretation and expression of meaning. While functional MRI investigations are able to consistently localize language and ToM to specific cortical regions, diffusion MRI investigations point to an inconsistent and sometimes overlapping set of white matter tracts associated with these two cognitive domains. To further examine the white matter tracts that may underlie these domains, we use a two-tensor tractography method to investigate the white matter microstructure of 809 participants from the Human Connectome Project. 20 association white matter tracts (10 in each hemisphere) are uniquely identified by leveraging a neuroanatomist-curated automated white matter tract atlas. The fractional anisotropy (FA), mean diffusivity (MD), and number of streamlines (NoS) are measured for each white matter tract. Performance on neuropsychological assessments of semantic memory (NIH Toolbox Picture Vocabulary Test, TPVT) and emotion perception (Penn Emotion Recognition Test, PERT) are used to measure critical subcomponents of the language and ToM networks, respectively. Regression models are constructed to examine how structural measurements of left and right white matter tracts influence performance across these two assessments. We find that semantic memory performance is influenced by the number of streamlines of the left superior longitudinal fasciculus III (SLF-III), and emotion perception performance is influenced by the number of streamlines of the right SLF-III. Additionally, we find that performance on both semantic memory & emotion perception is influenced by the FA of the left arcuate fasciculus (AF). The results point to multiple, overlapping white matter tracts that underlie the cognitive domains of language and ToM. Results are discussed in terms of hemispheric dominance and concordance with prior investigations.

YNIMG Journal 2021 Journal Article

Tractography dissection variability: What happens when 42 groups dissect 14 white matter bundles on the same dataset?

• Kurt G. Schilling
• François Rheault
• Laurent Petit
• Colin B. Hansen
• Vishwesh Nath
• Fang-Cheng Yeh
• Gabriel Girard
• Muhamed Barakovic

White matter bundle segmentation using diffusion MRI fiber tractography has become the method of choice to identify white matter fiber pathways in vivo in human brains. However, like other analyses of complex data, there is considerable variability in segmentation protocols and techniques. This can result in different reconstructions of the same intended white matter pathways, which directly affects tractography results, quantification, and interpretation. In this study, we aim to evaluate and quantify the variability that arises from different protocols for bundle segmentation. Through an open call to users of fiber tractography, including anatomists, clinicians, and algorithm developers, 42 independent teams were given processed sets of human whole-brain streamlines and asked to segment 14 white matter fascicles on six subjects. In total, we received 57 different bundle segmentation protocols, which enabled detailed volume-based and streamline-based analyses of agreement and disagreement among protocols for each fiber pathway. Results show that even when given the exact same sets of underlying streamlines, the variability across protocols for bundle segmentation is greater than all other sources of variability in the virtual dissection process, including variability within protocols and variability across subjects. In order to foster the use of tractography bundle dissection in routine clinical settings, and as a fundamental analytical tool, future endeavors must aim to resolve and reduce this heterogeneity. Although external validation is needed to verify the anatomical accuracy of bundle dissections, reducing heterogeneity is a step towards reproducible research and may be achieved through the use of standard nomenclature and definitions of white matter bundles and well-chosen constraints and decisions in the dissection process.

YNICL Journal 2020 Journal Article

Anatomical assessment of trigeminal nerve tractography using diffusion MRI: A comparison of acquisition b-values and single- and multi-fiber tracking strategies

• Guoqiang Xie
• Fan Zhang
• Laura Leung
• Michael A. Mooney
• Lorenz Epprecht
• Isaiah Norton
• Yogesh Rathi
• Ron Kikinis

BACKGROUND: and conventional diffusion tensor MRI (DTI) tractography methods. Though higher b-value acquisitions and multi-tensor tractography methods are known to be beneficial for tracking brain white matter fiber tracts, there have been no studies conducted to evaluate the performance of these advanced approaches on nerve tracking of the TGN, in particular on tracking different anatomical regions of the TGN. OBJECTIVE: We compare TGN tracking performance using dMRI data with different b-values, in combination with both single- and multi-tensor tractography methods. Our goal is to assess the advantages and limitations of these different strategies for identifying the anatomical regions of the TGN. METHODS: , using single-tensor (1T) and two-tensor (2T) unscented Kalman filter (UKF) tractography. This resulted in a total of six tracking strategies. The TGN was identified using an anatomical region-of-interest (ROI) selection approach. First, in a subset of the dataset we identified ROIs that provided good TGN tracking performance across all tracking strategies. Using these ROIs, the TGN was then tracked in all subjects using the six tracking strategies. An expert rater (GX) visually assessed and scored each TGN based on seven anatomical judgment criteria. These criteria included the presence of multiple expected anatomical segments of the TGN (true positive structures), specifically branch-like structures, cisternal portion, mesencephalic trigeminal tract, and spinal cord tract of the TGN. False positive criteria included the presence of any fibers entering the temporal lobe, the inferior cerebellar peduncle, or the middle cerebellar peduncle. Expert rating scores were analyzed to compare TGN tracking performance across the six tracking strategies. Intra- and inter-rater validation was performed to assess the reliability of the expert TGN rating result. RESULTS: The TGN was selected using two anatomical ROIs (Meckel's Cave and cisternal portion of the TGN). The two-tensor tractography method had significantly better performance on identifying true positive structures, while generating more false positive streamlines in comparison to the single-tensor tractography method. TGN tracking performance was significantly different across the three b-values for almost all structures studied. Tracking performance was reported in terms of the percentage of subjects achieving each anatomical rating criterion. Tracking of the cisternal portion and branching structure of the TGN was generally successful, with the highest performance of over 98% using two-tensor tractography and b = 1000 or b = 2000. However, tracking the smaller mesencephalic and spinal cord tracts of the TGN was quite challenging (highest performance of 37.5% and 57.07%, using two-tensor tractography with b = 1000 and b = 2000, respectively). False positive connections to the temporal lobe (over 38% of subjects for all strategies) and cerebellar peduncles (100% of subjects for all strategies) were prevalent. High joint probability of agreement was obtained in the inter-rater (on average 83%) and intra-rater validation (on average 90%), showing a highly reliable expert rating result. CONCLUSIONS: Overall, the results of the study suggest that researchers and clinicians may benefit from tailoring their acquisition and tracking methodology to the specific anatomical portion of the TGN that is of the greatest interest. For example, tracking of branching structures and TGN-T2 overlap can be best achieved with a two-tensor model and an acquisition using b = 1000 or b = 2000. In general, b = 1000 and b = 2000 acquisitions provided the best-rated tracking results. Further research is needed to improve both sensitivity and specificity of the depiction of the TGN anatomy using dMRI.

YNIMG Journal 2020 Journal Article

Creation of a novel trigeminal tractography atlas for automated trigeminal nerve identification

• Fan Zhang
• Guoqiang Xie
• Laura Leung
• Michael A. Mooney
• Lorenz Epprecht
• Isaiah Norton
• Yogesh Rathi
• Ron Kikinis

Diffusion MRI (dMRI) tractography has been successfully used to study the trigeminal nerves (TGNs) in many clinical and research applications. Currently, identification of the TGN in tractography data requires expert nerve selection using manually drawn regions of interest (ROIs), which is prone to inter-observer variability, time-consuming and carries high clinical and labor costs. To overcome these issues, we propose to create a novel anatomically curated TGN tractography atlas that enables automated identification of the TGN from dMRI tractography. In this paper, we first illustrate the creation of a trigeminal tractography atlas. Leveraging a well-established computational pipeline and expert neuroanatomical knowledge, we generate a data-driven TGN fiber clustering atlas using tractography data from 50 subjects from the Human Connectome Project. Then, we demonstrate the application of the proposed atlas for automated TGN identification in new subjects, without relying on expert ROI placement. Quantitative and visual experiments are performed with comparison to expert TGN identification using dMRI data from two different acquisition sites. We show highly comparable results between the automatically and manually identified TGNs in terms of spatial overlap and visualization, while our proposed method has several advantages. First, our method performs automated TGN identification, and thus it provides an efficient tool to reduce expert labor costs and inter-operator bias relative to expert manual selection. Second, our method is robust to potential imaging artifacts and/or noise that can prevent successful manual ROI placement for TGN selection and hence yields a higher successful TGN identification rate.

YNIMG Journal 2020 Journal Article

Cross-scanner and cross-protocol multi-shell diffusion MRI data harmonization: Algorithms and results

• Lipeng Ning
• Elisenda Bonet-Carne
• Francesco Grussu
• Farshid Sepehrband
• Enrico Kaden
• Jelle Veraart
• Stefano B. Blumberg
• Can Son Khoo

Cross-scanner and cross-protocol variability of diffusion magnetic resonance imaging (dMRI) data are known to be major obstacles in multi-site clinical studies since they limit the ability to aggregate dMRI data and derived measures. Computational algorithms that harmonize the data and minimize such variability are critical to reliably combine datasets acquired from different scanners and/or protocols, thus improving the statistical power and sensitivity of multi-site studies. Different computational approaches have been proposed to harmonize diffusion MRI data or remove scanner-specific differences. To date, these methods have mostly been developed for or evaluated on single b-value diffusion MRI data. In this work, we present the evaluation results of 19 algorithms that are developed to harmonize the cross-scanner and cross-protocol variability of multi-shell diffusion MRI using a benchmark database. The proposed algorithms rely on various signal representation approaches and computational tools, such as rotational invariant spherical harmonics, deep neural networks and hybrid biophysical and statistical approaches. The benchmark database consists of data acquired from the same subjects on two scanners with different maximum gradient strength (80 and 300 ​mT/m) and with two protocols. We evaluated the performance of these algorithms for mapping multi-shell diffusion MRI data across scanners and across protocols using several state-of-the-art imaging measures. The results show that data harmonization algorithms can reduce the cross-scanner and cross-protocol variabilities to a similar level as scan-rescan variability using the same scanner and protocol. In particular, the LinearRISH algorithm based on adaptive linear mapping of rotational invariant spherical harmonics features yields the lowest variability for our data in predicting the fractional anisotropy (FA), mean diffusivity (MD), mean kurtosis (MK) and the rotationally invariant spherical harmonic (RISH) features. But other algorithms, such as DIAMOND, SHResNet, DIQT, CMResNet show further improvement in harmonizing the return-to-origin probability (RTOP). The performance of different approaches provides useful guidelines on data harmonization in future multi-site studies.

YNICL Journal 2020 Journal Article

Mild traumatic brain injury impacts associations between limbic system microstructure and post-traumatic stress disorder symptomatology

• Valerie J. Sydnor
• Sylvain Bouix
• Ofer Pasternak
• Elisabeth Hartl
• Laura Levin-Gleba
• Benjamin Reid
• Yorghos Tripodis
• Jeffrey P. Guenette

BACKGROUND: Post-traumatic stress disorder (PTSD) is a psychiatric disorder that afflicts many individuals, yet the neuropathological mechanisms that contribute to this disorder remain to be fully determined. Moreover, it is unclear how exposure to mild traumatic brain injury (mTBI), a condition that is often comorbid with PTSD, particularly among military personnel, affects the clinical and neurological presentation of PTSD. To address these issues, the present study explores relationships between PTSD symptom severity and the microstructure of limbic and paralimbic gray matter brain regions, as well as the impact of mTBI comorbidity on these relationships. METHODS: Structural and diffusion MRI data were acquired from 102 male veterans who were diagnosed with current PTSD. Diffusion data were analyzed with free-water imaging to quantify average CSF-corrected fractional anisotropy (FA) and mean diffusivity (MD) in 18 limbic and paralimbic gray matter regions. Associations between PTSD symptom severity and regional average dMRI measures were examined with repeated measures linear mixed models. Associations were studied separately in veterans with PTSD only, and in veterans with PTSD and a history of military mTBI. RESULTS: Analyses revealed that in the PTSD only cohort, more severe symptoms were associated with higher FA in the right amygdala-hippocampus complex, lower FA in the right cingulate cortex, and lower MD in the left medial orbitofrontal cortex. In the PTSD and mTBI cohort, more severe PTSD symptoms were associated with higher FA bilaterally in the amygdala-hippocampus complex, with higher FA bilaterally in the nucleus accumbens, with lower FA bilaterally in the cingulate cortex, and with higher MD in the right amygdala-hippocampus complex. CONCLUSIONS: These findings suggest that the microstructure of limbic and paralimbic brain regions may influence PTSD symptomatology. Further, given the additional associations observed between microstructure and symptom severity in veterans with head trauma, we speculate that mTBI may exacerbate the impact of brain microstructure on PTSD symptoms, especially within regions of the brain known to be vulnerable to chronic stress. A heightened sensitivity to the microstructural environment of the brain could partially explain why individuals with PTSD and mTBI comorbidity experience more severe symptoms and poorer illness prognoses than those without a history of brain injury. The relevance of these microstructural findings to the conceptualization of PTSD as being a disorder of stress-induced neuronal connectivity loss is discussed.

YNICL Journal 2019 Journal Article

Abnormalities in gray matter microstructure in young adults with 22q11.2 deletion syndrome

• Zora Kikinis
• Nikos Makris
• Valerie J. Sydnor
• Sylvain Bouix
• Ofer Pasternak
• Ioana L. Coman
• Kevin M. Antshel
• Wanda Fremont

BACKGROUND: 22q11.2 Deletion Syndrome (22q11DS) is a genetic, neurodevelopmental disorder characterized by a chromosomal deletion and a distinct cognitive profile. Although abnormalities in the macrostructure of the cortex have been identified in individuals with 22q11DS, it is not known if there are additional microstructural changes in gray matter regions in this syndrome, and/or if such microstructural changes are associated with cognitive functioning. METHODS: This study employed a novel diffusion MRI measure, the Heterogeneity of Fractional Anisotropy (HFA), to examine variability in the microstructural organization of the cortex in healthy young adults (N = 30) and those with 22q11DS (N = 56). Diffusion MRI, structural MRI, clinical and cognitive data were acquired. RESULTS: Compared to controls, individuals with 22q11DS evinced increased HFA in cortical association (p = .003, d = 0.86) and paralimbic (p < .0001, d = 1.2) brain areas, whereas no significant differences were found between the two groups in primary cortical brain areas. Additionally, increased HFA of the right paralimbic area was associated with poorer performance on tests of response inhibition, i.e., the Stroop Test (rho = -0.37 p = .005) and the Gordon Diagnostic System Vigilance Commission (rho = -0.41 p = .002) in the 22q11DS group. No significant correlations were found between HFA and cognitive abilities in the healthy control group. CONCLUSIONS: These findings suggest that cortical microstructural disorganization may be a neural correlate of response inhibition in individuals with 22q11DS. Given that the migration pattern of neural crest cells is disrupted at the time of early brain development in 22q11DS, we hypothesize that these neural alterations may be neurodevelopmental in origin, and reflect cortical dysfunction associated with cognitive deficits.

YNIMG Journal 2019 Journal Article

Limits to anatomical accuracy of diffusion tractography using modern approaches

• Kurt G. Schilling
• Vishwesh Nath
• Colin Hansen
• Prasanna Parvathaneni
• Justin Blaber
• Yurui Gao
• Peter Neher
• Dogu Baran Aydogan

Diffusion MRI fiber tractography is widely used to probe the structural connectivity of the brain, with a range of applications in both clinical and basic neuroscience. Despite widespread use, tractography has well-known pitfalls that limits the anatomical accuracy of this technique. Numerous modern methods have been developed to address these shortcomings through advances in acquisition, modeling, and computation. To test whether these advances improve tractography accuracy, we organized the 3-D Validation of Tractography with Experimental MRI (3D-VoTEM) challenge at the ISBI 2018 conference. We made available three unique independent tractography validation datasets – a physical phantom and two ex vivo brain specimens - resulting in 176 distinct submissions from 9 research groups. By comparing results over a wide range of fiber complexities and algorithmic strategies, this challenge provides a more comprehensive assessment of tractography's inherent limitations than has been reported previously. The central results were consistent across all sub-challenges in that, despite advances in tractography methods, the anatomical accuracy of tractography has not dramatically improved in recent years. Taken together, our results independently confirm findings from decades of tractography validation studies, demonstrate inherent limitations in reconstructing white matter pathways using diffusion MRI data alone, and highlight the need for alternative or combinatorial strategies to accurately map the fiber pathways of the brain.

YNIMG Journal 2019 Journal Article

Retrospective harmonization of multi-site diffusion MRI data acquired with different acquisition parameters

• Suheyla Cetin Karayumak
• Sylvain Bouix
• Lipeng Ning
• Anthony James
• Tim Crow
• Martha Shenton
• Marek Kubicki
• Yogesh Rathi

A joint and integrated analysis of multi-site diffusion MRI (dMRI) datasets can dramatically increase the statistical power of neuroimaging studies and enable comparative studies pertaining to several brain disorders. However, dMRI data sets acquired on multiple scanners cannot be naively pooled for joint analysis due to scanner specific nonlinear effects as well as differences in acquisition parameters. Consequently, for joint analysis, the dMRI data has to be harmonized, which involves removing scanner-specific differences from the raw dMRI signal. In this work, we propose a dMRI harmonization method that is capable of removing scanner-specific effects, while accounting for minor differences in acquisition parameters such as b-value, spatial resolution and number of gradient directions. We validate our algorithm on dMRI data acquired from two sites: Philadelphia Neurodevelopmental Cohort (PNC) with 800 healthy adolescents (ages 8–22 years) and Brigham and Women's Hospital (BWH) with 70 healthy subjects (ages 14–54 years). In particular, we show that gender and age-related maturation differences in different age groups are preserved after harmonization, as measured using effect sizes (small, medium and large), irrespective of the test sample size. Since we use matched control subjects from different scanners to estimate scanner-specific effects, our goal in this work is also to determine the minimum number of well-matched subjects needed from each site to achieve best harmonization results. Our results indicate that at-least 16 to 18 well-matched healthy controls from each site are needed to reliably capture scanner related differences. The proposed method can thus be used for retrospective harmonization of raw dMRI data across sites despite differences in acquisition parameters, while preserving inter-subject anatomical variability.

YNIMG Journal 2018 Journal Article

An anatomically curated fiber clustering white matter atlas for consistent white matter tract parcellation across the lifespan

• Fan Zhang
• Ye Wu
• Isaiah Norton
• Laura Rigolo
• Yogesh Rathi
• Nikos Makris
• Lauren J. O'Donnell

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

YNICL Journal 2018 Journal Article

Diagnostic value of structural and diffusion imaging measures in schizophrenia

• Jungsun Lee
• Myong-Wuk Chon
• Harin Kim
• Yogesh Rathi
• Sylvain Bouix
• Martha E. Shenton
• Marek Kubicki

Objectives: Many studies have attempted to discriminate patients with schizophrenia from healthy controls by machine learning using structural or functional MRI. We included both structural and diffusion MRI (dMRI) and performed random forest (RF) and support vector machine (SVM) in this study. Methods: We evaluated the performance of classifying schizophrenia using RF method and SVM with 504 features (volume and/or fractional anisotropy and trace) from 184 brain regions. We enrolled 47 patients and 23 age- and sex-matched healthy controls and resampled our data into a balanced dataset using a Synthetic Minority Oversampling Technique method. We randomly permuted the classification of all participants as a patient or healthy control 100 times and ran the RF and SVM with leave one out cross validation for each permutation. We then compared the sensitivity and specificity of the original dataset and the permuted dataset. Results: Classification using RF with 504 features showed a significantly higher rate of performance compared to classification by chance: sensitivity (87.6% vs. 47.0%) and specificity (95.9 vs. 48.4%) performed by RF, sensitivity (89.5% vs. 48.0%) and specificity (94.5% vs. 47.1%) performed by SVM. Conclusions: Machine learning using RF and SVM with both volume and diffusion measures can discriminate patients with schizophrenia with a high degree of performance. Further replications are required.

YNIMG Journal 2018 Journal Article

Investigation into local white matter abnormality in emotional processing and sensorimotor areas using an automatically annotated fiber clustering in major depressive disorder

• Ye Wu
• Fan Zhang
• Nikos Makris
• Yuping Ning
• Isaiah Norton
• Shenglin She
• Hongjun Peng
• Yogesh Rathi

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

YNIMG Journal 2018 Journal Article

Suprathreshold fiber cluster statistics: Leveraging white matter geometry to enhance tractography statistical analysis

• Fan Zhang
• Weining Wu
• Lipeng Ning
• Gloria McAnulty
• Deborah Waber
• Borjan Gagoski
• Kiera Sarill
• Hesham M. Hamoda

This work presents a suprathreshold fiber cluster (STFC) method that leverages the whole brain fiber geometry to enhance statistical group difference analyses. The proposed method consists of 1) a well-established study-specific data-driven tractography parcellation to obtain white matter tract parcels and 2) a newly proposed nonparametric, permutation-test-based STFC method to identify significant differences between study populations. The basic idea of our method is that a white matter parcel's neighborhood (nearby parcels with similar white matter anatomy) can support the parcel's statistical significance when correcting for multiple comparisons. We propose an adaptive parcel neighborhood strategy to allow suprathreshold fiber cluster formation that is robust to anatomically varying inter-parcel distances. The method is demonstrated by application to a multi-shell diffusion MRI dataset from 59 individuals, including 30 attention deficit hyperactivity disorder patients and 29 healthy controls. Evaluations are conducted using both synthetic and in-vivo data. The results indicate that the STFC method gives greater sensitivity in finding group differences in white matter tract parcels compared to several traditional multiple comparison correction methods.

YNIMG Journal 2018 Journal Article

Whole brain white matter connectivity analysis using machine learning: An application to autism

• Fan Zhang
• Peter Savadjiev
• Weidong Cai
• Yang Song
• Yogesh Rathi
• Birkan Tunç
• Drew Parker
• Tina Kapur

In this paper, we propose an automated white matter connectivity analysis method for machine learning classification and characterization of white matter abnormality via identification of discriminative fiber tracts. The proposed method uses diffusion MRI tractography and a data-driven approach to find fiber clusters corresponding to subdivisions of the white matter anatomy. Features extracted from each fiber cluster describe its diffusion properties and are used for machine learning. The method is demonstrated by application to a pediatric neuroimaging dataset from 149 individuals, including 70 children with autism spectrum disorder (ASD) and 79 typically developing controls (TDC). A classification accuracy of 78. 33% is achieved in this cross-validation study. We investigate the discriminative diffusion features based on a two-tensor fiber tracking model. We observe that the mean fractional anisotropy from the second tensor (associated with crossing fibers) is most affected in ASD. We also find that local along-tract (central cores and endpoint regions) differences between ASD and TDC are helpful in differentiating the two groups. These altered diffusion properties in ASD are associated with multiple robustly discriminative fiber clusters, which belong to several major white matter tracts including the corpus callosum, arcuate fasciculus, uncinate fasciculus and aslant tract; and the white matter structures related to the cerebellum, brain stem, and ventral diencephalon. These discriminative fiber clusters, a small part of the whole brain tractography, represent the white matter connections that could be most affected in ASD. Our results indicate the potential of a machine learning pipeline based on white matter fiber clustering.

YNICL Journal 2017 Journal Article

Automated white matter fiber tract identification in patients with brain tumors

• Lauren J. O’Donnell
• Yannick Suter
• Laura Rigolo
• Pegah Kahali
• Fan Zhang
• Isaiah Norton
• Angela Albi
• Olutayo Olubiyi

We propose a method for the automated identification of key white matter fiber tracts for neurosurgical planning, and we apply the method in a retrospective study of 18 consecutive neurosurgical patients with brain tumors. Our method is designed to be relatively robust to challenges in neurosurgical tractography, which include peritumoral edema, displacement, and mass effect caused by mass lesions. The proposed method has two parts. First, we learn a data-driven white matter parcellation or fiber cluster atlas using groupwise registration and spectral clustering of multi-fiber tractography from healthy controls. Key fiber tract clusters are identified in the atlas. Next, patient-specific fiber tracts are automatically identified using tractography-based registration to the atlas and spectral embedding of patient tractography. Results indicate good generalization of the data-driven atlas to patients: 80% of the 800 fiber clusters were identified in all 18 patients, and 94% of the 800 fiber clusters were found in 16 or more of the 18 patients. Automated subject-specific tract identification was evaluated by quantitative comparison to subject-specific motor and language functional MRI, focusing on the arcuate fasciculus (language) and corticospinal tracts (motor), which were identified in all patients. Results indicate good colocalization: 89 of 95, or 94%, of patient-specific language and motor activations were intersected by the corresponding identified tract. All patient-specific activations were within 3mm of the corresponding language or motor tract. Overall, our results indicate the potential of an automated method for identifying fiber tracts of interest for neurosurgical planning, even in patients with mass lesions.

YNICL Journal 2017 Journal Article

Performance of unscented Kalman filter tractography in edema: Analysis of the two-tensor model

• Ruizhi Liao
• Lipeng Ning
• Zhenrui Chen
• Laura Rigolo
• Shun Gong
• Ofer Pasternak
• Alexandra J. Golby
• Yogesh Rathi

Diffusion MRI tractography is increasingly used in pre-operative neurosurgical planning to visualize critical fiber tracts. However, a major challenge for conventional tractography, especially in patients with brain tumors, is tracing fiber tracts that are affected by vasogenic edema, which increases water content in the tissue and lowers diffusion anisotropy. One strategy for improving fiber tracking is to use a tractography method that is more sensitive than the traditional single-tensor streamline tractography. We performed experiments to assess the performance of two-tensor unscented Kalman filter (UKF) tractography in edema. UKF tractography fits a diffusion model to the data during fiber tracking, taking advantage of prior information from the previous step along the fiber. We studied UKF performance in a synthetic diffusion MRI digital phantom with simulated edema and in retrospective data from two neurosurgical patients with edema affecting the arcuate fasciculus and corticospinal tracts. We compared the performance of several tractography methods including traditional streamline, UKF single-tensor, and UKF two-tensor. To provide practical guidance on how the UKF method could be employed, we evaluated the impact of using various seed regions both inside and outside the edematous regions, as well as the impact of parameter settings on the tractography sensitivity. We quantified the sensitivity of different methods by measuring the percentage of the patient-specific fMRI activation that was reached by the tractography. We expected that diffusion anisotropy threshold parameters, as well as the inclusion of a free water model, would significantly influence the reconstruction of edematous WM fiber tracts, because edema increases water content in the tissue and lowers anisotropy. Contrary to our initial expectations, varying the fractional anisotropy threshold and including a free water model did not affect the UKF two-tensor tractography output appreciably in these two patient datasets. The most effective parameter for increasing tracking sensitivity was the generalized anisotropy (GA) threshold, which increased the length of tracked fibers when reduced to 0.075. In addition, the most effective seeding strategy was seeding in the whole brain or in a large region outside of the edema. Overall, the main contribution of this study is to provide insight into how UKF tractography can work, using a two-tensor model, to begin to address the challenge of fiber tract reconstruction in edematous regions near brain tumors.

YNIMG Journal 2017 Journal Article

Precise Inference and Characterization of Structural Organization (PICASO) of tissue from molecular diffusion

• Lipeng Ning
• Evren Özarslan
• Carl-Fredrik Westin
• Yogesh Rathi

Inferring the microstructure of complex media from the diffusive motion of molecules is a challenging problem in diffusion physics. In this paper, we introduce a novel representation of diffusion MRI (dMRI) signal from tissue with spatially-varying diffusivity using a diffusion disturbance function. This disturbance function contains information about the (intra-voxel) spatial fluctuations in diffusivity due to restrictions, hindrances and tissue heterogeneity of the underlying tissue substrate. We derive the short- and long-range disturbance coefficients from this disturbance function to characterize the tissue structure and organization. Moreover, we provide an exact relation between the disturbance coefficients and the time-varying moments of the diffusion propagator, as well as their relation to specific tissue microstructural information such as the intra-axonal volume fraction and the apparent axon radius. The proposed approach is quite general and can model dMRI signal for any type of gradient sequence (rectangular, oscillating, etc.) without using the Gaussian phase approximation. The relevance of the proposed PICASO model is explored using Monte-Carlo simulations and in-vivo dMRI data. The results show that the estimated disturbance coefficients can distinguish different types of microstructural organization of axons.

YNIMG Journal 2016 Journal Article

A joint compressed-sensing and super-resolution approach for very high-resolution diffusion imaging

• Lipeng Ning
• Kawin Setsompop
• Oleg Michailovich
• Nikos Makris
• Martha E. Shenton
• Carl-Fredrik Westin
• Yogesh Rathi

Diffusion MRI (dMRI) can provide invaluable information about the structure of different tissue types in the brain. Standard dMRI acquisitions facilitate a proper analysis (e. g. tracing) of medium-to-large white matter bundles. However, smaller fiber bundles connecting very small cortical or sub-cortical regions cannot be traced accurately in images with large voxel sizes. Yet, the ability to trace such fiber bundles is critical for several applications such as deep brain stimulation and neurosurgery. In this work, we propose a novel acquisition and reconstruction scheme for obtaining high spatial resolution dMRI images using multiple low resolution (LR) images, which is effective in reducing acquisition time while improving the signal-to-noise ratio (SNR). The proposed method called compressed-sensing super resolution reconstruction (CS-SRR), uses multiple overlapping thick-slice dMRI volumes that are under-sampled in q-space to reconstruct diffusion signal with complex orientations. The proposed method combines the twin concepts of compressed sensing and super-resolution to model the diffusion signal (at a given b-value) in a basis of spherical ridgelets with total-variation (TV) regularization to account for signal correlation in neighboring voxels. A computationally efficient algorithm based on the alternating direction method of multipliers (ADMM) is introduced for solving the CS-SRR problem. The performance of the proposed method is quantitatively evaluated on several in-vivo human data sets including a true SRR scenario. Our experimental results demonstrate that the proposed method can be used for reconstructing sub-millimeter super resolution dMRI data with very good data fidelity in clinically feasible acquisition time.

YNIMG Journal 2016 Journal Article

Corrigendum to “A joint compressed-sensing and super-resolution approach for very high-resolution diffusion imaging”

• Lipeng Ning
• Kawin Setsompop
• Oleg Michailovich
• Nikos Makris
• Martha E. Shenton
• Carl-Fredrik Westin
• Yogesh Rathi

YNICL Journal 2015 Journal Article

Reconstruction of the arcuate fasciculus for surgical planning in the setting of peritumoral edema using two-tensor unscented Kalman filter tractography

• Zhenrui Chen
• Yanmei Tie
• Olutayo Olubiyi
• Laura Rigolo
• Alireza Mehrtash
• Isaiah Norton
• Ofer Pasternak
• Yogesh Rathi

BACKGROUND: Diffusion imaging tractography is increasingly used to trace critical fiber tracts in brain tumor patients to reduce the risk of post-operative neurological deficit. However, the effects of peritumoral edema pose a challenge to conventional tractography using the standard diffusion tensor model. The aim of this study was to present a novel technique using a two-tensor unscented Kalman filter (UKF) algorithm to track the arcuate fasciculus (AF) in brain tumor patients with peritumoral edema. METHODS: Ten right-handed patients with left-sided brain tumors in the vicinity of language-related cortex and evidence of significant peritumoral edema were retrospectively selected for the study. All patients underwent 3-Tesla magnetic resonance imaging (MRI) including a diffusion-weighted dataset with 31 directions. Fiber tractography was performed using both single-tensor streamline and two-tensor UKF tractography. A two-regions-of-interest approach was applied to perform the delineation of the AF. Results from the two different tractography algorithms were compared visually and quantitatively. RESULTS: Using single-tensor streamline tractography, the AF appeared disrupted in four patients and contained few fibers in the remaining six patients. Two-tensor UKF tractography delineated an AF that traversed edematous brain areas in all patients. The volume of the AF was significantly larger on two-tensor UKF than on single-tensor streamline tractography (p < 0.01). CONCLUSIONS: Two-tensor UKF tractography provides the ability to trace a larger volume AF than single-tensor streamline tractography in the setting of peritumoral edema in brain tumor patients.

AIIM Journal 2015 Journal Article

Sparse deconvolution of higher order tensor for fiber orientation distribution estimation

• Yuanjing Feng
• Ye Wu
• Yogesh Rathi
• Carl-Fredrik Westin

Purpose Higher order tensor (HOT) imaging approaches based on the spherical deconvolution framework have attracted much interest for their effectiveness in estimating fiber orientation distribution (FOD). However, sparse regularization techniques are still needed to obtain stable FOD in solving the deconvolution problem, particularly in very high orders. Our goal is to adequately characterize the actual sparsity lying in the FOD domain to develop accurate estimation approach for fiber orientation in HOT framework. Materials and methods We propose a sparse HOT regularization model by enforcing the sparse constraint directly on the representation of FOD instead of imposing it on coefficients of basis function. Then, we incorporate both the stabilizing effect of the l 2 penalty and the sparsity encouraging effect of the l 1 penalty in the sparse model to adequately characterize the actual sparsity lying in the FOD domain. Furthermore, a weighted regularization scheme is developed to iteratively solve the deconvolution problem. The deconvolution technique is compared against existing methods using l 2 or l 1 regularizer and tested on synthetic data and real human brain. Results Experiments were conducted on synthetic data and real human brain data. The synthetic experimental results indicate that crossing fibers are more easily detected and the angular resolution limit is improved by our method by approximately 20°–30° compared to existing HOT method. The detection accuracy is considerably improved compared with that of spherical deconvolution approaches using the l 2 regularizer and the reweighted l 1 scheme. Conclusions Results of testing the deconvolution technique demonstrate that it allows HOTs to obtain increasingly clean and sharp FOD, which in turn significantly increases the angular resolution of current HOT methods. With sparsity on FOD domain, this method efficiently improves the ability of HOT in resolving crossing fibers.