The swelling and shrinkage of biological tissues are modelled by a four-component mixture theory [J.M. Huyghe and J.D. Janssen, Int. J. Engng. Sci.35 (1997) 793–802; K. Malakpoor, E.F. Kaasschieter and J.M. Huyghe, Mathematical modelling and numerical solution of swelling of cartilaginous tissues. Part I: Modeling of incompressible charged porous media.ESAIM: M2AN41 (2007) 661–678]. This theory results in a coupled system of nonlinear parabolic differential equations together with an algebraic constraint for electroneutrality. In this model, it is desirable to obtain accurate approximations of the fluid flow and ions flow. Such accurate approximations can be determined by the mixed finite element method. The solid displacement, fluid and ions flow and electro-chemical potentials are taken as degrees of freedom. In this article the lowest-order mixed method is discussed. This results into a first-order nonlinear algebraic equation with an indefinite coefficient matrix. The hybridization technique is then used to reduce the list of degrees of freedom and to speed up the numerical computation. The mixed hybrid finite element method is then validated for small deformations using the analytical solutions for one-dimensional confined consolidation and swelling. Two-dimensional results are shown in a swelling cylindrical hydrogel sample.