Diffusion MRI (dMRI) streamline tractography is the gold-standard for in vivo estimation of white matter (WM) pathways in the brain. However, the high angular resolution dMRI acquisitions capable of fitting the microstructural models needed for tractography are often time-consuming and not routinely collected clinically, restricting the scope of tractography analyses. To address this limitation, we build on recent advances in deep learning which have demonstrated that streamline propagation can be learned from dMRI directly without traditional model fitting. Specifically, we propose learning the streamline propagator from T1w MRI to facilitate arbitrary tractography analyses when dMRI is unavailable. To do so, we present a novel convolutional-recurrent neural network (CoRNN) trained in a teacher-student framework that leverages T1w MRI, associated anatomical context, and streamline memory from data acquired for the Human Connectome Project. We characterize our approach under two common tractography paradigms, WM bundle analysis and structural connectomics, and find approximately a 5-15% difference between measures computed from streamlines generated with our approach and those generated using traditional dMRI tractography. When placed in the literature, these results suggest that the accuracy of WM measures computed from T1w MRI with our method is on the level of scan-rescan dMRI variability and raise an important question: is tractography truly a microstructural phenomenon, or has dMRI merely facilitated its discovery and implementation?

BackgroundWhile graph measures are used increasingly to characterize human connectomes, uncertainty remains in how to use these metrics in a quantitative and reproducible manner. Specifically, there is a lack of community consensus regarding the number of streamlines needed to generate connectomes.PurposeThe purpose was to define the relationship between streamline count and graph‐measure value, reproducibility, and repeatability.Study TypeRetrospective analysis of previously prospective study.PopulationTen healthy subjects, 70% female, aged 25.3 ± 5.9 years.Field Strength/SequenceA 3‐T, T1‐weighted sequences and diffusion‐weighted imaging (DWI) with two gradient strengths (b‐values = 1200 and 3000 sec/mm2, echo time [TE] = 68 msec, repetition time [TR] = 5.4 seconds, 120 slices, field of view = 188 mm2).AssessmentA total of 13 graph‐theory measures were derived for each subject by generating probabilistic whole‐brain tractography from DWI and mapping the structural connectivity to connectomes. The streamline count invariance from changes in mean, repeatability, and reproducibility were derived.Statistical TestsPaired t‐test with P value <0.05 was used to compare graph‐measure means with a reference, intraclass correlation coefficient (ICC) to measure repeatability, and concordance correlation coefficient (CCC) to measure reproducibility.ResultsModularity and global efficiency converged to their reference mean with ICC > 0.90 and CCC > 0.99. Edge count, small‐worldness, randomness, and average betweenness centrality converged to the reference mean, with ICC > 0.90 and CCC > 0.95. Assortativity and average participation coefficient converged with ICC > 0.75 and CCC > 0.90. Density, average node strength, average node degree, characteristic path length, average local efficiency, and average clustering coefficient did not converge, though had ICC > 0.90 and CCC > 0.99. For these measures, alternate definitions that converge a reference mean are provided.Data ConclusionModularity and global efficiency are streamline count invariant for greater than 6 million and 100,000 streamlines, respectively. Density, average node strength, average node degree, characteristic path length, average local efficiency, and average clustering coefficient were strongly dependent on streamline count.Evidence Level1.Technical EfficacyStage 1.

Complex graph theory measures of brain structural connectomes derived from diffusion weighted images (DWI) provide insight into the network structure of the brain. Further, as the number of available DWI datasets grows, so does the ability to investigate associations in these measures with major biological factors, like age. However, one key hurdle that remains is the presence of scanner effects that can arise from different DWI datasets and confound multisite analyses. Two common approaches to correct these effects are voxel-wise and feature-wise harmonization. However, it is still unclear how to best leverage them for graph-theory analysis of an aging population. Thus, there is a need to better characterize the impact of each harmonization method and their ability to preserve age related features. We investigate this by characterizing four complex graph theory measures (modularity, characteristic path length, global efficiency, and betweenness centrality) in 48 participants aged 55 to 86 from Baltimore Longitudinal Study of Aging (BLSA) and Vanderbilt Memory and Aging Project (VMAP) before and after voxel-and feature-wise harmonization with the Null Space Deep Network (NSDN) and ComBat, respectively. First, we characterize across dataset coefficients of variation (CoV) and find the combination of NSDN and ComBat causes the greatest reduction in CoV followed by ComBat alone then NSDN alone. Second, we reproduce published associations of modularity with age after correcting for other covariates with linear models. We find that harmonization with ComBat or ComBat and NSDN together improves the significance of existing age effects, reduces model residuals, and qualitatively reduces separation between datasets. These results reinforce the efficiency of statistical harmonization on the feature-level with ComBat and suggest that harmonization on the voxel-level is synergistic but may have reduced effect after running through the multiple layers of the connectomics pipeline. Thus, we conclude that featurewise harmonization improves statistical results, but the addition of biologically informed voxel-based harmonization offers further improvement.

Diffusion weighted magnetic resonance imaging (DW-MRI) captures tissue microarchitecture at a millimeter scale. With recent advantages in data sharing, large-scale multi-site DW-MRI datasets are being made available for multi-site studies. However, DW-MRI suffers from measurement variability (e.g., inter-and intra-site variability, hardware performance, and sequence design), which consequently yields inferior performance on multi-site and/or longitudinal diffusion studies. In this study, we propose a novel, deep learning-based method to harmonize DW-MRI signals for a more reproducible and robust estimation of microstructure. Our method introduces a data-driven scanner-invariant regularization scheme to model a more robust fiber orientation distribution function (FODF) estimation. We study the Human Connectome Project (HCP) young adults test-retest group as well as the MASiVar dataset (with inter-and intra-site scan/rescan data). The 8 th order spherical harmonics coefficients are employed as data representation. The results show that the proposed harmonization approach maintains higher angular correlation coefficients (ACC) with the ground truth signals (0.954 versus 0.942), while achieves higher consistency of FODF signals for intra-scanner data (0.891 versus 0.826), as compared with the baseline supervised deep learning scheme. Furthermore, the proposed data-driven framework is flexible and potentially applicable to a wider range of data harmonization problems in neuroimaging.