After this question popped up again, it seemed to me to scream out for a use of the Fourier Transform (FT).  I have decided to post this as an answer, since this approach is transparent and provides a systematic method to obtain all (distributional) solutions. 

A simple formal manipulation (to be justified below) shows that the FT $g$ of a solution satisfies the equation
$$(e^{-iy}+1+iy)g(y)=0.$$

In classical terms, this only provides the zero solution but in the sense of distributions it has many non-trivial ones, in fact, suitable combinations of $\delta$-functions with singularities at the zeros of the  function in brackets.  This easily translates into a precise description of all possible solutions of the original equation.

[Historical aside on the FT for distributions: The problem of extending the FT to distributions, which was motivated by applications in physics, was studied in detail at the middle of the last century. Laurent Schwartz considered the case of tempered distributions in his seminal monograph. But this can be done for many other situations. The basic starting point for such an extension is an initial setting where it is an isomorphism between two l.c. spaces
of test functions (that is, functions with good smoothness and/or growth properties). One then uses transposes to translate this into an isomorphism between corresponding spaces of distributions. Schwartz used the (symmetric) case of the smooth functions of rapid decrease but there are many non symmetric cases which are of interest.]

We require the FT’s for arbitrary distributions on the line and so  used the fact that the FT is an isomorphism between the smooth functions on the line with compact support and a suitable space of entire functions of exponential growth. (This is the Paley-Schwartz theorem, with details in Strichartz).  Dualising, we get a version of the FT which works for ALL distributions--it takes its values in a space of analytic functionals (a superspace of  the dual space of the  Fréchet space of entire functions).  This can be described explicitly, but for our purposes it suffices to know that it contains the delta functions on the complex plane, together with suitable combinations.