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 nondimensionality. 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. Natural convection has an outstanding applications in daily life. These applications are exist in petrochemical processes, cooling of electronic components, geothermal engineering, crystal growth processes, in the annular gap between the rotor and stator, thermal insulation system, food industry, growth of single silicon crystals, packed bed chemical reactors, grain storage installations, rotating systems, porous heat exchangers, fuel cells, solar ponds etc. Researchers paid extensive attention to work on convection. Venkatachalappa et al. 1 performed a study to analyze the role of rotation on the axisymmetric gravity driven complex flow in a cylindrical annulus whose side walls rotate about their axis with different angular velocities. They obtained the results for Grashof number, rotational speeds, Prandtl number, aspect ratio and compared the results with the existing data. Khan et al. 2 treated the movement in heating prevailing system of a differential type dispersion on an expanding medium using series solution. Sankar et al. 3 tested numerically the hydromagnetic field influence in axial or radial forms for natural convection of a low Prandtl number electrically conducting fluid in a vertical cylindrical annulus. Their outcomes showed that in shallow cavities the flow and heat transfer were suppressed sufficiently through an axial magnetic field and in tall cavities the radial magnetic field had an excellent output. Khan et al. 4 tested the thermal disorder, heat and mass transfer tiny dispersion movement with gyrotactic microorganisms in porous medium using heating wall information. Using the Brinkman-extended Darcy equation, Sankar et al. 5 investigated the natural convection flows in a vertical annulus filled with a fluid-saturated porous medium in which the inner wall was subjected to discrete heating, outer wall was subjected to isothermally at lower temperature and the adiabatic parts were the bottom and top walls including