This paper is devoted to the justification of the macroscopic, mean-field nutrient taxis system with doubly degenerate cross-diffusion proposed by Leyva et al. [39] to model the complex spatio-temporal dynamics exhibited by the bacterium B. subtilis during experiments run in vitro. This justification is based on a microscopic description of the movement of individual cells whose changes in velocity (in both speed and orientation) obey a velocity jump process [47], governed by a transport equation of Boltzmann type. For that purpose, the asymptotic method introduced by Hillen and Othmer [27,48] is applied, which consists of the computation of the leading order term in a regular Hilbert expansion for the solution to the transport equation, under an appropriate parabolic scaling and a first order perturbation of the turning rate of Schnitzer type [54]. The resulting parabolic limit equation at leading order for the bacterial cell density recovers the degenerate nonlinear cross diffusion term and the associated chemotactic drift appearing in the original system of equations. Although the bacterium B. subtilis is used as a prototype, the method and results apply in more generality.2010 Mathematics Subject Classification. 35K65, 92C17, 60J75.