Here's a way to approach this directly with complex analysis methods. By symmetry, it suffices to count zeroes in the upper half plane.
Observe that for any $k\in \mathbb{Z}$, along the line $\Im(z)=(2k+\frac{1}{2})\pi$ we have that $e^z$ is a positive real multiple of $i$. In particular, for any $t\in [0,1]$, the function $f_t(z)=e^z + (1-t)z+ti$ has no zero on that line for $k\geq 0$ (as then both $e^z$ and $(1-t)z+ti$ lie strictly in the upper half plane).
Now recall that $\frac{1}{2\pi i}\int_\gamma \frac{df_t}{f_t(z)}$ counts the number of zeroes of $f_t$ enclosed by the contour $\gamma$. If we let $\gamma$ be the boundary of the rectangle where $(2k+\frac{1}{2})\pi\leq \Im(z)\leq (2(k+1)+\frac{1}{2})\pi$ and $\Re(z) \in [-C_1,C_2]$, then by the above $f_t(z)$ has no zeroes on the horizontal boundary lines of this rectangle. But since $\frac{1}{2} \leq |(1-t)z+ti| \leq \max(|z|,1)$ we find $C_1$ and $C_2$ depending on $k$ ($C_1$ is roughly $\ln(2)$, and $C_2$ is roughly $\ln(2k\pi)$) for which $f_t(z)$ also has no zeroes on the vertical boundary lines of this rectangle. So $\frac{1}{2\pi i}\int_\gamma \frac{df_t}{f_t(z)}$ is an integer depending continuously on $t$, and so agrees for $t=0$ and $t=1$. For $t=0$ we get the number of zeroes in our rectangle for $f_0(z)=e^z+z$, while for $t=1$ we get $f_1(z)=e^z+i$. This has exactly one zero in our rectangle, $(2k+1)\pi i$. Since this argument works for all larger $C_1,C_2$ as well, it follows that $e^z+z$ has a unique zero in each of the horizontal strips $(2k+\frac{1}{2})\pi\leq \Im(z)\leq (2(k+1)+\frac{1}{2})\pi$, with real value bounded logarithmically by $k$. In particular, the absolute value of zeroes grows linearly.