Heat stress in dairy cattle alters lipid composition of milk

Heat stress in dairy cattle alters lipid composition of milk


Play all audios:


ABSTRACT Heat stress, potentially affecting both the health of animals and the yield and composition of milk, occurs frequently in tropical, sub-tropical and temperate regions. A simulated


acute heat stress experiment was conducted in controlled-climate chambers and milk samples collected before, during and after the heat challenge. Milk lipid composition, surveyed using


LC-MS, showed significant changes in triacylglycerol (TAG) and polar lipid profiles. Heat stress (temperature-humidity index up to 84) was associated with a reduction in TAG groups


containing short- and medium-chain fatty acids and a concomitant increase in those containing long-chain fatty acids. The abundance of five polar lipid classes including


phosphatidylethanolamine, phosphatidylserine, phosphatidylcholine, lysophosphatidylcholine and glucosylceramide, was found to be significantly reduced during heat stress.


Lysophosphatidylcholine, showing the greatest reduction in concentration, also displayed a differential response between heat tolerant and heat susceptible cows during heat stress. This


phospholipid could be used as a heat stress biomarker for dairy cattle. Changes in TAG profile caused by heat stress are expected to modify the physical properties of milk fat, whereas the


reduction of phospholipids may affect the nutritional value of milk. The results are discussed in relation to animal metabolism adaptation in the event of acute heat stress. SIMILAR CONTENT


BEING VIEWED BY OTHERS CARRY-OVER EFFECTS OF DRY PERIOD HEAT STRESS ON THE MAMMARY GLAND PROTEOME AND PHOSPHOPROTEOME IN THE SUBSEQUENT LACTATION OF DAIRY COWS Article Open access 22 April


2022 REVEALING THE DIVERSITY OF INTERNAL BODY TEMPERATURE AND PANTING RESPONSE FOR FEEDLOT CATTLE UNDER ENVIRONMENTAL THERMAL STRESS Article Open access 25 March 2023 A POSTBIOTIC FROM


_ASPERGILLUS ORYZAE_ ATTENUATES THE IMPACT OF HEAT STRESS IN ECTOTHERMIC AND ENDOTHERMIC ORGANISMS Article Open access 19 March 2021 INTRODUCTION The imposition of heat stress on domestic


production animals is commonplace in tropical, sub-tropical and temperate regions. The long term breeding of dairy cattle for increased milk production has led to higher metabolic heat


generation and therefore increased susceptibility to heat stress1, 2. The upper critical limit of the thermoneutral zone for dairy cattle is between 25 °C and 26 °C and temperature-humidity


index (THI) below 723. Heat stress not only affects the welfare and health of cows, but also their productivity, which in turn increases the herd management cost and reduces the


profitability of dairy farms. The frequency of high-temperature days is expected to increase with climate change4, 5, so mitigating the negative impact of heat stress on dairy production


will be an important task for the industry. Recent studies demonstrated that selection for heat-tolerant dairy cow genotypes is feasible and results in improvements in milk production and


feed intake during and after heat stress events5, 6. Under heat stress conditions, some physiological changes can be observed in dairy cattle, such as increased respiration rate, heart rate


and core body temperature1, 7, but the most prominent changes in relation to animal performance are the reduction in dry matter intake and milk yield2, 7,8,9. Heat dissipation associated


metabolic adaptation is an energy expensive process and is believed to be responsible for a proportion of the decline in milk production10. Apart from milk yield, the composition of milk as


influenced by heat stress has been investigated in numerous studies. Most studies found heat stress was associated with a decline in total protein content and total fat content8, 11,12,13.


However, some authors found no significant decrease in fat percentage for cows under heat stress14,15,16. Lipid is one of the main components of milk. The dominant fraction of milk fat is


TAG (about 98%) present in the form of fat globules17. In addition to being an energy source, TAG composition is implicated in human health and the property of dairy products18, 19. The


second most important fraction of milk fat is polar lipids, which are the main structural constituents of fat globule membrane and thus play a role of emulsifier ensuring the stability of


milk emulsion system20, 21. The principal classes of milk polar lipids include phosphatidylcholine (PC), phosphatidylethanolamine (PE), phosphatidylserine (PS), phosphatidylinositol (PI),


sphingomyelin (SM), lysophosphatidylcholine (LPC), lactosylceramide (LacCer) and glucosylceramide (GluCer)18. Besides their functional property, polar lipids in particular SM have other


beneficial effects on human health. For example, SM can reduce cholesterol absorption22; a recent study has also demonstrated a link between SM and infant cognitive development23. Despite


