Using a dual representation of lattice fermion models that is based on spin-charge transformation and fermionization of the original description, I derive an algorithm for diagrammatic Monte Carlo simulation of strongly correlated systems. This scheme allows eliminating large expansion parameters, as well as large corrections to the density matrix that generally prevent diagrammatic methods from being efficient in this regime. As an example, I compute the filling factor for the Hubbard model at infinite on-site repulsion and compare the results to controllable data obtained from numerical linked-cluster expansion. I find excellent agreement between the two methods, as well as rapid convergence of the diagrammatic series. I also report results for the momentum distribution and kinetic energy of the electrons.