Skip to main content
Download PDF

Transcriptome analysis reveals the mechanism of chronic heat stress on meat quality of broilers

Journal of Animal Science and Biotechnology volume 13, Article number: 110 (2022) Cite this article

Abstract

Background

Chronic heat stress has a negative impact on poultry meat quality. Although this has been extensively investigated, previous studies have primarily focused on metabolic alterations and oxidative stress in the pectoralis major (PM) muscle under chronic heat stress, and not all of the underlying molecular mechanisms are completely understood.

Methods

A total of 144 male Arbor Acres broilers (28 d old) were randomly allocated into 3 treatment groups: (1) the normal control (NC) group, with broilers raised at 22 °C and fed a basal diet; (2) the heat stress (HS) group, with birds raised at 32 °C and fed a basal diet; and (3) the pair-fed (PF) group, with birds raised at 22 °C and fed the amount of feed equal to the feed consumed on the previous day by the HS group. The experiment lasted for 14 d.

Results

Chronic heat stress decreased the average daily feed intake and average daily gain, increased feed:gain ratio (P < 0.05); and increased drip loss, cooking loss, shear force, hardness, and decreased pH, redness (a*); and springiness of PM muscle (P < 0.05). Furthermore, chronic heat stress decreased muscle fiber density, increased connective tissue, and led to intracellular vacuolation. The transcriptome analyses indicated that the effect of chronic heat stress on meat quality was not only related to metabolism and oxidative stress, but also to signal transduction, immune system, transport and catabolism, cell growth and death, and muscle structure.

Conclusions

Chronic heat stress has a negative impact on the growth performance, meat quality, and the PM muscle structure of broilers. Transcriptome analysis revealed a comprehensive understanding of the mechanism of the chronic heat stress-induced deterioration of broiler meat quality at the transcriptional level.

Introduction

Poultry meat has highly digestible proteins and low-fat content [1]. Modern broilers have a faster growth rate as a result of directional breeding in the past decades. Faster growth rate means greater metabolic activity and higher body heat production. Broilers are more susceptible to undesired high temperatures due to their thermoregulatory system (covered feathers and lack of sweat gland), and global warming renders broilers especially vulnerable to heat stress [2, 3].

Heat stress can cause sharp changes in the physiology and metabolism of birds, which negatively affect a variety of production parameters including feed intake and conversion efficiency, growth, egg production and quality, carcass traits, immune function, meat production, and quality [2, 4, 5]. In our previous study, we studied the metabolome changes caused by chronic heat stress and found that it caused a negative energy balance in broilers and rendering them unable to effectively oxidize fat for energy production, and instead resorting to protein decomposition [6]. Subsequently, our research team found that chronic heat stress impaired meat quality by altering redox status, energy metabolism, amino acid transport, glycolysis, intramuscular fat deposition, and protein synthesis [7, 8]. Although the negative impact of chronic heat stress on meat quality has been intensively studied, previous studies have primarily focused on metabolic changes and oxidative stress in the pectoralis major (PM) muscle under chronic heat stress, and the underlying molecular mechanisms are not completely understood.

The transcriptome is the complete set of transcripts and their quantity in cell or tissue. It can quantify the changing expression levels of each transcript during development and under different conditions [9]. Transcriptome analysis has been widely used to study the molecular mechanisms of tissue damage or metabolic changes caused by environmental stress and it can provide a definitive view of the role of the pathways and genes involved in stress [10, 11]. Therefore, in this study, the transcriptome analyses of broiler PM muscle were performed to identify the differentially expressed genes in broilers experiencing chronic heat stress and to identify the potential signal mechanisms involved. The results could comprehensively reveal the mechanism of inferior meat quality induced by chronic heat stress.

Materials and methods

Animals and experimental design

One hundred and forty-four chickens at 28 days of age (1460 ± 20 g) were randomly divided into three environmental chambers (Kooland Technology Co., Ltd., Beijing, China): normal control group (NC; 22 °C, ad libitum access to feed), heat stress group (HS; 32 °C, ad libitum access to feed) and pair-fed group (PF; 22 °C, fed the amount of feed equal to the feed intake of the previous day by the HS group). Each environmental chamber included six replicates (cages) (length × width × height: 1.2 m × 0.8 m × 0.45 m), and each cage included eight chickens. All birds had free access to water and received a common basal feed. The feed formula is shown in Table 1. The feed was provided in pellet form. The experiment lasted for 14 d.

Table 1 The compositions and nutrient levels of basal diet
Full size table

Daily feed intake was measured at 06:00 every day. The live weights of broilers were recorded on the 7th and 14th d of the experiment. Average daily gain (ADG), average daily feed intake (ADFI), and feed:gain ratio (F:G) were calculated at the end of the experiment.

Infrared thermal image analysis

At the 1st, 3rd, 7th, and 14th d of the experiment, the body surface temperature (BST) of the birds was measured at 10:00 by using a Fluke Ti400 thermal imaging camera (Fluke Corporation, Everett, WA, USA) (Fig. 1A, B, and C). The BST of each chicken was presented as the average temperature of 6 points of measurements.

Fig. 1
figure 1

The body surface temperature of broilers. A, B, C The Infrared thermal image (After 14 d of heat stress) from the NC, HS and PF groups. D Body surface temperature of broilers in three groups in different times. The result represents the mean value ± standard deviation (SD) (n = 6). NC, normal control group; HS, heat stress group; PF, pair-fed group.a,bDifferent letters represent significant differences (P < 0.05)

Full size image

Sample collection

At the 7th (35 days of age) and 14th (42 days of age) d of the experiment, 12 chickens (2 chickens per cage) were selected, electrically stunned, and slaughtered. The entire right PM muscle was stripped for meat quality determination. A sample of muscle (0.5 cm × 0.5 cm × 2 cm) was collected and fixed with 4% paraformaldehyde fixative solution. Additional PM muscle samples were collected in 3 cryovials and stored in liquid nitrogen.

Measurement of meat quality

pH

