Diffusion process models for evolution of neutral genes have the coalescent process underlying them. Models are reversible with transition functions having a diagonal expansion in orthogonal polynomial eigenfunctions of dimension greater than one, extending classical one-dimensional diffusion models with Beta stationary distribution and Jacobi polynomial expansions to models with Dirichlet or...