On the parabolic Harnack inequality for non-local diffusion equations