A pH/ORP meter (HI9125, HANNA Instrument, Padova, Italy) was used to measure the muscle pH at 45 min and 24 h postmortem. The results for each meat sample are the average of 3 measurements at different locations.

Meat color

A chroma meter (CR410, Konica Minolta Sensing Inc., Osaka, Japan) was used to measure meat color at 24 h postmortem. The results for each meat sample are the average of 3 measurements at different locations.

Drip loss

The drip losses were measured as previously described [6]. Briefly, PM muscle samples (1 cm × 1 cm × 3 cm) were placed into a valve bag stored at 4 °C for 24 h, the samples were wiped and weighed (initial weight and final weight). Then, the percentage of lost weight was calculated.

Cooking loss

At 24 h postmortem, PM muscle samples (3 cm × 5 cm × 5 cm) were weighed and placed into valve bags and cooked until the central temperature reached 70 °C in a 72 °C water bath. When the samples were chilled to room temperature, they were wiped and reweighed. Then the percentage of lost weight was calculated.

Shear force

The samples after cooking loss determination were cut into 3 strips (1 cm × 1 cm × 3 cm) parallel to the muscle fibers. Then, the shear force was measured 3 times using a digital meat tenderizer (C1LM3, Northeast Agricultural University, Harbin, China). The results for each meat sample are the average of 3 measurements at different locations.

Texture profile analysis

A cylindrical sample from cooked meat samples was taken with a sampler and then the textural characteristics were assessed by using a TMS-PRO texture analyzer (TA-XT plus, Stable Microsystem Ltd., Surrey, UK). The texture test parameters were as previously described [12].

Hematoxylin and eosin staining

After dehydrated in a graded series of ethanol, paraformaldehyde-fixed muscle samples were embedded in paraffin and cut into 6 µm thick sections, and stained with hematoxylin and eosin. Images were captured using a light microscope (Olympus BX50 microscope, Tokyo, Japan).

RNA extraction, sequencing and analysis

A total of 6 broilers in each group, were mixed in pairs to make 3 samples. Then, total RNA was extracted with 1 mL TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's instructions. RNA quality was assessed using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA). Sequencing libraries were constructed using an RNA sample preparation kit (Illumina, San Diego, CA, USA) according to the manufacturer's instructions. Then, the libraries were sequenced on the Illumina Hiseq Novaseq6000 Platform.

The fastp (version 0.18.0) was used to obtain high-quality, clean reads. The rRNA mapped reads were removed by using the Bowtie2 (version 2.2.8) tool. The software Stringtie and HISAT2 were used to align the clean reads to the reference genome.

The fragments per kilobase of transcript per million mapped fragments (FPKM) of each transcription region were calculated by using the RSEM software to quantify its expression abundance and variations. The R package gmodels [13] was used for principal component analysis (PCA).

DESeq2 software was used for differential expression analysis between two groups, and the genes with P < 0.05 and an absolute fold change >1.5 were considered differentially expressed. Then differentially expressed genes (DEGs) were carried on the GO functional and KEGG pathway enrichment analyses by using GOseq [14] and KOBAS [15], respectively. The P-values (P < 0.05) were used to determine the significantly enriched GO terms and KEGG pathways.

Real-time PCR analysis

Twelve genes were selected for qRT-PCR analysis to verify the reliability of RNA sequencing results (RNA-Seq). Specific primer pairs were designed by using the Primer3 software [16] (Table 2). Reverse transcription was performed by using HiScript III RT Super Mix (Vazyme Biotech Co., Ltd, Nanjing, China), according to the manufacturer's instructions. Real-time qRT-PCR was carried out on an ABI QuantStudio™ 5 (Applied Biosystems, Waltham, MA, USA) using ChamQ Universal SYBR qPCR Master Mix (Vazyme Biotech Co., Ltd, Nanjing, China), according to the manufacturer's instructions. The expressions of relative genes were calculated by using the 2−ΔΔCt method and β-actin was selected as reference gene for normalization.

Table 2 Primer sequences for real-time-quantitative-PCR analysis
Full size table

Statistical analysis

One-way ANOVA was used to analyze the data using SPSS Statistics (version 20.0; SPSS Inc., Chicago, IL, USA). The results are presented as the mean ± standard deviation (SD), and the significant difference was set at P < 0.05. The data of growth performance were analyzed with the cage as the experimental replicate and other indices were analyzed by the mean of two sampled chickens per cage as a replicate (n = 6).

Results

Skin temperature

As shown in Fig. 1, the BST of the NC and PF groups was stably maintained at 38—39 °C, whereas the BST of the HS group was over 40 °C, significantly higher than the NC and PF groups (P < 0.05).

Growth performance

Compared with the NC group, heat stress decreased the ADG and ADFI and increased the F:G ratio (P < 0.05, Table 3). Although the HS and PF groups had the same feed intake, the F:G ratio was higher in the HS group (P < 0.05, Table 3), the ADG of the PF group was significantly lower than that of the NC group after 14 d of the experiment (P < 0.05, Table 3). No significant difference was observed in the F:G ratio between NC and PF groups.

Table 3 Effects of chronic heat stress on the growth performance of broilers
Full size table

Meat quality

At the 7th d of the experiment, the pH24h of PM muscle in the HS group was lower than that in the NC and PF groups (P < 0.05, Table 4). The drip loss and shear force value of the PM muscle in the HS group were higher than those in the NC and PF groups (P < 0.05, Table 4). At the 14th d of the experiment, the pH24h and a* values of PM muscle in the HS group were lower; while the drip loss, cooking loss, and shear force were higher than those in the NC and PF groups (P < 0.05, Table 4).

Table 4 Effects of chronic heat stress on the meat qualities of pectoralis major (PM) muscle in broiler chickens
Full size table

Textural properties

As shown in Table 5, At the 7th d of the experiment, the PM muscle from the HS group had higher hardness and lower cohesiveness than those in the NC and PF groups (P < 0.05). Furthermore, the springiness of the PM muscle in the HS group was declined compared with that in the NC group (P < 0.05). At the 14th d of the experiment, the hardness of the PM muscle from the HS group increased, while springiness decreased (P < 0.05).

