Meixner polynomials in several variables satisfying bispectral difference equations

We construct a set $M_d$ whose points parametrize families of Meixner polynomials in $d$ variables. There is a natural bispectral involution $b$ on $M_d$ which corresponds to a symmetry between the variables and the degree indices of the polynomials. We define two sets of $d$ commuting partial difference operators diagonalized by the polynomials. One of the sets consists of difference operators acting on the variables of the polynomials and the other one on their degree indices, thus proving their bispectrality. The two sets of partial difference operators are naturally connected via the involution $b$.