Here's an attempt at a proof. I don't have a reference, or much intuition beyond this construction.
Consider first the case of a variable $X$ taking the value $a<0$ with probability $b/(|a|+b)$ and
taking the value $b>0$ with probability $|a|/(|a|+b).$ Pick random variables $Y,Z$ uniformly distributed on $[a,b],$ with $Y=Z+b$ when $Z<0$ and $Y=Z+a$ when $Z\geq 0.$ Then $Y-Z\stackrel{\mathrm d}= X$ as required.
For the general case we can take a mixture, by pairing up the events $X<0$ and $X>0$ in a suitable way.
Remark. We just need to write the law of $X$ as a weighted average of zero-mean probability distributions supported on two points. A slick way to do this is to use Choquet's theorem and a characterization of extreme points of the set of zero-mean variables, e.g. Winkler's "Extreme Points of Moment Sets". The remainder of this answer is just a more explicit version of this.
Let $\mu$ be the law of $X.$ Let $M=\tfrac 1 2 \mathbb E|X|.$
There exists a Borel function $f^+:[0,M]\to\mathbb R_{>0}$ such that the pushforward of Lebesgue measure $f^+_*\lambda$ is the measure $x d\mu(x).$
For example, the generalized inverse of the function $t\mapsto \int_0^t x d\mu(x).$
Similarly there exists a Borel function $f^-:[0,M]\to\mathbb R_{<0}$ such that $f^-_*\lambda$ is the measure $|x| d\mu(x).$
The pushforward of $(1/f^+(x))d\lambda(x)$ under $f^+$ is $\mu$ restricted to $\mathbb R_{>0}.$
And similarly the pushforward of $(1/|f^-(x)|)d\lambda(x)$ under $f^-$ is $\mu$ restricted to $\mathbb R_{<0}.$
Pick independent random variables $R,S,T$ such that
- $R$ is a Bernoulli random variable with $\mathbb P[R=0]=\mathbb P[X=0],$
- $S$ is a random variable on $[0,M]$ distributed according to the measure $((1/f^+(x))+(1/|f^-(x)|))d\lambda(x),$
- $T$ is uniformly distributed on $[0,1].$
Then define $Y$ and $Z$ by:
- wherever $R=0,$ we take $Y=Z=0,$ and skip the next steps.
- $Z$ is $T$ rescaled from $[0,1]$ to the interval $[f^-(S),f^+(S)].$
- $Y$ is $Z+f^+(S)$ when $Z<0$
- $Y$ is $Z+f^-(S)$ when $Z\geq 0.$
$Z$ has the same distribution as $Y.$
The random variable $Y-Z$ has the same probability of being zero as $X.$
For any Borel $U\subseteq\mathbb R_{>0},$ the probability that $f^+(S)\in U$ and $Y-Z=f^+(S)$ is $\int (1/f^+(x))1_{f^+(x)\in U}d\lambda(x)=\mu(U).$
So the laws of $Y-Z$ and $X$ are the same on $\mathbb R_{>0}$.
And similarly for $\mathbb R_{<0}.$
So this construction ensures $Y-Z\stackrel{\mathrm d}= X.$