Table 5 Effects of chronic heat stress on the textural properties of pectoralis major (PM) muscle in broiler chickens
Full size table

Muscle structure

The morphology of PM muscle in broilers is shown in Fig. 2. At the 7th d of the experiment, PM muscle fibers from the HS group were scattered, the connective tissue increased, and vacuoles appeared in some cells compared with the NC and PF groups. At the 14th d of the experiment, PM muscle damage in the HS group was further aggravated. At both time points, the morphology of the muscle fibers in the NC and PF groups were relatively intact and there was no damage. Heat stress had no significant effect on the fiber diameter of the PM muscle. Compared with the NC group, the muscle fiber density of the PM muscle in the HS and PF groups showed a decreasing trend (Table 6).

Fig. 2
figure 2

Hematoxylin and eosin staining of pectoralis major (PM) muscle sections. A, B, C The HE staining of the PM muscle (After 7 d of heat stress) from the NC, HS and PF groups. D, E, F The HE staining of the PM muscle (After 14 d of heat stress) from the NC, HS and PF groups. NC, normal control group; HS, heat stress group; PF, pair-fed group

Full size image
Table 6 Effects of chronic heat stress on the textural properties of pectoralis major (PM) muscle in broiler chickens
Full size table

RNA sequencing and differential expression analysis

From the PCA result (Fig. 3A) we found that the HS group had a unique transcriptional profile, and heat stress was the principal component (98.3%). In the Venn diagram (Fig. 3B), 8774, 9609, and 8554 genes were identified in the NC, HS, and PF groups, respectively. Compared with the NC group, 1873 differently expressed genes (DEGs) were identified in PM muscle from the HS group, including 1054 up-regulated genes and 369 down-regulated genes. Compared with the PF group, 2485 DEGs were identified in PM muscle from the HS group, including 2036 up-regulated genes and 449 down-regulated genes. Compared with the NC group, 668 DEGs were identified in PM muscle from the PF group, among which 213 genes were up-regulated and 445 genes down-regulated in the PF groups (Fig. 4A). Based on the Venn diagram analysis of DEGs (Fig. 4B), 472 unique DEGs were identified in HS vs. NC, 891 unique DEGs were identified in HS vs. PF and 250 unique DEGs were identified in NC vs. PF. We also found 1338 DEGs shared by HS vs. NC and HS vs. PF; and 99 DEGs shared by HS vs. NC, HS vs. PF, and NC vs. PF. To exclude the factor of feed intake and further investigated the effects of heat stress, we selected 1239 DEGs as the target DEGs, which were shared by HS vs. NC and HS vs. PF after removing the DEGs of NC vs. PF (Fig. 4B). Among the target DEGs, 1034 genes were up-regulated and 205 genes were down-regulated in the HS group compared with NC and PF groups. From these DEGs, we identified some genes that were potentially related to meat quality and listed in Additional file 1: Table S1.

Fig. 3
figure 3

Principal component analysis (PCA) and Wayne (VENN) analysis. A Principal component analysis (PCA) of the samples in three groups. B Venn diagram of the identified genes in pectoralis major (PM) muscle of three groups. NC, normal control group; HS, heat stress group; PF, pair-fed group

Full size image
Fig. 4
figure 4

Differentially expressed genes (DEGs) in HS vs. NC, HS vs. PF and NC vs. PF. A the numbers of DEGs in pectoralis major (PM) muscle after heat stress. B the Venn diagram of DEGs in pectoralis major (PM) muscle after heat stress. NC, normal control group; HS, heat stress group; PF, pair-fed group

Full size image

Heat stress significantly up-regulated genes related to lipid metabolism and synthesis processes such as acyl-CoA synthetase family members (ACSS1A, ACSS1B, ACSS2, AACS), acetyl-CoA carboxylase alpha (ACACA), carnitine palmitoyl transferase 2 (CPT2), acyl-CoA thioesterase 9 (ACOT9), diacylglycerol kinase zeta (DGKZ), glycerol-3-phosphate acyltransferase 3 (GPAT3), 3-hydroxyacyl-CoA dehydratase 3 (HACD3), lipoprotein lipase (LPL), lipase A (LIPA), lipin 2 (LPIN2), palmitoyl-protein thioesterase 1 (PPT1), and peroxisome proliferator-activated receptor gamma (PPARG). Glucose metabolic processes related genes such as fructose-bisphosphatase 2 (FBP2), bisphosphoglycerate mutase (BPGM), glucose-6-phosphate isomerase (GPI), glycogen synthase (GYS1), pyruvate dehydrogenase kinase 4 (PDK4), lactate dehydrogenase A (LDHA), probable D-lactate dehydrogenase mitochondrial isoform X2 (LDHD), pyruvate dehydrogenase phosphatase catalytic subunit 1 (PDP1), fructose-2,6-biphosphatase 3 (PFKFB3), phosphoglucomutase 1 (PGM1), phosphoglycerate mutase 1 (PGAM1), phosphoglucomutase 2 like 1 (PGM2L1), phosphoglycerate kinase 2 (PGK2), pyruvate kinase (PKM), succinate-CoA ligase ADP-forming beta subunit (SUCLA2) and triosephosphate isomerase 1 (TPI1) were down-regulated in PM muscle from HS group compared with that from the NC and PF groups. Aconitase 1 (ACO1), aldolase fructose-bisphosphate C (ALDOC), aconitate decarboxylase 1 (ACOD1), enolase family members (ENO1, and ENO2), fructose-bisphosphatase 1 (FBP1), aldolase fructose-bisphosphate C (GPD1L2), lactate dehydrogenase B (LDHB), phosphogluconate dehydrogenase (PGD), phosphoglucomutase 2 (PGM2), and solute carrier family 2 members (SLC2A14, and SLC2A6) were up-regulated in birds from the HS group. In the HS group, muscle structure related genes such as actin alpha1 (ACTA1), myotubularin related protein 7 (MTMR7), myosin family members (MYH1E, MYH1F, MYL1, MYLK2, MYO18B, MYOM1), tropomodulin 4 (TMOD4), troponinI 2 (TNNI2), troponin I3 (TNNT3), were down-regulated, while actin beta (ACTB), actinin alpha1 (ACTN1), actin filament associated protein 1 (AFAP1), ARP2 actin related protein 2 family members (ACTR2, ACTR3), anillin actin binding protein (ANLN), actin related protein 2/3 complex subunit family members (ARPC1B, ARPC2, ARPC3, ARPC5, ARPC5L), coactosin like F-actin binding protein 1 (COTL1), calponin 2 (CNN2), destrin actin depolymerizing factor (DSTN), Fascin actin-bundling protein 1(FSCN1), leiomodin 2 (LMOD2), myosin family members (MYH1B, MYO1E, MYO1F, MYO3AL, MYO7A), spire type actin nucleation factor 1 (SPIRE1), tropomodulin 1 (TMOD1) were up-regulated. Compared with the NC and PF groups, all genes related to apoptosis listed in Table 2 were up-regulated in birds from the HS group, except sestrin 1 (SESN1), solute carrier family 25 member 4 (SLC25A4), apoptosis-associated tyrosine kinase (AATK), heat shock protein 90 alpha family class B member 1 (HSP90AB1), mitofusin 2 (MFN2), and programmed cell death 2-like (PDCD2L).

