The effect of hydrodynamic mixing in bacterial populations due to bacterial chemotaxis is a well-described phenomenon known as bioconvection. Here we report the observation of buoyant plumes that result in hydrodynamic mixing, but in contrast to bioconvection the plumes form in the absence of bacterial motility. We propose that the buoyant flow originates from solute gradients created by bacterial metabolism, similar to solute-induced buoyant flow around growing protein crystals. In our experiments, metabolically-active non-motile Escherichia coli were layered along the bottom of flat-bottomed containers. The E. coli consumed glucose in the medium creating a lighter fluid beneath a heavier fluid. The situation is an example of Rayleigh–Taylor instability, in which a lighter fluid pushes on a heavier one. We developed a numerical model to study the effect of E. coli nutrient consumption and by-product excretion on extracellular solute gradients. The model solutions showed reduced-density fluid along the bottom of the fluid domain leading to buoyant plumes, which were qualitatively similar to the experimental plumes. We also used scaling analyses to study the dependence of plume formation on container size and cell size, and to investigate the effect of reduced gravity, such as the microgravity conditions encountered during spaceflight.