In this paper, we study a fully discrete weak Galerkin finite element method for solving parabolic integro-differential equations baesd on polygons/polyhedrons mesh of any shape. We prove the existence and uniqueness of the solutions of the fully discrete schemes of the method. Corresponding error estimates are derived. Numerical experiments are provided to verify the theoretical results.