The molecular mechanisms underlying heat stress tolerance in animals to high temperatures remain unclear. This experiment was performed with Sprague-Dawley rats housed at 22 ℃ (control group; CT), and 42 ℃ for 30 min (H30), 60 min (H60), and 120 min (H120). Transcriptomic analysis using RNA-Sequencing (RNA-Seq) from blood (CT and H120), liver (CT, H30, H60, and H120), and adrenal glands (CT, H30, H60, and H120) was performed. The differentially expressed mRNA isoforms (DEIs) were identified and annotated by the CLC Genomics Workbench. Biological process and metabolic pathway analyses were performed using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) database. A total of 225, 5,764, and 4,988 DEIs in the blood, liver, and adrenal glands were observed. Furthermore, the number of novel differentially expressed transcript lengths with annotated genes and the novel differentially expressed transcript with non-annotated genes were 136 and 8 in blood, 3,549 and 120 in the liver, as well as 3,078 and 220 in adrenal glands, respectively. A total of 35 genes were involved in the heat stress response, out of which Dnaja1, LOC680121, Chordc1, AABR07011951.1, Hsp90aa1, Hspa1b, Cdkn1a, Hmox1, Bag3 and Dnaja4 were commonly identified in the liver and adrenal glands, suggesting that these genes may regulate heat stress response through interactions between liver and adrenal glands. The results suggest that the identified mRNA isoforms could be considered as potential candidates for selecting mammals to improve thermotolerance.