Transcriptome analysis reveals that constant heat stress modifies the metabolism and structure of the porcine longissimus dorsi skeletal muscle

A systematic review on the trend of transcriptomic study in livestock: An effort to unwind the complexity of adaptation in a climate change environment

Heat stress has been proven to cause negative effects on livestock leading to lower productivity and economic value. Understanding how heat stress manifests within an animal's body is the first step in devising a heat stress mitigation strategy; transcriptomic studies are one of the methods used. Here, using a systematic literature review methodology, we examine the recent decade of transcriptomics' application to the study of livestock adaptation. We identified 152 studies that met our criteria for using transcriptome methods to heat stress adaptation and were published within the last ten years. Our analysis demonstrates the growing popularity and application of transcriptome approaches in the investigation of the response of ruminants, pigs, and poultry livestock to heat stress. Majority of the works was done in chicken and cattle using multiple organs as the sample, with qRT-PCR as the most employed technique. It has been established that a variety of biomarkers can be used to assess animals under heat stress, such as the HSPs, ILs, and TLRs. Although transcriptomics has lately been employed extensively to uncover the mechanism of heat adaptation, this adaptive feature's complex mechanism remains unclear, leaving many knowledge gaps for investigation. A more complex studies involving more various cell types, organs, or even model organisms using multi-omics approach could be the future research direction in understanding the heat stress effects on livestock better.

Applications of Next-Generation Sequencing Technologies and Statistical Tools in Identifying Pathways and Biomarkers for Heat Tolerance in Livestock

The climate change-associated abnormal weather patterns negatively influences the productivity and performance of farm animals. Heat stress is the major detrimental factor hampering production, causing substantial economic loss to the livestock industry. Therefore, it is important to identify heat-tolerant breeds that can survive and produce optimally in any given environment. To achieve this goal, a clearer understanding of the genetic differences and the underlying molecular mechanisms associated with climate change impacts and heat tolerance are a prerequisite. Adopting next-generation biotechnological and statistical tools like whole transcriptome analysis, whole metagenome sequencing, bisulphite sequencing, genome-wide association studies (GWAS), and selection signatures provides an opportunity to achieve this goal. Through these techniques, it is possible to identify permanent genetic markers for heat tolerance, and by incorporating those markers in marker-assisted breeding selection, it is possible to achieve the target of breeding for heat tolerance in livestock. This review gives an overview of the recent advancements in assessing heat tolerance in livestock using such 'omics' approaches and statistical models. The salient findings from this research highlighted several candidate biomarkers that have the potential to be incorporated into future heat-tolerance studies. Such approaches could revolutionise livestock production in the changing climate scenario and support the food demands of the growing human population.

Deciphering the molecular mechanisms of heat stress tolerance in goats: Insights from transcriptome and Gene Co-expression analysis

Climate change poses a significant threat to the sustainability of livestock production systems in developing countries, particularly impacting small ruminants like goats, which are highly susceptible to heat stress. This stressor not only reduces productivity but also undermines economic viability. This study aimed to delve into the molecular mechanisms underlying heat stress tolerance in goats by conducting a comprehensive transcriptome analysis of heat-tolerant (HT, n = 4) and heat-susceptible (HS, n = 6) Jamunapari goats. Physiological metrics, such as rectal temperature, respiratory rate, and heart rate, were meticulously monitored under extreme environmental conditions (Temperature Humidity Index >92) to effectively classify goats based on their distinct heat stress responses. Samples of blood were obtained, and peripheral blood mononuclear cells (PBMCs) were extracted for subsequent RNA extraction. RNA-Seq analysis revealed a sum of 734 differentially expressed genes (DEGs), comprising 251 upregulated and 483 downregulated genes in HT goats compared to their HS counterparts. The WGCNA revealed three key modules, darkorange (tolerance), paleturquoise (respiration rate), and darkmagenta (heart rate). Moreover, functional enrichment analysis revealed that DEGs within these modules played intricate roles in crucial biological processes and pathways, including mitochondrial function, cardiac function, immune response, genomic stability, and metabolic regulation. This research notably enhances our comprehension of the genetic underpinnings of thermo-tolerance in goats and provides invaluable guidance for formulating breeding strategies aimed at bolstering livestock resilience against the challenges of climate change.

Proteomic changes of the bovine blood plasma in response to heat stress in a tropically adapted cattle breed

Background

Identifying molecular mechanisms responsible for the response to heat stress is essential to increase production, reproduction, health, and welfare. This study aimed to identify early biological responses and potential biomarkers involved in the response to heat stress and animal's recovery in tropically adapted beef cattle through proteomic analysis of blood plasma.

Methods

Blood samples were collected from 14 Caracu males during the heat stress peak (HSP) and 16 h after it (heat stress recovery-HSR) assessed based on wet bulb globe temperature index and rectal temperature. Proteome was investigated by liquid chromatography-tandem mass spectrometry from plasma samples, and the differentially regulated proteins were evaluated by functional enrichment analysis using DAVID tool. The protein-protein interaction network was evaluated by STRING tool.

