We consider the quasilinear parabolic–parabolic Keller–Segel system { u t = ∇ ⋅ ( D ( u ) ∇ u ) − ∇ ⋅ ( S ( u ) ∇ v ) , x ∈ Ω , t > 0 , v t = Δ v − v + u , x ∈ Ω , t > 0 , under homogeneous Neumann boundary conditions in a convex smooth bounded domain Ω ⊂ R n with n ⩾ 1 . It is proved that if S ( u ) D ( u ) ⩽ c u α with α < 2 n and some constant c > 0 for all u > 1 , then the classical solutions to the above system are uniformly-in-time bounded, provided that D ( u ) satisfies some technical conditions such as algebraic upper and lower growth (resp. decay) estimates as u → ∞ . This boundedness result is optimal according to a recent result by the second author (Winkler, 2010 [27]), which says that if S ( u ) D ( u ) ⩾ c u α for u > 1 with c > 0 and some α > 2 n , n ⩾ 2 , then for each mass M > 0 there exist blow-up solutions with mass ∫ Ω u 0 = M . In addition, this paper also proves a general boundedness result for quasilinear non-uniformly parabolic equations by modifying the iterative technique of Moser–Alikakos (Alikakos, 1979 [1]).