Next Article in Journal
Inbreeding Alters the Gut Microbiota of the Banna Minipig
Next Article in Special Issue
The Impact of Temperature and Humidity on the Performance and Physiology of Laying Hens
Previous Article in Journal
Social Network Analysis of Agonistic Behaviour and Its Association with Economically Important Traits in Pigs
Previous Article in Special Issue
Association Analysis of Polymorphisms in the 5′ Flanking Region of the HSP70 Gene with Blood Biochemical Parameters of Lactating Holstein Cows under Heat and Cold Stress
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Heat Stress on Dairy Cow Performance and on Expression of Protein Metabolism Genes in Mammary Cells

Department of Agricultural, Food, Environmental and Animal Sciences, University of Udine, 33100 Udine, Italy
*
Author to whom correspondence should be addressed.
Animals 2020, 10(11), 2124; https://doi.org/10.3390/ani10112124
Submission received: 24 September 2020 / Revised: 15 October 2020 / Accepted: 20 October 2020 / Published: 16 November 2020
(This article belongs to the Special Issue Livestock and Heat)

Abstract

:

Simple Summary

Environmental temperatures are increasing, and consequent global warming also has negative effects on dairy cattle farms, which may result in reduced production and poorer milk quality. The protein content of casein, in particular, is important in influencing the coagulation properties of milk and, therefore, the production and quality of cheese. The aim of this study was to assess the effect of heat stress on animal performance and on the expression of selected genes involved in milk protein metabolism. Eight dairy cows were kept under thermoneutral conditions for 8 days. The same animals were then maintained under mild heat stress conditions for an additional 8 days. The results of this study revealed that mild heat stress reduced the feed intake and performance of dairy cows in terms of milk and protein yield, but not the expression of the target genes involved in milk protein metabolism, such as those coding for caseins.

Abstract

The aim of this study was to assess the effect of heat stress on dairy cow performance and on the expression of selected genes involved in milk protein metabolism. Eight Italian Holstein Friesian cows were kept under thermoneutral conditions (temperature–humidity index (THI) < 72, CON) for 8 days and under mild heat stress conditions (72 < THI < 78, HS) for an additional 8 days. The rectal temperature, feed intake, and milk yield were recorded during the last 3 days of the CON and HS periods. During the same time period, milk samples were collected to assess the composition and expression of selected genes involved in milk protein metabolism. Gene expression analyses were performed on somatic cells from milk, which are representative of mammary tissue. In terms of dairy cow performance, HS resulted in lower milk and protein yields and feed intake but higher rectal temperature than for CON (p < 0.05). Under HS, there were greater abundances of HSPA1A (p < 0.05) and BCL2 (p < 0.05), compared to CON, but similar levels of CSN2 (p > 0.05), CSN3 (p > 0.05), HSPA8 (p > 0.05), and STAT5B (p > 0.05) mRNA. Mild heat stress reduced the performance of dairy cows without affecting the expression of genes coding for caseins.

1. Introduction

Global warming and climate change are among the biggest issues facing the world, and their economic impact on dairy farming is a relevant issue. The reduction in milk yield related to heat stress has led to estimated losses of 5.4% of the monthly income of farmers during summer [1]. This impact is expected to significantly increase in the future, with the annual average land temperature in Europe expected to increase by 1.0–5.5 °C by the end of the century [2].
The temperature–humidity index (THI) [3] is widely used to predict heat stress events in cattle, and a value of 72 points was initially considered to be the threshold at which heat stress begins [4]. This threshold was later reduced to 68 points for dairy cows that produce more than 35 kg/d [5]. The content of true protein in milk is important in terms of determining the primary income of farmers and also the dairy industry due to its influence on milk coagulation properties and, hence, cheese production and quality [6]. Heat stress can have a detrimental effect not only on dairy cow performance but also on the milk content of true protein and particularly caseins [7]. These changes could be the result of the indirect effect of reduced dry matter intake and of the direct effect (tissue hyperthermia) on mammary synthesis [8]. There are many studies available which investigate the effect of heat stress on the performance of dairy cows [9,10]. West et al. [11] explained that every additional 1 °C in air temperature above the thermal neutral zone causes a 0.85 kg reduction in feed intake, and that heat stress can cause a considerable reduction in milk production. Fewer studies have characterized the effects of heat stress on milk protein composition or protein fractions, and the results are often contradictory [7,8]. The effects of heat stress on the expression of genes involved in the synthesis of milk proteins are far from being well understood, and few in vivo studies are available [12,13]. When measured in vitro, heat stress increased the expression of heat shock proteins (Hsp) and of genes that code for milk proteins [14,15].
In the past, many studies on gene expression in bovine mammary glands were performed on mammary tissue collected at slaughter or through mammary biopsies that were invasive, did not ensure animal welfare, required surgical procedures, or did not allow for repetitive sampling without injuring the mammary gland. With the aim of avoiding these problems, one study demonstrated that somatic cells from milk are an effective source of mammary transcripts [16]. To our best knowledge, only one study that used milk somatic cells for assessing the effect of heat stress on gene expression in mammary cells is available [17], and this study was performed on goats.
It has previously been hypothesized that heat stress could also have a detrimental effect on milk composition by reducing the expression of target genes in mammary glands. Therefore, the aim of this study was to assess the effect of heat stress on milk production and composition in dairy cows, and on the expression of selected genes involved in milk protein metabolism in milk somatic cells.

2. Materials and Methods

The study was conducted in accordance with EU Directive 2010/63/EU and Italian legislation (DL n. 26, 4 March 2014), and adhered to the rules of the University of Udine. No invasive procedures were applied, and the adopted procedures were routine. The ethical committee of the University of Udine approved the trial (prot. no. 4/2017).

2.1. Animals, Treatments and Sampling

