Next Article in Journal
Improving Terrestrial Squamate Surveys with Camera-Trap Programming and Hardware Modifications
Previous Article in Journal
Characterization of the Normal Portal and Hepatic Blood Flow of Adult Holstein-Friesian Cows
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Animals
Volume 9
Issue 6
10.3390/ani9060387
animals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptomic Analysis Provides Novel Insights into Heat Stress Responses in Sheep

by
Zengkui Lu
1,2,
Mingxing Chu
1,
Qing Li
1,
Meilin Jin
1,
Xiaojuan Fei
1,
Lin Ma
1,
Liping Zhang
2 and
Caihong Wei
1,*
1
Key Laboratory of Animal Genetics and Breeding and Reproduction of Ministry of Agriculture, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing 100193, China
2
College of Animal Science and Technology, Gansu Agricultural University, Lanzhou 730070, China
*
Author to whom correspondence should be addressed.
Animals 2019, 9(6), 387; https://doi.org/10.3390/ani9060387
Submission received: 18 April 2019 / Revised: 7 June 2019 / Accepted: 17 June 2019 / Published: 24 June 2019
(This article belongs to the Section Small Ruminants)
Browse Figures
Review Reports Versions Notes

Abstract

:

Simple Summary

The general increase in global temperatures has meant that heat stress has become an increasingly significant problem for sheep. This has both direct and indirect impact on their physiological functions, productivity, and health of sheep. Sheep generally live in high-temperature environments; however, the genes and pathways that play regulatory roles in the heat stress responses of sheep remain unclear. In this study, we applied RNA-Seq technology to analyze liver tissues of sheep from heat-stressed and control groups, and screened genes and pathways related to sheep heat stress. This work provides a theoretical foundation for the breeding and production of heat-resistant sheep.

Abstract

With the intensified and large-scale development of sheep husbandry and global warming, sheep heat stress has become an increasingly important issue. However, little is known about the molecular mechanisms related to sheep responses to heat stress. In this study, transcriptomic analysis of liver tissues of sheep in the presence and absence of heat stress was conducted, with the goal of identifying genes and pathways related to regulation when under such stress. After a comparison with the sheep reference genome, 440,226,436 clean reads were obtained from eight libraries. A p-value ≤ 0.05 and fold change ≥ 2 were taken as thresholds for categorizing differentially expressed genes, of which 1137 were identified. The accuracy and reliability of the RNA-Seq results were confirmed by qRT-PCR. The identified differentially expressed genes were significantly associated with 419 GO terms and 51 KEGG pathways, which suggested their participation in biological processes such as response to stress, immunoreaction, and fat metabolism. This study’s results provide a comprehensive overview of sheep heat stress-induced transcriptional expression patterns, laying a foundation for further analysis of the molecular mechanisms of sheep heat stress.

1. Introduction

Since its foundation in 1988, the Intergovernmental Panel on Climate Change has released five comprehensive climate reports and highlighted the indisputable fact that the world is warming. In the wild, air temperature has a major bearing on the ability of animals to thrive. Long-term high temperatures and excessively warm weather generate serious heat stress for animals. Heat stress directly or indirectly affects animal physiology, production, and health, among others [1,2,3], which is manifested by physiological features such as gasping for breath, accelerated respiration, elevated heart rate, rising rectal temperature, increasing in-vivo oxidative metabolism, and water and electrolyte disturbance [4,5,6,7]. Regarding reproduction, it can also cause reduced food intake, weight loss, and declining fertility [8,9,10], as well as declining immunity and disease resistance [11,12].
Transcriptomics refers to the study of all of the transcripts in cells, tissues, or living organisms at a specific developmental stage or under specific physiological conditions [13]. In contrast to the relatively stable genome, the transcriptome varies among developmental stages, physiological conditions, and upon exposure to external stimuli. The transcriptome, which has become a powerful tool for analyzing the relationship between genotype and phenotype, has been extensively applied to the analysis of complicated biological events.
To date, studies on heat stress in animals have mainly focused on dairy cow, broiler, and pig. These studies have mainly concentrated on features such as physiology and biochemistry, production performance, meat quality, and immunological functions. However, these studies were generally restricted to individual genes and pathways and could not characterize the complicated processes of animal responses to heat stress holistically. Moreover, very few studies have reported on heat stress in sheep. There is thus a need for analyses of the genes and pathways that exert regulatory effects when sheep are under heat stress. Against this background, Hu sheep were selected as experimental samples in this study, and subjected to RNA-Seq, bioinformatic, and molecular biology analyses to explore the molecular mechanisms of sheep heat stress, to seek genes and regulatory pathways related to sheep heat stress, and to provide a theoretical basis for sheep production and breeding practices that will help sheep producers deal with heat stress in their sheep.

2. Materials and Methods

2.1. Animals and Sample Collection

