The combinatorial identity goes as follows:
$$ \sum_{k=0}^\ell {n+k-1 \choose k} {n-k-1 \choose {\ell-k}} = {2n-1 \choose \ell } \ ,\ell \leqslant n-1 $$
Intuitively, the RHS counts all (0,1)-strings of length $(2n-1)$ which contain $\ell$ 0(s), and a string of this form must contain an $n-$th 1. Classifying the loci of the $n-$th 1 (from $n$ to $n+\ell$) yields the LHS.
However, is there a generating function that generates this identity? I find this intriguing. Thanks a lot if anybody could think about this.
Recall we can extend $\binom{n}k$ to negative $k$ with the convention $\binom{n}k=\frac{n(n-1)\cdots (n-k+1)}{k!}$, in which case we have $$ \binom{-n}k=(-1)^k\binom{n+k-1}k. $$ We can use this to write your equation as $$ \sum_{k=0}^\ell (-1)^k\binom{-n}k\cdot (-1)^{\ell-k}\binom{-n+\ell}{\ell-k}=(-1)^\ell\binom{-2n+\ell}{\ell},\tag{*} $$ Using Newton's binomial theorem $(1+x)^m=\sum_k \binom{m}k x^k$, valid for all $m\in \mathbb R$, $(*)$ is the $[x^\ell]$ cooeficient of the equation $$ \boxed{(1-x)^{-n}\cdot (1-x)^{-n+\ell}=(1-x)^{-2n+\ell}.} $$