Milk and dairy products obtained from small ruminants such as sheep and goats are gaining consumer interest for their nutritional value and for being a source of vitamins, fats, amino acids and probiotic compounds. In the Mediterranean area, sheep feeding managements are optimized to sustain milk production and the feeding is traditionally based on the extensive and semi-extensive exploitation of natural and/or cultivated pastures. Typically the chosen feeding system is based on environmental and climatic considerations as well as land availability and profitability factors. In Italy, and in particular in Sardinia, pasture feeding is anciently practiced thanks to the temperate climate and fertile soils factors that allow sustainable and competitive grass production. However, pasture availability is seasonal and often scarce, as in many Mediterranean areas, making supplementation necessary, at least at some times during the productive season (Molle et al., Reference Molle, Cabiddu, Decandia, Acciaro, Scanu, Addis, Fiori and Caredda2021). A drawback of the pasture-based feeding system is the seasonality of production (Scano, Reference Scano2019).
Milk composition is known to be dependent on animal breed, stage of lactation and, in particular, the feeding system. Grazing systems relying on semi-natural pasture vegetation and forage species and their phenological stage are known to strongly affect the fatty acid profile of milk (Addis et al., Reference Addis, Cabiddu, Pinna, Decandia, Piredda, Pirisi and Molle2005; Morand-Fehr et al., Reference Morand-Fehr, Fedele, Decandia and Le Frileux2007; Buccioni et al., Reference Buccioni, Decandia, Minieri, Molle and Cabiddu2012; Renna et al., Reference Renna, Ferlay, Lussiana, Bany, Graulet, Wyss, Enri, Battaglini and Coppa2020) and the functional and sensory characteristics of dairy products (Schönfeldt, Hall, and Smit, Reference Schönfeldt, Hall and Smit2012). The effects of different forage diets (grazing vs. hay) on the fatty acid profile and gene expression in the mammary gland of the Churra Tensina sheep breed were studied (Dervishi et al., Reference Dervishi, Joy, Sanz, Alvarez-Rodriguez, Molino and Calvo2012). These authors reported that the forage type affected the levels conjugated linoleic acid (CLA) and long-chain saturated fatty acid (LCFA) content, with higher percentages in grazing when compared to hay feeding animals.
The effects of different feeding systems cannot be studied solely from the analysis of milk gross compositional traits or fatty acid profile but rather by means of chemical fingerprinting obtained from hyphenated analytical techniques combined with chemometrics approaches. Lipidomics is an -omics science that allows the investigation of the structural and functional complexity of lipids in biological systems. Liquid chromatography and high-resolution mass spectrometry are the analytical approaches most used in lipid research. With this analytical approach and in particular with MS/MS experiments important structural elucidation of complex lipids such triacylglycerols (TG), phosphatidylcholines (PC), phosphatidylethanolamine (PE), phosphatidylserine (PS), ceramides and sphingomyelins (SM) can be obtained. Recently, the change of goat lipidome and compositional characteristics of milk from two feeding treatments, grazing vs. hay-feeding in confinement, were reported (Argov-Argaman et al., Reference Argov-Argaman, Glasser, Muklada, Hadaya, Mesilati-Stahy, Raz and Yan Landau2021). Milk fat and protein contents were higher in grazing goats. Except for saturated fatty acids, dietary management affected fatty acid composition, in particular polyunsaturated fatty acids content. Omega 6 and omega 3 fatty acids increased and omega6/omega3 ratio decreased in goat milk under grazing when compared with confined managements (Argov-Argaman et al., Reference Argov-Argaman, Glasser, Muklada, Hadaya, Mesilati-Stahy, Raz and Yan Landau2021).
In this work with a lipidomics approach we studied the impact of two different diets on the sheep milk lipidome. Animal diets consisted of 40% and 60% of ewe energy requirement derived from pasture integrated with hay and concentrate. Bulk milk samples were analyzed by a high-resolution mass spectrometric (UHPLC-QTOF-MS) platform. MS data were submitted to multivariate statistical analysis to differentiate the two diet regimens and to identify discriminant milk lipid compounds. MS/MS data were used for the annotation and confirmation purposes of complex lipids.
Materials and methods
Animals and diets
Forty-eight mature Sarda ewes were divided into two groups balanced for age, body weight, body condition score, lambing date and milk production and randomly assigned to the diets. The diets were either 40% (P40) or 60% (P60) of ewe energy requirement derived from pasture, with 4 and 6 h/day of grazing, respectively. Pasture was based on Persian clover, sulla and burr medic cultivated in mixture with Italian ryegrass. Diet supplementation consisted of maize and faba beans at 600 and 300 g/head day split into two meals at milking, and overnight 800 and 600 g/head day of a ryegrass or hay, for P40 and P60 respectively. After the morning milking the groups were at pasture while for the remaining time they were kept indoors in separate pens. Samples of bulk milk were collected for 10 weeks, from March to June 2020, for 3 d every two weeks for a total of 15 bulk milk samples for each group obtained by pooling milk from 24 ewes.
Sample preparation
150 μl of milk samples were transferred to Eppendorf tubes containing 10 μl of the internal mixture of standards (Splash, Lipidomics, Sigma Aldrich Milan, Italy) and added with 525 μl of methanol and 525 μl of MTBE. Samples were then centrifuged at 4000 rpm for 15 min. Subsequently, a second round of centrifuge was performed at 12 000 rpm for 5 min. Before transferring to autosampler vials, the supernatant was filtered through a 0.22 μm MS nylon syringe filter.
Lipidomics analysis
The supernatant of milk samples was analyzed with a LC-QTOF-MS coupled with an Agilent 1290 Infinity II LC system. An aliquot of 1.0 μl from each sample was injected in a Kinetex 5 μm EVO C18 100 A, 150 mm × 2.1 μm column (Agilent Technologies, Palo Alto, CA). The column was maintained at 50 °C at a flow rate of 0.4 ml/min. The mobile phase for positive ionization mode consisted of (A) 10 mm ammonium formate solution in 60% of milliQ water and 40% of acetonitrile and (B) 10 mm ammonium formate solution containing 90% of isopropanol and 10% of acetonitrile. In positive ionization mode, the chromatographic separation was obtained with the following gradient: initially 60% of A, then a linear decrease from 60% to 50% of A in 2 min then at 1% in 5 min staying at this percentage for 1.9 min and then brought back to the initial conditions in 1 min. The mobile phase for negative ionization mode differed only for the use of 10 mm ammonium acetate instead of ammonium formate (A and B). The analytical setup used was equipped with an Agilent jet stream technology source which was operated in both positive and negative ion modes with the following parameters: gas temperature, 200 °C; gas flow (nitrogen) 10 l/min; nebulizer gas (nitrogen), 50 psig; sheath gas temperature, 300 °C; sheath gas flow, 12 l/min; capillary voltage 3500 V for positive and 3000 V for negative; nozzle voltage 0 V; fragmentor 150 V; skimmer 65 V, octapole RF 7550 V; mass range, 40–1700 m/z; capillary voltage, 3.5 kV; collision energy 20 eV in positive and 25 eV in negative mode, mass precursor per cycle = 3; threshold for MS/MS 5000 counts.
Data analysis
The LCQTOF-MS data were uploaded to the web platform XCMS (Tautenhahn et al., Reference Tautenhahn, Patti, Rinehart and Siuzdak2012) and the resulting data matrix was submitted to multivariate statistical analysis as implemented in SIMCA-P+ software (version 14.1, Umetrics, Umeå, Sweden). Prior to this, QTOF-MS features were mean centered and scaled to unit variance column-wise. Principal component analysis, partial least-squares-discriminant analysis (PLS-DA) and its orthogonal variant (OPLS-DA) were performed (Manis et al., Reference Manis, Scano, Nudda, Carta and Pulina2021). The quality of the models was evaluated based on the cumulative parameters R2Y and Q2Y, the latter estimated by cross validation. The variable importance in projection (VIP) scores in the predictive component were analyzed and only those metabolites having VIP values >1 were considered as discriminants between the classes.
Results
Untargeted lipidomics
Gross composition of the P40 and P60 groups of bulk milk samples was studied and the results are reported in online Supplementary Table S1. P60 showed a higher milk yield, but lower proteins and fat content.
Data from XCMS in positive and negative ionization modes were submitted to multivariate statistical analysis. Principal component analysis allowed us to explore the data structure in a visual way to understand peculiar characteristics such as groupings of samples and interrelated variables. From the analysis of the derived score plots of both ionization modes no clusters of milk samples, based on different diets, were observed (online Supplementary Fig. S1). Attempts to differentiate milk samples based on the P40 and P60 groups were carried out by performing a PLS-DA for both positive and negative ionization modes. The P40 and P60 milk samples were correctly classified (statistical parameters: R2Y of 0.96 and Q2Y of 0.45) by the statistical model built using the positive ionization mode results. Unfortunately, the PLS-DA model of the negative ionization mode data was not able to correctly classify samples. Furthermore, an OPLS-DA of positive ionization results was performed with the aim to underline differences between the two classes. The corresponding score plot is shown in Fig. 1, and Table 1 reports the discriminant lipids and their VIP values.
OPLS-DA score plot of milk lipids acquired in the positive ionization mode (P40 green circles, P60 red boxes).
OPLS-DA discriminant lipid compounds with the variable importance in the projection (VIP) values.
After the analysis in the positive ionization mode the P40 group showed higher levels of the following complex lipids: Saturated fatty acid PC such as PC(28:0), PC(30:0), PC(34:0) were found upregulated along with the unsaturated fatty acid PE (PE(35:1), PE(36:1) PE(36:3)) and lyso-derivative compounds LPC(18:1) and LPE(18:1). The most highly discriminant compounds for the P40 group were: sphingomyelins SM(41:2) and SM(42:3) with a VIP value of 3.10 and 3.16 respectively. Among discriminants in the P40 group, the triacylglycerols class showed TG with a carbon number (CN) from TG (40:5) to TG (49:0) and, in detail, the TG (48:3) showed the highest VIP value.
Conversely, the P60 group showed more unsaturated PC such as PC(32:1), PC(33:1), PC(34:3), PC(36:1), PC(36:2), PC(36:3), PC(36:4). Interestingly, in the P60 group we found more TG with a high CN from TG (50:2) to TG (55:1) and a higher overall unsaturation degree (Table 1). The P60 was also characterized by highest level of the following SMs, SM (32:0), SM(36:1), SM(39:1) and SM(40:0).
Discussion
From the UHPLC-QTOF-MS analysis we detected several lipid classes, and with the OPLS-DA model we were able to annotate discriminant lipids that can be used as potential biomarkers to differentiate different nutritional regimens. From a physiological point of view, milk TG are synthetized from FA that are derived from de novo mammary gland synthesis, dietary lipids, and endogenous fat stores (adipose or hepatic lipids). In our model we noticed that the discriminants TG for the P40 group preferentially incorporated medium and short saturated chain FA (C4–C16) that are produced in the mammary gland by de novo synthesis. On the contrary, TGs of the P60 group incorporated longer chain fatty acid with more unsaturated FA. In the mammary gland, these TGs are synthetized from polyunsaturated FA, that originate from the diet or from lipid mobilization of the adipose tissue, or from the activity of stearoyl desaturase (Leroux et al., Reference Leroux, Le Provost, Petti, Bernard, Chilliard and Martin2003).
Furthermore, the same trend can be observed for PC discriminants. The P40 group is characterized by PC bearing medium chain saturated FA (i.e. lauric acid, myristic acid, and palmitic acid) thus suggesting an increase of the de novo synthesis of FA. Instead, the PC of the P60 group are characterized by unsaturated LCFA. Interestingly, discriminant lipids of the P40 group were lysophosphatidylcholine (18:1) and lysophosphatidylethanolamine (18:1) produced by the hydrolysis of the corresponding precursors PC and PE.
Several studies have shown that the challenges induced by biotic factors (such as changes in vegetation and heat) may be a source of animal stress, eliciting a greater degree of behavioral plasticity and adaptability in grazing animals (Villalba and Manteca, Reference Villalba and Manteca2019). The results obtained from our experiments show a downregulation in the P60 group of different lipid classes such as phosphocoline, lysophosphocoline and phosphoethanolamines. These results are in accordance with Liu's suggesting an intracellular stress condition of the animals under study that is reflected in the chemical composition of the milk. The condition of intracellular stress can inhibit the initial stages of the choline/phosphocholine biosynthesis process (Liu et al., Reference Liu, Ezernieks, Wang, Wanni Arachchillage, Garner, Wales, Cocks and Rochfort2017).
In conclusion, pairwise partial least square discriminant analysis of the lipid profile of ovine milk was carried out to to find discriminant lipids characteristic of lesser (P40) or greater (P60) herbage grazing. The two dietary groups were characterized by differences in triacylglycerols, phosphocholines and phosphatidylethanolamines levels. Discriminants of the P40 group were TG and PC containing in their backbone saturated medium chain FA thus suggesting greater de novo fatty synthesis in the mammary gland. On the other hand, the P60 group was characterized by TG and PC formed by unsaturated long chain FA originating from the diet or from lipid mobilization.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S0022029921000583
Acknowledgments
This work was supported by the Fondazione Banco di Sardegna under the project ‘Innovative antioxidant molecules for the food and health industry’ (CUP F71I17000180002) and the ‘Biomilkchina’ project funded by POR FESR Regione Sardegna. We acknowledge the CeSAR (Centro Servizi d'Ateneo per la Ricerca) of the University of Cagliari, Italy for the IM-QTOF-MS-MS experiments performed with an Agilent 6560.
