Hi Hugo,
I have been thinking about this but do not have a good insight, maybe one thing worth to try is used the masked array to code the data matrix (as in our examples), but passing a local RV to Variational Inference as coded for the latent variable. In the local RV you can have a small neural network that specifically dealt with the missing variables.