Connectomes from streamlines tractography: Assigning streamlines to brain parcellations is not trivial but highly consequential

When using diffusion MRI streamlines tractograms to construct structural connectomes, ideally, each streamline should connect exactly 2 regions-of-interest (i.e. network nodes)as defined by a given brain parcellation scheme. However, the ill-posed nature of termination criteria in many tractography algorithms can cause streamlines apparently being associated with zero, one, or more than two grey matter (GM)nodes; streamlines that terminate in white matter or cerebrospinal fluid may even end up being assigned to nodes if the definitions of these nodes are not strictly constrained to genuine GM areas, resulting in a misleading connectome in non-trivial ways. Based on both in-house MRI data and state-of-the-art data provided by the Human Connectome Project, this study investigates the actual influence of streamline-to-node assignment methods, and their interactions with fibre-tracking terminations and brain parcellations, on the construction of pairwise regional connectivity and subsequent connectomic measures. Our results show that the frequency of generating successful pairwise connectivity is heavily affected by the convoluted interactions between the applied strategies for connectome construction, and that minor changes in the mechanism can cause significant variations in the within- and between-module connectivity strengths as well as in the commonly-used graph theory metrics. Our data suggest that these fundamental processes should not be overlooked in structural connectomics research, and that improved data quality is not in itself sufficient to solve the underlying problems associated with assigning streamlines to brain nodes. We demonstrate that the application of advanced fibre-tracking techniques that are designed to correct for inaccuracies of track terminations with respect to anatomical information at the fibre-tracking stage is advantageous to the subsequent connectome construction process, in which pairs of parcellation nodes can be more robustly identified from streamline terminations via a suitable assignment mechanism.