All the experimental procedures mentioned in the present study were approved by the Science Research Department (in charge of animal welfare issue) of the Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (IAS-CAAS) (Beijing, China). Ethical approval on animal survival was given by the animal ethics committee of IAS-CAAS (No. IASCAAS-AE-03, 12 December 2016). Thirty-six purebred 5-month-old healthy Hu sheep (half females and half males with approximately equivalent weights. At this age, sheep of this breed have normally achieved puberty) were selected in August (heat-stressed group) and December 2017 (no heat stress, control group) and bred at Jiangsu Qianbao Sheep Industry Co., Ltd. (Yancheng, Jiangsu, China). Sheep in the heat-stressed group and the control group were fed the same food and had free access to water. Five temperature and humidity meters were hung 1.5 m above the pen, and temperature (Td) and relative humidity (RH) were recorded every 2 h from 08:00 to 20:00. Temperature-humidity index (THI) was used to evaluate the severity of heat stress [14]. When THI was continuously higher than 23.3 (severe heat stress) or continuously lower than 22.2 (no heat stress) for 7 d, four Hu sheep (half females and half males; HG1 and HG2 were heat-stressed males and HM1 and HM2 were heat-stressed females, while CG1 and CG2 were males in the control group and CM1 and CM2 were females in the control group) were randomly selected and slaughtered on an empty stomach the morning (sheep were fasted for 12 h before slaughter, but had free access to water), after which complete liver tissues were separated out, and some left hepatic lobe tissues were preserved in liquid nitrogen for RNA extraction.

2.2. Detection of Thyroid-Related Hormones in Sheep Serum

Blood was collected from the jugular vein on an empty stomach on the day of slaughtering the sheep. Five sheep were selected for this at random in heat stress and control groups, for which non-anticoagulant vacutainer tubes were used to collect 5 ml blood. Serum was then separated from it by centrifugation for 15 min at 3000 r/min. In accordance with the manufacturer’s instructions, radioimmunoassay was used to determine the T3 {Iodine [125I] 3,3’,5-Triiodothyronine(T3) Radioimmunoassay Kit, Beifang, Beijing, China} and T4 {Iodine [125I] Thyroxine Radioimmunoassay Kit, Beifang, Beijing, China} concentrations in sheep serum.

2.3. RNA Extraction and Quality Inspection

Total RNA in the liver tissues of the eight sheep was extracted using TRlzol reagent (Invitrogen, Carlsbad, CA, USA), in accordance with the manufacturer’s instructions. The integrity and concentration of RNA were determined by agarose gel electrophoresis and a Nanodrop 2000 (Thermo, Waltham, MA, USA). Finally, the RNA concentration was accurately quantified using Qubit 2.0 (Invitrogen, Carlsbad, CA, USA) and stored at −80 °C until subsequent analysis.

2.4. cDNA Library Construction and Sequencing

RNA samples that qualified through quality inspection were used for library construction in accordance with the instructions provided by Illumina (Illumina, San Diego, CA, USA). The specific steps were as follows: (1) Oligo-carrying magnetic beads were used for mRNA enrichment. (2) Fragmentation buffer was used to randomly break enriched mRNA into fragments of about 200 nucleotides in length. (3) The fragmented mRNA was used as a template; one-chain cDNA was synthesized through inverse transcription using a random primer. In the synthesis of second-chain cDNA, dTTP in dNTPs was replaced by dUTP. (4) AMPure XP beads were used to purify double-chain cDNA, and End Repair Mix was used for end filling and the addition of an A-tail and sequencing linkers. (5) USER enzyme was used to digest double-stranded cDNA so that the library only contained single-stranded cDNA. (6) PCR enrichment of the cDNA fulfilling the above criteria was conducted. (7) Qubit 2.0 and Agilent 2010 (Agilent, Palo Alto, CA, USA) were used for quality inspection of the concentration and fragment size of the established library. The Illumina Hiseq 2500 platform was used to conduct double-end sequencing of the eight established libraries.

2.5. Raw Data Preprocessing and Alignment

Raw reads were obtained after Illumina sequencing ended. Skewer software was used to dynamically remove the 3’ ends, linker sequences, and low-mass sequences of raw reads, after which clean reads were obtained. FastQC software was used to conduct quality control analysis of preprocessed data and determine indexes such as GC content, Q20, and Q30. STAR software was used to align obtained clean reads onto reference genomes of sheep (Oar_v4.0), where alignment parameters were set as —twopassMode Basic, —outSAMstrandField intronMotif, —alignSJstitchMismatchNmax 5, −1, 5, 5, and other parameters were set as the defaults. Those reads aligned to more than one (>1) genome were removed and the others were kept for subsequent analysis.

2.6. Identification and Analysis of Differentially Expressed Genes

StringTie software was used for counting the original sequences of known genes of the eight samples. Expression levels of known genes were expressed using fragments per kilobase of transcript per million fragments mapped (FPKM), where FPKM = Total fragments/Mapped reads (millions) × Exon length (KB). The genes were divided into eight intervals [(0.1–1), (1–5), (5–10), (10–20), (20–30), (30–40), (40–50), and (≥50)] according to their FPKM values, and the number of transcripts in each interval and its percentage relative to the total number of transcripts were calculated. It is generally considered that FPKM ≥ 0.1 indicates that the transcript is expressed. DESeq2 software was used to screen out genes that were differentially expressed between sheep in the heat-stressed group and those in the control group; in the screening, the thresholds of p-value ≤ 0.05 and fold change ≥ 2 were applied. Clustering analysis of the screened differentially expressed genes was subsequently carried out.

2.7. Functional Enrichment Analysis of Differentially Expressed Genes

The DAVID database was used in Gene Ontology (GO) functional enrichment analysis of the differentially expressed genes. All genes were placed on the background list; the differentially expressed genes were included in the candidate list, and hypergeometric distribution was used to calculate p values before multiple testing and Benjamini–Hochberg correction. GO items with a false discovery rate ≤ 0.05 were regarded as significantly enriched. The Kyoto Encyclopedia of Genes and Genomes (KEGG) database was used in the functional annotation and classification of genes associated with differential peaks.

