Effects of Supplementing Rumen-Protected Methionine and Lysine on Milk Performance and Oxidative Status of Dairy Ewes

There is limited information on the impact of dietary supplementation with separate rumen-protected (RP) amino acids (AA), or with their combination, on ewes’ oxidative status. Sixty ewes were divided into five groups; C: basal diet (control); M: basal diet + 6 g/ewe RP methionine; L: basal diet + 5 g/ewe RP lysine; LML: basal diet + 6 g methionine and 5 g lysine/ewe; and HML: basal diet + 12 g methionine + 5 g lysine/ewe. Milk’s fat content increased in RP-AA fed ewes, while that of protein in M and L only. In blood plasma, the malondialdehyde (MDA) content was reduced in the M, LML, and HML compared to C-fed ewes. An increase in glutathione transferase activity in the blood plasma of the M and LML compared to the C and HML-fed ewes were found. In milk, lower values of the ferric reducing ability of plasma (FRAP) in the LML and HML-fed ewes and of 2,2′-Azino-bis 3-ethylbenzthiazoline-6-sulfonic acid (ABTS) in LML only, were found. Lysine increased milk’s FRAP values and MDA content. Both L and HML diets increased milk’s protein carbonyls content. Methionine improves the organism’s oxidative status, without adversely affecting milk’s oxidative stability. Lysine dietary inclusion affects negatively the oxidative stability of milk.


Introduction
Proteins and amino acids (AA) constitute a historical scientific issue in ruminant nutrition that counts more than 100 years [1]. The AAs have been studied in high-yielding ruminants, as feed additives capable of fulfilling their protein requirements [2,3] and improve nitrogen efficiency resulting in lowering nutritional cost and minimizing the environmental burden [4]. Further to this, AAs such as methionine are involved in several cells' metabolic pathways including that of energy [5] since are principal generators of one-carbon compounds [6]. For this reason, AAs can improve energy supply [7] and consequently affect milk performance, chemical composition, and fatty acids (FA) profile [3]. Specifically, it has been reported that methionine is indirectly involved in long-chain FA (LCFA) transportation from blood to udder through the very-low-density lipoproteins (VLDL) [8]. Additionally, methionine has a pivotal role in the de novo short-and mediumchain FA synthesis in mammary cells, since it is a methyl donor for the transmethylation reactions in the biosynthesis of lipids [3]. Moreover, methionine and lysine are widely known as the first limiting AAs for milk protein synthesis in dairy ruminants as they constitute building blocks for caseins synthesis [7,9]. In this light, an increase in milk casein content has been related to improved cheese manufacturing assets in milk, including firmer curd formation and a shorter rennet coagulation interval [10].
Moreover, AAs are involved also, in cellular oxidative balance [6] since participate in taurine and glutathione (GSH) synthesis [11]. GSH is involved in both cellular detoxification through Glutathione transferase (GST) action and hydrogen peroxide neutralization via Glutathione peroxidase (GSH-Px), while Glutathione reductase (GR) is capable of regenerating the GSH by reducing its oxidized form (GSSG). Furthermore, methionine supplementation has been proven to increase the VLDL resulting in improved vitamin E circulation [12]. Hence, the dietary supplementation with rumen-protected amino acids (RP-AA) could suppress the detrimental effect of lipid peroxidation by-products such as malondialdehyde (MDA) [12].
Up to now, the dietary inclusion of rumen-protected methionine and lysine has been well documented in dairy cows mainly focusing on milk performance [13,14], immune response [13,15], oxidative status [16], and reproduction [13,17]. Recently, limiting AAs have been tested as immune function regulators during the early lactation of dairy sheep [18,19], while extensive studies have been carried out on small ruminants' milk performance [20]. However, to the best of our knowledge, scarce information exists about the optimum inclusion level and the synergetic effect of RP-AAs on milk chemical composition, and oxidative status of both organism and milk in ewes.
Considering the aforementioned, this study aimed to investigate the effect of dietary supplementation with RP methionine or lysine, as well as with a combination of these AAs in two different proportions, on (a) milk yield, chemical composition, fatty acid profile, (b) GSH-Px, GR, GST, catalase (CAT), and superoxide dismutase (SOD), activities in blood plasma and SOD, GR and lactoperoxidase (LPO) activities in milk and c) total antioxidant capacity and oxidative stress indicators [MDA and protein carbonyls (PCs)] in both blood plasma and milk of early lactating ewes.

