Block coordinate proximal gradient methods with variable Bregman functions for nonsmooth separable optimization