In this paper, we calculate the scalar $a_0(980)$-meson leading-twist wave function by using light-cone harmonic oscillator model (LCHO), where the model parameters are determined by fitting the $\xi$-moments $\langle\xi_{a_0}^n\rangle_\zeta$ of its light-cone distribution amplitudes. Then, the $a_0(980)$-meson leading-twist light-cone distribution amplitudes with three different scales $\zeta= (1.0, 2.0, 5.2)~{\rm GeV}$ are given. After constructing the relationship between $a_0(980)$-meson leading-twist parton distribution functions/valence quark distribution function and its LCHO wave function, we exhibit the $\mathpzc q^{a_0}(x,\zeta)$ and $x\mathpzc q^{a_0}(x,\zeta)$ with different scales. Furthermore, we also calculate the Mellin moments of the $a_0(980)$-meson's valence quark distribution function $\langle x^n \mathpzc q^{a_0}\rangle_\zeta$ with $n = (1,2,3)$, i.e. $\langle x \mathpzc q^{a_0}\rangle_{\zeta_5} = 0.027$, $\langle x^2 \mathpzc q^{a_0}\rangle_{\zeta_5} = 0.018$ and $\langle x^3 \mathpzc q^{a_0}\rangle_{\zeta_5} = 0.013$. Finally, the scale evolution for the ratio of the Mellin moments $\mathpzc x^n_{\,a_0}(\zeta,\zeta_k)$ are presented.