We study the function and show that the poly-Bernoulli numbers introduced in our previous paper are expressed as special values at negative arguments of certain combinations of these functions. As a consequence of our study, we obtain a series of relations among multiple zeta values.