Objective: To investigate whole brain white matter changes in multiple sclerosis (MS) by multiple diffusion indices, we examined patients with diffusion tensor imaging and utilized tract-based spatial statistics (TBSS) method to analyze the data. Methods: Forty-one relapsing-remitting multiple sclerosis (RRMS) patients and 41 age- and gender-matched normal controls were included in this study. Diffusion weighted images were acquired by employing a single-shot echo planar imaging sequence on a 1.5 T MR scanner. Voxel-wise analyses of multiple diffusion metrics, including fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD) and radial diffusivity (RD) were performed with TBSS. Results: The MS patients had significantly decreased FA (9.11%), increased MD (8.26%), AD (3.48%) and RD (13.17%) in their white matter skeletons compared with the controls. Through TBSS analyses, we found abnormal diffusion changes in widespread white matter regions in MS patients. Specifically, decreased FA, increased MD and increased RD were involved in whole-brain white matter, while several regions exhibited increased AD. Furthermore, white matter regions with significant correlations between the diffusion metrics and the clinical variables (the EDSS scores, disease durations and white matter lesion loads) in MS patients were identified. Conclusion: Widespread white matter abnormalities were observed in MS patients revealed by multiple diffusion metrics. The diffusion changes and correlations with clinical variables were mainly attributed to increased RD, implying the predominant role of RD in reflecting the subtle pathological changes in MS.