Accounting for epistasis improves genomic prediction of phenotypes with univariate and bivariate models across environments