If two multivariate normal random variables $X\sim\mathcal N(\mu_X,\Sigma_X)$ and $Y\sim\mathcal N(\mu_Y,\Sigma_Y)$ are independent (and of the same dimension), then their sum is still normal, and you can sum the mean and variance directly:
$$X+Y\sim\mathcal N(\mu_X+\mu_Y,\Sigma_X+\Sigma_Y).$$
The easiest way to see this is with the expression of the characteristic function, as can be found on the Wikipedia article for example:
$$\phi_X(u)=\exp(iu'\mu_X-u'\Sigma_Xu).$$
Since the characteristic function of the sum of two independent random variables is the product of their respective characteristic functions, then
\begin{eqnarray}
\phi_{X+Y}(u)
&=&
\exp(iu'\mu_X-u'\Sigma_Xu)\exp(iu'\mu_Y-u'\Sigma_Yu)\\
&=&
\exp(iu'\mu_X-u'\Sigma_Xu+iu'\mu_Y-u'\Sigma_Yu)\\
&=&
\exp(iu'(\mu_X+\mu_Y)-u'(\Sigma_X+\Sigma_Y)u).
\end{eqnarray}
Now back to your original question, it suffices to see that $-Y$ has the same variance as $Y$, but with opposite mean to conclude:
$$-Y\sim \mathcal N(-\mu_Y,\Sigma_Y),$$
so indeed
$$X-Y\sim\mathcal N(\mu_X-\mu_Y,\Sigma_X+\Sigma_Y).$$