We study a symmetric diffusion $X$ on $\mathbb{R}^d$ in divergence form in a stationary and ergodic environment, with measurable unbounded and degenerate coefficients. We prove a quenched local central limit theorem for $X$, under some moment conditions on the environment; the key tool is a local parabolic Harnack inequality obtained with Moser iteration technique.