We give a formula that expresses a sum of products of poly-Bernoulli numbers of negative index as a linear combination of poly-Bernoulli numbers. More generally, we show that if a two-variable formal power series satisfies a certain partial differential equation, then its coefficients satisfy this type of formula. As an appendix, we solve this partial differential equation.