Eight multiparous Italian Holstein Friesian dairy cows in the last stage of lactation and belonging to one commercial farm were considered and transferred to the experiment farm of the University of Udine (Azienda Agraria Universitaria Antonio Servadei) where the trial took place. At the beginning of the trial, the animals had 271 ± 3.7 (mean ± SE) days of milk production, a mean milk production of 14.7 ± 0.84 kg, and a mean milk somatic cell count (SCC) of 131,625 ± 18,643 cells/mL. All animals had an SCC lower than 200,000 cells/mL, which indicated that the mammary gland was free from infection or physiological stress [18,19]. The animals were housed in a tie-stall barn equipped with a fan cooling system and 4 mini data loggers (FT-102; Econorma SAS, Treviso, Italy) that monitored the environmental temperature and relative humidity every second, and recorded averages every thirty minutes. The THI was calculated [3] on the basis of these data considering the hourly average value measured by the probes. The dimensions of the barn were 25 × 10 m with an average height of 4 m. After an adaptation period of 3 weeks, the animals were kept under thermoneutral conditions (CON; THI < 72) for 8 days with a fan used to help maintain these conditions. The same animals were then kept under environmental barn conditions that corresponded to mild heat stress conditions (HS) (THI ranging from 72 to 78 [4]) for an additional 8 days. The average diurnal patterns of the THI during CON and HS conditions are reported in Figure 1. The cows were in good health and had regular veterinary inspections during the entire adaptation and experiment period. The animals were offered sorghum silage, 6 kg of dry matter (DM)—concentrate, 9.2 kg of DM—and hay ad libitum. The utilized concentrate consisted of maize (426 g/kg), soybean meal (215 g/kg), wheat bran (157 g/kg), sunflower meal (81 g/kg), wheat middlings (78 g/kg), and minerals and vitamins (44 g/kg), and this was fed to the cows twice daily during milking. The diet was formulated in order to meet the nutritional requirements of dairy cows in accordance with the Institut National de la Recherche Agronomique (INRA) standard [20]. Samples of hay, silage, and concentrate were collected every two days; dried at 65 °C in a forced draft oven for 48 h; and analyzed for dry matter (DM), ash, ether extract (EE), crude protein (CP) [21], and neutral detergent fiber (NDF) [22]. The energetic value of the feeds, expressed as net energy for lactation (NEL), was estimated following the INRA standard [20]. The ingredients and chemical composition of the feed offered to the dairy cows are reported in Table 1. The daily individual dry-matter intake (DMI) of hay, silage, and concentrate was determined by weighing the amounts of offered and refused feed. During the entire experiment period, the rectal temperature (RT) of the animals was measured three times a day (at 6:00, 14:00, and 22:00 h) with a thermometer (Digital thermometer, GIMA, Milan, Italy) inserted to a depth of 3 cm for at least 3 min. Individual daily milk yield was recorded, and individual milk samples were collected from all four quarters of every cow during morning milking for analyses of composition and gene expression. The milk intended for RNA analysis was filtered with sterile gauze, collected in sterile containers kept on ice, and immediately processed following sample collection. The individual values for the rectal temperature, milk variables (yield, composition, and protein fractions), feed intake, and gene expression obtained during the last 3 days of the CON and HS periods were averaged for each individual, and the result values were used in data analysis. The milk yield was also expressed as fat-corrected milk (FCM) [23].

2.2. Milk Composition and Protein Fractions

Samples obtained from 50 mL of milk were analyzed for fat, protein, and lactose according to ISO 9622:2013 using a MilkoScan FT6000 (FOSS Electric, Hillerød, Denmark) [24], and for somatic cell count (SCC), using a Foss-o-Matic (FOSS Electric, Hillerød, Denmark) according to ISO 1366-2:2008 [25].
Samples obtained from 50 mL of milk were analyzed for milk protein fractions by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) according to Laemmli [26] with modifications. Samples of defatted milk were prepared by adding an equal volume of double-concentrated (2×) sample buffer (0.125 mol/L Tris-HCl, pH 6.8, 4% SDS, 20% glycerol, 10% 2-mercaptoethanol, and 0.2% bromophenol blue), and heated at 98 °C for 5 min. Electrophoresis was performed in a Mighty Small II Mini Vertical Electrophoresis System (Hoefer, Holliston, MA, USA). A stacking gel was prepared with 4% (v/v) acrylamide-N,N′-methylenebisacrylamide (2% Bis Solution, Bio-Rad) at a 32:1 (w/w) ratio, solubilized in 0.125 mol/L Tris-HCl buffer, pH 6.8, with 0.1% (w/v) SDS. The running gel contained 13% acrylamide-N,N′-methylenebisacrylamide at a 36:1 (w/w) ratio, solubilized in 0.375 mol/L Tris-HCl, pH 8.8, with 0.1% (w/v) SDS; 0.1% (w/v) ammonium persulfate and 0.025% (v/v) N,N,N′,N′-tetramethylethylenediamine (TEMED) were added in both gel solutions to initiate polymerization. Electrophoresis was performed at a constant voltage (30 V) until the samples completely left the stacking gel, and then at a constant current of 15 mA (and 250 V) until the tracking dye reached the bottom of the running gel. After electrophoresis, the gels were placed under oscillation in a fixing solution (50% methanol and 10% acetic acid) for 1 h and subsequently placed in a staining solution (10% acetic acid and 0.25% (w/v) Coomassie Blue G-250) overnight. The gels were then destained in a 10% acetic acid solution. Protein fractions were identified on the basis of their electrophoretic behavior [27]. The band molecular weights were estimated by comparison with a molecular weight (MW) protein standard in the range of 10–250 kDa (Precision Plus Protein Dual Color Standards; Bio-Rad). The protein content of the bands from milk were estimated through comparison with bovine serum albumin (BSA) solutions of known concentrations ranging 0.25–2.5 µg/µL; the standard solutions was electrophoresed in a separate lane. Gel image acquisition was performed using the G:Box instrument and Gene Sys software, version 1.5.7.0 (Syngene, Cambridge, UK). The molecular weight estimation and protein quantification of the separated bands were performed with the GeneTools software, version 4.3.9.0 (Syngene).

2.3. Milk Somatic Cell Isolation

Two liters of milk per animal was collected in RNase-free tubes and processed for RNA isolation and Western blot analyses.
For the assessment of somatic cell (SC) gene expression, total RNA was first extracted using the method of Tudisco et al. [28] with modifications. Milk samples were divided into 500 mL tubes (four tubes per animal) and SCs were pelleted by centrifugation at 1500× g for 15 min at 4 °C. The cell pellets were gently washed 3–5 times with 5–10 mL of phosphate-buffered saline (PBS) buffer (pH 7.2) with 0.5 mmol/L ethylenediaminetetraacetic acid (EDTA), to eliminate casein and fat globules, and 0.1% diethylpyrocarbonate (DEPC). After washing, the pellets were gently resuspended in 5–10 mL of PBS–EDTA–DEPC solution, transferred into 50 mL tubes (pellets from two 500 mL tubes were combined into one 50 mL tube), and made up to 50 mL volume using PBS–EDTA–DEPC solution. The tubes were centrifuged at 1000× g for 10 min at 4 °C, and the supernatant was discarded. The cell pellets were resuspended with 100–200 µL PBS–EDTA–DEPC solution and were finally transferred to 2 mL tubes. A portion of the suspended cells was sent for RNA analysis, while the other was stored at −80 °C for protein immunodetection analysis.

2.4. RNA Analysis