Functional annotation and pathway enrichment analysis of the target DEGs

In order to identify the biological function of the target DEGs, they were classified according to GO functional enrichment (Fig. 5A). The most enriched terms in Biological Process were cellular process, biological regulation, metabolic process, regulation of biological process, and response to stimulus metabolic process. In Molecular Function the most enriched terms were binding and catalytic activity. In Cellular Component the most enriched terms were cell part, cell, organelle, membrane. Additionally, immune system process, cytokine production, cell adhesion, ERK1 and ERK2 cascade, cellular response to chemical stimulus, inflammatory response were significantly enriched (Fig. 5B).

Fig. 5
figure 5

GO terms enrichment of target DEGs. A The GO enrichment of the target DEGs. B Top 20 GO enrichment for the target DEGs. NC, normal control group; HS, heat stress group; PF, pair-fed group

Full size image

Furthermore, KEGG pathway analysis was used to identify the biological pathways of the target DEGs (Fig. 6A and B). The KEGG annotation results revealed that the three most enriched metabolic pathways were carbohydrate metabolism, lipid metabolism and amino acid metabolism. The immune system, endocrine system and nervous system were the three most enriched organismal systems. The top three enriched cellular processes pathways were transport and catabolism, cellular community-eukaryotes, and cell growth and death. The three most enriched genetic information processing pathways were folding, sorting and degradation, translation, replication and repair. The top three enriched environmental information processing pathways were signal transduction, signaling molecules and interaction, and membrane transport.

Fig. 6
figure 6

KEGG pathways enrichment of target DEGs. A The KEGG enrichment of the target DEGs. B Top 20 of KEGG pathways enrichment for the target DEGs. NC, normal control group; HS, heat stress group; PF, pair-fed group

Full size image

From the top 20 significantly enriched KEGG pathways (Fig. 6B) we found that metabolic pathways were the most enriched, including glycolysis/gluconeogenesis, pentose phosphate pathway, pyruvate metabolism, biosynthesis of amino acids, glyoxylate and dicarboxylate metabolism, carbon metabolism and sphingolipid metabolism. In addition, lysosome, phagosome, regulation of actin cytoskeleton, HIF-1 signaling pathway, cytokine-cytokine receptor interaction, and toll-like receptor signaling pathway were also significantly enriched.

qRT-PCR validation of differentially expressed genes

The relative mRNA level trend of the genes were consistent with the RNA sequencing data (Fig. 7), indicating that the transcriptomic of the current study were reliable.

Fig. 7
figure 7

qRT-PCR of DEGs compared with RNAseq. The data presented as the relative expression of both RPKM (n = 3) and qRT-PCR (n = 6) and data that are shown as the mean ± standard deviation (SD) values. NC, normal control group; HS, the heat stress group; PF, the pair-fed group. a,bDifferent letters represent significant differences (P < 0.05)

Full size image

Discussion

Broiler is a homeothermic animal that maintains a relatively constant body temperature through the thermoregulatory system by adjusting behavioral and physiological responses to increase or decrease body heat loss. When the ambient temperature is raised beyond the ability of animals to dissipate heat, the animals experience heat stress [17]. Body surface temperature is considered as an ideal physiological indicator of heat stress in animals [2, 18,19,20]. In the current research, BST was used to evaluate the heat stress model. Heat stress significantly elevated the BST from the first day of the experiment, indicating that the broilers were in a state of heat stress.

Previous studies reported that chronic heat stress reduces growth performance by impairing the gut integrity and appetite of broilers [2, 4]. The current study showed that chronic heat stress significantly decreased the ADFI and ADG, and increased the F:G. Although the HS and the PF groups had the same feed intake, the ADG and feed conversion ratio of birds from the HS group were lower than those from the PF group, indicating that the decrease of growth performance is caused by a combination of high temperature and feed intake decrease. In line with previous reports [6, 21]. Our previous research reported that chronic heat stress significantly increased L* values, drip loss, and cooking loss and decreased the pH of broiler PM muscle [7]. In the current study, we found that chronic heat stress increased the drip loss, cooking loss, and shear force value and decreased the pH24h and a* value of the PM muscle. It value of the PM muscle. It has been reported that all texture properties (hardness, elasticity, cohesion, and chewiness) increased after a cyclic high-temperature treatment [22]. Similarly, we found that chronic heat stress resulted in an increase in hardness, but the springiness of the PM muscle was decreased.

