The falling factorials are $(x)_n=x(x-1)\cdots(x-n+1)$ with $(x)_0:=1$. Define the forward shift $E$ and the discrete derivative $\delta=E-1$, respectively, by $$Ef(x)=f(x+1) \qquad \text{and} \qquad \delta f(x)=f(x+1)-f(x).$$ Let $\lambda=(\lambda_1,\dots,\lambda_k)\vdash n$ be a partition of a positive integer $n$. Now, define the operator $\Psi_{\lambda}$ to act on polynomials according to $$\Psi_{\lambda}=(E-1)^{-1}\prod_{j=1}^{k}(E^{\lambda_j}-1)$$ (where $k $ is the length of $\lambda$) so in particular denote $\Psi_{\lambda}(x):=\Psi_{\lambda}((x)_n)$.
For example, if $n=9$ and $\lambda=(4,2,2,1)$ then $(E-1)^{-1}(E^4-1)(E^2-1)^2(E-1)=(E^4-1)(E^2-1)^2=E^8-2E^6+2E^2-1$ and hence $\Psi_{\lambda}(x)=(x+8)_9-2(x+6)_9+2(x+2)_9-(x)_9$, or $$\Psi_{\lambda}(x)=8064x^6+141120x^4+213696x^2.$$ Experiments suggest to ask
is it true that $\Psi_{\lambda}(x)$ is always an even or an odd polynomial?