Modern biomedical and tribological systems are increasingly deploying combinations of nanofluids and bioconvecting microorganisms which enable improved control of thermal management. Motivated by these developments, in this study, a new mathematical model is developed for the combined nanofluid bioconvection axisymmetric squeezing flow between rotating circular plates (an important configuration in, for example, rotating bioreactors and lubrication systems). The Buongiorno two-component nanoscale model is deployed, and swimming gyrotactic microorganisms are considered which do not interact with the nanoparticles. Thermal radiation is also included, and a Rosseland diffusion flux approximation is utilized. Appropriate similarity transformations are implemented to transform the nonlinear, coupled partial differential conservation equations for mass, momentum, energy, nanoparticle species and motile microorganism species under suitable boundary conditions from a cylindrical coordinate system into a dimensionless nonlinear ordinary differential boundary value problem. An efficient scheme known as differential transform method (DTM) combined with Padé-approximations is then applied to solve the emerging nonlinear similarity equations. The impact of different non-dimensional parameters i.e. squeezing Reynolds number, rotational Reynolds number, Prandtl number, thermophoresis parameter, Brownian dynamics parameter, thermal radiation parameter, Schmidt number, bioconvection number and Péclet number on velocity, temperature, nanoparticle concentration and motile gyrotactic microorganism density number distributions is computed and visualized graphically. The torque effects on both plates, i.e. the lower and the upper plate, are also determined. From the graphical results, it is seen that momentum in the squeezing regime is suppressed clearly as the upper disk approaches the lower disk. This inhibits the axial flow and produces axial flow retardation. Similarly, by enhancing the value of squeezing Reynolds number, the tangential velocity distribution also decreases. More rigorous squeezing clearly therefore also inhibits tangential momentum development in the regime and leads to tangential flow deceleration. Tables are also provided for multiple values of flow parameters. The numerical values obtained by DTM-Padé computation show very good agreement with shooting quadrature. DTM-Padé is shown to be a precise and stable semi-numerical methodology for studying rotating multi-physical flow problems. Radiative heat transfer has an important influence on the transport characteristics. When radiation is neglected, different results are obtained. It is important therefore to include radiative flux in models of rotating bioreactors and squeezing lubrication dual disk damper technologies since high temperatures associated with radiative flux can impact significantly on combined nanofluid bioconvection which enables more accurate prediction of actual thermofluidic characteristics. Corrosion and surface degradation effects may therefore be mitigated in designs.