In the current study, the transcriptome of PM muscle was used to explore the mechanisms of inferior meat quality induced by chronic heat stress. The PCA result demonstrated that heat stress was the main factor responsible for the difference in the transcriptional profiles of the HS group (98.3%), indicating that chronic heat stress was the major factor causing the deterioration of growth performance and meat quality. Confirmed the previous hypothesis that heat stress reduces growth performance has another mechanism, besides the reduction in feed intake. In order to investigate the effects of chronic heat stress and eliminate the influence of feed intake, we deleted the DEGs between the NC and PF groups. Then, selected 1239 DEGs as the target DEGs from those remaining genes in the HS groups, which were significantly compared with the NC and PF groups. The target genes were related to metabolic process, response to stimulus, catalytic activity, cell part, cell growth and death, HIF-1 signaling pathway, lysosome, regulation of actin cytoskeleton, and sphingolipid metabolism. Among the target DEGs, the genes involved in cell part, membrane, and sphingolipid metabolism are related to cell membrane integrity and apoptosis under heat stress [23]. Evidence shows that damage to muscle membranes may lead to a reduction in meat quality [24], and we observed damaged muscle cells and loose muscle fibers in the HS group. This could be one of the reasons for the increased drip loss and cooking loss observed in this experiment.

Heat stress also affects the redox status and induces oxidative stress. The cytochrome P450 superfamily members (CYPs) are a group of heme-containing proteins that catalyze the activation of molecular oxygen, resulting in the generation of reactive oxygen species (ROS) [25, 26]. In tissues, excessive ROS accumulation can cause oxidative stress [27]. Due to their ability to induce ROS production, the CYPs could be considered an effective RNA biomarker for evaluating oxidative stress. Wang et al. reported that 42 d of NH3 exposure induced up-regulation of the CYPS and oxidative stress in the jejunum of broilers [28]. In this current study, RNA-seq analysis showed that the mRNA expression of CYPS (CYP1C1, CYP2J21, CYP2J23, CYP39A1, CYP4A22) were up-regulated in birds from the HS group (Additional file 1: Table S1). This is in line with our previous finding that chronic heat stress generates excess ROS, resulting to oxidative stress [7]. Noeman et al. showed that ROS induces respiratory chain oxidative damage which in turn could increase ROS generation and a vicious cycle ensues [29]. In the current study, the damage of the muscle fiber caused by heat stress was aggravated with time, demonstrating the occurrence of oxidative stress in the PM muscle induced by heat stress.

Chronic heat stress has been shown to cause lipid metabolism changes [30]. In the present study, some genes associated with fat synthesis were up-regulated, including ACACA, GPAT3, and PPARG. This might lead to the deposition of intramuscular fat, which is consistent with the previous study [7]. Zhao et al. indicated that heat stress significantly reduced the oxidation of fatty acids in skeletal muscle [31]. In this study, although the genes associated with lipolysis such as LPIN2, LPL, and LIPA were up-regulated while the carnitine palmitoyl transferase1 (CPT1), one of the rate-limiting enzymes of fat oxidation, was down-regulated in HS group compared with the NC (log2(FC): −0.35) and PF (log2(FC): −1.04) groups. These results may be the reason for the decrease in beta-oxidation in heat-stressed birds.

Heat stress also affects energy metabolism. In this study, we found that the expression of isocitrate dehydrogenase 3 (NAD ( +)) beta, a rate-limiting enzyme of the tricarboxylic acid cycle, was down-regulated of broilers experienced heat stress, which might result in a reduction in ATP production. These findings were consistent with previous results which showed that chronic heat stress damaged the mitochondrial structure and affected muscle energy metabolism [7, 32]. At the same time, we found that the glucose transporters SLC2A14, SLC2A6, and ATPase family members were up-regulated (Additional file 1: Table S1). This may be a negative feedback regulation of energy reduction. A series of reports have demonstrated that the rate of muscle glycogen breakdown is faster in heat-stressed chickens than that in normal birds [33, 34]. Similarly, we found that the expression of GYS1 in muscle was down-regulated in birds in the HS groups. These results imply that chronic heat stress affects meat quality by affecting muscle energy metabolism.

Several studies have reported that heat stress causes muscle damage. Chronic heat stress causes histopathological damage in broilers [35]. Sula et al. found that heat stress induces multifocal myocyte degeneration, necrosis, and ultimately rhabdomyolysis in the muscles of lambs [36]. Heat stress induces red cusk-eel skeletal muscle atrophy at the molecular level [37]. In the current study, the expression of the actin family members and related protein members, except ACTA1, were up-regulated in the HS group. The expression of myosin family members MYH1E, MYH1F, MYL1, MYLK2, and MYO18B were down-regulated, while MYH1B, MYO1E, MYO1F, MYO3AL, and MYO7A were up-regulated in the HS group. Moreover, the expression of TNNT3 and TNNI2 in the HS group were decreased. Actin, myosin, and troponin are structural proteins of muscle cells and they play a pivotal role in the mechanism of muscular contraction [38, 39]. The molecular architecture of the sarcomere as defined by the myosin filaments and their S-1 and S-2 units, the interaction with the actin filaments, and the boundaries formed by the Z-disks, will subsequently influence basic meat quality traits such as tenderness and water-holding capacity [40]. GO and KEGG analysis showed that regulation of actin cytoskeleton organization and ECM-receptor interaction were significantly enriched. In the present study, histological observations showed that some muscle fibers were degraded, and connective tissue increased in the HS group, demonstrating that chronic heat stress causes muscle damage. This may be related to the increased drip loss, cooking loss, shear force, and hardness observed in the HS group.

Several studies have demonstrated that heat stress affects apoptosis. Heat stress activates the FasL/Fas and TNF-α systems and induces apoptosis of follicular cells in laying hens [41]. Heat exposure induces apoptosis in mouse gastrocnemius muscle cells [42]. In addition, heat stress increased the apoptosis of porcine semitendinosus skeletal muscle cells [43]. Heat stress caused endoplasmic reticulum stress, and oxidative stress activated C/EBP homologous protein (CHOP) and subsequently activated cell apoptosis [44, 45]. Chen et al. found that heat stress reduced the viability and induced apoptosis of myoblasts [46]. In this study, many apoptosis-related genes including TNF receptor superfamily member (TNFRSF), transforming growth factor beta 1 (TGFB1), v-myc avian myelocytomatosis viral oncogene homolog (MYC), Janus kinase 1 (JAK1), and programmed cell death 1 (PDCD1) were up-regulated. In the present study, the histological observation found that vacuoles were observed inside the cell and some muscle fibers were degraded in the HS groups. These results demonstrate that chronic heat stress causes cell apoptosis. Our previous study showed that oxidative stress led to apoptosis and autophagy and resulted in a decline in meat quality by decreasing muscle pH24h and increasing shear force value [47]. This may be one of the reasons for the increased shear force in this study.

