Single-nucleus RNA-sequencing reveals NRF1/NFE2L1 as a key factor determining the thermogenesis and cellular heterogeneity and dynamics of brown adipose tissues in mice