Results

A total of 1,550 proteins were detected in both time points, of which 84 and 65 were downregulated and upregulated during HSR, respectively. Among the differentially regulated proteins with the highest absolute log-fold change values, those encoded by the GABBR1, EPHA2, DUSP5, MUC2, DGCR8, MAP2K7, ADRA1A, CXADR, TOPBP1, and NEB genes were highlighted as potential biomarkers because of their roles in response to heat stress. The functional enrichment analysis revealed that 65 Gene Ontology terms and 34 pathways were significant (P < 0.05). We highlighted those that could be associated with the response to heat stress, such as those related to the immune system, complement system, hemostasis, calcium, ECM-receptor interaction, and PI3K-Akt and MAPK signaling pathways. In addition, the protein-protein interaction network analysis revealed several complement and coagulation proteins and acute-phase proteins as important nodes based on their centrality and edges.

Conclusion

Identifying differentially regulated proteins and their relationship, as well as their roles in key pathways contribute to improve the knowledge of the mechanisms behind the response to heat stress in naturally adapted cattle breeds. In addition, proteins highlighted herein are potential biomarkers involved in the early response and recovery from heat stress in tropically adapted beef cattle.

Multi-tissue metabolic and transcriptomic responses to a short-term heat stress in swine

BACKGROUND

Heat stress (HS) is an increasing threat for pig production with a wide range of impacts. When submitted to high temperatures, pigs will use a variety of strategies to alleviate the effect of HS. While systemic adaptations are well known, tissue-specific changes remain poorly understood. In this study, thirty-two pigs were submitted to a 5-day HS at 32 °C.

RESULTS

Transcriptomic and metabolomic analyses were performed on several tissues. The results revealed differentially expressed genes and metabolites in different tissues. Specifically, 481, 1774, 71, 1572, 17, 164, and 169 genes were differentially expressed in muscle, adipose tissue, liver, blood, thyroid, pituitary, and adrenal glands, respectively. Regulatory glands (pituitary, thyroid, and adrenal) had a lower number of regulated genes, perhaps indicating an earlier sensitivity to HS. In addition, 7, 8, 2, and 8 metabolites were differentially produced in muscle, liver, plasma, and urine, respectively. The study also focused on the oxidative stress pathway in muscle and liver by performing a correlation analysis between genes and metabolites.

CONCLUSIONS

This study has identified various adaptation mechanisms in swine that enable them to cope with heat stress (HS). These mechanisms include a global decrease in energetic metabolism, as well as changes in metabolic precursors that are linked with protein and lipid catabolism and anabolism. Notably, the adaptation mechanisms differ significantly between regulatory (pituitary, thyroid and adrenal glands) and effector tissues (muscle, adipose tissue, liver and blood). Our findings provide new insights into the comprehension of HS adaptation mechanisms in swine.

Potential role of miR-8159-x in heat stress response in rainbow trout (Oncorhynchus mykiss)

Rainbow trout (Oncorhynchus mykiss) is a representative species of cold-water fish. Elevated temperatures during summer often result in significant high mortality rates. MicroRNAs (miRNAs) are class of small non-coding RNAs that play a crucial role as post-transcriptional regulators in various biological processes. Emerging evidence suggests that miRNAs are important regulators role during heat stress. Analyzing previously obtained miRNA-sequencing data, we observed substantial down regulation of miR-8159-x in the liver tissue of heat stressed rainbow trout. In this study, we conducted a dual luciferase reporter assay to validate that miR-8159-x target, a key gene involved in heat stress in rainbow trout. By examining the expression patterns of miR-8159-x and hsp90a1 in the liver tissue at 18 °C (CG) and 24 °C (HS) groups, we propose that miR-8159-x may negatively regulate hsp90a1. Furthermore, in vitro hepatocyte assay, transfection with miR-8159-x mimics significantly reduced the expression level of hsp90a1, whereas transfection with a miR-8159-x inhibitor yielded the opposite effect. Additionally, overexpression of miR-8159-x inhibited cell proliferation and induced apoptosis in normal rainbow trout hepatocytes. We further investigated the effects of miR-8159-x overexpression or inhibition on the mRNA and protein levels of the target gene hsp90a1 under heat stress conditions. In conclusion, our findings suggest that miR-8159-x participates in the biological response to heat stress by targeting hsp90a1. These results contribute to a better understanding of the molecular mechanisms underlying heat stress in rainbow trout and provide valuable insights for future research.

Competing endogenous RNA (ceRNA) in a non-model animal: Non-coding RNAs respond to heat stress in rainbow trout (Oncorhynchus mykiss) through ceRNA-regulated mechanisms

BACKGROUND

