and the Lagrangian density of this interaction (W- W+ and others) can comes from the:
[itex]\mathcal{L} = - \frac{1}{4} F_{\mu \nu}^b F^{\mu \nu}_b[/itex]
which for non-abelian gauge group theories will give you interactions of fields with themselves (not only the kinetic terms).
The interactions of bosons with the Higgs scalar come from the Higg's potential term which contains: [itex]D_{\mu} \phi D^{\mu} \phi^*[/itex], where the gauge boson fields exist in the covariant derivatives [itex]D_{\mu}[/itex] (as connections)