We will discuss harmonic functions with respect to Markov processes. Recent results for estimates of heat kernel for jump processes (or non-local operators) will be introduced.