An anisotropic subgrid-scale parameterization for large-eddy simulations of stratified turbulence