Blood Transcriptomic Analyses Reveal Functional Pathways Associated with Thermotolerance in Pregnant Ewes Exposed to Environmental Heat Stress

Environmental heat stress triggers a series of compensatory mechanisms in sheep that are dependent on their genetic regulation of thermotolerance. Our objective was to identify genes and regulatory pathways associated with thermotolerance in ewes exposed to heat stress. We performed next-generation RNA sequencing on blood collected from 16 pregnant ewes, which were grouped as tolerant and non-tolerant to heat stress according to a physiological indicator. Additional samples were collected to measure complete blood count. A total of 358 differentially expressed genes were identified after applying selection criteria. Gene expression analysis detected 46 GO terms and 52 KEGG functional pathways. The top-three signaling pathways were p53, RIG-I-like receptor and FoxO, which suggested gene participation in biological processes such as apoptosis, cell signaling and immune response to external stressors. Network analysis revealed ATM, ISG15, IRF7, MDM4, DHX58 and TGFβR1 as over-expressed genes with high regulatory potential. A co-expression network involving the immune-related genes ISG15, IRF7 and DXH58 was detected in lymphocytes and monocytes, which was consistent with hematological findings. In conclusion, transcriptomic analysis revealed a non-viral immune mechanism involving apoptosis, which is induced by external stressors and appears to play an important role in the molecular regulation of heat stress tolerance in ewes.

Molecular markers for thermo-tolerance are associated with reproductive and physiological traits in Pelibuey ewes raised in a semiarid environment

Pelibuey sheep exhibit reproductive activity through the year, but warm weather lowers their fertility and demonstrates physiological limitations of environmental heat stress. Single nucleotide polymorphisms (SNPs) associated with heat stress tolerance in sheep have been reported previously. The objective was to validate the association of seven thermo-tolerance SNP markers with reproductive and physiological traits in Pelibuey ewes raised in a semiarid region. Pelibuey ewes were assigned to a cool (January 1st^.- March 31st^.; n = 101) or warm (April 1st^.- August 31st^.; n = 104) experimental group. All ewes were exposed to fertile rams and assessed for pregnancy diagnosis 90 days later; lambing day was reported at birth. These data served to calculate the reproductive traits of services per conception, prolificacy, days to estrus, days to conception, conception rate and lambing rate. Rectal temperature, rump/leg skin temperature and respiratory rate were measured and reported as physiological traits. Blood samples were collected and processed to extract DNA, which was genotyped using the TaqMan allelic discrimination method and qPCR. A mixed effects statistical model was used to validate associations between SNP genotypes and phenotypic traits. The SNPs rs421873172, rs417581105 and rs407804467 were confirmed as markers associated with reproductive and physiological traits (P < 0.05), and these SNPs were in the genes PAM, STAT1 and FBXO11, respectively. Interestingly, these SNP markers resulted as predictors for the evaluated traits but only in ewes from the warm group, which indicated their association with heat-stress tolerance. An additive SNP effect was confirmed with the highest contribution (P < 0.01) of the SNP rs417581105 for the evaluated traits. Reproductive performance improved (P < 0.05) and physiological parameters decreased in ewes carrying favorable SNP genotypes. In conclusion, three thermo-tolerance SNP markers were associated with improved reproductive and physiological traits in a prospective population of heat-stressed ewes raised in a semiarid environment.

Genome-wide association study of a thermo-tolerance indicator in pregnant ewes exposed to an artificial heat-stressed environment

Environmental heat stress negatively influences sheep production in warm semi-arid regions. An animal's ability to tolerate warm weather is difficult to measure naturally due to environmental variability and genetic variation between animals. In this study we developed a thermo-tolerance indicator (TTI) to define heat stress tolerance in pregnant sheep in a controlled environment. Next, we performed a genome-wide association study (GWAS) to identify genomic regions and target genes associated with thermo-tolerance in sheep. Pregnant Columbia-Rambouillet crossbred ewes (n = 127) were heat-stressed inside a climate-controlled chamber for 57 days by increasing the temperature-humidity index to ≥30. Rectal temperature (RT) and feed intake (FI) data were collected daily and used for the predictive TTI analysis. After the tenth day of heat stress, the regression analyses revealed that FI was stable; however, when the ewe's RT exceeded 39.8 °C their FI was less than thermo-tolerant ewes. This average predicted temperature was used to classify each ewe as heat stress tolerant (≤39.8 °C) and non-heat stress tolerant (>39.8 °C). A GWAS analysis was performed and genomic regions were compared between heat stress tolerant and non-tolerant ewes. The single-marker genomic analysis detected 16 single nucleotide polymorphisms (SNP) associated with heat stress tolerance (P < 0.0001), whereas the multi-marker Bayesian analysis identified 8 overlapped 1-Mb chromosomal regions accounting for 11.39% of the genetic variation associated with tolerance to heat stress. Four intragenic SNP showed a remarkable contribution to thermo-tolerance, and these markers were within the genes FBXO11 (rs407804467), PHC3 (rs414179061), TSHR (rs418575898) and STAT1 (rs417581105). In conclusion, genomic regions harboring four intragenic SNP were associated with heat stress tolerance, and these candidate genes are proposed to influence heat tolerance in pregnant ewes subjected to an artificially induced warm climate. Moreover, these genetic markers could be suitable for use in further genetic selection programs in sheep managed in semi-arid regions.

Discovery and validation of candidate SNP markers associated to heat stress response in pregnant ewes managed inside a climate-controlled chamber

Sheep production in desert environments during summer is challenging due to heat stress which reduces feed intake, growth, and fertility. Despite warm conditions, some ewes are able to maintain a normal performance suggesting the existence of genetic bases underlying heat tolerance. Our objective was to discover and validate genetic markers associated with thermo-tolerance in pregnant ewes exposed to warm environmental conditions. Using a well-defined model laboratory of heat stress in sheep, pregnant Columbia-Rambouillet crossbred ewes (n = 100) were examined. Following acclimation to the laboratory at thermo-neutral conditions, heat stress was induced in ewes by increasing the temperature-humidity index in a control environmental chamber during mid-gestation. Feed intake, water consumption, and rectal temperature were recorded daily and used to establish the heat stress tolerance index (HSTI) for each ewe. Rectal temperature was a predictor (P < 0.05) of feed intake, and the regression coefficient was used to classify the HSTI. In a subset of 24 ewes, a genome-wide association study (GWAS) was performed using the Illumina OvineSNP50 BeadChip. Single-marker analysis detected 3 intragenic SNPs associated with HSTI (P value = 10^-5). Bayesian multi-marker approach discovered 26 chromosomal regions across the genome which accounted for 9.8% of the variation associated with HSTI. In an independent sheep population (n = 42), the three discovered SNPs were validated as molecular markers associated with thermo-tolerance phenotypic traits. These SNPs were located within the genes F13A1, PAM, and PRELID2. In conclusion, three SNPs appear to be novel molecular markers associated with heat stress tolerance in pregnant ewes providing new knowledge about genetic foundations of thermo-tolerance.