We apply the Virtual Element approach for the discretization of advection-diffusion-reaction problems with variable coefficients on polygonal grids. We analyze both the primal and the mixed formulation, and prove optimal convergence results in H1 and in L2. We provide numerical tests confirming the theoretical results.