We consider Kolmogorov operator $-\nabla \cdot a \cdot \nabla + b \cdot \nabla$ with measurable uniformly elliptic matrix $a$ and prove Gaussian lower and upper bounds on its heat kernel under minimal assumptions on the vector field $b$ and its divergence ${\rm div\,}b$. More precisely, we prove: (1) Gaussian lower bound, provided that ${\rm div\,}b \geq 0$, and $b$ is in the class of form-bounded vector fields (containing e.g.\,the class $L^d$, the weak $L^d$ class, as well as some vector fields that are not even in $L_{\rm loc}^{2+\varepsilon}$, $\varepsilon>0$); in these assumptions, the Gaussian upper bound is in general invalid; (2) Gaussian upper bound, provided that $b$ is form-bounded, and the positive part of ${\rm div\,}b$ is in the Kato class; in these assumptions, the Gaussian lower bound is in general invalid; (3) Gaussian upper and lower bounds, provided that $b$ is form-bounded, ${\rm div\,}b$ is in the Kato class; (4) A priori Gaussian upper and lower bounds, provided that $b$ is in a large class containing the class of form-bounded vector fields, ${\rm div\,}b$ is in the Kato class.