Let $\sigma>0$ and $$Q(x,\;\cdot\;):=\mathcal N(x,\sigma^2)\;\;\;\text{for }x\in\mathbb R.$$ Note that $Q$ is a Markov kernel on $(\mathbb R,\mathcal B(\mathbb R))$.
Now, let $X,Y$ be real-valued random variables on a common probability space $(\Omega,\mathcal A,\operatorname P)$. Let $X_\ast\operatorname P$ and $Y_\ast\operatorname P$ denote the distribution of $X$ and $Y$, respectively, and asume that $$Y_\ast\operatorname P=(X_\ast\operatorname P)Q,$$ where the right-hand side denotes the composition of $X_\ast\operatorname P$ and $Q$.
Are we able to show that $(Y-X)_\ast\operatorname P=\mathcal N(0,\sigma^2)$?
Intuitively, if $X=x\in\mathbb R$, then $Y\sim\mathcal N(x,\sigma^2)$ and it's a well-known fact that $Y-x\sim\mathcal N(0,\sigma^2)$.
From the mentioned property of the normal distribution, we know that $$\int 1_B(y-x)\:\mathcal N(x,\sigma^2)({\rm d}y)=\mathcal N(0,\sigma^2)(B)\;\;\;\text{for all }(x,B)\in\mathbb R\times\mathcal B(\mathbb R)\tag1.$$ Now, it's easy to see that $$\operatorname P\left[(X,Y)\in\;\cdot\;\right]=X_\ast\operatorname P\otimes\:Q\tag2$$ and hence \begin{equation}\begin{split}\operatorname P\left[Y-X\in B\right]&=\int\operatorname P\left[X\in{\rm d}x\right]\int Q(x,{\rm d}y)1_B(y-x)\\&=\int\operatorname P\left[X\in{\rm d}x\right]\int\mathcal N(0,\sigma^2)({\rm d}y)1_B(y)=\mathcal N(0,\sigma^2)(B)\end{split}\tag3\end{equation} for all $B\in\mathcal B(\mathbb R)$.