We consider a quasilinear parabolic system which has the structure of Patlak-Keller-Segel model of chemotaxis and contains a class of models with degenerate diffusion. A cell population is described in terms of volume fraction or density. In the latter case, it is assumed that there is a threshold value which the density of cells cannot exceed. Existence and uniqueness of solutions to the corresponding initial-boundary value problem and existence of space inhomogeneous stationary solutions are discussed. In the 1D case a classification of stationary solutions for some model example is provided.