Rainbow trout (Oncorhynchus mykiss) is a typical cold-water fish. With global warming and extreme heat, high summer temperatures are the biggest threat to rainbow trout farming. Rainbow trout initiate stress defense mechanisms in response to thermal stimuli, and competing endogenous RNA (ceRNA) regulation of target genes (mRNAs) mediated by non-coding RNAs (microRNAs [miRNAs], long non-coding RNAs) may be the main strategy for responding to thermal stimuli and enhancing adaptation.

RESULTS

We screened the LOC110485411-novel-m0007-5p-hsp90ab1 ceRNA relationship pairs for affect heat stress in rainbow trout and validated their targeting relationships and functions based on preliminary high-throughput sequencing analysis results. The transfection of exogenous novel-m0007-5p mimics and inhibitors into primary rainbow trout hepatocytes effectively bound and inhibited the target genes hsp90ab1 and LOC110485411 without significant effects on hepatocyte viability, proliferation, and apoptosis. The inhibitory effect of novel-m0007-5p overexpression on hsp90ab1 and LOC110485411 under heat stress was time-efficient. Similarly, small interfering RNAs (siRNAs) affected hsp90ab1 mRNA expression by silencing LOC110485411 expression time-efficiently.

CONCLUSIONS

In conclusion, we found that in rainbow trout, LOC110485411 and hsp90ab1 can bind competitively to novel-m0007-5p via 'sponge adsorption' and that interference with LOC110485411 affects hsp90ab1 expression. These results provide potential for anti-stress drug screening in rainbow trout.

Transcriptome Analyses Reveal Essential Roles of Alternative Splicing Regulation in Heat-Stressed Holstein Cows

Heat stress (HS) severely impacts the productivity and welfare of dairy cows. Investigating the biological mechanisms underlying HS response is crucial for developing effective mitigation and breeding strategies. Therefore, we evaluated the changes in milk yield, physiological indicators, blood biochemical parameters, and alternative splicing (AS) patterns of lactating Holstein cows during thermoneutral (TN, N = 19) and heat-stress (HS, N = 17) conditions. There was a significant (p < 0.05) decline in milk yield as physiological indicators increased after exposure to natural HS conditions. The levels of eight out of 13 biochemical parameters of HS were also significantly altered in the presence of HS (p < 0.05). These results demonstrate that HS negatively influences various biological processes of Holstein cows. Furthermore, we investigated AS events based on the RNA-seq data from blood samples. With HS, five common types of AS events were generally increased by 6.7−38.9%. A total of 3470 AS events corresponding to 3143 unique genes were differentially alternatively spliced (DSGs) (p-adjusted < 0.05) between TN and HS groups. The functional annotation results show that the majority of DSGs are involved in mRNA splicing and spliceosomal complex, followed by enrichment in immune and metabolic processes. Eighty-seven out of 645 differentially expressed genes (DEGs) (fold change ≥ 1.5 and false discovery rate < 0.05) overlapped with DSGs. Further analyses showed that 20 of these genes were significantly enriched for the RNA splicing, RNA binding, and RNA transport. Among them, two genes (RBM25 and LUC7L3) had strong interrelation and co-expression pattern with other genes and were identified as candidate genes potentially associated with HS responses in dairy cows. In summary, AS plays a crucial role in changing the transcriptome diversity of heat-stress-related genes in multiple biological pathways and provides a different regulation mechanism from DEGs.

Identification and functional prediction of long non-coding RNAs related to skeletal muscle development in Duroc pigs

Objective

The growth of pigs involves multiple regulatory mechanisms, and modern molecular breeding techniques can be used to understand the skeletal muscle growth and development to promote the selection process of pigs. This study aims to explore candidate lncRNAs and mRNAs related to skeletal muscle growth and development among Duroc pigs with different average daily gain (ADG).

Methods

A total of 8 pigs were selected and divided into two groups: H group (high-ADG) and L group (low-ADG). And followed by whole transcriptome sequencing to identify differentially expressed (DE) lncRNAs and mRNAs.

Results

In RNA-seq, 703 DE mRNAs (263 up-regulated and 440 down-regulated) and 74 DE lncRNAs (45 up-regulated and 29 down-regulated) were identified. In addition, 1,418 Transcription factors (TFs) were found. Compared with mRNAs, lncRNAs had fewer exons, shorter transcript length and open reading frame length. DE mRNAs and DE lncRNAs can form 417 lncRNA-mRNA pairs (antisense, cis and trans). DE mRNAs and target genes of lncRNAs were enriched in cellular processes, biological regulation, and regulation of biological processes. In addition, quantitative trait locus (QTL) analysis was used to detect the functions of DE mRNAs and lncRNAs, the most of DE mRNAs and target genes of lncRNAs were enriched in QTLs related to growth traits and skeletal muscle development. In single-nucleotide polymorphism/insertion-deletion (SNP/INDEL) analysis, 1,081,182 SNP and 131,721 INDEL were found, and transition was more than transversion. Over 60% of percentage were skipped exon events among alternative splicing events.

Conclusion

