If we denote the Bernoulli numbers by $B_n$, then $$ {a^{k+1}(a^{2k} - 1) B_{2k} \over 2k}\in \mathbb{Z}$$ for any $a\in \mathbb{Z}$ and $k\in\mathbb{N}$. This fact is often called the Sylvester–Lipschitz theorem, although I think it was known earlier; see for example A Note on Bernoulli Numbers by I. Sh. Slavutskii, J. Number Theory 53 (1995), 309–310. If $a=2$ then (up to a factor of $\pm 2^{k-1}$) there is a combinatorial interpretation of the above formula in terms of alternating permutations.
Is there a combinatorial interpretation for other values of $a$?