Diets and Experimental Design
This study continued the analytical approach initiated in previous research works [18,19]. The study was conducted with respect to the guidelines of the European Union Directive on the defense of animals used for scientific purposes (EU 63/2010; Council of the European Union 2010). Briefly, sixty, two-to three-year-old dairy ewes of pure Chios breed, at 50 days in milk (DIM) were grouped into 5 homogenous subgroups (n = 12 ewes/treatment) based on their fat corrected milk (FCM 6% ) yields (2.20 ± 0.39 kg/day), ages, and body weights (BW; 63 ± 6 kg). Each group was allocated to one of the following five groups: C: basal diet (control); M: basal diet + 6 g/ewe RP methionine (MetaSmart™); L: basal diet + 5 g/ewe RP lysine (LysiGEM™); LML: basal diet + 6 g/ewe MetaSmart™ + 5 g/ewe LysiGEM™; and HML: basal diet + 12 g/ewe MetaSmart™ + 5 g/ewe LysiGEM™. MetaSmart™ and LysiGEM™ were incorporated into the concentrates. The basal diet consisted of 1.5 kg of concentrate, 1.1 kg of alfalfa hay, and 0.6 kg of wheat straw/ewe and fulfilled the average nutritional requirements of each dietary group. The forages (alfalfa hay and wheat straw) were provided separately from the concentrates.
The diets were formulated by NDS Professional software (Ver. 3.9.7.11, Rumen Sas, Reggio Emilia, Italy) using the equations for small ruminants [21]. Animals were fed on a group basis, considering their average energy and nutritional requirements, in order the experimental design to represent the typical commercial farm feeding management and the results to have practical implications for small ruminants. Parturitions in small ruminants are in fact gathered in a very short interval which means that the animals have similar requirements and consequently feed intakes. The available feeding space was higher than the one recommended for adult housed ewes (0.33 m per animal) considering in order to favor simultaneous access and lower competitive interactions at the feeder among animals. Forage was provided with the concentrate in three equal portions after milking while no concentrate was offered in the milking parlor. Diet consumption was being recorded on daily basis ensuring no refusals occurred.
The RP-AA contained 60% isopropyl ester of hydroxyl analog of methionine, and 68% L-Lysine monohydrochloride, respectively. Therefore, the mean daily intake of ewes that were fed the commercial products of RP methionine and lysine were 3.6 and 3.4 g respectively. The diet composition is presented in Table 1, while further justifications about amino acids chosen levels and ratios are described by [18]. The whole experimental period lasted 75 days.

Sample Collection
The ewes were milked three times per day (06:00, 13:00, and 21:00) with a milking machine. Individual milk samples (n = 180; 12 ewes/group × 5 groups × 3 sampling times) were collected on the 25th, 50th, and 75th experimental day and used for milk chemical composition, fatty acid profile, and oxidative status analyses. Milk yield was recorded at the same experimental days taking into account the three milked quantities, while each of the aforementioned individual milk samples was performed by the mixture of 5% of the milk volume obtained by three milkings aiming to ensure the highest reliability.
Individual blood samples were also taken on the 25th, 50th, and 75th experimental day from the jugular vein of each ewe (n = 180) after milking prior to access on feeds. Approximately, 10 mL of whole blood were immediately transferred to heparin-containing tubes (170 units heparin; BD Vacutainer, Plymouth, UK) and stored in an icebox (Thomas Scientific, Swedesboro, NJ, USA) until their transfer to the Laboratory of Nutritional Physiology and Feeding. Then, the blood samples were centrifuged (SL16R, Thermo Fisher Scientific, Waltham, MA, USA) at 2500 rpm for 15 min at 4 • C to separate plasma from the cells.
Milk samples for chemical composition were analyzed on the collection day, while milk and blood plasma samples were stored at −80 • C, prior to fatty acids and oxidative status analyses.

Feed Analyses
Feed chemical composition is presented in Table 1. The analytical procedures have been previously described by Tsiplakou et al. [19].