The results showed that different ADG among Duroc pigs with the same diet maybe due to the DE mRNAs and DE lncRNAs related to skeletal muscle growth and development.

RNA-seq Analysis Reveals Alternative Splicing Under Heat Stress in Rainbow Trout (Oncorhynchus mykiss)

Rainbow trout (Oncorhynchus mykiss) is one of the most economically important cold-water farmed species in the world, and transcriptomic studies in response to heat stress have been conducted and will be studied in depth. Alternative splicing (AS), a post-transcriptional regulatory process that regulates gene expression and increases proteomic diversity, is still poorly understood in rainbow trout under heat stress. In the present study, 18,623 alternative splicing events were identified from 9936 genes using RNA transcriptome sequencing technology (RNA-Seq) and genomic information. A total of 2731 differential alternative splicing (DAS) events were found among 2179 differentially expressed genes (DEGs). Gene ontology analysis revealed that the DEGs were mainly enriched in cellular metabolic process, cell part, and organic cyclic compound binding under heat stress. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis displayed that the DEGs were enriched for 39 pathways, and some key pathways, such as lysine degradation, are involved in the regulation of heat stress in liver tissues of rainbow trout. The results were validated by qRT-PCR, confirming reliability of our bioinformatics analysis.

Selection signatures for heat tolerance in Brazilian horse breeds

Since domestication, horse breeds have adapted to their environments and differentiated from one another. This paper uses two methods to detect selection signatures in 23 horse breeds, eight of which are Brazilian (610 animals), both cold-blooded and warm-blooded, from temperate and tropical regions. These animals were genotyped using the GGP Equine BeadChip and we analysed the data by Principal Component Analysis (PCA). The samples were separated into groups based on their geographical area of origin and PCA results studied. The genomic regions under selection were detected by hapFLK and PCAdapt methodologies, identifying six regions under selection with at least one Brazilian horse breed. These regions contain genes associated with heat tolerance, skin colour, body size, energy production/metabolism, genes involved in protein degradation/turnover/DNA repair, genes reducing the impact of oxidative stress/cellular repair, and transcriptional regulation. This work confirmed LCORL and NCAPG gene regions in previous studies associated with body size on Equine Chromosome Autosome 3 (ECA3). On the same ECA3, a region implicating genes linked to coat colour was identified, also previously related to heat stress. Regions with genes coding heat shock proteins were found on ECA1 and 2, and many candidate genes for oxidation-reduction which are a natural response to heat stress. However, a larger sample size and whole-genome SNPs are needed to understand better and identify new candidate regions as well as their functional relation with heat tolerance.

Environmental stress and livestock productivity in hot-humid tropics: Alleviation and future perspectives

Tropical environments are characterized by persistently high temperature and relative humidity and the harsh environmental conditions pose a serious limitation on the optimal performance of the animals raised in this region. Heat stress causes deleterious effects on welfare, immunology and physiology of farm animals with a resultant impact on their productivity as the use of body resources is re-organized and the metabolic priorities of animals shift away from production, growth, health and reproduction. It is imperative to understand the mechanisms involved in the thermoregulation of animals under tropical conditions in order to develop appropriate strategies for their improvement. This review focuses on the available data on the increasing global temperature and the adverse impact of tropical conditions on animals' adaptive mechanism affected during thermal stress on production performance, intestinal and ileal microbiome, physiological responses, antioxidant system, metabolic responses, cellular and molecular response, adaptive mechanism strategies to heat stress and also strategies to palliate environmental stress on livestock under humid tropical conditions including environmental manipulation, genetic opportunity, epigenetic and feeding modification. Overall, the present review has identified the disturbance in the physiological indices of tropical livestock and the need for concerted efforts in ameliorating the adverse impacts of high ambient temperature aggravated by high humidity on livestock in tropical environments. Further research is needed on genotype-by-environment interaction on the thermotolerance of different livestock species in the tropics.

Production, reproduction and some adaptation characteristics of Boran cattle breed under changing climate: A systematic review and meta-analysis

Introduction

Climate change affects livestock production and productivity, which could threaten livestock-based food security in pastoral and agro-pastoral production systems of the tropics and sub-tropics. Boran cattle breed is one of the hardiest Zebu cattle reared by Borana Oromo pastoralists for milk and meat production. However, there is limited comprensive information on production, reproduction and adaption traits of the Boran cattle in Ethiopia. Thus, this paper aims to compile the main production, reproduction and some adaptation traits of Boran cattle based on systematic review and meta-analysis of peer reviewed published and unpublished literature.

Methodology

A combination of systematic review and meta-analysis based on PRISMA guideline was employed. Accordingly, out of 646 recorded articles identified through database searching, 64 were found to be eligible for production, reproduction and adaptation characteristics of the Boran cattle, 28 articles were included in qualitative systematic review while 36 articles were used for quantitative meta-analysis.

Result

