Bioconvection observed in a culture of motile micro-organisms was analyzed numerically. The governing equations are the Navier-Stokes equations with the Boussinesq approximation and a diffusion equation for the motile micro-organism. A transition from a static condition to periodic oscillation was observed according to the increase of the Rayleigh number. It was found that the system of bioconvection could be led into chaotic conditions via a single-frequency oscillatory behavior to a sequence of period-doubling bifurcations by increasing the Rayleigh number, which is analogous to Bénard convection.