We present a unifying algebraic representation for orthogonal polynomial sequences both in the univariate and multivariate case. This representation gives new formulae for orthogonal polynomials arising from probability and alternative proof of the Favard’ theorem. Moreover orthogonal polynomials are connected via cumulants to invariant theory of binary forms.