Antioxidant Enzymes Activities and Oxidative Status Indicators
The assays for antioxidant enzyme activities, oxidative stress indicators, and the total antioxidant capacity were performed using a UV/Vis spectrophotometer (GENESYS 180, Thermo Fisher Scientific, Waltham, MA, USA) as previously described [27]. The GSTs activities were recorded by monitoring the conjunction of GSH to 1-chloro-2,4-dinitrobenzene (CDNT) at 340 nm. CAT activity was performed using a commercial spectrophotometric kit (Catalase Assay Kit; CAT100, Sigma-Aldrich, St. Louis, MO, USA). GSH-Px activity was assayed according to Paglia and Valentine [28]. GR activity was performed by measuring the reduction of oxidized glutathione (GSSG) to reduce glutathione in presence of nicoti-namide adenine dinucleotide phosphate (NADPH) at 340 nm. SOD activity was recorded by monitoring the inhibition of cytochrome c oxidation at 550 nm. LPO activity in milk was performed by monitoring the oxidation of 2,2 -Azino-bis (3-ethylbenzthiazoline-6-sulfonic acid) ABTS presence of hydrogen peroxide at 340 nm. MDA was measured according to Nielsen et al. [29] with some modifications described by Mavrommatis et al. [30]. The protein carbonyls (PC) were assayed according to the method of Patsoukis et al. [31]. The ABTS [32,33] and the ferric reducing ability of plasma (FRAP) [34] assays were used to assess the total antioxidant capacity.

Statistics
The dataset was evaluated in SPSS.IBM software (Version 20.0. Armonk, NY: IBM Corp.) and the results are depicted as mean ± standard error of means (SEM). The effect of dietary treatment was assessed by performing a GLM for repeated measures analysis of variance. The dietary treatments (T) (T = C, M, L, LML, and HML) were defined as the fixed factor and the sampling time (S) as the repeated measure, while their interactions (T × S) were also assessed, according to the following model: where Υijk is the dependent variable, µ the overall mean, Ti the effect of dietary treatment (i = 5; C, M, L, LML, and HML), Sj the effect of sampling time (j = 3; 25th, 45th, and 75th day), Ak the animal's random effect, (T × S)ij the interaction between dietary treatments and sampling time, and eijk the residual error. A total of 180 observations (12 ewes × 5 dietary groups × 3 sampling times) were emerged by 60 experimental units (ewes) and used for each investigated parameter (except for body weight; two sampling times = 120 observations). Posthoc analysis was applied when appropriate using Tukey's multiple range test. For all tests, the significance level was set at p = 0.05. Simplifying the visualization of the results, GraphPad Prism 6.0 (2012) (GraphPad, San Diego, CA, USA) depicted interleaved bars. Furthermore, a discriminant analysis was also performed (variables were entered independent together) on fatty acids data to establish those variables capable of distinguishing and classifying samples among the five dietary groups. Wilk's lambda(λ) criterion was used for selecting discriminant variables. Thirty-five variables (FA and group values) were entered to develop a model to discriminate the one-hundred 80 samples. Moreover, correlations between energy corrected milk yield, blood urea and b-hydroxybutyrate concentrations (data were obtained by our previous study [19]) and oxidative status of both blood plasma and milk were explored using Spearman's correlation coefficients and are presented as a 3-color heat map graph.

Animals and Milk Performance
The dietary treatments did not affect the milk yield, fat-(FCM), and energy corrected milk (ECM), while a significant reduction on these parameters throughout the experimental period was observed (p < 0.001) ( Table 2). Milk fat content was significantly increased (p < 0.001) in RP-AAs fed ewes, with the highest increase being reported in the HML group (Table 2). Milk protein content increased (p = 0.002) in the M and L compared to the C-fed ewes. The total milk casein content was significantly increased in the L compared to the C-fed ewes (p = 0.039). Total solids no-fat of milk increased (p = 0.011) in the L and HML compared to the C-fed ewes, while the total solids increased (p < 0.001) in the milk of supplemented ewes. Furthermore, significant interactions (treatment × sampling time) were observed in milk's protein and lactose content (Table 2, Figure S1).  Figure 1 depicts a discriminant plot of the five dietary treatments (C, M, L, LML, and HML) considering the data throughout the experimental period. The proportions of the samples that were correctly classified were 70.6%. Wilks' (λ) was reported at 0.113 for Function 1 (p < 0.001) and 0.307 for Function 2 (p < 0.001), while the proportions of C 4:0 , C 14:1 , C 20:0 , C 13:0 , C 6:0 , C 18:3n3 , C 15:0 and SCFA were the variables that contributed the most. The milk samples of the HML group were located in the lowest right-hand corner of the plot, far away from those of the C, M, and L groups. On the contrary, samples from the LML group were clustered between the aforementioned groups showing a moderate overlapping. Overall, the milk FAs were not considerably altered by RP-AAs supplementation. Confirming this, Table 3 presents minor changes in milk individual FAs. More specifically, the proportions of C 4:0 and C 6:0 FA were significantly decreased (p < 0.001) in the milk of the HML compared to the C-, M-, L-, and LML-fed ewes, while also significant (p < 0.001) interactions between dietary treatment and sampling time were unveiled (Table 3, Figure S2). Hence, the proportion of SCFA was also significantly decreased (p = 0.001) in the milk of HML compared to the C-fed ewes. Milk's C 14:1 proportion, showed a significant (p = 0.013) increase in M-fed ewes, while a decline (p = 0.013) was observed in those consuming the LML compared to the C diet. Interestingly, the C 14:1 proportion decreased (p < 0.001) after the 25th experimental day. These fluctuations induced significant (p < 0.001) interactions between dietary treatment and sampling time (Table 3, Figure S2). The proportions of C 18:0 and LCFA tended to increase (p = 0.066 and p = 0.099 respectively), while that of C 18:1 t11 tended to decrease (p = 0.068) in the milk of HML compared to the C-fed ewes. The cis-9, trans-11 C 18:2/ trans-11 C 18:1 ∆ −9 desaturase index, enhanced significantly in the milk of L in comparison to C-fed ewes (Table 3).

