Hi good question, we are working on a more automatic system to make working with fermions easier. However, it might be the end of this summer at best before it's put in and fully tested.
In the meantime, using swap gates with fermionic operators should work well as long as you define fermionic swap gates. As you may know, these are like usual swap gates, except they include a minus sign when the swap operation swaps two sites that are both in fermion-parity-odd states.
If you do proceed with this route, you could test your code by using the AutoMPO system to help you out, because AutoMPO does handle fermion signs automatically and the toExpH function lets you approximately exponentiate a Hamiltonian defined via AutoMPO.