In this work, the shape evolution and the possible shape coexistence are explored in the odd-$A$ Ne isotopes in the framework of the multidimensionally constrained relativistic-mean-field~(MDC-RMF)~model. By introducing the $s_\Lambda$ and $p_{\Lambda}$ hyperons, the impurity effects to the nuclear shape, energy, size, and density distribution are investigated. For the $NN$ interaction, the PK1 parameter set was adopted and for the $\Lambda N$ interaction, the PK1-Y1 parameter set was taken. By blocking the unpaired odd neutron in different orbitals around the Fermi surface, the nuclear ground state along with the low-lying excited states are determined. The potential energy curves~(PECs), quadrupole deformations, nuclear r.m.s.~radii, binding energies, and the density distributions for the core nuclei as well as the corresponding hypernuclei are analyzed. By examining the PECs, great possibilities of the shape coexistences in $^{27,29}$Ne and a triple shape coexistence in $^{31}$Ne are found. For the impurity effects of $\Lambda$ hyperons, as behaved in the even-even Ne hypernuclear isotopes, the $s_{\Lambda}$ hyperon shows obvious shrinkage effect, which reduces the nuclear size and drives the nuclear shape to be more spherical; the $p_{\Lambda}$ hyperon occupying the $1/2^-[110]$ orbital is prolate and makes the nuclear shape more prolate; and the $p_{\Lambda}$ hyperon on the $3/2^-[101]$ orbital displays oblate shape and drives the nuclei to be more oblate.