The Hubbard model is described by nearest-neighbor hopping $H_t = -t\sum [c^{\dagger}_{\sigma,i}c_{\sigma,j}+\rm{h.c.} ]$, which preserves the spin-orientation, plus an on-site Coulomb replusions $H_U = U \sum n_{i,\uparrow}n_{i,\downarrow}$. Both of these terms are invariant to a global rotation of the spins. The first term will still describe each flavor of spins hopping from site to site, and the other will still be $U$ if there are two electrons on the site and zero otherwise. Therefore, a global rotation of the spins, which is described by $SU(2)$, leaves the Hamiltonian invariant. This helps us as we can characterize the states with total spin and total $S_z$, which are conserved.
In order to expand to an $SU(N)$ symmetry, we need to characterize the particles with higher degrees of freedom, such that their symmetry is not spin-rotations of $SU(2)$ but a general $SU(N)$.
Edit after question in comments:
For $SU(N)$ it's not enough to have larger spins. An example of how to generate $SU(N)$ might be along the following lines - instead of just one type of electrons let's say we have $N$-types of electrons, each coming from a different source. For example, in the 1d Hubbard model, we can proximitize $N$ wires of spinless electrons. So now instead of labeling each by the spin index we label it by the index $\nu=1,\ldots,N$ that tells us at which wire the electron is. We can now define "rotations" in the basis of these wires, which simply signify basis change. These are rotations in $N$-dimensional space, that will be characterized by the symmetry group $SU(N)$. Depending on the type of Hamiltonian and which terms we allow, it might be symmetric under this rotations (that is - it will not care in which basis we write our electrons). We usually call this extra degree of freedom "flavor".