Thank you so much @Nicholas_Bokulich! My last idea was to use the log2FC (between my 2 time points) and look at the treatment effect with ANCOM. Not appropriate neither?
I am working on the q2-gneiss and will give a try as rhapsody!
It is difficult to find the best way to analyze these data because actually I am equally interested by the time AND the treatment effect! 2 different questions but I would love to address both!