Discrete orthogonal polynomials and difference equations of several variables

The goal of this work is to characterize all second order difference operators of several variables that have discrete orthogonal polynomials as eigenfunctions. Under some mild assumptions, we give a complete solution of the problem.