In this talk, I will present recent progress in the study of heat kernels and parabolic Harnack inequalities for symmetric Markov processes that have both diffusive and jumping parts on general metric measure spaces. Under general volume doubling condition and some mild assumptions on the scaling functions, we establish stability results for two-sided estimates for heat kernels in terms of the jumping kernels, the generalized capacity inequalities, and Poincare inequalities. Stable characterizations of the associated parabolic Harnack inequalities will also be given. Our results hold on spaces even when the underlying spaces have walk dimensions are larger than 2.