Not all voxels are created equal: reducing estimation bias in regional NODDI metrics using tissue-weighted means

Abstract

Neurite orientation dispersion and density imaging (NODDI) estimates microstructural properties of brain tissue relating to the organisation and processing capacity of neurites, which are essential elements for neuronal communication. Descriptive statistics of NODDI tissue metrics are commonly analysed in regions-of-interest (ROI) to identify brain-phenotype associations. Here, the conventional method to calculate the ROI mean weights all voxels equally. However, this produces biased estimates in the presence of CSF partial volume. This study introduces the tissue-weighted mean, which calculates the mean NODDI metric across the tissue within an ROI, utilising the tissue fraction estimate from NODDI to reduce estimation bias. We demonstrate the proposed mean in a study of white matter abnormalities in young onset Alzheimers disease (YOAD). Results show the conventional mean induces significant bias that correlates with CSF partial volume, primarily affecting periventricular regions and more so in YOAD subjects than in healthy controls. Due to the differential extent of bias between healthy controls and YOAD subjects, the conventional mean under- or over-estimated the effect size for group differences in many ROIs. This demonstrates the importance of using the correct estimation procedure when inferring group differences in studies where the extent of CSF partial volume differs between groups. These findings are robust across different acquisition and processing conditions. Bias persists in ROIs at higher image resolution, as demonstrated using data obtained from the third phase of the Alzheimers disease neuroimaging initiative (ADNI); and when performing ROI analysis in template space. This suggests that conventional ROI means of NODDI metrics are biased estimates under most contemporary experimental conditions, the correction of which requires the proposed tissue-weighted mean. The tissue-weighted mean produces accurate estimates of ROI means and group differences when ROIs contain voxels with CSF partial volume. In addition to NODDI, the technique can be applied to other multi-compartment models that account for CSF partial volume, such as the free water elimination method. We expect the technique to help generate new insights into normal and abnormal variation in tissue microstructure of regions typically confounded by CSF partial volume, such as those in individuals with larger ventricles due to atrophy associated with neurodegenerative disease.