The Boran cattle breed has the ability to survive, produce and reproduce under high ambient temperature, utilize low quality forage resources, and resist water shortage or long watering intervals and tick infestations. The review revealed that the breed employs various adaptation responses (morphological, physiological, biochemical, metabolic, cellular and molecular responses) to cope with harsh environmental conditions including climate change, rangeland degradation, seasonal feed and water shortages and high incidences of tick infestations. The meta-analysis using a random-effects model allowed provision of pooled estimates of heritability and genetic correlations for reproduction and production traits, which could be used to solve genetic prediction equations under a population level in purebred Boran cattle. In addition, heritability and genetic-correlation estimates found in the present study suggest that there is high genetic variability for most traits in Boran cattle, and that genetic progress is possible for all studied traits in this breed.

Conclusion

The Boran cattle breed has the ability to survive, produce and reproduce under high ambient temperature, utilize low quality forage resources, and resist water shortage or long watering intervals and tick infestations. However, currently there are several challenges such as recurrent droughts, pasture deterioration and lack of systematic selection and breeding programs that play to undermine the realization of the potential of the breed. Thus, we recommend systematic selection for enhancing the reproductive and production performances without compromising the adaptation traits of the breed coupled with improved management of rangelands.

Characterization and Comparative Transcriptomic Analysis of Skeletal Muscle in Pekin Duck at Different Growth Stages Using RNA-Seq

Simple Summary

Skeletal muscle is an important tissue and its development is strictly regulated by genes. In this study, in order to understand the muscle-related gene expression in Pekin duck, RNA-seq was performed to analyze and compare skeletal muscle at different growth stages. Alternative splicing, single nucleotide polymorphisms and insertion–deletions were detected, and 299 novel genes were discovered. MYL4, IGF2BP1, CSRP3, SPP1, KLHL31, LAMB2, LAMA2, ITGB1 and OPN played crucial roles in skeletal muscle development. Oxidative phosphorylation, ECM-receptor interaction, focal adhesion, carbon metabolism, and biosynthesis of amino acids participated in the regulation of skeletal muscle development in Pekin duck. This study provides an important reference for revealing the developmental mechanisms of pectoral and leg muscles in duck.

Abstract

Skeletal muscle, accounting for approximately 50% of body weight, is the largest and most important tissue. In this study, the gene expression profiles and pathways in skeletal muscle of Pekin duck were investigated and compared at embryonic day 17, 21, and 27 and postnatally at 6 months of age. An average of 49,555,936 reads in each sample was obtained from the transcriptome libraries. Over 70.0% of alternative splicing (AS) in each sample was mainly alternative 5′ first exon (transcription start site)—the first exon splicing (TSS) and alternative 3′ last exon (transcription terminal site)—the last exon splicing (TTS), indicating that TSS and TTS were the most common AS event in Pekin ducks, and these AS events were closely related to the regulation of muscle development at different growth stages. The results provided a valuable genomic resource for selective breeding and functional studies of genes. A total of 299 novel genes with ≥2 exons were obtained. There were 294 to 2806 differentially expressed genes (DEGs) in each pairwise comparison of Pekin duck. Notably, 90 DEGs in breast muscle and 9 DEGs in leg muscle were co-expressed at all developmental points. DEGs were validated by qPCR analysis, which confirmed the tendency of the expression. DEGs related to muscle development were involved in biological processes such as “endodermal cell differentiation”, “muscle cell cellular homeostasis”, “skeletal muscle tissue growth” and “skeletal muscle cell differentiation”, and were involved in pathways such as oxidative phosphorylation, ECM-receptor (extracellular matrix receptor) interaction, focal adhesion, carbon metabolism, and biosynthesis of amino acids. Some DEGs, including MYL4, IGF2BP1, CSRP3, SPP1 and KLHL31, as well as LAMB2, LAMA2, ITGB1 and OPN, played crucial roles in muscle growth and development. This study provides valuable information about the expression profile of mRNAs and pathways from duck skeletal muscle at different growth stages, and further functional study of these mRNAs and pathways could provide new ideas for studying the molecular networks of growth and development in duck skeletal muscle.

Differential expression of miRNAs and related mRNAs during heat stress in buffalo heifers

