I have reduced solving this question to proving the following identity, for $n, \ell \ge 0$: $$ (n-2\ell+1)^{n-1} \binom{n}{\ell-1} = \\ \frac{1}{2} \sum_{n_1+n_2=n-1}\left[ (n_1+1)^{n_1-1} (n_2-2\ell+1)^{n_2} \binom{n_2}{\ell-1} + \\ \ell \sum_{\ell_1+\ell_2=\ell} \left[(n_1-2\ell_1+1)^{n_1} \binom{n_1}{\ell_1-1} \frac{1}{\ell_1} (n_2-2\ell_2+1)^{n_2} \binom{n_2}{\ell_2-1} \frac{1}{\ell_2}\right] + (n_2+1)^{n_2-1} (n_1-2\ell+1)^{n_1} \binom{n_1}{\ell-1} \right] $$ where $n_1,n_2$ and non-negative integers and $\ell_1, \ell_2$ are positve integers.
I have checked this computationally for small values of $n$ and $\ell$, and it seems to be good.
It is not hypergeometric, but I think it may be "Abel-type" in the sense of this paper.
How could I prove an identity like this?