When 1 < γ < n, we study the regularity of W 1,γ -solutions of the anisotropic PDEBy focusing on the γ-homogeneity, we show that nonnegative subsolutions and supersolutions u are suitably bounded to ensure that full solutions are locally bounded. When u is nonnegative, this also implies a (weak) Harnack inequality. In the case when γ = 2, this includes divergence form uniformly elliptic PDEs in the special case when ρ(x, •) is assumed to be uniformly strongly convex. Analogous to the elliptic setting, when f, F ≡ 0, this ensures a Bernstein-type theorem for W 1,γ -solutions in R n .