While diffusion MRI is typically used to estimate microstructural properties of tissue in volumetric elements (voxels), more specificity can be obtained by separately modelling the properties of individual fibre populations within a voxel. In the context of cross-subjects modelling, these fixel-based analyses are usually performed in two stages. Crossing fibre modelling is first performed in each subject to produce fixels, and these are subsequently modelled across subjects following registration and fibre population reassignment. Here, we introduce a new hierarchical framework for fitting crossing fibre models to diffusion MRI data in a population of subjects. This hierarchical setup guarantees that the crossing fibres are consistent by construction and, therefore, comparable across subjects. We propose an expectation-maximisation algorithm to fit the model, which can scale to large number of subjects. This approach produces a template for crossing fibre populations in the white matter which can be used to estimate fibre-specific parameters that are consistent across subjects, hence providing data that are by construction suitable for fixel-based statistical analyses.

Anatomical magnetic resonance imaging (MRI) templates of the brain are essential to group-level analyses and image processing pipelines, as they provide a reference space for spatial normalisation. While it has become common for studies to acquire multimodal MRI data, many templates are still limited to one type of modality, usually either scalar or tensor based. Aligning each modality in isolation does not take full advantage of the available complementary information, such as strong contrast between tissue types in structural images, or axonal organisation in the white matter in diffusion tensor images. Most existing strategies for multimodal template construction either do not use all modalities of interest to inform the template construction process, or do not use them in a unified framework. Here, we present multimodal, cross-sectional templates constructed from UK Biobank data: the Oxford-MultiModal-1 (OMM-1) template and age-dependent templates for each year of life between 45 and 81 years. All templates are fully unbiased to represent the average shape of the populations they were constructed from, and internally consistent through jointly informing the template construction process with T1-weighted (T1), T2-weighted fluid-attenuated inversion recovery (T2-FLAIR), and diffusion tensor imaging (DTI) data. The OMM-1 template was constructed with a multiresolution, iterative approach using 240 individuals in the 50–55-year age range. The age-dependent templates were estimated using a Gaussian process, which describes the change in average brain shape with age in 37,330 individuals. All templates show excellent contrast and alignment within and between modalities. The global brain shape and size are not preconditioned on existing templates, although maximal possible compatibility with MNI-152 space was maintained through rigid alignment. We showed benefits in registration accuracy across two datasets (UK Biobank and HCP), when using the OMM-1 as the template compared with FSL’s MNI-152 template, and found that the use of age-dependent templates further improved accuracy to a small but detectable extent. All templates are publicly available and can be used as a new reference space for uni- or multimodal spatial alignment.

We present MMORF—FSL’s MultiMOdal Registration Framework—a newly released nonlinear image registration tool designed primarily for application to magnetic resonance imaging (MRI) images of the brain. MMORF is capable of simultaneously optimising both displacement and rotational transformations within a single registration framework by leveraging rich information from multiple scalar and tensor modalities. The regularisation employed in MMORF promotes local rigidity in the deformation, and we have previously demonstrated how this effectively controls both shape and size distortion, leading to more biologically plausible warps. The performance of MMORF is benchmarked against three established nonlinear registration methods—FNIRT, ANTs, and DR-TAMAS—across four domains: FreeSurfer label overlap, diffusion tensor imaging (DTI) similarity, task-fMRI cluster mass, and distortion. The evaluation is based on 100 unrelated subjects from the Human Connectome Project (HCP) dataset registered to the Oxford-MultiModal-1 (OMM-1) multimodal template via either the T1w contrast alone or in combination with a DTI/DTI-derived contrast. Results show that MMORF is the most consistently high-performing method across all domains—both in terms of accuracy and levels of distortion. MMORF is available as part of FSL, and its inputs and outputs are fully compatible with existing workflows. We believe that MMORF will be a valuable tool for the neuroimaging community, regardless of the domain of any downstream analysis, providing state-of-the-art registration performance that integrates into the rich and widely adopted suite of analysis tools in FSL.