If $(u,v)$ are isothermal coordinates, then $u$ and $v$ are harmonic functions with respect to the Laplace-Beltrami operator on your Riemannian manifold, that is $\Delta u=\Delta v=0$.
Now, the equation $\Delta f=0$ characterizes the stationary states for the heat equation. The level curves for a harmonic function are therefore the isothermal curves for some heat distribution.
It follows that if $(u,v)$ is an isothermal system of coordinates, then the level curves, _i.e._ , the coordinate lines, are isothermal curves.