  In this paper we introduce three combinatorial models for symmetrized
poly-Bernoulli numbers. Based on our models we derive generalizations of some
identities for poly-Bernoulli numbers. Finally, we set open questions and
directions of further studies.
