High-dimensional support union recovery in multivariate regression