  The strong $D^{*}_sD_{s}\eta^{(\prime)}$ and $B^{*}_sB_{s}\eta^{(\prime)}$
vertices are studied and the relevant couplings are calculated in the context
of the light-cone QCD sum rule method with twist-4 accuracy by including the
next-to-leading order corrections. In the analysis, both the quark and gluon
components of the $\eta$ and $\eta^{\prime}$ mesons and the axial anomaly
corrected higher twist distributions are included.
