A pointwise gradient bound for weak solutions to Dirichlet problem for quasilinear elliptic equations −div(A(x, ∇u)) = µ is established via Wolff type potentials. It is worthwhile to note that the model case of A here is the non-degenerate p-Laplacian operator. The central objective is to extend the pointwise regularity results in [Q.-H. Nguyen, N. C. Phuc, Pointwise gradient estimates for a class of singular quasilinear equations with measure data, J. Funct. Anal. 278(5) (2020), 108391] to the very singular case 1 < p ≤ 3n−2 2n−1 , where the data µ on right-hand side is assumed belonging to some classes that close to L 1 . Moreover, a global pointwise estimate for gradient of weak solutions to such problem is also obtained under the additional assumption that Ω is sufficiently flat in the Reifenberg sense.