Buongiorno's nanofluid model is followed to study the bioconvection in two stretchable rotating disks with entropy generation. Similarity transformations are used to handle the problem equations for non-dimensionality. For the simulation of the modeled equations, Homotopy Analysis Method is applied. The biothermal system is explored for all the embedded parameters whose effects are shown through different graphs. There exists interesting results due to the effects of different parameters on different profiles. Radial velocity decreases with increasing stretching and magnetic field parameters. Temperature increases with Brownian motion and thermophoresis parameters. Nanoparticles concentration decreases on increasing Lewis number and thermophoresis parameter while motile gyrotactic microorganisms profile increases with increasing Lewis and Peclet numbers. Convergence of the solution is found and good agreement is obtained when the results are compared with published work.