the importance of milk lipids in human health and dairy product property, compared with total fat content, how lipid composition is affected by heat stress has attracted much less attention.


Through correlation analysis between meteorological data and fatty acid (FA) traits predicted by mid-infrared spectra, it was found that the increase of THI was associated with a decrease


in content of short-chain and medium-chain FA, and an increase in that of long-chain FA12. A similar observation was also made for milk from heat stressed cows in an earlier study16.


However, these results were not always conclusive, since the effect of heat stress was often confounded with that of different feeding patterns across seasons. As for the effect of heat


stress on TAG profile and polar lipid composition, no information was available. We hypothesised that heat stress would modify the profile of both TAG and polar lipids of milk through


changing FA pool and animal metabolism. This alteration in lipid composition may have an impact on calf nutrition, human nutrition as well as the processing and quality of dairy products. A


simulated heat-stress experiment was conducted in controlled-climate chambers recently with the aim to validate the genomic selection approach in heat tolerance improvement of dairy cattle5.


In addition to the monitoring of milk yield, dry matter intake, respiration rate and body temperature, these authors also determined the milk composition (lactose, protein and fat %) using


infrared-based technique5. As part of a collaborative project, an aliquot of milk samples from the same experiment was brought to our laboratory and subjected to detailed lipid composition


analysis. We describe here the methodologies for detailed lipidomic analysis using LC-MS as well as the lipid composition change in heat stress conditions for both heat tolerant and heat


susceptible cows; 58 most abundant TAG groups and eight classes of polar lipids were surveyed at the molecular species level. RESULTS AND DISCUSSION EFFECT OF HEAT STRESS ON TAG PROFILE OF


MILK An exhaustive LC separation combined with the use of MS/MS spectrum-based automated structural assignment software (LipidSearch) has been adopted for lipid species identification in


this study. If a TAG group is defined as a series of species having the same total acyl carbon number (CN) and the same number of total double bonds (DB), 94 TAG groups comprising over 400


species (positional isomers not counted) have been identified in raw milk. Among them, 58 most abundant TAG groups were chosen for detailed study; the chemical formula, accurate mass as well


as the FA composition of the main isomers (based on abundance) of these TAG groups are summarised in Table 1. Even using two LC columns and a 90-min gradient elution, complete


chromatographic separation of all isomer species within most TAG groups was not achieved. For example, in the case of TAG 26:0 group, among the five species identified by LipidSearch, the


three main isomers were not resolved completely (Fig. S1B, Supporting Information), so reliable quantification at the species (isomer) level is not currently feasible. Consequently, our


survey on TAG abundance was conducted at the group level using a short LC method, which allowed all main isomers to be eluted as a single peak and enabled also a higher throughput sample


analysis (Fig. S1A, Supporting Information). A preliminary survey using samples collected at five time points from one cohort of six cows showed that the overall TAG profile, as revealed by


PCA, shifted noticeably from the baseline after 2 days’ heat challenge and a complete separation from the baseline was observed after 4 days’ heat challenge, whereas no remarkable difference


in TAG profile was found between pre-heat challenge (baseline) and post-heat challenge (recovery) samples (Fig. 1). As a result, further analysis on all the 5 cohorts of 30 cows was focused


only on the comparison between baseline and D4 heat challenge samples. Due to the large difference in abundance across the 58 TAG groups, the effect of heat stress on the abundance of each


TAG group is presented as the abundance ratio between D4 heat stress and baseline (control) samples (Fig. 2). It appears that the effect of heat stress on the relative abundance of TAG


groups varies mostly with their CN. The heat challenge caused a significant reduction of TAG groups containing 26–48 CN regardless of the level of unsaturation, whereas it increased the


level of those with 51–56 CN, irrespective of the number of BD (Fig. 2). The TAG groups showing the greatest reduction are TAG 28:0 and TAG 30:0 (by about 40% for both), and those showing


the greatest increase are TAG 54:3 and TAG 54:2 (by about 50% for both). For all TAG groups with 49–50 CN, heat stress had little impact. It is interesting to note that three TAG groups


containing a FA with odd number of carbons (TAG 45:2, TAG 45:1 and TAG 45:0) showed a significant reduction in heat stressed cows (Fig. 2). At the group level, TAG 28:0 sustained the


greatest reduction by heat stress, but whether all isomer species of this group are affected at the same magnitude remains a question, given that each TAG group contains on average 4–5


isomer species. The long LC separation method allowed us to resolve chromatographically 6 species from this TAG group with different FA composition, and also revealed that one of the minor


