We establish Harnack inequality and shift Harnack inequality for stochastic differential equation driven by $G$-Brownian motion. As applications, the uniqueness of invariant linear expectations and estimates on the $sup$-kernel are investigated, where the $sup$-kernel is introduced in this paper for the first time.