Disentangling the variability of the superficial white matter organization using regional-tractogram-based population stratification

•Study of the short superficial white matter connections and their relationship with the variability of cortical morphology.•Method for inferring fine grained regional atlases of the short bundles, from dMRI data, using a stratification strategy, based on a distance measure between tractography subsets.•New distance measures: (i) fibers angle-based distance and (ii) tractography subsets fiber-based distance.•A greater number of bundles is found in cortical variable areas when using a stratification strategy. Each variation of the cortical folding pattern implies a particular rearrangement of the geometry of the fibers of the underlying white matter. While this rearrangement only impacts the ends of the long pathways, it may affect most of the trajectory of the short bundles. Therefore, mapping the short fibers of the human brain using diffusion-based tractography requires a dedicated strategy to overcome the variability of the folding patterns. In this paper, we propose a fiber-based stratification strategy splitting the population into homogeneous groups for disentangling the superficial white matter bundle organization. This strategy introduces a new refined fiber distance which includes angular considerations for inferring fine-grained atlases of the short bundles surrounding a specific sulcus and a subtractogram distance that quantifies the similitude between fiber sets of two different subjects. The stratification splits the population into groups with similar regional fiber organization using manifold learning. We first successfully test the hypothesis that the main source of variability of the regional fiber organization is the variability of the regional folding pattern. Then, in each group, we proceed with the automatic identification of the most stable bundles, at a higher granularity level than what can be achieved with the non-stratified whole population, enabling the disentanglement of the very variable configuration of the short fibers. Finally, the method searches for bundle correspondence across groups to build a population level atlas. As a proof of concept, the atlas refinement achieved by this strategy is illustrated for the fibers that surround the central sulcus and the superior temporal sulcus using the HCP dataset.