Oxidative Status
Figures 2 and 3 present the total antioxidant capacity, oxidative stress indicators, and antioxidant enzyme activities in blood plasma and milk, respectively. In blood plasma, the malondialdehyde (MDA) content was significantly reduced (p = 0.011) in methionine (M, LML, and HML) compared to the un-supplemented fed ewes. A significant increase in GST activity (p = 0.007) in the blood plasma of M and LML compared to the C and HML-fed ewes were found. These variations also induced significant (p < 0.05) interactions between dietary treatment and sampling time (Table S1, Figure S3). In blood plasma, the highest SOD and GST activities, and total antioxidant capacity, (measured by ABTS assay), were observed on the 75th experimental day (Table S1). The total antioxidant capacity of milk, measured with FRAP assay, was significantly declined in LML and HML compared to the C, M and L-fed ewes (p = 0.002), while the opposite was observed when the animals consumed the L diet. Significantly lower, ABTS values in the milk of LML compared to the C, M, and L-fed ewes were found. Lysine inclusion in ewes' diets caused a significant increase (p = 0.048) in milk MDA content. The PCs concentration was elevated (p = 0.009) in the milk of L and HML compared to the C-and M-fed ewes. The SOD activity tended (p = 0.087) to increase in the milk of L and LML compared to the C-fed ewes. The sampling time affected also, some of the oxidative parameters in ewes' milk. More specifically, the lowest SOD and LPO activities, and total antioxidant capacity, (measured by FRAP and ABTS assays), were observed at the 75th experimental day while the opposite was found for the PCs content in ewes' milk (Table S1). These alterations induced significant (p < 0.05) interactions between dietary treatment and sampling time for GR and SOD activities, MDA content, and FRAP and ABTS values (Table S1, Figures S3 and S4).
In Figure 4, a Spearman correlation depicts the regression between milk and blood antioxidant indexes, energy, and nitrogen balance metabolites (b-hydroxybutyric acid and urea respectively), and energy corrected milk yield based on the results given by Tsiplakou et al. [19]. Interestingly, ECM yield was negatively correlated with blood GST, GSH-Px, and SOD activities and ABTS values, while a significant positive correlation was observed with milk ABTS values and LPO activity, as well as GR and CAT activities in blood plasma. B-HBA concentration was negatively correlated with GSH-Px activity in blood, MDA levels in milk, and PCs concentration in both milk and blood plasma. Urea in blood was negatively correlated with GR activities in both milk and blood plasma.  In Figure 4, a Spearman correlation depicts the regression between milk and blood antioxidant indexes, energy, and nitrogen balance metabolites (b-hydroxybutyric acid and urea respectively), and energy corrected milk yield based on the results given by Tsiplakou et al. [19]. Interestingly, ECM yield was negatively correlated with blood GST,

