Identification of functional features underlying heat stress response in Sprague–Dawley rats using mixed linear models