2.8. Verification of Differentially Expressed Genes

The RNA-Seq results were verified using qRT-PCR. A cDNA Synthesis Kit (TransStart Green, Beijing, China) was used for RNA reverse transcription. Oligo 7 software was used for cross-exon design of all primers and specific detection in NCBI. The qPCR SuperMix (TransStart Green, Beijing, China) and LightCycler 480 apparatus (Roche, Basel, Switzerland) were used for quantitative real-time PCR (qRT-PCR). Each sample was measured in triplicate to ensure the accuracy of the quantification. The relative expression of target genes was calculated by the 2-ΔΔCt method using β-actin (ACTB) as an internal reference gene. Information on the primers used in this study is listed in Supplementary Table S1.

2.9. Statistical Analysis

All data reported are expressed as mean ± SE. Student’s t-test was carried out using SPSS software for statistical analysis of the data. A p value of < 0.05 was considered to be statistically significant.

3. Results

3.1. Influence of Heat Stress on Serum T3 and T4 Levels of Sheep

During heat stress, THI was ≥ 26.18, indicating that the sheep were exposed to extremely severe heat stress; in the control group, THI ≤ 11.82, indicating that the sheep experienced no heat stress (Supplementary Figure S1) [14]. T3 and T4 act in regulating the body’s metabolic heat production. In this study, the levels of T3 and T4 in the serum of heat-stressed and control sheep were determined and they were found to be extremely highly expressed in the control group sheep (Table 1).

3.2. Summary of Sheep Pen THI and Sequencing Data

The Illumina Hiseq 2500 platform was used to conduct paired-end sequencing of eight established libraries; the results of data quality control are shown in Table 2. A total of 513,487,656 raw reads were generated in the eight libraries, 440,226,436 clean reads were left after quality control, and the average proportion of clean reads was 99.06%. The GC contents of the eight samples were within 49.82%–51.25%, and were thus in accordance with base composition rules. Q20 ≥ 98.10% and Q30 ≥ 94.40%. Clean reads were aligned to the sheep reference genome using STAR software, and over 93.25% of the reads could be accurately aligned with a high matching rate. Approximately 7.05%–12.37% of these clean reads had multiple aligned positions, 81.99%–86.73% of them had single aligned positions, and 3.66%–6.76% of them were not aligned to the sheep reference genome. Those reads that aligned on the sheep reference genome at a single position were used for further bioinformatic analysis.

3.3. Analysis of Overall Gene Expression Levels

Through a comparison of the FPKM values of all genes, the overall gene expression levels of each sample could be evaluated. As shown in Figure 1, the overall gene expression levels of the eight samples were similar. To further evaluate the overall gene expression levels, the genes were divided into eight intervals (0.1–1, 1–5, 5–10, 10–20, 20–30, 30–40, 40–50, and ≥50) according to their FPKM values, and the number of transcripts in each interval and its percentage relative to the total number of transcripts were calculated. There were small percentage changes of various expression quantities in the eight samples, and the number of transcripts gradually decreased as the FPKM value increased (Supplementary Table S2). Overall, most transcripts were expressed at low levels in sheep liver in the heat-stressed group and control group (Supplementary Table S2).

3.4. Screening and Clustering Analysis of Differentially Expressed Genes

DESeq2 software was used to screen out sheep genes differentially expressed between the heat-stressed group and the control group; 1,137 genes fulfilling both criteria of p-value ≤ 0.05 and fold change ≥ 2 were identified (Figure 2A and Supplementary Table S3). Compared with the levels in the control group,
798 genes in sheep in the heat-stressed group presented upregulated expression, while 339 genes were downregulated (Figure 2B). To obtain a deeper understanding of the gene expression patterns of sheep in the heat-stressed group and control group, clustering analysis of these differentially expressed genes was conducted. The results indicated that the differentially expressed genes of sheep in the heat-stressed group and control group were clustered into a single class (Figure 2C). In addition, genes with upregulated differential expression and those with downregulated differential expression were clustered into one class. We also analyzed the effects of heat stress on gene expression in sheep of different genders. 4,203 genes were found to be differentially expressed between heat-stressed rams and rams from the control group (Supplementary Table S4), while 1,369 were differentially expressed between heat-stressed ewes and control group ewes (Supplementary Table S5).

3.5. Functional Enrichment Analysis of Differentially Expressed Genes

Sheep genes that were differentially expressed between the heat-stressed group and the control group, were significantly enriched in 419 GO terms, including 357 biological processes, 27 cellular components, and 35 molecular functions (Supplementary Table S6). These GO terms are associated with stress response, immunoreaction, and fat metabolism, among others (Figure 3A). KEGG pathway analysis found that these differentially expressed genes were significantly enriched in 51 pathways (Supplementary Table S7). Through a further analysis, many of these pathways were shown to be related to fat metabolism, such as regulation of lipolysis in adipocytes, MAPK (mitogen-activated protein kinase), and PI3K-AKt (phosphatidylinositol 3-kinase and protein kinase B) (Figure 3B). They were also particularly associated with signaling pathways regulating stress reactions, including TNF (tumor necrosis factor) and Rap1 (Figure 3B). Besides, genes differentially expressed in heat-stressed rams and control group rams were particularly associated with the process of fat metabolism and stress response (Supplementary Tables S8 and S9). The results for genes differentially expressed between heat-stressed ewes and control group ewes were also particularly associated with the process of fat metabolism and stress response (Supplementary Tables S10 and S11).

