In the Derivation of CAVI, it rewrites the ELBO using iterated expectation. Then, it becomes the negative KL divengence between $q_j(z_j)$ and $q^*_j(z_j)$.

[Variational Inference: A Review for Statisticians][1]


  [1]: https://arxiv.org/abs/1601.00670