In this work we propose a new multivariate pure jump model. We fully characterize a multivariate Lévy process with finite- and infinite-activity components in positive and negative jumps. This process generalizes the variance gamma process, featuring a ‘stochastic volatility’ effect due to Poisson randomized intensities of positive and negative gamma jumps. Linear and nonlinear dependence is introduced, without restrictions on marginal properties, separately on both positive and negative jumps and on both finite- and infinite-activity jumps. Such a new approach provides greater flexibility in calibrating nonlinear dependence than in other comparable Lévy models in the literature. The model is very tractable and a straightforward multivariate simulation procedure is available. An empirical analysis shows an almost perfect fit of option prices across a span of moneyness and maturities and a very accurate multivariate fit of stock returns in terms of both linear and nonlinear dependence. A sensitivity analysis of multi-asset option prices emphasizes the importance of the proposed new approach for modeling dependence.