  In this paper, we establish sharp two-sided estimates for transition
densities of a large class of subordinate Markov processes. As applications, we
show that the parabolic Harnack inequality and H\"older regularity hold for
parabolic functions of such processes, and derive sharp two-sided Green
function estimates.