species TAG 4:0/6:0/18:0 was not reduced at the same level as other species by heat stress (Fig. S2, Supporting Information). It is well established that most of the C4:0 to C14:0 and almost


half of the C16:0 FA in milk are synthesized _de novo_ in the mammary gland, whereas the rest of the C16:0 and approximately all long-chain FA originate from blood lipids24, 25. If we


classify the major FA into 3 categories, namely short chain FA (SCFA, C4-C10), medium-chain FA (MCFA, C12-C16) and long-chain FA (LCFA, ≥C18), the main isomer species of the 58 TAG groups


feature a number of configurations, _i.e_. S-S-M, S-M-M, S-M-L, S-L-L, M-L-L and L-L-L (Table 1). TAG groups that showed a significant reduction after heat stress belong mainly to the first


3 types of configurations, which are composed of predominantly SCFA and MCFA, whereas those that are induced by heat stress are from the last 2 types of configurations containing mainly


LCFA. Global FA composition analysis of milk samples from three cohorts (18 animals) by GC-MS revealed that the level of C4:0-C15:0 was indeed significantly reduced, whereas that of LCFA


especially C18:0, C18:1 and C18:2 significantly increased by the heat challenge (Table 2). These results are in agreement with the data reported by Lacetera _et al_.16 and Hammami _et


al_.12, who described a decrease in content of SCFA and MCFA and an increase in that of LCFA from heat-stressed cows. It is interesting to mention that the ratios of C14:1 to C14:0, C16:1 to


C16:0 and C18:1 to C18:0 were very close between baseline and D4 heat stress samples, implying the activity of Δ9-desaturase was not affected by heat stress. As _de novo_ synthesis of


C4-C14 FA is from acetate and butyrate, which are generated in the rumen by fermentation of feed components17, heat stress-induced reduction of C4:0-C14:0 FA-rich TAG groups may partly


result from lowered dry matter intake as recorded by Garner _et al_.5 in this same experiment, which led to suppressed ruminal fermentation and consequently reduced supply of acetate and


butyrate. However, diminished anabolic activities including synthesis of fatty acids in mammary gland due to partitioning of energy in thermoregulation-related processes in response to heat


stress may also be responsible1, 4. So the reduction of TAG groups with 26–48 CN can be attributed to the reduced synthesis of _de novo_ FA under heat stress. This could also explain the


reduction of TAG 45:2, TAG 45:1 and TAG 45:0 which contain C15:0, a FA thought to be synthesised by the bacteria in the rumen26. By contrast, half of C16:0 and all of LCFA are not


synthesised in the mammary gland, but originate from dietary lipids and from lipolysis of adipose tissue TAG25, so the supply of these FA may not be affected or even be induced by heat


stress. Consequently, LCFA-rich TAG species (>50 CN) show an increase in relative proportions. TAG composition is known to affect the property of dairy product. For example, the spread


ability of butter produced in winter is much lower than that of summer butter, due to the higher level of unsaturated FA in summer milk27, 28. A combination of reduced proportion of


SCFA-rich TAG and an increased proportion of LCFA-rich TAG is expected to increase the viscosity and melting point of milk fat produced by heat stressed cows, which may have impact on the


processing and property of dairy products. The 5 cohorts of 30 cows used in our study were equally distributed for heat sensitivity, _i.e_. 15 were heat tolerant (HT) and 15 were heat


susceptible (HS). When the TAG profile of these two types of cows was compared, no significant difference was found for any of the TAG groups in D4 stress samples, nor in baseline samples


(results not shown). So it appears that the heat sensitivity trait in dairy cows is not associated with the TAG composition of milk fat. EFFECT OF HEAT STRESS ON POLAR LIPID CONTENT Eight


classes of polar lipids were surveyed in this study. Over 100 species were identified by LipidSearch based on MS/MS spectrum, of which 58 most abundant species were selected for relative


quantification; the formula, accurate mass and FA composition of these species are given in Table 3. The absolute content of the 8 classes of polar lipids in raw milk was described in our


previous report29. In this study, only relative quantification was performed at the species level (based on peak area of each species) and at the class level (based on the sum of all species


within the same class). To illustrate the effect of heat stress on the abundance of each polar lipid class, the abundance ratio between D4 heat stress and baseline samples is again


presented. At the class level, a 4-day heat challenge had no significant influence on the abundance of PI, SM and LacCer, but reduced that of PS, PE, PC, LPC and GluCer (Fig. 3). While only


a slight decrease was observed for PS (9%) and PE (11%), a substantial reduction was recorded for PC (17%) and GluCer (21%). However, the most remarkable change in polar lipid content


