In this work, Thermal Large Eddy Simulation (TLES) is performed to study the behavior of weakly compressible Newtonian fluids with anisotropic temperature-dependent viscosity in forced convection turbulent flow. A systematic analysis of variable-viscosity effects, isolated from gravity, with relevance to industrial cooling/heating applications is being carried out. A LES of a planar channel flow with significant heat transfer at a low Mach number was performed to study effects of fluid property variation on the near-wall turbulence structure. In this flow configuration the top wall is maintained at a higher temperature (Thot) than the bottom wall (Tcold). The temperature ratio (Rθ = Thot/Tcold) is fixed at 1.01, 2 and 3 to study the effects of property variations at low Mach number. Results indicate that average and turbulent fields undergo significant changes. Compared with isothermal flow with constant viscosity, we observe that turbulence is enhanced in the cold side of the channel, characterized by locally lower viscosity whereas a decrease of turbulent kinetic energy is found at the hot wall. The turbulent structures near the cold wall are very short and densely populated vortices but near the hot wall there seems to be a long streaky structure or large elongated vortices. Spectral study reveals that turbulence is completely suppressed at the hot side of the channel at a large temperature ratio because no inertial zone is obtained (i.e. index of Kolmogorov scaling law is zero) from the spectra in these region.