The purpose of this study was to analyze microstructural alterations in cerebral gray matter using non-Gaussian diffusion kurtosis imaging (DKI) in neuromyelitis optica spectrum disorder (NMOSD) patients with optic neuritis (NMOSD-ON). DKI was performed in 14 NMOSD-ON patients and 22 normal controls (NCs). DKI-derived metrics, including mean kurtosis (MK), radial kurtosis (RK), axial kurtosis (AK), fractional anisotropy (FA), and mean diffusivity (MD), were voxel-wisely compared by two-sample t-tests with gaussian random field (GRF) correction between the two groups. The correlations between altered DKI metrics and clinical features were analyzed. Compared with NCs, NMOSD-ON patients showed significantly decreased MK and RK both in the left inferior temporal gyrus (ITG), and decreased AK in the bilateral calcarine (CAL). While increased MD in the left fusiform gyrus (FFG), right CAL, and right hippocampus (HIP)/parahippocampal gyrus (PHG) were found. Furthermore, correlation analysis showed that mean deviation was negatively correlated with AK values of bilateral CAL and positively correlated with MD values of right CAL (q < 0.05, false discovery rate (FDR) corrected). For NMOSD-ON patients, microstructural abnormalities in the occipital visual cortex are correlated with clinical disability. These findings may provide complementary information to understand the neuropathological mechanisms underlying the impairments of cerebral gray matter in NMOSD-ON.