brought about by heat stress was detected with LPC, for which a dramatic reduction (52%) was observed after 4 days’ heat stress treatment (Fig. 3). It is worth mentioning that the abundance


of all the five species within this class was dramatically reduced (results not shown). A detailed survey of LPC level was conducted for all the 5 time points of the experiment. Indeed, the


level of LPC was significantly decreased (by 34%) after 2 days’ heat challenge, which dropped further after 4 days’ treatment and then rose to a normal level after 5 days’ recovery (Fig. 


4A). In addition to the dramatic decrease after 2 days’ heat treatment, LPC displayed a differential response between HT and HS cows in the event of heat stress. Figure 4B shows that the LPC


level is not significantly different between HT and HS cows at the baseline point, but the heat induced reduction of this polar lipid class is much less with HT as compared to HS animals


regardless of the duration of heat challenge. This LPC changing pattern is similar to the results reported by Garner _et al_.5, who observed a lower body temperature, a higher dry matter


intake and a higher milk yield of HT cows as compared to HS cows during the heat challenge. Consequently, LPC can be regarded as a heat stress biomarker for dairy cattle. In this regard, it


would be interesting to investigate the LPC level after a short period (for example, a few hours) of heat stress to further validate this finding. The effect of heat stress on polar lipid


composition has not been extensively investigated. A recent study by Tian _et al_.30 using untargeted metabolic profiling approach identified a large number of potential heat stress


biomarkers including some lipid species, but one single LPC species included in that list was found to be induced in heat stressed cows. This discrepancy may be attributable to different


experimental designs, ambient conditions (baseline level) _versus_ controlled environment (heat stress) in our experiment as compared to summer samples _versus_ spring samples in their


study. Compared with TAG, polar lipids are minor component in milk fat (<1%). However, our study revealed that one polar lipid class LPC is much more sensitive to heat stress than TAG as


a whole, since maximum reduction observed on TAG 28:0 and TAG 30:0 was 40% after a 4-day heat challenge, as compared to over 50% reduction in the case of LPC. While the heat-induced decrease


of some TAG groups appears to be caused by a decline in _de novo_ FA, this same mechanism solely cannot explain the reduction of LPC which contains only C16 and LCFA in the five species


surveyed (Table 3). Some up-stream steps in the biosynthesis pathways, such as the formation of choline/phosphocholine may be inhibited under heat stress. Clearly, further investigation is


needed to understand the underlying mechanisms. It is also worth noting that the simulated heat challenge in the current experiment is of moderate magnitude (maximum temperature 33 °C and


THI up to 84). How milk lipid composition changes in more severe heat stress conditions remains to be determined. We have found that 5 out of the 8 polar lipid classes were affected by heat


stress, but the overall level of polar lipids is not expected to be remarkably reduced, since LPC and GluCer are minor classes, whereas PS, PE and PC were only slightly and moderately


suppressed. In addition, our previous study demonstrated that the major classes of polar lipids (PC, PE, PS, PI and SM) in milk were correlated with each other29. Such an inter-class


correlation was not affected by heat stress. As an example, a rather similar correlation level was observed between PE and PC in D4 stress samples as compared to the baseline samples (Fig. 


S3, Supporting Information). It should be pointed out that our simulated heat stress in a controlled climate chamber can be considered as an acute heat stress, because it was conducted in


late winter/early spring with an ambient temperature below 20 °C and THI below 70. The effect of such an acute heat stress on animal metabolism may be different from the seasonal long term


heat stress experienced by animals in summer, as in the latter case heat load imposed on dairy cows is gradual but prolonged. A parallel study on milk lipid composition as influenced by a


chronic heat stress is still lacking. CONCLUSION In conclusion, heat stress could alter the TAG profile of milk, which is featured by a decrease of TAG groups containing predominantly SCFA


to MCFA and a concomitant increase of those containing mostly LCFA. The heat challenge also significantly reduced the level of some polar lipid classes especially LPC, which appears to be a


lipid marker for heat stress in dairy cattle. While the TAG profile was modified in a similar manner for both HT and HS cows, a higher level of LPC was detected for HT cows as compared to HS


cows during the heat stress. MATERIALS AND METHODS COWS AND EXPERIMENTAL DESIGN The experiment received animal ethics approval from the Agricultural Research and Extension Animal Ethics


Committee of the Department of Economic Development, Jobs, Transports and Resources, Victoria, Australia. All methods were performed by approved staff members in accordance with the relevant


standard operating procedures approved by the above mentioned ethics committee. Detailed information on animals and experimental design were described by Garner _et al_.5. In brief, 48


