Large sieve estimate for multivariate polynomial moduli and applications