Comprehensive RNA-Seq Profiling Reveals Temporal and Tissue-Specific Changes in Gene Expression in Sprague–Dawley Rats as Response to Heat Stress Challenges
Open Access
- 9 April 2021
- journal article
- research article
- Published by Frontiers Media SA in Frontiers in Genetics
- Vol. 12
Abstract
Understanding heat stress physiology and identifying reliable biomarkers are paramount for developing effective management and mitigation strategies. However, little is known about the molecular mechanisms underlying thermal tolerance in animals. In an experimental model of Sprague–Dawley rats subjected to temperatures of 22 ± 1°C (control group; CT) and 42°C for 30 min (H30), 60 min (H60), and 120 min (H120), RNA-sequencing (RNA-Seq) assays were performed for blood (CT and H120), liver (CT, H30, H60, and H120), and adrenal glands (CT, H30, H60, and H120). A total of 53, 1,310, and 1,501 differentially expressed genes (DEGs) were significantly identified in the blood (P< 0.05 and |fold change (FC)| >2), liver (P< 0.01, false discovery rate (FDR)–adjustedP= 0.05 and |FC| >2) and adrenal glands (P< 0.01, FDR-adjustedP= 0.05 and |FC| >2), respectively. Of these, four DEGs, namelyJunb, P4ha1, Chordc1, andRT1-Bb, were shared among the three tissues in CT vs. H120 comparison. Functional enrichment analyses of the DEGs identified in the blood (CT vs. H120) revealed 12 biological processes (BPs) and 25 metabolic pathways significantly enriched (FDR = 0.05). In the liver, 133 BPs and three metabolic pathways were significantly detected by comparing CT vs. H30, H60, and H120. Furthermore, 237 BPs were significantly (FDR = 0.05) enriched in the adrenal glands, and no shared metabolic pathways were detected among the different heat-stressed groups of rats. Five and four expression patterns (P< 0.05) were uncovered by 73 and 91 shared DEGs in the liver and adrenal glands, respectively, over the different comparisons. Among these, 69 and 73 genes, respectively, were proposed as candidates for regulating heat stress response in rats. Finally, together with genome-wide association study (GWAS) results in cattle and phenome-wide association studies (PheWAS) analysis in humans, five genes (Slco1b2, Clu, Arntl, Fads1, andNpas2) were considered as being associated with heat stress response across mammal species. The datasets and findings of this study will contribute to a better understanding of heat stress response in mammals and to the development of effective approaches to mitigate heat stress response in livestock through breeding.Keywords
This publication has 104 references indexed in Scilit:
- Role of the circadian clock gene Per2 in adaptation to cold temperatureMolecular Metabolism, 2013
- Segregation of Regulatory Polymorphisms with Effects on the Gluteus Medius Transcriptome in a Purebred Pig PopulationPLOS ONE, 2012
- Synchronization of Circadian Per2 Rhythms and HSF1-BMAL1:CLOCK Interaction in Mouse Fibroblasts after Short-Term Heat Shock PulsePLOS ONE, 2011
- Phylogenetic-based propagation of functional annotations within the Gene Ontology consortiumBriefings in Bioinformatics, 2011
- HSFs and regulation of Hsp70.1 (Hspa1b) in oocytes and preimplantation embryos: new insights brought by transgenic and knockout mouse modelsCell Stress and Chaperones, 2011
- Intracisternal Administration of NR2 Subunit Antagonists Attenuates the Nociceptive Behavior and p-p38 MAPK Expression Produced by Compression of the Trigeminal Nerve RootMolecular Pain, 2011
- Endoplasmic reticulum stress inhibits collagen synthesis independent of collagen-modifying enzymes in different chondrocyte populations and dermal fibroblastsBiochemistry and Cell Biology, 2010
- BEDTools: a flexible suite of utilities for comparing genomic featuresBioinformatics, 2010
- edgeR: a Bioconductor package for differential expression analysis of digital gene expression dataBioinformatics, 2009
- Mapping and quantifying mammalian transcriptomes by RNA-SeqNature Methods, 2008