Gaussian estimates for degenerate diffusion