RNA extraction was performed using TRI Reagent (Sigma-Aldrich, St. Louis, MI, USA), and deoxyribonuclease (DNase) treatment was carried out using a Precision DNase kit (PrimerDesign, Southampton, UK) according to the manufacturer’s protocols. The concentration and purity of the extracted RNA were assessed using a NanoDrop One spectrophotometer (Thermo Scientific, Waltham, MA, USA). RNA integrity and the absence of genomic DNA were assessed by agarose gel electrophoresis through comparison with MW standards to highlight the presence of only distinct and clear 18S and 28S rRNA bands. Complementary DNA (cDNA) was prepared using the iScript cDNA Synthesis kit (Bio-Rad, Hercules, CA, USA) according to the manufacturer’s protocol. Each 20 μL reaction contained 4 μL reverse transcription reaction mix, 1 μL iScript reverse transcriptase, the appropriate volume of RNA solution required for 50 ng/μL final RNA concentration, and the corresponding volume of nuclease-free water to make up the final volume. The mixture was incubated for 5 min at 25 °C, 60 min at 42 °C, and 5 min at 85 °C before being cooled on ice. Qualitative polymerase chain reaction (PCR) was carried out on a sample pool to validate primer pair specificity for all target genes (β-casein, CSN2; k-casein, CSN3; heat shock protein family A (Hsp70) member 1A, HSPA1A; heat shock protein family A (Hsp70) member 8, HSPA8; B-cell lymphoma 2, BCL2; mammalian target of rapamycin, MTOR; signal transducer and activator of transcription 5B, STAT5B; and reference genes (Table 2)). PCR was performed using the Bio-Rad CFX96 system (Bio-Rad) on a reaction volume of 20 μL containing 0.3 μL of each forward and reverse primer (0.3 μmol/L), 10 μL SsoAdvanced Universal SYBR Green Supermix (Bio-Rad), 8.4 μL sterile water, and 1 μL cDNA. The amplification conditions were one cycle of 3 min at 95 °C; 40 PCR cycles of 15 s at 95 °C, 30 s at 60 °C, and 30 s at 72 °C; and 1 min at 95 °C, followed by a melt curve of 55–95 °C with 0.5 °C increments every 5 s. The occurrence of nonspecific amplification and primer–dimer formation was excluded by checking the melt curves. The amplicon length was verified by agarose gel electrophoresis, comparing with MW standards (Table 2). Quantitative PCR (qPCR) was performed using the same instrument, reagents, and volumes as those for qualitative PCR. Each sample was analyzed in triplicate, and relative gene expression was calculated according to the efficiency-corrected 2−ΔΔCt method [29,30]. β-Actin (ACTB), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), and ribosomal protein large P0 (RPLP0) were initially considered as reference genes. ACTB, being the most stable gene, was used for the normalization of the qPCR data. The data are presented as fold-change ratios using the CON group as a reference. The primer efficiency was calculated using the standard curve obtained by the amplification of the serial dilution of the pooled cDNA [29] (Table 2).

2.5. Immunodetection of Bcl2 and Hsp70 in Somatic Cells

The expression levels of two proteins of interest, Bcl2 and Hsp70, were assessed by SDS-PAGE followed by Western blotting. SC samples were prepared by adding 0.25 vol of 5× sample buffer (0.5 mol/L Tris-HCl, pH 6.8, 10% SDS, 50% glycerol, 25% 2-mercaptoethanol, and 0.2% bromophenol blue) and heating at 98 °C for 5 min. SDS-PAGE was performed as described above. The proteins were then transferred to a nitrocellulose membrane (Hybond ECL; Amersham Pharmacia Biotech, Little Chalfont, UK) using blotting apparatus (Novex X Cell II Blot Module—X Cell Sure Lock Electrophoresis Cell; Life Technologies, Carlsbad, CA, USA). The transfer was performed at 30 V (and 150 mA) for 45 min. After staining with Ponceau S solution (0.1% (w/v) Ponceau S, 5% acetic acid) to verify transfer efficiency, the nitrocellulose sheets were saturated with 5% (w/v) nonfat dry milk in washing buffer (1× PBS, 0.1% Tween 20) for 2 h at room temperature under oscillation. After blocking, the nitrocellulose sheets were washed with washing buffer, cut at 37 kDa in the case of Bcl2 detection and at 50 kDa in the case of Hsp70-Hspa1a detection, and then incubated with the primary specific antibodies at 4 °C overnight under oscillation. For Bcl2 immunodetection, the upper membrane part was probed with antibodies against Actin (anti-Actin, Cat. N. A2066; Sigma-Aldrich; 42 kDa) at a concentration of 1:2000, while the lower part was probed with antibodies against Bcl2 (anti-Bcl2, N. SAB4500003; Sigma-Aldrich; 25 kDa), diluted 1:500. For Hsp70 immunodetection, the upper membrane part was probed with antibodies against Hsp70 (Anti-Hspa1a, N. SAB2107600; Sigma-Aldrich; 70 kDa) diluted 1:500, while the lower part was probed with antibodies against Actin diluted 1:2000. All blots were then rinsed three times with washing buffer and incubated for 2 h at room temperature with the peroxidase-coupled secondary antibody anti-rabbit-IgG (N. SAB3700934; Sigma Aldrich) diluted 1:5000. The membranes were then washed three times in washing buffer and developed with ECL Star (Enhanced Chemiluminescent Substrate; EuroClone, Pero, Mi, Italy), according to the manufacturer’s protocol. The acquisition of blot images was performed with the G:Box instrument and Syngene software (Syngene, Cambridge, UK). The quantification of the proteins was performed with the GeneTools software (Syngene), measuring the band intensity and normalizing it to that of Actin.

2.6. Statistical Analysis

Statistical analysis was performed using R software v. 3.4.0 [38]. The normality of data distribution was tested using the Shapiro–Wilk test. Differences in the means for all the variables between CON and HS were assessed using a paired-sample t-test. In cases where the assumptions of this test were not satisfied, the Wilcoxon signed-rank test was used. Values are reported as mean ± SE; a p value less than 0.05 was considered to indicate significant differences, while a p value less than 0.10 was considered as a tendency towards statistical significance.

3. Results and Discussion