primiparous Holstein-Friesian cows were used in this experiment; 24 cows were genomically selected to be heat tolerant (HT) and 24 cows heat susceptible (HS). At the beginning of the


experiment, mean DIM was 67 for HT cows and 68 for HS cows and mean body weight was 477 kg for HT cows and 488 kg for HS cows. The total duration of the experiment was 25 consecutive days


including a 7 day baseline measurement period outside in ambient conditions, a 4 day heat challenge in controlled-climate chambers, and a 14 day recovery period outside in ambient


conditions. The 48 cows were introduced to controlled climate chambers in 8 cohorts. Temperature and relative humidity inside the controlled-climate chambers was varied to mimic the diurnal


patterns in heat load imposed on dairy cows that occur during heat-wave events in southern Australia. The conditions in the controlled-climate chambers were designed to remain above THI 72


and not exceed THI 84 to impose a moderate level of heat stress. The climatic conditions programmed into the control system were 25 °C and 60% RH (THI 74) between 6 pm and 6 am, 30 °C and


50% RH (THI 80) between 6 am and 12 noon, and 33 °C and 50% RH (THI 84) between 12 noon and 6 pm. The 12 hour light and 12 hour dark cycle was controlled manually. The animal diet during the


experimental period was described in Garner _et al_.5. Milk samples were collected at five time points during the experiment period, on day 3 of the baseline period (baseline), on day 2 and


4 of the heat challenge in controlled-climate chambers (D2 stress and D4 stress), and at day 5 and 12 of the recovery period (D5 recovery and D12 recovery). Cows were milked twice daily


(6:00 am and 3:00 pm) and the samples from afternoon milking were analysed for lipid composition. The milk samples were transported to the laboratory on ice and stored in −80 °C. A subset of


samples from 5 randomly selected cohorts totalling 30 cows (15 HT and 15 HS) were analysed to determine the effect of heat stress on lipid composition of raw milk. CHEMICALS AND REAGENTS


One TAG species (TAG tri-20:1) used as internal standard for TAG analysis was purchased from Sigma-Aldrich. One PS species (PS 17:0/17:0) used as internal standard for polar lipid analysis


was purchased from Avanti Lipids. Solvents used for lipid extraction and mobile phase preparation were of chromatographic grade and were from Merck (methanol, butanol and acetonitrile) and


Sigma-Aldrich (chloroform and isopropanol). Ammonium formate, used as mobile phase additive, was of analytical grade (Sigma-Aldrich). LIPID EXTRACTION FROM MILK FOR LC-MS ANALYSIS Raw milk


samples were diluted by adding 2 volumes of Milli-Q water. Lipid extraction from the diluted milk samples was conducted using the one phase method recently developed31. Briefly, one mL of


lipid extraction mix (butanol/methanol/chloroform at a 3:5:4 ratio) was added to 100 µL of diluted milk. The mixture was shaken by vortex for 20 s, sonicated for 20 min and then centrifuged


for 15 min (15000 _g_). The supernatant was transferred to an injection vial and analysed directly by LC-MS. LC-MS METHOD Chromatographic separation for TAG identification was achieved using


two Poroshell 120 EC-C18 columns (150 × 4.6 mm, 2.7 µm, Agilent Technologies) connected in series on an Agilent 1290 Infinity HPLC system. The column compartment was maintained at 40 °C and


the auto-sampler at 12 °C. The mobile phase was composed of acetonitrile/water (60:40, v/v) containing 10 mM ammonium formate (A) and acetonitrile/isopropanol (10:90, v/v) containing 10 mM


ammonium formate (B). The flow rate was 0.5 mL/min with a gradient elution of 60 to 100% B over 90 min. The injection volume was 4 µL. Chromatographic separation for TAG quantification was


conducted using an Acquity UPLC HSS T3 column (100 × 2.1 mm, 1.8 µm, Waters) on the same Agilent HPLC system. The column compartment was maintained at 50 °C and the auto-sampler at 12 °C.


The mobile phase was composed of acetonitrile/water (60:40, v/v) containing 10 mM ammonium formate (A) and acetonitrile/isopropanol (10:90, v/v) containing 10 mM ammonium formate (B). The


flow rate was 0.28 mL/min with a gradient elution of 20 to 100% B over 20 min. The injection volume was 2 µL. Chromatographic separation for polar lipid identification and quantification was


performed using a Luna HILIC column (250 × 4.6 mm, 5 µm, Waters) on the aforementioned HPLC system. The column compartment was maintained at 30 °C and the auto-sampler at 12 °C. The mobile


