In this paper, we mainly study two spectral Galerkin approximation schemes of the FitzHugh-Nagumo (FHN) neuron model for the nerve impulse propagation, which contains a small parameter perturbation and strong nonlinearity. To this end, we consider two kinds of temporal discretization by a semi-implicit Euler method and a Crank-Nicolson method, and the spatial discretization by the spectral Galerkin method for the perturbed FHN neuron model. In particular, in order to obtain the second-order approximation in time for the numerical scheme, we modify the nonlinear term f(un) in the Crank-Nicolson spectral Galerkin (CNSG) scheme. We derive the error estimations with optimal convergence rates of the numerical solutions by the first-order semi-implicit spectral Galerkin (SISG) approximation scheme and the second-order modified CNSG (MCNSG) approximation scheme for the perturbed FHN neuron model. Finally, some numerical examples of one-dimensional and two-dimensional nonlinear FHN models with periodic boundary conditions are carried out to verify that the approximate solutions of the proposed SISG and MCNSG schemes satisfy the proved theoretical results.