A registration method for improving quantitative assessment in probabilistic diffusion tractography

Diffusion MRI-based probabilistic tractography is a powerful tool for non-invasively investigating normal brain architecture and alterations in structural connectivity associated with disease states. Both voxelwise and regionof-interest methods of analysis are capable of integrating population differences in tract amplitude (streamline count or density), given proper alignment of the tracts of interest. However, quantification of tract differences (between groups, or longitudinally within individuals) has been hampered by two related features of white matter. First, it is unknown to what extent healthy individuals differ in the precise location of white matter tracts, and to what extent experimental factors influence perceived tract location. Second, white matter lacks the gross neuroanatomical features (e.g., gyri, histological subtyping) that make parcellation of grey matter plausible - determining where tracts should lie within larger white matter structures is difficult. Accurately quantifying tractographic connectivity between individuals is thus inherently linked to the difficulty of identifying and aligning precise tract location. Tractography is often utilized to study neurological diseases in which the precise structural and connectivity abnormalities are unknown, underscoring the importance of accounting for individual differences in tract location when evaluating the strength of structural connectivity. We set out to quantify spatial variance in tracts aligned through a standard, whole-brain registration method, and to assess the impact of location mismatch on groupwise assessments of tract amplitude. We then developed a method for tract alignment that enhances the existing standard whole brain registration, and then tested whether this method improved the reliability of groupwise contrasts. Specifically, we conducted seed-based probabilistic diffusion tractography from primary motor, supplementary motor, and visual cortices, projecting through the corpus callosum. Streamline counts decreased rapidly with movement from the tract center (-35% per millimeter); tract misalignment of a few millimeters caused substantial compromise of amplitude comparisons. Alignment of tracts peak-to-peak is essential for accurate amplitude comparisons. However, for all transcallosal tracts registered through the whole-brain method, the mean separation distance between an individual subject's tract and the average tract (3.2 mm) precluded accurate comparison: at this separation, tract amplitudes were reduced by 74% from peak value. In contrast, alignment of subcortical tracts (thalamo-putaminal, pallido-rubral) was substantially better than alignment for cortical tracts; whole-brain registration was sufficient for these subcortical tracts. We demonstrated that location mismatches in cortical tractography were sufficient to produce false positive and false negative amplitude estimates in both groupwise and longitudinal comparisons. We then showed that our new tract alignment method substantially reduced location mismatch and improved both reliability and statistical power of subsequent quantitative comparisons.