Abstract The main objective of this article is to investigate sub-grid scale turbulence–radiation interaction (SGS TRI) effects on SGS turbulence kinetic energy (TKE) fluctuations and varying thermophysical properties in a partially premixed combustion system for a laboratory-piloted methane/air flame. The large eddy simulation approach is employed to simulate the turbulence of the compressible reactive flow. SGS quantities, including turbulent stress and fluxes of enthalpy and species in the sub-grid scale, are computed using the standard Smagorinsky–Lilly model. The radiative transfer equation is modeled by applying the spherical harmonic P1 approximation by considering the radiative heat source related to the SGS TRI contribution. Optically thin fluctuation approximation is utilized to simplify the radiative absorption term. A chemical reaction mechanism comprising 41 steps and 16 species is applied to model methane–air mixture combustion. Diffusion flamelet-generated manifolds are employed to govern the species transport equation. About 87% of TKE is resolved by applying the finest grid consisting of 1,822,580 cells. Impacts of SGS TRI on the spatially filtered density, eddy viscosity, SGS velocity and TKE, overall radiative emission, RMS temperature fluctuations, and nitrogen monoxide (NO) formation are studied. The results reveal that considering SGS TRI in the simulation leads to remarkable discrepancies, particularly in SGS velocity and TKE by 6.70% and 7.40%, respectively. Meanwhile, SGS density and eddy viscosity deviate negligibly in the presence of SGS TRI. Also, the filtered mass fraction of NO reduces up to 17.54% on average by considering TRI.