The aim of this paper is to establish regularity for weak solutions to the nondiagonal quasilinear degenerate elliptic systems related to Hörmander's vector fields, where the coefficients are bounded with vanishing mean oscillation. We first prove L p (p ≥ 2) estimates for gradients of weak solutions by using a priori estimates and a known reverse Hölder inequality, and consider regularity to the corresponding nondiagonal homogeneous degenerate elliptic systems. Then we get higher Morrey and Campanato estimates for gradients of weak solutions to original systems and Hölder estimates for weak solutions.Regularity for solutions to elliptic systems in Euclidean spaces has been studied by many authors and a lot of important conclusions were got. Campanato in [2] obtained gradient estimates for weak solutions to linear elliptic system with discontinuous coefficients.For related articles, we quote [1,14] and the references therein. Huang in [18] derived Morrey estimates for uniformly elliptic systems by applying Campanato's technique. Zheng and Feng in [28] established Hölder estimates for weak solutions to quasilinear elliptic systems by reverse Hölder inequality and Dirichlet growth theorem, where the coefficients belong to L ∞ ∩ V MO, and low terms satisfy controlled