An automated pipeline for mouse brain morphometric MRI
Marco Pagani1,2, Mario Damiano1, Alberto Galbusera1, Sotirios A Tsaftaris3,4, and Alessandro Gozzi1

1Istituto Italiano di Tecnologia, Rovereto, Italy, 2Center for Mind and Brain Sciences, Rovereto, Italy, 3IMT - Institute for Advanced Studies, Lucca, Italy, 4Institute of Digital Communications, School of Engineering, The University of Edinburgh, Edinburgh, United Kingdom

Synopsis

We provide a detailed description of registration-based procedures for voxel based morphometry, cortical thickness estimation and automated anatomical labelling of the mouse brain. To illustrate our procedures, we described their application to quantify morphological differences in two inbred mouse strains characterised by different social behaviour. We show that our approach can reliably detect both focal and large scale gray matter alterations using complementary readouts. The operational workflows described here are expected to help the implementation of rodent morphoanatomical methods by non-expert users, and promote the use of these tools across the preclinical neuroimaging community.

PURPOSE

Brain morphometric mapping in genetically modified mouse lines may permit to relate genetic mutations to clinically relevant endophenotypes without the complexity of genetic heterogeneity and the uncontrolled impact of gene-gene and gene-environment interactions typical of human populations1. However, substantial differences in the anatomical organisation of human and rodent brain prevent a straightforward extension of clinical neuroimaging tools to mouse brain imaging and most published studies lack a detailed description of the complex workflow required to process, analyse and quantify structural MRI alterations. To address these issues, here we provide a detailed methodological description of our semi-automated operational workflow for voxel based morphometry, cortical thickness estimation and automated anatomical mapping of the mouse brain. To show the capabilities of our methods based on the ANTs toolkit2, we describe its application to quantify morphological alterations in socially-impaired BTBR T+Itpr3tf/J mice with respect to normo social C57BL/6J controls3-4.

METHODS

Preprocessing. An automatic registration-based approach has been devised to automate skull-stripping and avoid error-prone manual segmentation, and two independent bias corrections are executed to correct for intensity non-uniformity within the same tissue type (Fig.1). An iterative five step multi-scale alignment process is performed using diffeomorphic registration to construct a study based template from skull stripped and bias corrected wild type mice and establish a common reference space for all the subsequent analyses. Automated anatomical labelling. Registration based labelling is performed to assess volumetric anatomical differences in gross morphology and relies on non-linear registrations of two labelled neuroanatomical reference MRI atlas via the study based template. Cortical thickness. A diffeomorphic registration based cortical thickness measure is used to establish a continuous one-to-one correspondence between the inner and outer cortical surfaces and estimate a voxelwise distance between the two interfaces. VBM. Individual images are registered to the study based template and segmented to calculate tissue probability maps (Fig.2). The separation of the different tissues is improved by initializing the process with the probability maps of the study based template previously segmented. The Jacobian determinants of the deformation field were extracted and applied to modulate the grey matter probability maps calculated during the segmentation. This expedient permits the analysis of grey matter probability maps in terms of local volumetric variation instead of tissue density. The resulting modulated grey matter probability maps were then smoothed using a Gaussian kernel with a sigma of three voxel width. Statistics. Standard non-parametric Monte Carlo test with 5000 random permutations was performed coupled with threshold-free cluster enhancement and multiple comparisons correction using a cluster-based threshold of 0.01.

RESULTS

Cross-strain volumetric analysis highlighted the presence of a general reduction in cortical and subcortical volume in BTBR mice with respect to controls. Further voxelwise investigations using whole-brain VBM revealed widespread and bilateral reductions in GM volume across dorsofrontal, cingulate, retrosplenial, occipital and parietal cortex in BTBR (Fig.3). In good agreement with the results of automated anatomical labelling and VBM mapping, a widespread reduction in mean cortical thickness was observed in BTBR mice (Fig.4). Importantly, manual measurements of exemplificative regions were highly correlated with cortical thickness estimates (Fig.5).

DISCUSSION

The vast majority of current morphometric studies on the mouse brain lacks a detailed description of the complex workflow required to process and analyse different readouts. We provide a detailed description of methods to quantify morphoanatomical alterations in socially-impaired BTBR mice. We show that these approaches can reliably detect both focal and large scale gray matter alterations using complementary readouts. Importantly, our VBM workflow can be straightforwardly extended to perform tensor based morphometry, where Jacobian maps can be used to localise inter-group differences in the local volumes of brain structures at the voxel level. Morphological alterations comparable to those we mapped have been recently described in BTBR mice by others5, thus permitting an empirical cross-laboratory validation of our findings.

CONCLUSION

The use of complementary semi-automated MRI morphometric measurers can help to pinpoint the pathological bases of brain morphometric changes of neuropathological origin. The detailed operational workflow of the present work is expected to help the implementation of fine-grained morphoanatomical mapping methods by non-expert users, promote future comparative study, and, above all, facilitate forward and back translation of MRI preclinical and clinical research findings.

Acknowledgements

No acknowledgement found.

References

1. Nestler EJ, Hyman SE. Animal models of neuropsychiatric disorders. Nat Neurosci, 2010; 13:1161-9.

2. Avants BB, Tustison N, Song G. Advanced normalization tools (ANTS). Insight J, 2009

3. Dodero L, Damiano M, Galbusera A et al. Neuroimaging evidence of major morpho-anatomical and functional abnormalities in the BTBR T+TF/J mouse model of autism. PLoS One, 2013; 8

4. Squillace M, Dodero L, Federici M et al. Dysfunctional dopaminergic neurotransmission in asocial BTBR mice. Transl Psychiatry, 2014; 4:e427.

5. Ellegood J, Babineau BA, Henkelman RM et al. Neuroanatomical analysis of the BTBR mouse model of autism using magnetic resonance imaging and diffusion tensor imaging. Neuroimage, 2013; 70:288-300.

Figures

Preprocessing results. The original subject image (a) is bias corrected before (b) and after (e) skull stripping (d). Note the improved bias field correction after skull stripping (f) with respect to the bias correction prior skull stripping (c), especially in the ventral part of the brain and in the ventricles

Study based template and tissue segmentation. A study based template of the B6 mice population obtained using the iterative diffeomorphic registration process and its corresponding tissue segmentation (a). The template is segmented using 6 different tissue classes which are used as a-priori information for individual estimation of gray matter in VBM. The different tissue classes are combined to obtain gray matter (b) and non gray matter components (c, white matter, plus ventricles and CSF).

VBM. Differences in local gray matter volume (GMV) are assessed combining gray matter probability maps and local Jacobian determinants. Statistical comparison showed widespread and bilateral reductions in GMV across dorsofrontal, cingulate, retrosplenial, occipital and parietal cortex as well as in subcortical structures in BTBR. VBM highlighted also small foci of increased GMV in the olfactory bulbs, in the medial prefrontal and insular cortex, in the amygdala and in the dorsal hippocampus

Cortical thickness estimation. Three-dimensional rendering views of average cortical thickness in BTBR and B6 mice (a). Statistical comparison showed significant cortical thickness thinning in parietal, temporal and peri-hippocampal cortex of BTBR mice. Increased thickness was observed in medial prefrontal and anterior insular regions of this strain (b)

Correlation plot between DiReCT outputs and manual measurements of cortical thickness. Secondary motor (M2), secondary somatosensory (S2) and auditory cortex (Au) were chosen as representative areas to validate our cortical thickness methodology. Representative measures from DiReCT and manual estimates are reported for selected cortical regions. A correlation plot of manual and automatic measurements highlighted an excellent correspondence between the readouts in terms of Pearson’s correlation



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
3376