Introduction: Brain viscoelasticity as assessed by magnetic resonance elastography (MRE) has been discussed as a promising surrogate of microstructural alterations due to neurodegenerative processes. Existing studies indicate that multiple sclerosis (MS) is associated with a global reduction in brain stiffness. However, no study to date systematically investigated the MS-related characteristics of brain viscoelasticity separately in normal-appearing white matter (NAWM), deep gray matter (DGM) and T2-hyperintense white matter (WM) lesions. Methods: 70 MS patients and 42 healthy volunteers underwent whole-cerebral MRE using a stimulated echo sequence (DENSE) with a low-frequency mechanical excitation at 20 Hertz. The magnitude |G∗| (Pa) and phase angle φ (rad) of the complex shear modulus G∗ were reconstructed by multifrequency dual elasto-visco (MDEV) inversion and related to structural imaging and clinical parameters. Results: We observed φ in the thalamus to be higher by 4.3 % in patients relative to healthy controls (1.11 ± 0.07 vs. 1.06 ± 0.07, p