High-dimensional multivariate forecasting with low-rank Gaussian Copula Processes