phase was composed of 5 mM aqueous ammonium formate (A) and acetonitrile +0.1% formic acid (B). The flow rate was 0.6 mL/min with a gradient elution of 2 to 21% A over 25 min. The injection


volume was 5 µL. The detection of lipids was by LTQ-Orbitrap Elite mass spectrometer (Thermo Scientific) operated in electrospray ionization positive (for analysis of TAG and most polar


lipid classes) or negative (for analysis of PI) Fourier transform mode. The resolution was set to 60,000 for both positive and negative modes. Identification of lipid species present in milk


was performed based on accurate mass of parent ions (±5 ppm) and product ions (±10 ppm) as well as top5 MS/MS spectra (CE 35) using LipidSearch software (version 4.1, Thermo Scientific)32,


followed by manual verification. Selected lipid species was quantified at a relative scale using peak area of parent ions after normalization by the internal standard. FATTY ACID PROFILING


BY GC-MS For FA composition analysis of milk fat, total lipid was extracted by the method of Bligh and Dyer33. After removal of chloroform under a stream of nitrogen, the transesterification


of the extracted lipids was carried out by adding acidic methanol (6% H2SO4) and heating at 80 °C for 3 h34. The released FA methyl esters (FAMEs) were extracted by hexane and analysed by


GC-MS. The separation of FAMEs was achieved by a BPX-70 column (50 m × 0.22 mm ID, 0.25 µm film thickness, SGE Analytical Science) with a constant flow of 1.0 mL/min helium as carrier gas


and the following oven temperature program: 120 °C to 245 °C ramping at 3 °C/min, with a total run time of 42 min. The injection volume was 1 µL in split mode (1:40). The detection was by an


Agilent 7000 GC/MS Triple Quad with the following settings: scanning mass range of 40–550 amu, transfer line temperature of 240 °C, source temperature of 280 °C, and quad temperature of 150


 °C. A standard mix (C4-C24, Supelco) containing 37 FAMEs was used to provide standard curve for quantification. STATISTICAL ANALYSIS OF DATA All lipid content data were subjected to ANOVA


(XLSTAT, Microsoft Excel); where significant differences were found between treatments, a Tukey’s HSD test was conducted for pairwise comparisons. Principal component analysis (PCA) for


unsupervised classification of samples was performed with MATLAB R2014a (MathWorks, Natick, MA) utilising PLS Toolbox (Eigenvector Research, Manson, WA). REFERENCES * Kadzere, C. T., Murphy,


M. R., Silanikove, N. & Maltz, E. Heat stress in lactating dairy cows: a review. _Livestock Production Science_ 77, 59–91 (2002). Article  Google Scholar  * West, J. W. Effects of


heat-stress on production in dairy cattle. _J Dairy Sci_ 86, 2131–44 (2003). Article  CAS  PubMed  Google Scholar  * Berman, A. _et al_. Upper critical temperatures and forced ventilation


effects for high-yielding dairy cows in a subtropical climate. _J Dairy Sci_ 68, 1488–95 (1985). Article  CAS  PubMed  Google Scholar  * Das, R. _et al_. Impact of heat stress on health and


performance of dairy animals: A review. _Vet World_ 9, 260–8 (2016). Article  ADS  PubMed  PubMed Central  Google Scholar  * Garner, J. B. _et al_. Genomic Selection Improves Heat Tolerance


in Dairy Cattle. _Sci Rep_ 6, 34114 (2016). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  * Nguyen, T. T., Bowman, P. J., Haile-Mariam, M., Pryce, J. E. & Hayes, B. J.


Genomic selection for tolerance to heat stress in Australian dairy cattle. _J Dairy Sci_ 99, 2849–62 (2016). Article  CAS  PubMed  Google Scholar  * Wheelock, J. B., Rhoads, R. P., Vanbaale,


M. J., Sanders, S. R. & Baumgard, L. H. Effects of heat stress on energetic metabolism in lactating Holstein cows. _J Dairy Sci_ 93, 644–55 (2010). Article  CAS  PubMed  Google Scholar


  * Bouraoui, R., Lahmar, M., Majdoub, A., Djemali, Mn. & Belyea, R. The relationship of temperature-humidity index with milk production of dairy cows in a Mediterranean climate. _Anim.


Res._ 51, 479–491 (2002). Article  Google Scholar  * Rhoads, M. L. _et al_. Effects of heat stress and plane of nutrition on lactating Holstein cows: I. Production, metabolism, and aspects