3.6. Validation of RNA-Seq Data by qRT-PCR

To verify the sequencing results of RNA-Seq, five genes with high expression and five with low expression in the heat-stressed group in this study were selected for confirmatory analysis by qRT-PCR. As shown in Figure 4, although their expression levels differed to a certain degree, their expression patterns were identical, thus verifying that the RNA-Seq results were accurate and reliable.

4. Discussion

Heat stress impacts animal production and endangers animal welfare, causing an annual economic loss of over US$1.2 billion in animal husbandry globally [15]. In conventional breeding methods, the resistance of animals to heat can only be measured by determining some key production indexes, but the genetic changes of animals in association with environmental changes cannot be determined. In contrast, high-throughput sequencing can be used to shed light on animal heat resistance and genes associated with sensitivity to heat, and it has been proven to be a feasible approach in studies of animals such as cow, pig, and chicken [16,17,18]. Therefore, starting from transcriptomics, this study was aimed at identifying sheep heat resistance-related genes and pathways. Four sheep under heat stress and four without such stress were selected in this study, from which 1,137 differentially expressed genes were screened out using RNA-Seq. qRT-PCR results revealed expression patterns of 10 genes that were identical to those in the RNA-Seq results, proving that the sequencing results were accurate and reliable. These differentially expressed genes were mainly related to stress response, energy metabolism, and immunity. Further functional enrichment analysis of these differentially expressed genes was carried out, which showed that most genes were particularly associated with signaling pathways such as stress response, fat metabolism, and immunity.

4.1. T3 and T4

T3 and T4 are hormones secreted from the thyroid gland and its peripheral tissues. T3 is a type of thyroid hormone that works on anabolism and catabolism. The reduction in secretion of T3 is considered to reflect an animal’s efforts to reduce its body heat production and accumulation as well as to maintain its body heat balance [1,19]. T4 also plays an important role in the adaptation of animals to environmental changes. It stimulates intracellular oxygen consumption and heat production, thereby increasing basal metabolic rate, increasing glucose utilization, and altering lipid metabolism [19,20]. In this study, it was found that the concentrations of T3 and T4 in the blood of heat-stressed sheep were significantly lower than those in the sheep of the control group, which was consistent with the results in studies of other sheep breeds [7]. A decrease in the concentrations of T3 and T4 indicates that heat stress reduces the metabolic activity of the sheep, so as to inhibit the body from producing more heat. In addition, the decrease in thyroxine concentration may also be related to a decrease in blood glucose concentration, which is closely related to energy metabolism during heat stress.

4.2. Regulation of Body Temperature

5-Hydroxytryptamine receptor 4 (HTR4) presented significantly high expression in the heat-stressed sheep in this study, while HTR1B was significantly more highly expressed in sheep in the control group. It has been indicated that 5-HT plays a significant role in the process of regulating animal body temperature [21]. When animals are in a high-temperature environment, the 5-HT neuroendocrine system is activated and combines with 5-HTR to activate downstream cyclic adenosine monophosphate (cAMP) and cyclic guanosine monophosphate (cGMP) signaling pathways, which enhance neural activity sensitive to heat so as to inhibit heat production [22]. Fewer studies have been performed on HTR4 with regard to body temperature regulation, but it has been shown in animals to exert important regulatory effects in gastrointestinal sensitivity, memory, and food intake, among others [23]. Under heat stress, sheep exhibited reduced food intake and rumination activities, which further reduced the amount of heat generated by metabolism so as to maintain a stable body temperature [6,9]. During this period, HTR4 might play a significant role. HTR1B is also involved in body temperature regulation, and was shown to be activated to reduce the body temperature of rats [24]. Here, HTR1B presented high expression in sheep in the control group, suggesting its importance in regulating body temperature.

4.3. Regulation of Stress Reactions

Under heat stress, animals need to undertake a series of nonspecific responses to maintain homeostasis. The functional enrichment analysis of differentially expressed genes in sheep in the heat-stressed group and the control group showed the significant enrichment of genes in stress-related pathways, such as Rap1, MAPK, and PI3K-Akt. Rap1 is regulated by a wide range of external stimuli as a molecular switch when stress occurs [25]. MAPK is one of the main signaling pathways in mammals and can also be activated by external stimuli, but the activity of the MAPK signaling pathway is stimulated by Rap1 in multiple cells [26,27]. In a similar way, the PI3K-Akt pathway is activated by external stimuli of many cellular types and can involve a cascade reaction with the MAPK signaling pathway [28,29]. In summary, these stress-related signaling pathways can form a complicated cascade reaction to respond to heat stress in sheep, and reduce the harm brought about by such stress.

4.4. Regulation of Energy Metabolism

Under heat stress, sheep reduce their food intake in order to reduce metabolic heat production. This study showed that differentially expressed genes in the heat-stressed group were particularly associated with fat metabolism-related pathways. Besides regulating stress reactions, MAPK and PI3K-Akt signaling pathways exert important effects on the fat metabolic process [30,31]. In addition, the cAMP signaling pathway also plays a significant role in the hepatic fat metabolic process and can activate protein kinase A (PKA) to promote β-oxidation of hepatic fat [32]. Genes associated with these signaling pathways were also particularly common among the differentially expressed ones, such as NPR1 (natriuretic peptide receptor 1), ANGPT2 (angiopoietin 2), and SLC13A5 (sodium-dependent citrate transporter member 5); this study indicated that all of these genes participate in the regulation of fat metabolism [33,34,35]. Sheep needed to mobilize a large amount of energy to respond to the heat stress, which enhanced catabolism and reduced anabolism correspondingly. Although genes associated with some glycometabolism-related signaling pathways were also commonly identified in this study, compared with fat, saccharide energy content was low with easy hydration. Therefore, fat mobilization in sheep under heat stress could help them adapt to this condition.

