Multivariate orthogonal polynomials and operator theory