The gradient flow exact renormalization group (GFERG) is an exact renormalization group motivated by the Yang-Mills gradient flow and its salient feature is a manifest gauge invariance. We generalize this GFERG, originally formulated for the pure Yang-Mills theory, to vector-like gauge theories containing fermion fields keeping the manifest gauge invariance. For the chiral symmetry, we have two options: One possible formulation preserves the conventional form of the chiral symmetry and the other simpler formulation realizes the chiral symmetry in a modified form ร la Ginsparg-Wilson. We work out a gauge-invariant local Wilson action in QED to the lowest nontrivial order of perturbation theory. This Wilson action reproduces the correct axial anomaly in D = 2.