Harmonic functions for a class of integro-differential operators