This paper aims to study the asymptotic behaviour of the fundamental solutions (heat kernels) of non-local (partial and pseudo differential) equations with fractional operators in time and space. In particular, we obtain exact asymptotic formulas for the heat kernels of time-changed Brownian motions and Cauchy processes. As an application, we obtain exact asymptotic formulas for the fundamental solutions to the $n$-dimensional fractional heat equations in both time and space \begin{gather*} \frac{\partial^β}{\partial t^β}u(t,x) = -(-Δ_x)^γu(t,x), \quad β,γ\in(0,1). \end{gather*}