The present experiment was aimed to study differential expression of miRNAs and related mRNAs during heat stress (HS) in buffalo heifers. Twelve Murrah buffalo heifers aged between 1.5 and 2.0 years, weighting between 250 and 300 Kg were randomly assigned into two equal groups. The animals were kept in the psychrometric chamber under Thermo-neutral (TN; THI = 72) and HS (THI = 87-90) conditions for 6 h every day between 1000 and 1600 h for 21 days. The blood sampling was done at 1500 h on 15th day of the experiment and physiological parameters viz. pulse rate (PR), respiratory rate (RR) and rectal temperature (RT) were recorded at 1500 h on day -5, -3, -1, 0, +1, +3, +5 with respect to blood sampling. PBMCs were used for extraction of miRNAs and total RNA; and first strand cDNA was synthesized. qPCR was performed for relative gene expression studies. Physiological, hematological (erythrocytic indices), biochemical (triglycerides, urea, ALT, AST, LDH), redox (SOD, ROS) and endocrine parameters (T₄) altered significantly (P < 0.05) during HS as compared to TN. Out of eight targeted miRNAs only four were expressed in buffalo heifers. The relative expression of bta-mir-142, bta-mir-1248 and bta-mir-2332 was significantly (P < 0.05) up-regulated whereas expression of bta-mir-2478 was significantly (P < 0.05) down-regulated during HS as compared to TN. The relative expression of the predicted target genes i.e. HSF1, HSP60, HSP70, HSPA8 and HSP90 were significantly (P < 0.05) up-regulated whereas HSF4 expression was significantly (P < 0.05) down-regulated during HS as compared to TN. It can be concluded that a THI of 87-90 could lead to a moderate HS in buffalo heifers. Differential expression studies of miRNAs and related mRNAs in present study deciphers the role of miRNAs in the heat tolerance in buffalo heifers.

A novel missense mutation within KRT75 gene strongly affects heat stress in Chinese cattle

The KRT75 gene (Keratin 75) is a member of the type II epithelial α-keratin gene family which plays a key role in hair and nail formation. And the coat conformation affects heat tolarence in mammals. Therefore, the aim of this study was to identify novel single nucleotide polymorphisms (SNPs) of the KRT75 gene and further evaluate its relation to heat stress in Chinese cattle. A missense mutation (NC_037332.1: g.1052 T > C) of the bovine KRT75 was identified using the Bovine Genome Variation Database (BGVD). The g.1052 T > C variant was then genotyped in 519 individuals of 22 cattle breeds. Further analyses showed that the frequency of T allele in Chinese indigenous cattle breeds gradually diminished from northern groups to southern groups, whereas the frequency of C allele displayed a contrary patternl. Simultaneously, the frequency of the CC and CT genotype for southern groups was much greater than that of the TT genotype. Additionally, association analysis showed the genotypes were remarkably associated with mean annual temperature (T), relative humidity (RH) and temperature humidity index (THI) (P < 0.01). Our results demonstrated that the KRT75 gene might be a candidate gene associated with the heat stress in Chinese cattle.

Omics Application in Animal Science-A Special Emphasis on Stress Response and Damaging Behaviour in Pigs

Increasing stress resilience of livestock is important for ethical and profitable meat and dairy production. Susceptibility to stress can entail damaging behaviours, a common problem in pig production. Breeding animals with increased stress resilience is difficult for various reasons. First, studies on neuroendocrine and behavioural stress responses in farm animals are scarce, as it is difficult to record adequate phenotypes under field conditions. Second, damaging behaviours and stress susceptibility are complex traits, and their biology is not yet well understood. Dissecting complex traits into biologically better defined, heritable and easily measurable proxy traits and developing biomarkers will facilitate recording these traits in large numbers. High-throughput molecular technologies (“omics”) study the entirety of molecules and their interactions in a single analysis step. They can help to decipher the contributions of different physiological systems and identify candidate molecules that are representative of different physiological pathways. Here, we provide a general overview of different omics approaches and we give examples of how these techniques could be applied to discover biomarkers. We discuss the genetic dissection of the stress response by different omics techniques and we provide examples and outline potential applications of omics tools to understand and prevent outbreaks of damaging behaviours.

Genome-Wide Transcriptome and Metabolome Analyses Provide Novel Insights and Suggest a Sex-Specific Response to Heat Stress in Pigs

Heat stress (HS) negatively impacts pig production and swine health. Therefore, to understand the genetic and metabolic responses of pigs to HS, we used RNA-Seq and high resolution magic angle spinning (HR-MAS) NMR analyses to compare the transcriptomes and metabolomes of Duroc pigs (n = 6, 3 barrows and 3 gilts) exposed to heat stress (33 °C and 60% RH) with a control group (25 °C and 60% RH). HS resulted in the differential expression of 552 (236 up, 316 down) and 879 (540 up, 339 down) genes and significant enrichment of 30 and 31 plasma metabolites in female and male pigs, respectively. Apoptosis, response to heat, Toll-like receptor signaling and oxidative stress were enriched among the up-regulated genes, while negative regulation of the immune response, ATP synthesis and the ribosomal pathway were enriched among down-regulated genes. Twelve and ten metabolic pathways were found to be enriched (among them, four metabolic pathways, including arginine and proline metabolism, and three metabolic pathways, including pantothenate and CoA biosynthesis), overlapping between the transcriptome and metabolome analyses in the female and male group respectively. The limited overlap between pathways enriched with differentially expressed genes and enriched plasma metabolites between the sexes suggests a sex-specific response to HS in pigs.

Heat Stress Impairs the Physiological Responses and Regulates Genes Coding for Extracellular Exosomal Proteins in Rat