Conclusions

In conclusion heat stress has a negative impact on growth performance and meat quality. The transcriptome analyses indicated that the effect of chronic heat stress on meat quality was not only related to metabolism and oxidative stress, but also to signal transduction, immune system, transport and catabolism, cell growth and death and muscle structure. The results provide a comprehensive understanding of the mechanism of chronic heat stress-induced deterioration of broiler meat quality at the transcriptional level.

Availability of data and materials

All data generated or analyzed during this study are available from the corresponding authors on reasonable request.

Abbreviations

NC:

Normal control

HS:

Heat stress

PF:

Pair-fed

PM:

Pectoralis major

ADG:

Average daily gain

ADFI:

Average daily feed intake

F:G:

Feed:gain ratio

BST:

Body surface temperature

References

  1. Marangoni F, Corsello G, Cricelli C, Ferrara N, Ghiselli A, Lucchin L, et al. Role of poultry meat in a balanced diet aimed at maintaining health and wellbeing: an Italian consensus document. Food Nutr Res. 2015;59:27606. https://doi.org/10.3402/fnr.v59.27606.

    Article  PubMed  Google Scholar 

  2. He X, Lu Z, Ma B, Zhang L, Li JL, Jiang Y, et al. Effects of chronic heat exposure on growth performance, intestinal epithelial histology, appetite-related hormones and genes expression in broilers. J Sci Food Agric. 2018;98(12):4471–8. https://doi.org/10.1002/jsfa.8971.

    Article  CAS  PubMed  Google Scholar 

  3. Lu Z, He XF, Ma BB, Zhang L, Li JL, Jiang Y, et al. The alleviative effects and related mechanisms of taurine supplementation on growth performance and carcass characteristics in broilers exposed to chronic heat stress. Poult Sci. 2019;98(2):878–86. https://doi.org/10.3382/ps/pey433.

    Article  CAS  PubMed  Google Scholar 

  4. He X, Lu Z, Ma B, Zhang L, Li JL, Jiang Y, et al. Chronic heat stress alters hypothalamus integrity, the serum indexes and attenuates expressions of hypothalamic appetite genes in broilers. J Therm Biol. 2019;81:110–7. https://doi.org/10.1016/j.jtherbio.2019.02.025.

    Article  CAS  PubMed  Google Scholar 

  5. Vandana GD, Sejian V, Lees AM, Pragna P, Silpa MV, Maloney SK. Heat stress and poultry production: impact and amelioration. Int J Biometeorol. 2021;65(2):163–79. https://doi.org/10.1007/s00484-020-02023-7.

    Article  CAS  PubMed  Google Scholar 

  6. Lu Z, He X, Ma B, Zhang L, Li JL, Jiang Y, et al. Serum metabolomics study of nutrient metabolic variations in chronic heat-stressed broilers. Br J Nutr. 2018;119(7):771–81. https://doi.org/10.1017/S0007114518000247.

    Article  CAS  PubMed  Google Scholar 

  7. Lu Z, He X, Ma B, Zhang L, Li JL, Jiang Y, et al. Chronic heat stress impairs the quality of breast-muscle meat in broilers by affecting redox status and energy-substance metabolism. J Agric Food Chem. 2017;65(51):11251–8. https://doi.org/10.1021/acs.jafc.7b04428.

    Article  CAS  PubMed  Google Scholar 

  8. Ma B, He X, Lu Z, Zhang L, Li JL, Jiang Y, et al. Chronic heat stress affects muscle hypertrophy, muscle protein synthesis and uptake of amino acid in broilers via insulin like growth factor-mammalian target of rapamycin signal pathway. Poult Sci. 2018;97(12):4150–8. https://doi.org/10.3382/ps/pey291.

    Article  CAS  PubMed  Google Scholar 

  9. Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10(1):57–63. https://doi.org/10.1038/nrg2484.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Perini F, Cendron F, Rovelli G, Castellini C, Cassandro M, Lasagna E. Emerging genetic tools to investigate molecular pathways related to heat stress in chickens: a review. Animals (Basel). 2021;11(1):46. https://doi.org/10.3390/ani11010046.

  11. Wang H, Zeng X, Zhang X, Liu H, Xing H. Ammonia exposure induces oxidative stress and inflammation by destroying the microtubule structures and the balance of solute carriers in the trachea of pigs. Ecotoxicol Environ Saf. 2021;212:111974. https://doi.org/10.1016/j.ecoenv.2021.111974.

    Article  CAS  PubMed  Google Scholar 

  12. Gao T, Li JL, Zhang L, Jiang Y, Ma R, Song L, et al. Effect of different tumbling marination treatments on the quality characteristics of prepared pork chops. Asian-Australas J Anim Sci. 2015;28(2):260–7. https://doi.org/10.5713/ajas.14.0511.

    Article  PubMed  PubMed Central  Google Scholar 

  13. The r project for statistical computing. https://www.r-project.org/. Accessed 25 May 2022

  14. Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11(2):R14. https://doi.org/10.1186/gb-2010-11-2-r14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S, et al. Kobas 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011;39(Suppl_2):W316–22. https://doi.org/10.1093/nar/gkr483.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Primer3web version 4.1.0. http://primer3.ut.ee/. Accessed 25 May 2022.

  17. Bernabucci U, Lacetera N, Baumgard LH, Rhoads RP, Ronchi B, Nardone A. Metabolic and hormonal acclimation to heat stress in domesticated ruminants. Animal. 2010;4(7):1167–83. https://doi.org/10.1017/S175173111000090X.

    Article  CAS  PubMed  Google Scholar 

  18. Candido M, Tinoco I, Albino L, Freitas L, Santos TC, Cecon PR, et al. Effects of heat stress on pullet cloacal and body temperature. Poult Sci. 2020;99(5):2469–77. https://doi.org/10.1016/j.psj.2019.11.062.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Yan G, Li H, Zhao W, Shi Z. Evaluation of thermal indices based on their relationships with some physiological responses of housed lactating cows under heat stress. Int J Biometeorol. 2020;64(12):2077–91. https://doi.org/10.1007/s00484-020-01999-6.

    Article  PubMed  Google Scholar 

  20. Jacob BT, Leon P, Rodney H, Matthew ST, Camie S. Pilot study predicting core body temperatures in hot work environments using thermal imagery. J Chem Health Saf. 2019;26(6):75–83. https://doi.org/10.1016/j.jchas.2019.08.001.

    Article  Google Scholar 

  21. Lu Q, Wen J, Zhang H. Effect of chronic heat exposure on fat deposition and meat quality in two genetic types of chicken. Poult Sci. 2007;86(6):1059–64. https://doi.org/10.1093/ps/86.6.1059.

    Article  CAS  PubMed  Google Scholar 

  22. Abudabos AM, Suliman GM, Al-Owaimer AN, Sulaiman A, Alharthi AS. Effects of Nano Emulsified vegetable oil and betaine on growth traits and meat characteristics of broiler chickens reared under cyclic heat stress. Animals (Basel). 2021;11(7):1911. https://doi.org/10.3390/ani11071911.

    Article  Google Scholar 

  23. Lingwood D, Simons K. Lipid rafts as a membrane-organizing principle. Science. 2010;327(5961):46–50. https://doi.org/10.1126/science.1174621.

    Article  CAS  PubMed  Google Scholar 

  24. Scheffler TL, Gerrard DE. Mechanisms controlling pork quality development: The biochemistry controlling postmortem energy metabolism. Meat Sci. 2007;77(1):7–16. https://doi.org/10.1016/j.meatsci.2007.04.024.

    Article  CAS  PubMed  Google Scholar 

  25. Werck-Reichhart D, Feyereisen R. Cytochromes p450: a success story. Genome Biol. 2000;1(6):REVIEWS3003. https://doi.org/10.1186/gb-2000-1-6-reviews3003.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Leung TM, Nieto N. CYP2E1 and oxidant stress in alcoholic and non-alcoholic fatty liver disease. J Hepatol. 2013;58(2):395–8. https://doi.org/10.1016/j.jhep.2012.08.018.

    Article  CAS  PubMed  Google Scholar 

  27. Xing T, Gao F, Tume RK, Zhou G, Xu X. Stress effects on meat quality: a mechanistic perspective. Compr Rev Food Sci Food Saf. 2019;18(2):380–401. https://doi.org/10.1111/1541-4337.12417.

    Article  PubMed  Google Scholar 

  28. Wang S, Li X, Wang W, Zhang H, Xu S. Application of transcriptome analysis: Oxidative stress, inflammation and microtubule activity disorder caused by ammonia exposure may be the primary factors of intestinal microvilli deficiency in chicken. Sci Total Environ. 2019;696:134035. https://doi.org/10.1016/j.scitotenv.2019.134035.

    Article  CAS  PubMed  Google Scholar 

  29. Noeman SA, Hamooda HE, Baalash AA. Biochemical study of oxidative stress markers in the liver, kidney and heart of high fat diet induced obesity in rats. Diabetol Metab Syndr. 2011;3(1):17. https://doi.org/10.1186/1758-5996-3-17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Belhadj SI, Najar T, Ghram A, Abdrrabba M. Heat stress effects on livestock: molecular, cellular and metabolic aspects, a review. J Anim Physiol Anim Nutr (Berl). 2016;100(3):401–12. https://doi.org/10.1111/jpn.12379.

    Article  Google Scholar 

  31. Zhao L, McMillan RP, Xie G, Giridhar S, Baumgard LH, El-Kadi S, et al. Heat stress decreases metabolic flexibility in skeletal muscle of growing pigs. Am J Physiol Regul Integr Comp Physiol. 2018;315(6):R1096–106. https://doi.org/10.1152/ajpregu.00404.2017.

    Article  CAS  PubMed  Google Scholar 

  32. Ma X, Wang L, Shi Z, Chen W, Yang X, Hu Y, et al. Mechanism of continuous high temperature affecting growth performance, meat quality, and muscle biochemical properties of finishing pigs. Genes Nutr. 2019;14:23. https://doi.org/10.1186/s12263-019-0643-9.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Zaboli G, Huang X, Feng X, Ahn DU. How can heat stress affect chicken meat quality? - A review. Poult Sci. 2019;98(3):1551–6. https://doi.org/10.3382/ps/pey399.

    Article  CAS  PubMed  Google Scholar 

  34. Zhang ZY, Jia GQ, Zuo JJ, Zhang Y, Lei J, Ren L, et al. Effects of constant and cyclic heat stress on muscle metabolism and meat quality of broiler breast fillet and thigh meat. Poult Sci. 2012;91(11):2931–7. https://doi.org/10.3382/ps.2012-02255.

    Article  CAS  PubMed  Google Scholar 

  35. Xu Y, Lai X, Li Z, Zhang X, Luo Q. Effect of chronic heat stress on some physiological and immunological parameters in different breed of broilers. Poult Sci. 2018;97(11):4073–82. https://doi.org/10.3382/ps/pey256.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Sula MJ, Winslow CM, Boileau MJ, Barker LD, Panciera RJ. Heat-related injury in lambs. J Vet Diagn Invest. 2012;24(4):772–6. https://doi.org/10.1177/1040638712445772.

    Article  PubMed  Google Scholar 

  37. Dettleff P, Zuloaga R, Fuentes M, Gonzalez P, Aedo J, Estrada JM, et al. Physiological and molecular responses to thermal stress in red cusk-eel (Genypterus chilensis) juveniles reveals atrophy and oxidative damage in skeletal muscle. J Therm Biol. 2020;94:102750. https://doi.org/10.1016/j.jtherbio.2020.102750.

    Article  CAS  PubMed  Google Scholar 

  38. Squire J. Special Issue: The actin-myosin interaction in muscle: Background and overview. Int J Mol Sci. 2019;20(22):5715. https://doi.org/10.3390/ijms20225715.

    Article  CAS  PubMed Central  Google Scholar 

  39. Tobacman LS. Troponin Revealed: Uncovering the structure of the thin filament on-off switch in striated muscle. Biophys J. 2021;120(1):1–9. https://doi.org/10.1016/j.bpj.2020.11.014.

    Article  CAS  PubMed  Google Scholar 

  40. Ertbjerg P, Puolanne E. Muscle structure, sarcomere length and influences on meat quality: A review. Meat Sci. 2017;132:139–52. https://doi.org/10.1016/j.meatsci.2017.04.261.

    Article  CAS  PubMed  Google Scholar 

  41. Li GM, Liu LP, Yin B, Liu YY, Dong WW, Gong S, et al. Heat stress decreases egg production of laying hens by inducing apoptosis of follicular cells via activating the FasL/Fas and TNF-α systems. Poult Sci. 2020;99(11):6084–93. https://doi.org/10.1016/j.psj.2020.07.024.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Chen Y, Yu T. Mouse liver is more resistant than skeletal muscle to heat-induced apoptosis. Cell Stress Chaperones. 2021;26(1):275–81. https://doi.org/10.1007/s12192-020-01163-4.

    Article  PubMed  Google Scholar 

  43. Ganesan S, Pearce SC, Gabler NK, Baumgard LH, Rhoads RP, Selsby JT. Short-term heat stress results in increased apoptotic signaling and autophagy in oxidative skeletal muscle in Sus scrofa. J Therm Biol. 2018;72:73–80. https://doi.org/10.1016/j.jtherbio.2018.01.003.

    Article  CAS  PubMed  Google Scholar 

  44. Sharma S, Chaudhary P, Sandhir R, Bharadwaj A, Gupta RK, Khatri R, et al. Heat-induced endoplasmic reticulum stress in soleus and gastrocnemius muscles and differential response to UPR pathway in rats. Cell Stress Chaperones. 2021;26(2):323–39. https://doi.org/10.1007/s12192-020-01178-x.

    Article  CAS  PubMed  Google Scholar 

  45. Cui Y, Zhou X, Chen L, Tang Z, Mo F, Li XC, et al. Crosstalk between endoplasmic reticulum stress and oxidative stress in heat exposure-induced apoptosis is dependent on the ATF4-CHOP-CHAC1 signal pathway in IPEC-J2 cells. J Agric Food Chem. 2021;69(51):15495–511. https://doi.org/10.1021/acs.jafc.1c03361.

    Article  CAS  PubMed  Google Scholar 

  46. Chen KL, Wang Y, Lin ZP, Li HX. The protective effect of rosmarinic acid on myotube formation during myoblast differentiation under heat stress. In Vitro Cell Dev Biol Anim. 2020;56(8):635–41. https://doi.org/10.1007/s11626-020-00498-7.

    Article  CAS  PubMed  Google Scholar 

  47. Chen X, Zhang L, Li JL, Gao F, Zhou G. Hydrogen peroxide-induced change in meat quality of the breast muscle of broilers is mediated by ros generation, apoptosis, and autophagy in the nf-kappab signal pathway. J Agric Food Chem. 2017;65(19):3986–94. https://doi.org/10.1021/acs.jafc.7b01267.

    Article  CAS  PubMed  Google Scholar 