Heat stress increased RT (p < 0.01; Table 3) by 0.5 °C on average. An increase in RT therefore indicates the effectiveness of heat stress [39]. Heat stress is the condition in which dairy cows are not able to completely dissipate the heat load that is derived from their metabolism and the environment [11]; consequently, rectal temperature can increase. In agreement with the results of the present study, Bouraoui et al. [40] observed an increase of 0.5 °C (from 38.4 to 38.9 °C) in the average daily RT of low-producing Holstein Friesian dairy cows when the THI increased from 68 to 78 points. Ting et al. [41] reported rectal temperature values of 39.1 and 38.3 °C in animals maintained for 4 days at a THI greater than 76 and in animals not in conditions of heat stress, respectively. Moreover, RT is significantly correlated with THI, and even variations of less than 1 °C can negatively affect the milk yield and DMI of dairy cattle [42,43]. Karimi et al. [39] explained that the effect of heat stress on DMI is mediated by RT increases. Indeed, HS resulted in lower total DMI than for CON (p < 0.05; Table 4), by about 14%. However, this significant reduction involved the forage fraction, observed for both silage (23%, p < 0.05; Table 4) and hay (30%, p < 0.05; Table 4) but not concentrate (p > 0.05). Garner et al. [44] explained that a reduction in ingestion is an adaptive mechanism for high temperatures that allows the animal to reduce metabolic heat production derived from rumen fermentation. However, the production of endogenous heat due to rumen fermentation is influenced by feed characteristics. Generally, a diet rich in fiber can increase the production of endogenous heat compared to a diet rich in concentrates. Indeed, Coppock and West [45] found that by increasing the hay in diets from 50% to 100%, the conversion efficiency of the metabolizable energy in milk was reduced from 65% to 34%. In agreement with our results, Cowley et al. [8], in considering dairy cows after 7 days of heat stress (a THI of approximately 78), observed a 14% reduction in ingestion. Bouraoui et al. [40] observed a 10% reduction in DMI in heat-stressed dairy cows. Under HS, lower milk and FCM yields were observed compared to those under CON (p < 0.05; Table 3), by about 18% and 13%, respectively. Polsky et al. [46] explained that the reduction in milk production as a result of high THI values can be due to both a reduction in DMI and an increase in the maintenance needs of the animal, which can reach +25%. Many authors found a reduction in milk production in heat-stressed dairy cows. Pragna et al. [47] reviewed that heat stress could lead to a reduction in milk production of up to 30%–40%, and the magnitude of this reduction was positively related to the intensity of the heat stress and milk yield [48]. Ting et al. [41] observed a 19% reduction in milk yield in high-producing dairy cows that had undergone 4 days of heat stress.
Considering milk composition, heat stress tended to increase the content of protein (p < 0.10) and fat (p < 0.10), but not the content of lactose (p > 0.10; Table 3). At the same time, HS reduced the daily protein yield (p < 0.05) without affecting the daily fat yield (p > 0.05; Table 3). The review by Tao et al. [49] concluded that the effect of heat stress on milk fat and protein content is largely inconsistent. In the present study, the tendency toward higher percentage of protein and fat in the milk of heat-stressed cows could be due to the reduction in milk yield and subsequent concentration of protein and fat in addition to possibly greater non-protein nitrogen contents in the milk produced from cows under HS rather than CON. Similar to the present study, Cowley et al. [8] observed that dairy cows under heat stress had reduced daily protein yields. The SCC tended to be higher under HS than CON (p < 0.10; Table 3). Bernabucci et al. [7] also showed a tendency for SCC increase in milk collected in summer in comparison to milk collected in the winter or spring season, whereas Bouraoui et al. [40] observed that the SCC increased significantly as a result of heat stress. Total milk casein (p > 0.05) as well as β-casein, κ-casein, and whey protein (p > 0.05) were not affected by heat stress. Conversely, αs-casein tended to be lower in milk from cows under HS than CON (p < 0.10; Table 3). The average values of casein and whey protein concentrations were similar to those reported by Bernabucci et al. [7] during summer, while the milk casein reported by Cowley et al. [8] was slightly higher at 27.5 g/L. The average values of αs-casein and κ-casein were higher and lower, respectively, than those reported by Cowley et al. [8] by 50.1% and 11.0%, respectively. Barlowska et al. [50] reported that the relationship between αS-casein and κ-casein can vary. The proportion of κ-casein observed in the present study is comparable to that reported by Islam et al. [51] (7%), who explained that the proportion of k-casein to whole casein is highly variable in relation to many factors, such as genetic polymorphism, the stage of lactation, and diet. Conflicting results about the effect of heat stress on milk protein fractions are reported in the literature. Bernabucci et al. [7] reported reduced milk casein and αS-casein concentration in the summer season, though κ-casein was increased without any effect on β-casein when expressed as a percentage of total casein. Cowley et al. [8] reported that heat stress reduced milk casein and αs2-casein concentrations, with no effect on κ-casein, but milk whey protein concentrations increased. More recently, Ma et al. [52] failed to detect differences in casein fractions (αS-, β-, and κ-casein) in dairy cows exposed to heat stress (82.4 for the THI).
The effect of heat stress on the mRNA abundance of genes involved in milk protein metabolism in milk somatic cells is reported in Table 5. Under HS, there was significantly higher mRNA expression of the HSPA1A (p < 0.05) and BCL2 (p < 0.05) genes and a tendency for higher expression of MTOR (p < 0.10) gene than under CON. An effect of heat stress was not found on the mRNA expression of the CSN2, CSN3, HSPA8, and STAT5B genes (p > 0.05). The significant effects of heat stress on BCL2 and HSPA1A expression were also supported by the results of Western blot analysis (Table 6). Indeed, the protein expression of BCL2 (p < 0.05) and HSPA1A (p < 0.05) genes was higher under HS than CON. The HSPA8 gene codes for Hsc70 (known as heat shock cognate 71 kDa protein), which is 86% identical to Hsp70 encoded by HSPA1A [53]. The expression of HSPA8 was not affected by heat stress; indeed, Hsc70 is constitutively expressed, having a key role in the cell in protein folding and unfolding, and it can be considered as a housekeeping gene [53,54]. As expected, Hsp70 was highly expressed during heat stress. Its role is to protect the cell and bind proteins to prevent their aggregation, and it has an antiapoptotic function [55]. In agreement with the results of the present study, Salama et al. [56] and Hu et al. [14] observed an increase in HSPA1A gene expression and the mRNA abundance of Hsp70 in mammary cells incubated at 42 °C in comparison with cells incubated at 36/38 °C. Considering other genes involved in apoptosis, heat stress increased the expression of BCL2, which is an antiapoptotic gene. In their vitro study, Hu et al. [57] observed that under heat stress, cells increased the expression of antiapoptotic factors to reduce heat damage. Under HS, there was tendency for higher MTOR mRNA expression than under CON. Kaufman et al. [58] observed that heat stress reduced MTOR signaling pathway activity in mammary cells in vitro. Conversely, Salama et al. [56] failed to find an effect of heat stress on the mRNA abundance of MTOR in bovine mammary cells in vitro. On the other hand, Chou et al. [59] observed that MTOR has a key role in the synthesis of Hsp; in particular, a reduction in MTOR is accompanied by a reduction in the synthesis of Hsp70; in this way, cells become more sensitive to heat shock. The induction of the gene expression of STAT5 is crucial for the expression of the gene coding for β-casein, CSN2 [33,60]. For this reason, an effect of heat stress on the expression of the CSN2 gene was not found. Moreover, the expression of the gene coding for κ-casein, CSN3, was not affected by heat stress, in agreement with the results for κ-casein in milk.

4. Conclusions

