Polynomial reduction, designed first for hypergeometric terms, can be used to automatically prove and generate new hypergeometric identities from old ones. In this paper, we extend the reduction method to holonomic sequences. As applications, we describe an algorithmic way to prove and generate new multi-summation identities. Especially we present new families of $π$-series involving Domb numbers and Franel numbers.