Animals' Performance
In agreement with our findings, the dietary supplementation with RP methionine at 0.2% DM [35] or 2.75 g/animal/day [3] did not affect sheep milk yield. In dairy cows, the dietary supplementation with RP methionine (8, 16, or 24 g/day) increased milk production and its protein content [36]. Additionally, the dietary supplementation with a combination of 18.2 g methionine and 11.7 g lysine daily, increased cow's milk yield, protein, and caseins content [37], while with RP methionine only (0.17% of DM), enhanced the ECM yield [38]. Discrepancies among studies concerning the effect of dietary supplementation with methionine on milk yield have been reported in a meta-analysis study [39]. On the other hand, it is well documented the improvement of milk protein content by the dietary inclusion of methionine and/or lysine [37], since both of them, are two limiting amino acids for milk protein synthesis [40]. The enhancement of both milk protein and casein content by the RP-AA administration is very important for the dairy industry since affects milk coagulation properties [41].
In this study, the fat content increased in the milk of supplemented ewes. However, the highest fat content was reported in HML-fed ewes indicating that not only the dietary supplementation level of RP-AAs but also the ratio between them determines their effectiveness on milk performance. In agreement with our findings, an increase in milk fat content was observed in RP methionine-fed sheep [42] and goats [43]. On the contrary, some other studies reported no considerable alterations in the milk fat content of goats [44,45] and cows [46][47][48] supplemented with RP methionine. In the mammary gland, AAs can be used as an energy supply as well [49,50]. Indeed, a transcriptomic study by Bionaz et al. [51], indicated the dependence of the tricarbonic acid cycle (TCA) by the AAs as an

Animals' Performance
In agreement with our findings, the dietary supplementation with RP methionine at 0.2% DM [35] or 2.75 g/animal/day [3] did not affect sheep milk yield. In dairy cows, the dietary supplementation with RP methionine (8, 16, or 24 g/day) increased milk production and its protein content [36]. Additionally, the dietary supplementation with a combination of 18.2 g methionine and 11.7 g lysine daily, increased cow's milk yield, protein, and caseins content [37], while with RP methionine only (0.17% of DM), enhanced the ECM yield [38]. Discrepancies among studies concerning the effect of dietary supplementation with methionine on milk yield have been reported in a meta-analysis study [39]. On the other hand, it is well documented the improvement of milk protein content by the dietary inclusion of methionine and/or lysine [37], since both of them, are two limiting amino acids for milk protein synthesis [40]. The enhancement of both milk protein and casein content by the RP-AA administration is very important for the dairy industry since affects milk coagulation properties [41].
In this study, the fat content increased in the milk of supplemented ewes. However, the highest fat content was reported in HML-fed ewes indicating that not only the dietary supplementation level of RP-AAs but also the ratio between them determines their effectiveness on milk performance. In agreement with our findings, an increase in milk fat content was observed in RP methionine-fed sheep [42] and goats [43]. On the contrary, some other studies reported no considerable alterations in the milk fat content of goats [44,45] and cows [46][47][48] supplemented with RP methionine. In the mammary gland, AAs can be used as an energy supply as well [49,50]. Indeed, a transcriptomic study by Bionaz et al. [51], indicated the dependence of the tricarbonic acid cycle (TCA) by the AAs as an energy source. Additionally, it was already reported that citrate, an intermediate in the TCA cycle, has a pivotal function in fat synthesis by reducing equivalents in the form of NADPH [52,53].

Fatty Acid Profile
The dietary supplementation with rumen-protected AA had a negligible impact on ewes' milk FA profile. In accordance with our previous work, the proportion of C 14:1 increased in M-fed ewes [3]. Additionally, the dietary supplementation with methionine and lysine reduced the proportions of C 4:0 and C 6:0 FA, and consequently the short-chain fatty acids, in ewes' milk, while an opposite trend was observed for the C 18:0 and LCFA [54]. It is well-justified that low-density lipoprotein, a main source of medium-and LCFA for milk fat synthesis, is strongly dependent on methionine supply through s-adenosylmethionine (SAM) substrate [54]. Considering the tendency for lower proportions of C 18:1 t11 and higher of its biohydrogenation product (C 18:0 ) in the higher AA supplementation level (HML) it could be assumed that the increasing levels of nitrogenous compounds within the rumen could upsurge both lipolysis and biohydrogenation function. Remarkably, Giallongo et al. [55] also reported a positive linkage between rumen digestible protein levels and C18-fatty acids biohydrogenation within the rumen. However, further studies are required to validate this hypothesis.
Finally, the highest proportions (but no statistically significant in all groups), of LCFA, C 18:0, and cis-9 C 18:1 in the milk of AA s fed ewes may be related to the highest fat content which was found in these groups since these fatty acids are able to regulate the fluidity of milk fat globules and consequently enhance milk fat content [26].

