In this paper, we establish the asymptotic validity and analyse the finite sample performance of a simple bootstrap procedure for constructing multi-step multivariate forecast densities in the context of non-Gaussian unrestricted VAR models. This bootstrap procedure avoids the backward representation, and, as a consequence, can be used to obtain multivariate forecast densities in, for example, VARMA or VAR-GARCH models. In the context of bivariate stationary VAR(p) models, we show that its finite sample properties are comparable to those of alternatives based on the backward representation. The bootstrap procedure is also implemented in a VAR-DCC model which lacks a backward representation. Finally, joint forecast densities of US quarterly ...