I explored different packages (such as mgcv, MCMCglmm, glmmADMB, etc) but they all have limitations. Either they don't allow for zero-inflated, zero-altered distributions or they don't allow for temporal auto-correlation through functions as corAR() or corARMA().