In the present study, dairy cows in the last stage of lactation that were under mild heat stress clearly had reduced daily milk, feed intake, and protein yield. However, little difference was observed in the casein and whey proteins of milk. The expression of selected genes involved in milk protein metabolism in mammary cells was characterized in milk somatic cells, which allow for repetitive sampling without injuring the mammary gland. From this point of view, mRNA and Western blot analyses showed that mammary cells under mild heat stress had increased expression of the antiapoptotic genes HSPA1A and BCL2, probably involved in the protection from heat damage, and the mRNA abundance of genes coding for caseins was not affected. The results of the present study expand our knowledge of the animal’s responses to heat stress and thus lay the foundations for improving or introducing breeding techniques that can improve dairy cow welfare and production.

Author Contributions

Conceptualization, M.C., E.S., G.L., F.D.B., and E.P.; investigation, M.C., A.R., V.F., F.D.B., and E.P.; formal analysis, E.S., G.L., A.R., and F.D.B.; data curation, M.C., G.L., E.S., and E.P.; methodology, E.S. and G.L.; funding acquisition, M.C. and E.P. All authors were involved in developing, writing, and commenting on the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by was funded by the Start-up 2018 project, Department of Agricultural, Food, Environmental and Animal Sciences, University of Udine, and by the Regional Agency for Rural Development of Friuli Venezia Giulia Region (ERSA).

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsor had no role in the design of this study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the final decision to publish the results.

