Multiple commutation relations in the models with gl(2|1) symmetry

We consider quantum integrable models with $\mathfrak{gl}(2|1)$ symmetry. We derive a set of multiple commutation relations between the monodromy matrix entries. These multiple commutation relations allow us to obtain different representations for Bethe vectors.