The aim of this work is to show how to renormalize the nucleon-nucleon interaction at next-to-next-to-leading order using a. systematic subtractive renormalization approach with multiple subtractions. As an example, we calculate the phase shifts for the partial waves with total angular momentum J = 2. The intermediate driving terms at each recursive step as well as the renormalized T-matrix are also shown. We conclude that our method is reliable for singular potentials such as the two-pion exchange and derivative contact interactions.