We show that the generating function of the symmetric group with respect to five particular statistics gives rise to an exponential Riordan array, whose inverse is the coefficient array of the associated orthogonal polynomials. This also provides us with an LDU factorization of the Hankel matrix of the associated moments.