References

  1. Hempel, S.; Menz, C.; Pinto, S.; Galán, E.; Janke, D.; Estellés, F.; Müschner-Siemens, T.; Wang, X.; Heinicke, J.; Zhang, G.; et al. Heat stress risk in European dairy cattle husbandry under different climate change scenarios—Uncertainties and potential impacts. Earth Syst. Dynam. 2019, 10, 859–884. [Google Scholar] [CrossRef] [Green Version]
  2. European Environmental Agency. Global and European Temperature. Available online: https://www.eea.europa.eu/data-and-maps/indicators/global-and-european-temperature-9/assessment (accessed on 26 August 2020).
  3. Yousef, M.K.; Johnson, H.D. Endocrine system and thermal environment. In Stress Physiology in Livestock; Yousef, M.K., Ed.; CRC Press: Boca Raton, FL, USA, 1985; pp. 133–142. [Google Scholar]
  4. Armstrong, D.V. Heat stress interaction with shade and cooling. J. Dairy Sci. 1994, 77, 2044–2050. [Google Scholar] [CrossRef]
  5. Collier, R.J.; Laun, W.H.; Rungruang, S.; Zimbleman, R.B. Quantifying Heat Stress and Its Impact on Metabolism and Performance, Proceedings of the Florida Ruminant Nutrition Symposium, Gainesville, FL, USA, 31 January–1 February 2012; University of Florida: Gainesville, FL, USA, 2012; pp. 74–83. [Google Scholar]
  6. Wedholm, A.; Larsen, L.B.; Lindmark-Månsson, H.; Karlsson, A.H.; Andrén, A. Effect of protein composition on the cheese-making properties of milk from individual dairy cows. J. Dairy Sci. 2006, 89, 3296–3305. [Google Scholar] [CrossRef]
  7. Bernabucci, U.; Barisicò, L.; Morera, P.; Dipasquale, D.; Vitali, A.; Piccioli Cappelli, F.; Calamari, L. Effect of summer season on milk protein fractions in Holstein cows. J. Dairy Sci. 2015, 98, 1815–1827. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Cowley, F.C.; Barber, D.G.; Houlihan, A.V.; Poppi, D.P. Immediate and residual effects of heat stress and restricted intake on milk protein and casein composition and energy metabolism. J. Dairy Sci. 2015, 98, 2356–2368. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Bernabucci, U.; Biffani, S.; Buggioni, L.; Vitali, A.; Lacetera, N.; Nardone, A. The effects of heat stress in Italian Holstein dairy cattle. J. Dairy Sci. 2014, 97, 471–786. [Google Scholar] [CrossRef] [PubMed]
  10. Das, R.; Sailo, L.; Verma, N.; Bharti, P.; Saikia, J.; Imtiwati, I.; Kumar, R. Impact of heat stress on health and performance of dairy animals: A review. Vet. World 2016, 9, 260–268. [Google Scholar] [CrossRef] [Green Version]
  11. West, J.W. Effects of heat-stress on production in dairy cattle. J. Dairy Sci. 2003, 86, 2131–2144. [Google Scholar] [CrossRef]
  12. Dado-Senn, B.; Skibiel, A.L.; Fabris, T.F.; Zhang, Y.; Dahl, G.E.; Peñagaricano, F.; Laporta, J. RNA-Seq reveals novel genes and pathways involved in bovine mammary involution during the dry period and under environmental heat stress. Sci. Rep. 2018, 8, 11096. [Google Scholar] [CrossRef]
  13. Gao, S.T.; Ma, L.; Zhou, Z.; Zhou, Z.K.; Baumgard, L.H.; Jiang, D.; Bionaz, M.; Bu, D.P. Heat stress negatively affects the transcriptome related to overall metabolism and milk protein synthesis in mammary tissue of mid-lactating dairy cows. Physiol. Genom. 2019, 51, 400–409. [Google Scholar] [CrossRef]
  14. Hu, H.; Zhang, Y.; Zheng, N.; Cheng, J.; Wang, J. The effect of heat stress on gene expression and synthesis of heat-shock and milk proteins in bovine mammary epithelial cells. Anim. Sci. J. 2016, 87, 84–91. [Google Scholar] [CrossRef] [PubMed]
  15. Kapila, N.; Sharma, A.; Kishore, A.; Sodhi, M.; Tripathi, P.K.; Mohanty, A.K.; Mukesh, M. Impact of heat stress on cellular and transcriptional adaptation of mammary epithelial cells in riverine buffalo (Bubalus Bubalis). PLoS ONE 2016, 11, e0157237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Cánovas, A.; Rincón, G.; Bevilacqua, C.; Islas-Trejo, A.; Brenaut, P.; Hovey, R.C.; Boutinaud, M.; Morgenthaler, C.; Van Klompenberg, M.K.; Martin, P.; et al. Comparison of five different RNA sources to examine the lactating bovine mammary gland transcriptome using RNA-Sequencing. Sci. Rep. 2015, 4, 5297. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Salama, A.A.K.; Caja, G.; Hamzaoui, S.; Badaoui, B.; Castro-Costa, A.; Façanha, D.E.; Guilhermino, M.M.; Bozzi, R. Different levels of response to heat stress in dairy goats. Small Rumin. Res. 2014, 121, 73–79. [Google Scholar] [CrossRef]
  18. Bradley, A.; Green, M. Use and interpretation of somatic cell count data in dairy cows. Practice 2005, 27, 310–315. [Google Scholar] [CrossRef]
  19. Dohoo, I.R.; Leslie, K.E. Evaluation of changes in somatic cell counts as indicators of new intramammary infections. Prev. Vet. Med. 1991, 10, 225–238. [Google Scholar] [CrossRef]
  20. Agabriel, J. Alimentation de Bovins, Ovins et Caprins; Quae: Versailles, France, 2010. [Google Scholar]
  21. AOAC International. Official Methods of Analysis; AOAC International: Arlington, VA, USA, 2016.
  22. Goering, H.K.; Van Soest, P.J. Forage Fiber Analyses (Apparatus, Reagents, Procedures, and Some Applications). In Agricultural Handbook No.379; USDA Agricultural Research Service: Washington, DC, USA, 1970; pp. 387–598. [Google Scholar]
  23. Gaines, W.L. The Energy Basis of Measuring Milk Yield in Dairy Cows. Bulletin No. 308; Univeristy of Illinois: Urbana, IL, USA, 1928. [Google Scholar]
  24. ISO. Milk and Liquid Milk Products. Guidelines for the Application of Mid-Infrared Spectroscopy (ISO 9622:2013); International Organization for Standardization: Geneve, Switzerland, 2013. [Google Scholar]
  25. ISO. Milk- Enumeration of Somatic Cells—Part 2. Guidance on the Operation of Fluoro-Opto-Electronic Counters (ISO 1366-2:2008); International Organization for Standardization: Geneve, Switzerland, 2008. [Google Scholar]
  26. Laemmli, U.K. Cleavage of structural proteins during the assembly of the head of bacteriophage T4. Nature 1970, 227, 680–685. [Google Scholar] [CrossRef]
  27. Giacinti, G.; Basiricò, L.; Ronchi, B.; Bernabucci, U. Lactoferrin concentration in buffalo milk. Ital. J. Anim. Sci. 2013, 12, e23. [Google Scholar] [CrossRef] [Green Version]
  28. Tudisco, R.; Grossi, M.; Calabrò, S.; Cutrignelli, M.I.; Musco, N.; Addi, L.; Infascelli, F. Influence of pasture on goat milk fatty acids and Stearoyl-CoAdesaturase expression in milk somatic cells. Small Rumin. Res. 2014, 122, 38–43. [Google Scholar] [CrossRef]
  29. Pfaffl, M.W. A new mathematical model for relative quantification in realtime RT-PCR. Nucleic Acids Res. 2001, 29, 2002–2007. [Google Scholar] [CrossRef]
  30. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  31. Boutinaud, M.; Ben Chedly, M.H.; Delamaire, E.; Guinard-Flament, J. Milking and feed restriction regulate transcripts of mammary epithelial cells purified from milk. J. Dairy Sci. 2008, 91, 988–998. [Google Scholar] [CrossRef] [PubMed]
  32. Li, C.; Cai, W.; Zhou, C.; Yin, H.; Zhang, Z.; Loor, J.J.; Sun, D.; Zhang, Q.; Liu, J.; Zhang, S. RNA-Seq reveals 10 novel promising candidate genes affecting milk protein concentration in the Chinese Holstein population. Sci. Rep. 2016, 6, 26813. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Yang, J.X.; Wang, C.H.; Xu, Q.B.; Zhao, F.Q.; Liu, J.X.; Liu, H.Y. Methionyl-methionine promoted milk protein synthesis in bovine mammary gland explants by enhancing intracellular substrate availability and activating JAK2-STAT5 and mTOR-mediated signaling pathways. J. Nutr. 2015, 145, 1748–1753. [Google Scholar] [CrossRef] [PubMed]
  34. Li, S.S.; Loor, J.J.; Liu, H.Y.; Liu, L.; Hosseini, A.; Zhao, W.S.; Liu, J.X. Optimal ratios of essential amino acids stimulate β-casein synthesis via activation of the mammalian target of rapamycin signaling pathway in MAC-T cells and bovine mammary tissue explants. J. Dairy Sci. 2017, 100, 6676–6688. [Google Scholar] [CrossRef] [PubMed]
  35. Duckett, S.K.; Pratt, S.L.; Pavan, E. Corn oil or corn grain supplementation to steers grazing endophyte-free tall fescue. II. Effects on subcutaneous fatty acid content and lipogenic gene expression. J. Anim. Sci 2009, 87, 1120–1128. [Google Scholar] [CrossRef] [Green Version]
  36. Bernabucci, U.; Basiricò, L.; Lacetera, N.; Morera, P.; Ronchi, B.; Accorsi, P.A.; Seren, E.; Nardone, A. Photoperiod affects gene expression of leptin and leptin receptors in adipose tissue from lactating dairy cows. J. Dairy Sci. 2006, 89, 4678–4686. [Google Scholar] [CrossRef]
  37. Wang, Y.H.; Bower, N.I.; Reverter, A.; Tan, S.H.; De Jager, N.; Wang, R.; McWilliam, S.M.; Café, L.M.; Greenwood, P.L.; Lehnert, S.A. Gene expression patterns during intramuscular fat development in cattle. J. Anim. Sci. 2009, 87, 119–130. [Google Scholar] [CrossRef]
  38. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2015; Available online: https://www.R-project.org (accessed on 24 May 2017).
  39. Karimi, M.T.; Ghorbani, G.R.; Kargar, S.; Drackley, J.K. Late-gestation heat stress abatement on performance and behaviour of Holstein dairy cows. J. Dairy Sci. 2015, 98, 6865–6875. [Google Scholar] [CrossRef] [Green Version]
  40. Bouraoui, R.; Lahmar, M.; Majdoub, A.; Djemali, M.; Belyea, R. The relationship of temperature-humidity index with milk production of dairy cows in a mediterranean climate. Anim. Res. 2002, 51, 479–491. [Google Scholar] [CrossRef] [Green Version]
  41. Ting, H.; Zheng, N.; Wang, W.; Cheng, J.; Li, S.; Zhang, Y.; Wang, J. Integrated metabolomics study of the milk of heat-stressed lactating dairy cows. Sci. Rep. 2016, 6, 24208. [Google Scholar]
  42. Johnson, H.D.; Ragsdale, A.C.; Berry, I.L.; Shanklin, M.D. Temperature-Humidity Effects Including Influence of Acclimation in Fed and Water Consumption of Holstein Cattle, Research Bulletin, No. 846; University of Missouri: Columbia, NY, USA, 1963. [Google Scholar]
  43. Liu, J.; Li, L.; Chen, X.; Lu, Y.; Wang, D. Effects of heat stress on body temperature, milk production, and reproduction in dairy cows: A novel idea for monitoring and evaluation of heat stress—A review. Asian Australas. J. Anim. Sci. 2019, 32, 1332–1339. [Google Scholar] [CrossRef] [PubMed]
  44. Garner, J.B.; Douglas, M.L.; Williams, S.R.O.; Wales, W.J.; Marett, L.C.; Nguyen, T.T.T.; Reich, C.M.; Hayes, B.J. Genomic selection improves heat tolerance in dairy cattle. Sci. Rep. 2016, 6, 34114. [Google Scholar] [CrossRef]
  45. Coppock, C.E.; West, J.W. Nutritional adjustment to reduce heat stress in lactating dairy cows. In Proceedings of the Georgia Nutrition Conference for the Feed Industry, Atlanta, GA, USA, 19–21 February 1986; pp. 19–62. [Google Scholar]
  46. Polsky, L.; Marina, A.G.; von Keyserlingk, G. Effects of heat stress on dairy cattle welfare. J. Dairy Sci. 2017, 100, 8645–8657. [Google Scholar] [CrossRef] [Green Version]
  47. Pragna, P.; Archana, P.R.; Aleena, J.; Sejian, V.; Krishnan, G.; Bagath, M.; Manimaran, A.; Beena, V.; Kurien, E.K.; Varma, G.; et al. Heat stress and dairy cow: Impact on both milk yield and composition. Int. J. Dairy Sci. 2016, 12, 1–11. [Google Scholar] [CrossRef] [Green Version]
  48. Kanca, H. Milk production and composition in ruminants under heat stress. In Nutrients in Dairy and Their Implication on Health and Disease; Watson, R.R., Collier, R.J., Preedy, V.R., Eds.; Academic Press: Cambridge, MS, USA, 2017. [Google Scholar]
  49. Tao, S.; Orellana, R.M.; Weng, X.; Marins, T.N.; Dahl, G.E.; Bernard, J.K. The influences of heat stress on bovine mammary gland function. J. Dairy Sci. 2018, 101, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Barlowska, J.; Szwajkowska, M.; Litwińczuk, Z.; Krol, J. Nutritional value and technological suitability of milk from various animal species used for dairy production. Compr. Rev. Food Sci. Food Saf. 2011, 10, 291–302. [Google Scholar] [CrossRef]
  51. Islam, M.A.; Alam, M.K.; Islam, M.N.; Khan, M.A.S.; Ekeberg, D.; Rukke, E.O.; Vegarud, G.E. Principal milk components in buffalo, holstein cross, indigenous cattle and red chittagong cattle from Bangladesh. Asian Australas. J. Anim. 2014, 27, 886–897. [Google Scholar] [CrossRef] [Green Version]
  52. Ma, L.; Yang, Y.; Zhao, X.; Wang, F.; Gao, S.; Bu, D. Heat stress induces proteomic changes in the liver and mammary tissue of dairy cows independent of feed intake: An iTRAQ study. PLoS ONE 2019, 14, e0209182. [Google Scholar] [CrossRef]
  53. Daugaard, M.; Rohde, M.; Jäättelä, M. The heat shock protein 70 family: Highly homologous proteins with overlapping and distinct functions. FEBS Lett. 2007, 581, 3702–3710. [Google Scholar] [CrossRef] [Green Version]
  54. Quintanilla, R.A.; Opazo, F.C.; Pallo, S.P.; Chesser, A.S.; Johnson, G.V.W. Stimulation of tau degradation. In Developing Therapeutics for Alzheimer’s Disease: Progress and Challenges; Wolfe, M., Ed.; Academic Press: Cambridge, MA, USA, 2016. [Google Scholar]
  55. Hassan, F.; Nawaz, A.; Rehman, M.S.; Ali, M.A.; Dilshad, S.M.R.; Yang, C. Prospects of HSP70 as a genetic marker for thermo-tolerance and immuno-modulation in animals under climate change scenario. Anim. Nutr. 2019, 5, 340–350. [Google Scholar] [CrossRef]
  56. Salama, A.A.K.; Duque, M.; Wang, L.; Shahzed, K.; Olivera, M.; Loor, J.J. Enhanced supply of methionine or arginine alters mechanistic target of rapamycin signaling proteins, messenger RNA, and microRNA abundance in heat-stressed bovine mammary epithelial cells in vitro. J. Dairy Sci. 2019, 102, 2469–2480. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Hu, H.; Wang, J.; Gao, H.; Li, S.; Zhang, Y.; Zheng, N. Heat-induced apoptosis and gene expression in bovine mammary epithelial cells. Anim. Prod. Sci. 2016, 56, 918–926. [Google Scholar] [CrossRef]
  58. Kaufman, J.D.; Kassube, K.R.; Almeida, R.A.; Ríus, A.G. High incubation temperature in bovine mammary epithelial cells reduced the activity of the mTOR signaling pathway. J. Dairy Sci. 2018, 101, 7478–7486. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Chou, S.-D.; Prince, T.; Gong, J.; Calderwood, S.K. mTOR Is Essential for the proteotoxic stress response, HSF1 activation and Heat Shock Protein synthesis. PLoS ONE 2012, 7, e39679. [Google Scholar] [CrossRef]
  60. Bionaz, M.; Hurley, W.; Loor, J.J. Milk protein synthesis in the lactating mammary gland: Insight from transcriptomics analysis. Milk Protein 2012, 11, 285–324. [Google Scholar]