of circulating somatotropin. _J Dairy Sci_ 92, 1986–97 (2009). Article  CAS  PubMed  Google Scholar  * Bernabucci, U., Bani, P., Ronchi, B., Lacetera, N. & Nardone, A. Influence of


short- and long-term exposure to a hot environment on rumen passage rate and diet digestibility by Friesian heifers. _J Dairy Sci_ 82, 967–73 (1999). Article  CAS  PubMed  Google Scholar  *


Bernabucci, U. _et al_. Effect of summer season on milk protein fractions in Holstein cows. _J Dairy Sci_ 98, 1815–27 (2015). Article  CAS  PubMed  Google Scholar  * Hammami, H. _et al_.


Genetic analysis of heat stress effects on yield traits, udder health, and fatty acids of Walloon Holstein cows. _J Dairy Sci_ 98, 4956–68 (2015). Article  CAS  PubMed  Google Scholar  *


Hill, D. L. & Wall, E. Dairy cattle in a temperate climate: the effects of weather on milk yield and composition depend on management. _Animal_ 9, 138–49 (2015). Article  CAS  PubMed 


Google Scholar  * Roman-Ponce, H., Thatcher, W. W., Buffington, D. E., Wilcox, C. J. & Van Horn, H. H. Physiological and Production Responses of Dairy Cattle to a Shade Structure in a


Subtropical Environment. _J Dairy Sci_ 60, 424–430 (1977). Article  Google Scholar  * Knapp, D. M. & Grummer, R. R. Response of lactating dairy cows to fat supplementation during heat


stress. _J Dairy Sci_ 74, 2573–9 (1991). Article  CAS  PubMed  Google Scholar  * Lacetera, N., Bernabucci, U., Ronchi, B. & Nardone, A. Physiological and productive consequences of heat


stress: the case of dairy ruminants. In _Interactions between climate and animal production_ (eds. Lacetera, N. _et al._) 45–59 (Wageningen Academic Publishers, 2003). * Mansson, H. L. Fatty


acids in bovine milk fat. _Food Nutr Res_ 52 (2008). * Jensen, R. G. The composition of bovine milk lipids: January 1995 to December 2000. _J Dairy Sci_ 85, 295–350 (2002). * Palmquist, D.


L. Milk Fat: Origin of Fatty Acids and Influence of Nutritional Factors Thereon. In _Advanced Dairy Chemistry Volume 2 Lipids_ (eds Fox, P. F. & McSweeney, P. L. H.) 43–92 (Springer US,


Boston, MA, 2006). * Fong, B. Y., Norris, C. S. & MacGibbon, A. K. H. Protein and lipid composition of bovine milk-fat-globule membrane. _International Dairy Journal_ 17, 275–288 (2007).


Article  CAS  Google Scholar  * Sánchez-Juanes, F., Alonso, J. M., Zancada, L. & Hueso, P. Distribution and fatty acid content of phospholipids from bovine milk and bovine milk fat


globule membranes. _International Dairy Journal_ 19, 273–278 (2009). Article  Google Scholar  * Rombaut, R. & Dewettinck, K. Properties, analysis and purification of milk polar lipids.


_Int Dairy J_ 16, 1362–1373 (2006). Article  CAS  Google Scholar  * Tanaka, K. _et al_. The pilot study: sphingomyelin-fortified milk has a positive association with the neurobehavioural


development of very low birth weight infants during infancy, randomized control trial. _Brain Dev_ 35, 45–52 (2013). Article  CAS  PubMed  Google Scholar  * Chilliard, Y., Ferlay, A.,


Mansbridge, R. M. & Doreau, M. Ruminant milk fat plasticity: nutritional control of saturated, polyunsaturated, trans and conjugated fatty acids. _Ann. Zootech._ 49, 181–205 (2000).


Article  CAS  Google Scholar  * Parodi, P. W. Milk fat in human nutrition. _Australian J. of Dairy Technol_ 59 (2004). * German, J. B. & Dillard, C. J. Composition, structure and


absorption of milk lipids: a source of energy, fat-soluble nutrients and bioactive molecules. _Crit Rev Food Sci Nutr_ 46, 57–92 (2006). Article  CAS  PubMed  Google Scholar  * Collomb, M.


_et al_. Seasonal variation in the fatty acid composition of milk supplied to dairies in the mountain regions of Switzerland. _Dairy Science & Technology_ 88, 631–647 (2008). Article 


CAS  Google Scholar  * Ozcan, T., Yaslioglu, E., Kilic, l & Simsek, E. The influence of the season and milking time on the properties and the fatty acid composition of the milk in