Download references

Funding

This study was supported by the National Natural Science Foundation of China (32072780, 31872374), the National Key Research and Development Program of China (2016YFD0500501, 2018YFD0500405), and the Earmarked Fund for Jiangsu Agricultural Industry Technology System (JATS [2021] 459).

Author information

Authors and Affiliations

  1. College of Animal Science and Technology, Key Laboratory of Animal Origin Food Production and Safety Guarantee of Jiangsu Province, Jiangsu Collaborative Innovation Center of Meat Production and Processing, Quality and Safety Control, Nanjing Agricultural University, No. 1 Weigang, Weigang, 210095, Nanjing, China

    Zhen Liu, Yingsen Liu, Tong Xing, Lin Zhang & Feng Gao

  2. Institute of Agri-Products Processing, Jiangsu Academy of Agricultural Sciences, Nanjing, 210014, People’s Republic of China

    Jiaolong Li

  3. School of Food Science and Pharmaceutical Engineering, Nanjing Normal University, Nanjing, 210023, People’s Republic of China

    Yun Jiang

Authors
  1. Zhen Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Yingsen Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Tong Xing
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Jiaolong Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Lin Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Yun Jiang
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Feng Gao
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

ZL: Conceptualization, Methodology, Investigation, Formal Analysis, Writing-Original Draft. YSL: Investigation, Data Curation. TX: Investigation, Methodology, Project administration, Supervision. JLL: Writing-Reviewing and Editing. LZ: Writing-Reviewing and Editing. YJ: Resources, Supervision. FG: Conceptualization, Funding Acquisition, Resources, Supervision, Writing -Review & Editing. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Feng Gao.

Ethics declarations

Ethics approval and consent to participate

Animal testing procedures were approved by the Animal Care and Use Committee of Nanjing Agricultural University, Nanjing, China (SYXK (Su)2017–0007).

Consent for publication

All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no financial interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

 Table S1. List of some differentially expressed genes related to meat quality.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Liu, Z., Liu, Y., Xing, T. et al. Transcriptome analysis reveals the mechanism of chronic heat stress on meat quality of broilers. J Animal Sci Biotechnol 13, 110 (2022). https://doi.org/10.1186/s40104-022-00759-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40104-022-00759-3

Keywords

Journal of Animal Science and Biotechnology

ISSN: 2049-1891