4.5. Regulation of Immunoreactions

Moderate heat stress has the effect of promoting immune system activity in animals so that they can better adapt to the environment. However, severe heat stress has a negative influence on the immune systems of animals and results in deterioration of their disease resistance. IL1R (interleukin type 1 receptor) has two subtypes, IL1R1 and IL1R2, both of which exert important effects on immune and inflammatory reactions [36]. HSPA2 belongs to the HSP70 (heat stress protein 70) family. It has been shown that HSP70 can increase the secretion of IL-1β from cells and facilitate the transport of antigens to T cells [37,38,39]. In this study, it was found that IL1R1, IL1R2, and HSPA2 presented high expression in heat-stressed sheep, suggesting that these genes play significant roles in the immunoregulatory process of sheep.

5. Conclusions

Upon RNA-Seq analysis of eight libraries of heat-stressed and control groups, 1,137 differentially expressed genes were screened out. These genes are mainly involved in biological processes such as stimulus-response reaction, immunoreaction, and fat metabolism. These candidate genes and pathways should deepen our understanding of the molecular mechanisms involved in heat stress in sheep, and provide a theoretical basis for improved sheep production and breeding of heat-tolerant varieties.

Availability of Supporting Data

All the RNA-Sequencing data used in this study has been deposited in the Sequence Read Archive (SRA) public databases under BioProject (PRJNA531367).

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-2615/9/6/387/s1, Figure S1: THI of sheep 1 week before slaughter, Table S1: Primers used in this study for qRT-PCR, Table S2: Number of genes in different FPKM intervals, Table S3: Differentially expressed genes identified in heat-stressed and control sheep groups by RNA-Seq, Table S4: Differentially expressed genes identified in heat-stressed rams and control group rams, Table S5: Differentially expressed genes identified in heat-stressed ewes and control group ewes, Table S6: GO analysis of differentially expressed genes in heat-stressed and control sheep groups, Table S7: KEGG analysis of differentially expressed genes in heat-stressed and control sheep groups. Table S8: GO analysis of differentially expressed genes in heat-stressed rams and control group rams, Table S9: KEGG analysis of differentially expressed genes in heat-stressed rams and control group rams, Table S10: GO analysis of differentially expressed genes in heat-stressed ewes and control group ewes, Table S11: KEGG analysis of differentially expressed genes in heat-stressed ewes and control group ewes.

Author Contributions

Conceptualization, Z.L. and C.W.; formal analysis, Z.L., C.W., and M.C.; investigation, Z.L., M.J., X.F., and L.M.; resources, L.Z. and Q.L.; writing—original draft preparation, Z.L. and C.W.; writing—review and editing, Z.L. and C.W.; project administration, Z.L.

Funding

This research was supported by the National Natural Science Foundation of China (No. 31672380).

Acknowledgments

We thank Jiangsu Qianbao Sheep Industry Co., Ltd. (Yancheng, Jiangsu, China) for providing the experimental animals.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ACTBβ-actin
ANGPT2angiopoietin 2
cAMPcyclic adenosine monophosphate
cGMPcyclic guanosinc monophosphate
FPKMfragments per kilobase of transcript per million fragments mapped
GOgene ontology
HSP70heat stress protein 70
HTR45-Hydroxytryptamine receptor 4
IL1Rinterleukinl type 1 receptor
KEGGkyoto encyclopedia of genes and genomes
MAPKmitogen-activated protein kinase
NPR1natriuretic peptide receptor 1
PI3K-AKtphosphatidylinositol 3-kinase and protein kinase B
PKAprotein kinase A
qRT-PCRquantitative real time PCR
RHrelative humidity
SLC13A5sodium-dependent citrate transporter member 5
Tdtemperature
THItemperature-humidity index

