We establish the boundedness of solutions of reaction-diffusion systems with qua-dratic (in fact slightly super-quadratic) reaction terms that satisfy a natural entropy dissipation property, in any space dimension N > 2. This bound imply the C ∞-regularity of the solutions. This result extends the theory which was restricted to the two-dimensional case. The proof heavily uses De Giorgi's iteration scheme, which allows us to obtain local estimates. The arguments relies on duality reasonings in order to obtain new estimates on the total mass of the system, both in L (N +1)/N norm and in a suitable weak norm. The latter uses C α regularization properties for parabolic equations.