Figure 1. Average diurnal patterns of temperature–humidity index (THI) for control and heat stress conditions.
Figure 1. Average diurnal patterns of temperature–humidity index (THI) for control and heat stress conditions.
Animals 10 02124 g001
Table 1. Ingredients and chemical composition of feed offered to dairy cows.
Table 1. Ingredients and chemical composition of feed offered to dairy cows.
ItemHaySorghum SilageConcentrate
Ingredient offered (kg DM/d)Ad libitum69.2
Chemical composition
  DM (%)90.427.587.6
  Ash (% DM)7.46.108.5
  EE (% DM)1.62.44.8
  CP (% DM)7.66.920.5
  NDF (% DM)71.258.823.8
  NEL (MJ/kg DM)4.75.87.6
DM: dry matter; EE: ether extract; CP: crude protein; NDF: neutral detergent fiber; NEL: net energy for lactation.
Table 2. Primer sequences, amplification products, and amplification efficiency for the target genes and reference genes.
Table 2. Primer sequences, amplification products, and amplification efficiency for the target genes and reference genes.
GenePrimers Sequence, 5′ to 3′Amplicon Length (bp)Accession Number
(NCBI; GenBank)
Efficiency
(%)
R2
CSN2 1F: CCCTAACAGCCTCCCACA
R: AGCCATAGCCTCCTTCAC
112KC993858.1101.10.996
CSN3 2F: TGCAATGATGAAGAGTTTTTTCCTAG
R: GATTGGGATATATTTGGCTATTTTGT
150NM_174294.297.10.999
HSPA1A 3F: AACATGAAGAGCGCCGTGGAGG
R: GTTACACACCTGCTCCAGCTCC
171NM_203322.3101.00.988
HSPA8 4F: CGAATCATCAATGAGCCAACTG
R: TGCCACCCCCTAAATCAAAG
100NM_174345.4102.00.999
BCL2 5F: TGTGGATGACCGAGTACCTGAA
R: GACAGCCAGGAGAAATCAAACAG
124NM_001166486.1100.60.997
MTOR 6F: CGTTCCTCTCAACATGGACACA
R: AGCTTCTCCGCGTCTTTACAA
102XM_002694043.596.30.988
STAT5B 7F: GCCAACAATGGTACTTCTCCG
R: TGTGTGACCAGTCGCAGCTC
101NM_174617.498.50.997
ACTB 8F: CTCTTCCAGCCTTCCTTCCT
R: GGGCAGTGATCTCTTTCTGC
177NM_173979.3100.40.999
GAPDH 9F: TCATCCCTGCTTCTACTGGC
R: CCTGCTTCACCACCTTCTTG
176NM_001034034100.70.998
RPLP0 10F: CAACCCCGAAGTGCTTGACAT
R: AGGCAGATGGATCAGCCA
226NM_001012682.1100.20.999
CSN2: β-casein; CSN3: κ-casein; HSPA1A: heat shock protein 70 isoform A1A; HSP8: heat shock protein 70 isoform A8; BCL2: B-cell lymphoma 2; MTOR: mammalian target of rapamycin; STAT5B: signal transducer and activator of transcription 5B; ACTB: β-Actin; GAPDH: glyceraldehyde-3-phosphate dehydrogenase; RPL0: ribosomal protein large P0; F: forward; R: reverse. References: 1 Hu et al. [14]; 2,5 Boutinaud et al. [31]; 3 Kapila et al. [15]; 4 Li et al. [32]; 6 Yang et al. [33]; 7 Li et al. [34]; 8 Duckett et al. [35]; 9 Bernabucci et al. [36]; 10 Wang et al. [37].
Table 3. Effect of heat stress on rectal temperature, milk yield, milk composition, and protein fractions of dairy cows (n = 8).
Table 3. Effect of heat stress on rectal temperature, milk yield, milk composition, and protein fractions of dairy cows (n = 8).
ItemTreatment GroupSEp-Value
CONHS
Rectal temperature (°C)38.639.10.1160.002
Milk yield (kg/d)14.612.00.750.009
FCM (kg/d)14.412.60.6070.020
Milk composition
Protein (%)3.273.370.0430.057
Fat (%)3.834.330.2620.097
Lactose (%)4.444.500.0990.529
SCC × 1000 (cells/mL)12819735.40.091
Milk total casein (%)2.202.140.6100.416
αs-casein (% total casein)60.7759.530.5720.066
β-casein (% total casein)32.9234.170.6690.105
κ-casein (% total casein)6.306.300.4480.998
Milk whey protein (%)0.970.980.0210.734
Protein yield (g/d)48140023.90.012
Fat yield (g/d)56851827.00.111
CON: thermoneutral conditions (temperature–humidity index < 72); HS: heat stress conditions (temperature–humidity index > 72); SE: standard error of the difference; FCM: fat-corrected milk; SCC: somatic cell count.
Table 4. Effect of heat stress on daily intake of dairy cows (n = 8).
Table 4. Effect of heat stress on daily intake of dairy cows (n = 8).
ItemTreatment GroupSEp-Value
CONHS
Concentrate (kg DM)8.618.300.2380.641
Forage (kg DM)7.665.670.6130.014
  Sorghum silage (kg DM)4.393.370.2700.007
  Hay (kg DM)3.272.290.3970.044