Oxidative Status
A significant decline in the MDA content in the blood plasma of methionine-fed ewes (M, LML, HML) was observed accompanied by an increase in GST activity of Mand LML-fed animals. In accordance with our findings, a significant enhancement of GST activity in the intestine and hepatopancreas of methionine-fed juvenile Jian carp was found accompanied by a significant reduction in the MDA content of the aforementioned tissues [56]. Moreover, the dietary supplementation with methionine increases GST activity and glutathione concentration in ewes blood plasma [3] and cows' liver [57], respectively. A significant increase of GST activity in rats' hepatocytes was also reported when were incubated with L-methionine and L-cysteine simultaneously in vitro [58]. Methionine forms SAM, a precursor for cysteine production [59], which is the limiting AA in the de novo synthesis of GSH [60]. High GSH concentrations protect against cellular damage and tissue degeneration [61]. However, it should be mentioned that the highest inclusion level of methionine in combination with lysine (HML) did not modify the GST activity in ewes' blood plasma. Similarly, the lowest GSH content in turkeys' liver was found when their diets were supplemented with the highest methionine content (higher than that recommended by the National Research Council; NRC) [62]. Thus, special attention should be given to the supplemented quantities and in the synergistic effects (ratio) of AAs in animal diets not only to address their optimal inclusion levels but also to avoid physiological imbalance and consequently adverse effects (oxidative stress).
Milk appears to be an optimal environment for MDA production due to its high-fat content. This was reported in cows [63] and lysine-fed ewes (L, LML, HML) of our study. In vitro studies reported that although MDA can react with a broad variety of AA, the epsilon amino group of lysine is the main target [64,65]. Additionally, proteins compared to the free amino acids are more easily damaged by MDA under physiological conditions, as has been indicated in vitro [64,66]. The significant increase in milk PC content of L and HML-fed ewes supports this physiological mechanism. Besides, PCs are formed not only from oxidative deamination of AAs side chains by a direct attack of oxidants but also from products of lipids oxidation [67,68] as indicated by the positive correlation between them which was found in this study. The numerical increase in milk PCs in the M and LML groups may be attributed to methionine inclusion. Without narrowing out the involvement of other AA, methionine can also trigger to some extent the MDA concentration [64,66]. It could be assumed that not only the type of the AA but their mode of action (ratio), as already mentioned in blood plasma, seems to affect differently this metabolic pathway as well. Indeed, the dietary supplementation with these two AAs simultaneously caused the most intense decline in the total antioxidant capacity (FRAP and ABTS) of milk. This indicates that the dietary supplementation with high AA levels, which are involved in cells' detoxification mechanisms, might have a negative impact on the antioxidant status of milk as previously discussed for the HML-fed ewes' organism (blood).

Conclusions
Supplementing ewes' diet with methionine (6 g/animal/day) not only improved the organism's oxidative status, but also milk's chemical composition without affecting the oxidative stability of the latter. Lysine dietary inclusion (5 g/animal/day) induced oxidative stress in ewes' milk being more intense when combined with methionine at different ratios. More research is needed to define the optimum inclusion levels and the combination (ratio) of RP-AA in ewes' diets, in order to achieve a positive, synchronize response of the antioxidant defense system of both organism and milk.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/antiox10050654/s1, Figure S1. The mean individual protein and lactose content of ewes' milk were fed with the five diets (C, M, L, LML, and HML) illustrated as connected superimposed symbols (interaction between the means of five dietary treatments and time response). Figure S2. The mean individual fatty acids (FA) (% of total FA) of ewes' milk were fed with the five diets (C, M, L, LML, and HML) illustrated as connected superimposed symbols (interaction between the means of five dietary treatments and time response). Figure S3. The mean individual antioxidant enzyme activities and oxidative indices of ewes' blood plasma were fed with the five diets (C, M, L, LML, and HML) illustrated as connected superimposed symbols (interaction between the means of five dietary treatments and time response). Figure S4. The mean individual antioxidant enzyme activities and oxidative indices of ewes' milk were fed with the five diets (C, M, L, LML, and HML) illustrated as connected superimposed symbols (interaction between the means of five dietary treatments and time response). Table S1. Means of total antioxidant capacity, oxidative status biomarkers, and enzyme activities (Units/mL), in the blood plasma and milk of ewes fed the five diets (C, M, L, LML, and HML) within the three sampling times.