Biological transformation of volatile organic compounds is one of the key factors that influence contaminant-plume evolution and thus natural attenuation. In this study we investigate the effect of biological transformation on the transport of contaminants in the aqueous and gaseous phases. The analysis includes the study of the effect of density-driven advection of contaminants in the gaseous phase on multiphase and multispecies flow, fate and transport modeling in the subsurface. Trichloroethylene (TCE) and its two byproducts, dichloroethylene and vinyl chloride, are analyzed as the target contaminants. Our results indicate that density-driven advection of the gaseous phase, which is initiated by evaporation of TCE as a nonaqueous phase liquid, increases the downward and also the lateral migration of TCE within the unsaturated zone. This process also influences the location of high-concentration zones of the byproducts of TCE in the unsaturated and the saturated zones. Biotransformation of TCE contributes to the reduction of dissolved TCE plume development as expected. The daughter byproducts, which are introduced into the subsurface system, show distinct transport patterns as they are affected by their independent degradation kinetics and density-driven advection. These observations, which are based on our simulation results for biotransformation and transport of TCE and its byproducts, are useful in evaluating the natural attenuation processes, its potential health hazards and also the evaluation of potential plume development at contaminated sites.