References

  1. Pantoja, M.H.A.; Esteves, S.N.; Jacinto, M.A.C.; Pezzopane, J.R.M.; Paz, C.C.P.; Silva, J.; Lourenco Junior, J.B.; Brandao, F.Z.; Moura, A.B.B.; Romanello, N.; et al. Thermoregulation of male sheep of indigenous or exotic breeds in a tropical environment. J. Therm. Biol. 2017, 69, 302–310. [Google Scholar] [CrossRef] [PubMed]
  2. Vilas Boas Ribeiro, B.P.; Lanferdini, E.; Palencia, J.Y.P.; Lemes, M.A.G.; Teixeira de Abreu, M.L.; de Souza Cantarelli, V.; Ferreira, R.A. Heat negatively affects lactating swine: A meta-analysis. J. Therm. Biol. 2018, 74, 325–330. [Google Scholar] [CrossRef]
  3. Cui, Y.; Hao, Y.; Li, J.; Bao, W.; Li, G.; Gao, Y.; Gu, X. Chronic heat stress induces immune response, oxidative stress response, and apoptosis of finishing pig liver: A proteomic approach. Int. J. Mol. Sci. 2016, 17, 393. [Google Scholar] [CrossRef]
  4. Daramola, J.O.; Adeloye, A.A. Physiological adaptation to the humid tropics with special reference to the West African Dwarf (WAD) goat. Trop. Anim. Health Prod. 2009, 41, 1005–1016. [Google Scholar] [CrossRef] [PubMed]
  5. Marai, I.F.M.; Habeeb, A.A.M.; Gad, A.E. Reproductive traits of female rabbits as affected by heat stress and lighting regime under subtropical conditions of Egypt. Anim. Sci. 2004, 78, 119–127. [Google Scholar] [CrossRef]
  6. Li, F.K.; Yang, Y.; Jenna, K.; Xia, C.H.; Lv, S.J.; Wei, W.H. Effect of heat stress on the behavioral and physiological patterns of Small-tail Han sheep housed indoors. Trop. Anim. Health Prod. 2018, 50, 1893–1901. [Google Scholar] [CrossRef]
  7. Rathwa, S.D.; Vasava, A.A.; Pathan, M.M.; Madhira, S.P.; Patel, Y.G.; Pande, A.M. Effect of season on physiological, biochemical, hormonal, and oxidative stress parameters of indigenous sheep. Vet. World 2017, 10, 650–654. [Google Scholar] [CrossRef] [Green Version]
  8. Das, R.; Sailo, L.; Verma, N.; Bharti, P.; Saikia, J.; Imtiwati; Kumar, R. Impact of heat stress on health and performance of dairy animals: A review. Vet. World 2016, 9, 260–268. [Google Scholar] [CrossRef]
  9. De, K.; Kumar, D.; Saxena, V.K.; Thirumurugan, P.; Naqvi, S.M.K. Effect of high ambient temperature on behavior of sheep under semi-arid tropical environment. Int. J. Biometeorol. 2017, 61, 1269–1277. [Google Scholar] [CrossRef]
  10. Narayan, E.; Sawyer, G.; Parisella, S. Faecal glucocorticoid metabolites and body temperature in Australian merino ewes (Ovis aries) during summer artificial insemination (AI) program. PLoS ONE 2018, 13, e0191961. [Google Scholar] [CrossRef]
  11. Dahl, G.E.; Tao, S.; Monteiro, A.P.A. Effects of late-gestation heat stress on immunity and performance of calves. J. Dairy Sci. 2016, 99, 3193–3198. [Google Scholar] [CrossRef]
  12. Saelao, P.; Wang, Y.; Gallardo, R.A.; Lamont, S.J.; Dekkers, J.M.; Kelly, T.; Zhou, H. Novel insights into the host immune response of chicken Harderian gland tissue during Newcastle disease virus infection and heat treatment. BMC Vet. Res. 2018, 14, 280. [Google Scholar] [CrossRef] [PubMed]
  13. Wang, Z.; Gerstein, M.; Snyder, M. RNA-Seq: A revolutionary tool for transcriptomics. Nat. Rev. Genet. 2009, 10, 57–63. [Google Scholar] [CrossRef] [PubMed]
  14. Marai, I.F.M.; El-Darawany, A.A.; Fadiel, A.; Abdel-Hafez, M.A.M. Physiological traits as affected by heat stress in sheep—A review. Small Rumin. Res. 2007, 71, 1–12. [Google Scholar] [CrossRef]
  15. Baumgard, L.H.; Rhoads, R.P., Jr. Effects of heat stress on postabsorptive metabolism and energetics. Annu. Rev. Anim. Biosci. 2013, 1, 311–337. [Google Scholar] [CrossRef] [PubMed]
  16. Salleh, S.M.; Mazzoni, G.; Lovendahl, P.; Kadarmideen, H.N. Gene co-expression networks from RNA sequencing of dairy cattle identifies genes and pathways affecting feed efficiency. BMC Bioinform. 2018, 19, 513. [Google Scholar] [CrossRef]
  17. Horodyska, J.; Wimmers, K.; Reyer, H.; Trakooljul, N.; Mullen, A.M.; Lawlor, P.G.; Hamill, R.M. RNA-seq of muscle from pigs divergent in feed efficiency and product quality identifies differences in immune response, growth, and macronutrient and connective tissue metabolism. BMC Genom. 2018, 19, 791. [Google Scholar] [CrossRef]
  18. Khan, S.; Wu, S.B.; Roberts, J. RNA-sequencing analysis of shell gland shows differences in gene expression profile at two time-points of eggshell formation in laying chickens. BMC Genom. 2019, 20, 89. [Google Scholar] [CrossRef]
  19. Koluman, N.; Daskiran, I. Effects of ventilation of the sheep house on heat stress, growth and thyroid hormones of lambs. Trop. Anim. Health Prod. 2011, 43, 1123–1127. [Google Scholar] [CrossRef]
  20. Todini, L.; Malfatti, A.; Valbonesi, A.; Trabalza-Marinucci, M.; Debenedetti, A. Plasma total T3 and T4 concentrations in goats at different physiological stages, as affected by the energy intake. Small Rumin. Res. 2007, 68, 285–290. [Google Scholar] [CrossRef]
  21. Cox, B.; Lee, T.F. Further evidence for a physiological role for hypothalamic dopamine in thermoregulation in the rat. J. Physiol. 1980, 300, 7–17. [Google Scholar] [CrossRef] [PubMed]
  22. Boulant, J.A. Role of the preoptic-anterior hypothalamus in thermoregulation and fever. Clin. Infect. Dis. 2000, 31 (Suppl. 5), S157–S161. [Google Scholar] [CrossRef] [PubMed]
  23. Lucas, G. Serotonin receptors, type 4: A new hope? Curr. Drug Targets 2009, 10, 1085–1095. [Google Scholar] [CrossRef] [PubMed]
  24. Jin, Z.Y.; Ku, B.S.; Yao, H.Y.; Li, Z.H.; Rao, Y. Relationship between the function of 5-HT1A,2 receptors and waking-sleeping stages. Chin. Pharmacol. Bull. 2000, 16, 208–211. [Google Scholar]
  25. Frische, E.W.; Zwartkruis, F.J. Rap1, a mercenary among the Ras-like GTPases. Dev. Biol. 2010, 340, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Wu, C.; Ramirez, A.; Cui, B.; Ding, J.; Delcroix, J.D.; Valletta, J.S.; Liu, J.J.; Yang, Y.; Chu, S.; Mobley, W.C. A functional dynein-microtubule network is required for NGF signaling through the Rap1/MAPK pathway. Traffic 2007, 8, 1503–1520. [Google Scholar] [CrossRef] [PubMed]
  27. Kam, C.Y.; Dubash, A.D.; Magistrati, E.; Polo, S.; Satchell, K.J.F.; Sheikh, F.; Lampe, P.D.; Green, K.J. Desmoplakin maintains gap junctions by inhibiting Ras/MAPK and lysosomal degradation of connexin-43. J. Cell Biol. 2018, 217, 3219–3235. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Yang, Z.Z.; Tschopp, O.; Baudry, A.; Dummler, B.; Hynx, D.; Hemmings, B.A. Physiological functions of protein kinase B/Akt. Biochem. Soc. Trans. 2004, 32, 350–354. [Google Scholar] [CrossRef] [PubMed]
  29. Borrie, S.C.; Brems, H.; Legius, E.; Bagni, C. Cognitive dysfunctions in intellectual disabilities: The contributions of the Ras-MAPK and PI3K-AKT-mTOR pathways. Annu. Rev. Genom. Hum. Genet. 2017, 18, 115–142. [Google Scholar] [CrossRef]
  30. Lee, H.W.; Rhee, D.K.; Kim, B.O.; Pyo, S. Inhibitory effect of sinigrin on adipocyte differentiation in 3T3-L1 cells: Involvement of AMPK and MAPK pathways. Biomed. Pharmacother. 2018, 102, 670–680. [Google Scholar] [CrossRef]
  31. Wu, M.; Liu, D.; Zeng, R.; Xian, T.; Lu, Y.; Zeng, G.; Sun, Z.; Huang, B.; Huang, Q. Epigallocatechin-3-gallate inhibits adipogenesis through down-regulation of PPARgamma and FAS expression mediated by PI3K-AKT signaling in 3T3-L1 cells. Eur. J. Pharmacol. 2017, 795, 134–142. [Google Scholar] [CrossRef] [PubMed]
  32. Reverte-Salisa, L.; Sanyal, A.; Pfeifer, A. Role of cAMP and cGMP signaling in brown fat. In Handbook of Experimental Pharmacology; Springer: Cham, Switzerland, 2018; Volume 251, pp. 161–182. [Google Scholar]
  33. Huang, H.Y.; Zhao, G.P.; Liu, R.R.; Li, Q.H.; Zheng, M.Q.; Li, S.F.; Liang, Z.; Zhao, Z.H.; Wen, J. Brain natriuretic peptide stimulates lipid metabolism through its receptor NPR1 and the glycerolipid metabolism pathway in chicken adipocytes. Biochemistry 2015, 54, 6622–6630. [Google Scholar] [CrossRef] [PubMed]
  34. Rudnicki, M.; Abdifarkosh, G.; Rezvan, O.; Nwadozi, E.; Roudier, E.; Haas, T.L. Female mice have higher angiogenesis in perigonadal adipose tissue than males in response to high-fat diet. Front. Physiol. 2018, 9, 1452. [Google Scholar] [CrossRef] [PubMed]
  35. Wang, L.Y.; Jiang, J.; Ma, H.M. Molecular characterization, tissue expression profile, and SNP analysis of porcine SLC13A5. Genet. Mol. Res. 2015, 14, 16090–16101. [Google Scholar] [CrossRef] [PubMed]
  36. Balasubramanian, S.P.; Azmy, I.A.; Higham, S.E.; Wilson, A.G.; Cross, S.S.; Cox, A.; Brown, N.J.; Reed, M.W. Interleukin gene polymorphisms and breast cancer: A case control study and systematic literature review. BMC Cancer 2006, 6, 188. [Google Scholar] [CrossRef] [PubMed]
  37. Zininga, T.; Ramatsui, L.; Shonhai, A. Heat shock proteins as immunomodulants. Molecules 2018, 23, 2846. [Google Scholar] [CrossRef] [PubMed]
  38. Binder, R.J. Heat-shock protein-based vaccines for cancer and infectious disease. Expert Rev. Vaccines 2008, 7, 383–393. [Google Scholar] [CrossRef]
  39. Shevtsov, M.; Multhoff, G. Heat shock protein-peptide and HSP-based immunotherapies for the treatment of cancer. Front. Immunol. 2016, 7, 171. [Google Scholar] [CrossRef]