Total (kg DM)16.2613.970.7710.021
CON: thermoneutral conditions (temperature–humidity index < 72); HS: dairy heat stress conditions (temperature–humidity index > 72); SE: standard error of the difference; DM: dry matter.
Table 5. Fold change ratio for relative mRNA expression of genes involved in milk protein metabolism (n = 8) in dairy cows 1.
Table 5. Fold change ratio for relative mRNA expression of genes involved in milk protein metabolism (n = 8) in dairy cows 1.
GeneTreatment GroupSEp-Value
CONHS
CSN212.350.4950.215
CSN314.680.9040.164
HSPA811.110.2780.461
HSPA1A11.280.3620.020
BCL211.660.4030.016
MTOR13.000.5420.059
STAT5B12.430.5070.129
CON: thermoneutral conditions (temperature–humidity index < 72); HS: heat stress conditions (temperature–humidity index > 72); SE: standard error of the mean; CSN2: β-casein; CSN3: κ-casein; HSPA1A: heat shock protein 70 isoform A1A; HSP8: heat shock protein 70 isoform A8; BCL2: B-cell lymphoma 2; MTOR: mammalian target of rapamycin; STAT5B: signal transducer and activator of transcription 5B. 1 β-Actin was used for data normalization.
Table 6. Fold change ratio, relative to Actin, of protein expression of BCL2 and HSPA1A genes in dairy cows (n = 8).
Table 6. Fold change ratio, relative to Actin, of protein expression of BCL2 and HSPA1A genes in dairy cows (n = 8).
GeneTreatment GroupSEp-Value
CONHS
BCL20.1000.1330.00500.007
HSPA1A0.1700.3040.04820.049
CON: thermoneutral conditions (temperature–humidity index < 72); HS: heat stress conditions (temperature–humidity index > 72); SE: standard error of the difference; BCL2: B-cell lymphoma 2; HSPA1A: heat shock protein family A (Hsp70) member 1A.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Corazzin, M.; Saccà, E.; Lippe, G.; Romanzin, A.; Foletto, V.; Da Borso, F.; Piasentier, E. Effect of Heat Stress on Dairy Cow Performance and on Expression of Protein Metabolism Genes in Mammary Cells. Animals 2020, 10, 2124. https://doi.org/10.3390/ani10112124

AMA Style

Corazzin M, Saccà E, Lippe G, Romanzin A, Foletto V, Da Borso F, Piasentier E. Effect of Heat Stress on Dairy Cow Performance and on Expression of Protein Metabolism Genes in Mammary Cells. Animals. 2020; 10(11):2124. https://doi.org/10.3390/ani10112124

Chicago/Turabian Style

Corazzin, Mirco, Elena Saccà, Giovanna Lippe, Alberto Romanzin, Vinicius Foletto, Francesco Da Borso, and Edi Piasentier. 2020. "Effect of Heat Stress on Dairy Cow Performance and on Expression of Protein Metabolism Genes in Mammary Cells" Animals 10, no. 11: 2124. https://doi.org/10.3390/ani10112124

APA Style

Corazzin, M., Saccà, E., Lippe, G., Romanzin, A., Foletto, V., Da Borso, F., & Piasentier, E. (2020). Effect of Heat Stress on Dairy Cow Performance and on Expression of Protein Metabolism Genes in Mammary Cells. Animals, 10(11), 2124. https://doi.org/10.3390/ani10112124

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

Article Metrics

Back to TopTop