The Lambert $W$ function has the property that
$$ W'(x) = \frac{W(x)}{x[1+W(x)]}, $$
and using this one can show that its Taylor expansion about $x=a$ has the form
$$ W(x) = W(a) + \sum_{n=1}^{\infty} \frac{P_n[W(a)]}{n!\,[1+W(a)]^{2n-1}} \left(\frac{W(a)}{a}\right)^n (x-a)^n, $$
where $(P_n(x))_{n\in\mathbb{N}}$ is a sequence of polynomials defined recursively by $P_1(x) = 1$ and
$$ \tag{$*$} P_{n+1}(x) = \Bigl[1-(n+1)(3+x)\Bigr] P_n(x) + (1+x) P_n'(x). $$
Calculating the first 12 polynomials in the sequence I see that the constant term $P_n(0)$ is nonzero. I would like to know that $P_n(0) \neq 0$ for all $n$, but I don't see how to show it.
It appears that the recurrence $(*)$ can be solved "explicitly" as a mess of sums and products, but this seems like trading one problem for a bigger one.