Animals 09 00387 g001 550
Figure 1. Distribution of gene expression levels. CG1 and CG2 were males in the control group and CM1 and CM2 were females in the control group, while HG1 and HG2 were heat-stressed males and HM1 and HM2 were heat-stressed females.
Figure 1. Distribution of gene expression levels. CG1 and CG2 were males in the control group and CM1 and CM2 were females in the control group, while HG1 and HG2 were heat-stressed males and HM1 and HM2 were heat-stressed females.
Animals 09 00387 g001
Animals 09 00387 g002 550
Figure 2. Sheep genes identified to be differentially expressed between heat-stressed and control groups. (A) Venn diagram showing the differentially expressed genes. (B) Volcanic plot of the differentially expressed genes. (C) Heat map of the differentially expressed gene clustering analysis. HGHM were heat-stressed males and females, while CGCM were males and females in the control group.
Figure 2. Sheep genes identified to be differentially expressed between heat-stressed and control groups. (A) Venn diagram showing the differentially expressed genes. (B) Volcanic plot of the differentially expressed genes. (C) Heat map of the differentially expressed gene clustering analysis. HGHM were heat-stressed males and females, while CGCM were males and females in the control group.
Animals 09 00387 g002
Animals 09 00387 g003 550
Figure 3. Functional analysis of sheep genes differentially expressed between the heat-stressed and control groups. (A) GO annotation of the differentially expressed genes. (B) KEGG enrichment analysis of the differentially expressed genes.
Figure 3. Functional analysis of sheep genes differentially expressed between the heat-stressed and control groups. (A) GO annotation of the differentially expressed genes. (B) KEGG enrichment analysis of the differentially expressed genes.
Animals 09 00387 g003
Animals 09 00387 g004 550
Figure 4. Validation of differentially expressed genes by qRT-PCR.
Figure 4. Validation of differentially expressed genes by qRT-PCR.
Animals 09 00387 g004
Table
Table 1. Influence of heat stress on serum T3 and T4 concentration of sheep.
Table 1. Influence of heat stress on serum T3 and T4 concentration of sheep.
ItemSexNumberControl GroupHeat Stress Group
T3 (ng/mL)♂53.92 ± 0.041 A1.24 ± 0.020 B
♀53.89 ± 0.045 A1.28 ± 0.016 B
T4 (ng/mL)♂582.92 ± 1.55 A41.89 ± 1.23 B
♀576.05 ± 1.17 A42.73 ± 0.74 B
In the same line, values with different upper-case letters are highly significantly different (p < 0.01).
Table
Table 2. Summary of sequenced RNA-Seq data.
Table 2. Summary of sequenced RNA-Seq data.
ItemsCG1CG2CM1CM2HG1HG2HM1HM2
Raw reads69,008,62839,367,26269,141,27668,952,96678,420,39651,776,41076,392,67460,428,044
Clean reads68,367,48438,903,46868,478,22468,463,51277,762,86251,308,88275,471,93859,933,578
Clean ratio (%)99.0798.8299.0499.2999.1699.198.7999.18
GC content (%)49.8250.0951.1950.9251.1851.0450.8151.25
Q20 (%)98.6598.698.598.798.6598.598.198.6
Q30 (%)95.895.7595.3595.995.895.3594.495.55
Total mapped65,863,115 (96.34%)36,273,376 (93.24%)65,257,697 (95.30%)65,546,575 (95.74%)73,440,220 (94.44%)48,208,876 (93.96%)71,214,523 (94.36%)56,581,847 (94.41%)
Multiple mapped6,565,779 (9.60%)2,743,494 (7.05%)6,088,334 (8.89%)6,757,881 (9.87%)7,483,120 (9.62%)6,034,386 (11.76%)9,336,634 (12.37%)5,753,481 (9.60%)
Unique mapped59,297,336 (86.73%)33,529,882 (86.19%)59,169,363 (86.41%)58,788,694 (85.87%)65,957,100 (84.82%)42,174,490 (82.20%)61,877,889 (81.99%)50,828,366 (84.81%)
Unmapped2,504,369 (3.66%)2,630,092 (6.76%)3,220,527 (4.70%)2,916,937 (4.26%)4,322,642 (5.56%)3,100,006 (6.04%)4,257,415 (5.64%)3,351,731 (5.59%)
CG1 and CG2 were males in the control group and CM1 and CM2 were females in the control group, while HG1 and HG2 were heat-stressed males and HM1 and HM2 were heat-stressed females.