Heat stress (HS) is challenging in humans and animals as it is a complicated regulatory mechanism. This prompted us to characterize the physiological and molecular responses of a HS-animal model. In this study, a rat model system was developed by using three temperature treatments (40 ℃, 42 ℃, and 43 ℃) and sixteen biochemical indicators in blood at 42 ℃ for 30 min (H30), 60 min (H60), and 120 min (H120). In addition, transcriptomic profiling was carried out in H120-rats’ blood, liver, and adrenal gland samples for detection of the genes of interest. Our findings demonstrated that the adrenocorticotropic hormone, catalase, prolactin, growth hormone, and lactic acid have significant spatiotemporal variation in the H120-rats as compared with the control. Furthermore, through transcriptomic screening, we documented a high ratio of differentially expressed genes (DEGs) in adrenal glands, liver, and blood, respectively. Among them, Nup153, Plxnb2, Stx7, Hspa9, Chordc1, Pde4d, Gm2α, and Rnf125 were associated with the regulation of HS and immune response processes. Notably, 36 and 314 of DEGs in blood and adrenal glands were detected in the composition of the extracellular exosome, respectively. Furthermore, the correlation analysis between gene transcripts and biochemical indicator levels identified the Lgals3, S1006, Fn1,F2, and Kng1l1 as key candidate genes for HS encoding extracellular exosomal proteins. On the basis of our results, it was concluded that the current rat model provides a molecular basis for future research in HS resistance in humans and livestock.

Transcriptome Analysis Reveals Potential Regulatory Genes Related to Heat Tolerance in Holstein Dairy Cattle

Heat stress affects the physiology and production performance of Chinese Holstein dairy cows. As such, the selection of heat tolerance in cows and elucidating its underlying mechanisms are vital to the dairy industry. This study aimed to investigate the heat tolerance associated genes and molecular mechanisms in Chinese Holstein dairy cows using a high-throughput sequencing approach and bioinformatics analysis. Heat-induced physiological indicators and milk yield changes were assessed to determine heat tolerance levels in Chinese Holstein dairy cows by Principal Component Analysis method following Membership Function Value Analysis. Results indicated that rectal temperature (RT), respiratory rate (RR), and decline in milk production were significantly lower (p < 0.05) in heat tolerant (HT) cows while plasma levels of heat shock protein (HSP: HSP70, HSP90), and cortisol were significantly higher (p < 0.05) when compared to non-heat tolerant (NHT) Chinese Holstein dairy cows. By applying RNA-Seq analysis, we identified 200 (81 down-regulated and 119 up-regulated) significantly (|log2fold change| ≥ 1.4 and p ≤ 0.05) differentially expressed genes (DEGs) in HT versus NHT Chinese Holstein dairy cows. In addition, 14 of which were involved in protein–protein interaction (PPI) network. Importantly, several hub genes (OAS2, MX2, IFIT5 and TGFB2) were significantly enriched in immune effector process. These findings might be helpful to expedite the understanding for the mechanism of heat tolerance in Chinese Holstein dairy cows.

Heat Stress-Responsive Transcriptome Analysis in the Liver Tissue of Hu Sheep

Heat stress has a severe effect on animal health and can reduce the productivity and reproductive efficiency; it is therefore necessary to explore the molecular mechanism involved in heat stress response, which is helpful for the cultivation of an animal breed with resistance to heat stress. However, little research about heat stress-responsive molecular analysis has been reported in sheep. Therefore, in this study, RNA sequencing (RNA-Seq) was used to investigate the transcriptome profiling in the liver of Hu sheep with and without heat stress. In total, we detected 520 and 22 differentially expressed mRNAs and lncRNAs, respectively. The differentially expressed mRNAs were mainly associated with metabolic processes, the regulation of biosynthetic processes, and the regulation of glucocorticoid; additionally, they were significantly enriched in the heat stress related pathways, including the carbon metabolism, the PPAR signaling pathway, and vitamin digestion and absorption. The co-located differentially expressed lncRNA Lnc_001782 might positively influence the expression of the corresponding genes APOA4 and APOA5, exerting co-regulative effects on the liver function. Thus, we made the hypothesis that Lnc_001782, APOA4 and APOA5 might function synergistically to regulate the anti-heat stress ability in Hu sheep. This study provides a catalog of Hu sheep liver mRNAs and lncRNAs, and will contribute to a better understanding of the molecular mechanism underlying heat stress responses.

Identifying Hub Genes for Heat Tolerance in Water Buffalo (Bubalus bubalis) Using Transcriptome Data

