Calculation of NMR chemical shifts at second-order many-body perturbation theory using gauge-including atomic orbitals