different dairy cattle farms. _Mljekarstvo_ 64, 9–17 (2015). Article  Google Scholar  * Liu, Z., Moate, P., Cocks, B. & Rochfort, S. Comprehensive polar lipid identification and


quantification in milk by liquid chromatography-mass spectrometry. _J Chromatogr B Analyt Technol Biomed Life Sci_ 978–979, 95–102 (2015). Article  PubMed  Google Scholar  * Tian, H. _et


al_. Integrated Metabolomics Study of the Milk of Heat-stressed Lactating Dairy Cows. _Sci Rep_ 6, 24208 (2016). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  * Liu, Z.,


Rochfort, S. & Cocks, B. G. Optimization of a single phase method for lipid extraction from milk. _J Chromatogr A_ 1458, 145–9 (2016). Article  CAS  PubMed  Google Scholar  * Yamada, T.


_et al_. Development of a lipid profiling system using reverse-phase liquid chromatography coupled to high-resolution mass spectrometry with rapid polarity switching and an automated lipid


identification software. _J Chromatogr A_ 1292, 211–8 (2013). Article  CAS  PubMed  Google Scholar  * Bligh, E. G. & Dyer, W. J. A rapid method of total lipid extraction and


purification. _Can J Biochem Physiol_ 37, 911–7 (1959). Article  CAS  PubMed  Google Scholar  * Taylor, A. G. & Savage, C. Fatty acid composition of New Zealand green-lipped mussels,


Perna canaliculus: Implications for harvesting for n-3 extracts. _Aquaculture_ 261, 430–439 (2006). Article  CAS  Google Scholar  Download references ACKNOWLEDGEMENTS Thanks to Ben Hayes,


Coralie Reich and Thuy Nguyen for assisting in establishing the animal resources used in this study. AUTHOR INFORMATION AUTHORS AND AFFILIATIONS * Biosciences Research, Agriculture Victoria,


AgriBio, 5 Ring Road, Bundoora, Victoria, 3083, Australia Z. Liu, V. Ezernieks, J. Wang & N. Wanni Arachchillage * Farming Systems Research, Agriculture Victoria, Ellinbank Centre, 1301


Hazeldean Rd, Ellinbank, Victoria, 3821, Australia J. B. Garner & W. J. Wales * School of Applied Systems Biology, La Trobe University, Bundoora, Victoria, 3083, Australia B. G. Cocks 


& S. Rochfort Authors * Z. Liu View author publications You can also search for this author inPubMed Google Scholar * V. Ezernieks View author publications You can also search for this


author inPubMed Google Scholar * J. Wang View author publications You can also search for this author inPubMed Google Scholar * N. Wanni Arachchillage View author publications You can also


search for this author inPubMed Google Scholar * J. B. Garner View author publications You can also search for this author inPubMed Google Scholar * W. J. Wales View author publications You


can also search for this author inPubMed Google Scholar * B. G. Cocks View author publications You can also search for this author inPubMed Google Scholar * S. Rochfort View author


publications You can also search for this author inPubMed Google Scholar CONTRIBUTIONS Z.L. performed all LC-MS analysis and wrote the manuscript. V.E. conducted the fatty acid analysis by


GC-MS. J.W. and N.W.A. contributed to the sample preparation and data analysis. J.B.G. and W.J.W. managed all animal experimental work. B.G.C. and S.R. supervised the study and were


responsible for the project development and funding acquisition. CORRESPONDING AUTHOR Correspondence to Z. Liu. ETHICS DECLARATIONS COMPETING INTERESTS The authors declare that they have no


competing interests. ADDITIONAL INFORMATION PUBLISHER'S NOTE: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


ELECTRONIC SUPPLEMENTARY MATERIAL SUPPLEMENTARY INFORMATION: HEAT STRESS IN DAIRY CATTLE ALTERS LIPID COMPOSITION OF MILK 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 license, and indicate if changes were made. The images or other third party material in this article are


included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license 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


license, visit http://creativecommons.org/licenses/by/4.0/. Reprints and permissions ABOUT THIS ARTICLE CITE THIS ARTICLE Liu, Z., Ezernieks, V., Wang, J. _et al._ Heat Stress in Dairy


Cattle Alters Lipid Composition of Milk. _Sci Rep_ 7, 961 (2017). https://doi.org/10.1038/s41598-017-01120-9 Download citation * Received: 06 December 2016 * Accepted: 22 March 2017 *


Published: 19 April 2017 * DOI: https://doi.org/10.1038/s41598-017-01120-9 SHARE THIS ARTICLE Anyone you share the following link with will be able to read this content: Get shareable link


Sorry, a shareable link is not currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt content-sharing initiative