Heat stress has a detrimental effect on the physiological and production performance of buffaloes. Elucidating the underlying mechanisms of heat stress is challenging, therefore identifying candidate genes is urgent and necessary. We evaluated the response of buffaloes (n = 30) to heat stress using the physiological parameters, ELISA indexes, and hematological parameters. We then performed mRNA and microRNA (miRNA) expression profiles analysis between heat tolerant (HT, n = 4) and non-heat tolerant (NHT, n = 4) buffaloes, as well as the specific modules, significant genes, and miRNAs related to the heat tolerance identified using the weighted gene co-expression network analysis (WGCNA). The results indicated that the buffaloes in HT had a significantly lower rectal temperature (RT) and respiratory rate (RR) and displayed a higher plasma heat shock protein (HSP70 and HSP90) and cortisol (COR) levels than those of NHT buffaloes. Differentially expressed analysis revealed a total of 753 differentially expressed genes (DEGs) and 16 differentially expressed miRNAs (DEmiRNAs) were identified between HT and NHT. Using the WGCNA analysis, these DEGs assigned into 5 modules, 4 of which were significantly correlation with the heat stress indexes. Interestingly, 158 DEGs associated with heat tolerance in the turquoise module were identified, 35 of which were found within the protein-protein interaction network. Several hub genes (IL18RAP, IL6R, CCR1, PPBP, IL1B, and IL1R1) were identified that significantly enriched in the Cytokine-cytokine receptor interaction. The findings may help further elucidate the underlying mechanisms of heat tolerance in buffaloes.

Genetic control of longissimus dorsi muscle gene expression variation and joint analysis with phenotypic quantitative trait loci in pigs

Background

Economically important growth and meat quality traits in pigs are controlled by cascading molecular events occurring during development and continuing throughout the conversion of muscle to meat. However, little is known about the genes and molecular mechanisms involved in this process. Evaluating transcriptomic profiles of skeletal muscle during the initial steps leading to the conversion of muscle to meat can identify key regulators of polygenic phenotypes. In addition, mapping transcript abundance through genome-wide association analysis using high-density marker genotypes allows identification of genomic regions that control gene expression, referred to as expression quantitative trait loci (eQTL). In this study, we perform eQTL analyses to identify potential candidate genes and molecular markers regulating growth and meat quality traits in pigs.

Results

Messenger RNA transcripts obtained with RNA-seq of longissimus dorsi muscle from 168 F2 animals from a Duroc x Pietrain pig resource population were used to estimate gene expression variation subject to genetic control by mapping eQTL. A total of 339 eQTL were mapped (FDR ≤ 0.01) with 191 exhibiting local-acting regulation. Joint analysis of eQTL with phenotypic QTL (pQTL) segregating in our population revealed 16 genes significantly associated with 21 pQTL for meat quality, carcass composition and growth traits. Ten of these pQTL were for meat quality phenotypes that co-localized with one eQTL on SSC2 (8.8-Mb region) and 11 eQTL on SSC15 (121-Mb region). Biological processes identified for co-localized eQTL genes include calcium signaling (FERM, MRLN, PKP2 and CHRNA9), energy metabolism (SUCLG2 and PFKFB3) and redox hemostasis (NQO1 and CEP128), and results support an important role for activation of the PI3K-Akt-mTOR signaling pathway during the initial conversion of muscle to meat.

Conclusion

Co-localization of eQTL with pQTL identified molecular markers significantly associated with both economically important phenotypes and gene transcript abundance. This study reveals candidate genes contributing to variation in pig production traits, and provides new knowledge regarding the genetic architecture of meat quality phenotypes.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5386-2) contains supplementary material, which is available to authorized users.

Proteomic changes of the porcine skeletal muscle in response to chronic heat stress

BACKGROUND

Heat stress (HS) has an adverse effect on meat quality; however, the underlying molecular mechanisms altering meat quality due to muscle responses to stress remain unclear. Sixteen castrated male crossbreeds between Landrace × Yorkshire sows and Duroc boars (79.00 ± 1.50 kg body weight) were exposed to either thermal neutral (22 °C, n = 8) or HS (30 °C, n = 8) conditions for 3 weeks. Subsequently, the longissimus dorsi (LD) muscle of all pigs was assayed for meat quality parameters and proteome analysis.

RESULTS

HS decreased post mortem (24 h) pH and intramuscular fat, changed ultimate L*, a* and b* values and increased drip loss and shear force. Proteome analysis of the LD was conducted by two-dimensional gel electrophoresis and mass spectrometry. A total of 23 differentially expressed proteins were identified, of which three were verified by western blotting analysis. The identified proteins were involved in six types of biological process: carbohydrate metabolism, myofibrillar and cytoskeleton structure, stress response, antioxidant and detoxification, calcium binding and cellular apoptosis. Interestingly, HS induced higher levels of heat shock protein, antioxidants and calcium binding proteins, which are involved in the mechanisms of defense and homeostasis.

CONCLUSION

The results indicate that HS-induced changes in the expression of myofibrillar proteins, glucose and energy metabolism-related proteins, heat shock protein and antioxidant enzymes might, at least partly, contribute to increase in meat tenderness. These findings will provide the foundation for developing future mitigating solutions and preventative therapies to reduce the detrimental effects of chronic HS on muscle function, metabolism and meat quality. © 2017 Society of Chemical Industry.