
Heat stress in dairy cattle alters lipid composition of milk
- Select a language for the TTS:
- UK English Female
- UK English Male
- US English Female
- US English Male
- Australian Female
- Australian Male
- Language selected: (auto detect) - EN
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