Tree-based variational inference for Poisson log-normal models