This paper studies the Sobolev regularity estimates of weak solutions of a class of singular quasilinear elliptic problems of the form u t − div[A(x, t, u, ∇u)] = div[F] with homogeneous Dirichlet boundary conditions over bounded spatial domains. Our main focus is on the case that the vector coefficients A are discontinuous and singular in (x, t)-variables, and dependent on the solution u. Global and interior weighted W 1,p (Ω, ω)-regularity estimates are established for weak solutions of these equations, where ω is a weight function in some Muckenhoupt class of weights. The results obtained are even new for linear equations, and for ω = 1, because of the singularity of the coefficients in (x, t)-variables.