Share and Cite

MDPI and ACS Style

Lu, Z.; Chu, M.; Li, Q.; Jin, M.; Fei, X.; Ma, L.; Zhang, L.; Wei, C. Transcriptomic Analysis Provides Novel Insights into Heat Stress Responses in Sheep. Animals 2019, 9, 387. https://doi.org/10.3390/ani9060387

AMA Style

Lu Z, Chu M, Li Q, Jin M, Fei X, Ma L, Zhang L, Wei C. Transcriptomic Analysis Provides Novel Insights into Heat Stress Responses in Sheep. Animals. 2019; 9(6):387. https://doi.org/10.3390/ani9060387

Chicago/Turabian Style

Lu, Zengkui, Mingxing Chu, Qing Li, Meilin Jin, Xiaojuan Fei, Lin Ma, Liping Zhang, and Caihong Wei. 2019. "Transcriptomic Analysis Provides Novel Insights into Heat Stress Responses in Sheep" Animals 9, no. 6: 387. https://doi.org/10.3390/ani9060387

APA Style

Lu, Z., Chu, M., Li, Q., Jin, M., Fei, X., Ma, L., Zhang, L., & Wei, C. (2019). Transcriptomic Analysis Provides Novel Insights into Heat Stress Responses in Sheep. Animals, 9(6), 387. https://doi.org/10.3390/ani9060387

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop