Glycerol-based formulations increases ovulation rate in sheep similarly to eCG but red blood cell indices may be affected

Background: The objective of this study was to investigate the metabolic and osmotic effects of different doses of glycerol or a glycerol – propylene mixture in Sarda sheep with the aim to identify those able to beneficially modify ewe’s metabolic status without harmful changes in red blood cell (RBC) indices. Thereafter, the selected doses were tested for their effects on ewe’s ovarian activity during an induced follicular phase and compared to the effects of a hormonal treatment with equine chorionic gonadotrophin (eCG). Results: Glycerol was administered alone (G groups: 90% glycerol and 10% water; % v/v) or in combination with propylene glycol (M groups: 70% glycerol, 20% propylene glycol, 10% water; % v/v). Treatments were formulated to provide 100%, 75%, 50% and 25% of the amount of energy supplied in previous experiments. Obtained results showed that the formulations G75 and M75 (22.5% and 18.2% on DM basis, respectively) induce metabolic changes comparable to those induced by M100. The latter dose has been already evaluated for its effects on sheep metabolism and reproductive performance. However, with these high doses, plasma osmolality increased significantly, and RBC indices showed significant alterations. The low dose groups (G25 and M25, 8.6% and 6.9% on DM basis, respectively) did not show any alterations in plasma osmolality and RBC indices, but the metabolic milieu differed markedly from that of M100. Between the medium dose groups, M50 (12.9% on DM basis) showed a more comparable milieu to M100 than G50 (15.9% on DM basis) and no RBC alterations. Therefore, M75, G75 and M50 doses were tested for their effect on ovarian functions and proved to be equally effective as eCG. Conclusion: The results of the present study evidenced an alteration of RBC indices, and possibly of their functions, as a side effect of glycerol administration at high doses in the diet of ewes. Therefore, protocols foreseeing the administration and RBC indices among the different groups before and 2 hrs after glucogenic formulations administration. Pearson’s correlation analyses showed that the moles of glycerol and propylene glycol permeability ions osmotic fragility These results are confirmed by another study in which the glycerol- dependent shape alterations (increased MCV and decreased MCHC) were also associated with an electrolyte shift (34). Glycerol is indeed a polar molecule which is known to alter the ionic strength and the dielectric constant of aqueous solutions (46). In sheep, incubation of RBCs in glycerol media causes haemolysis due to a gradual increase of the intracellular glycerol concentration (47). We can thus speculate that the RBC shape alterations observed in the high dose groups of this investigation may be caused by the glycerol permeating cell membrane following a concentration gradient (32) generated by the significant rise in glycerol circulating concentrations consequent to the dietary administration. Shape alterations of sheep RBCs reported in this trial are similar to those described in human glycerolized RBCs (34,45), and might thus be associated with erythrocyte impaired functionality. In addition, when the average size of RBC’s increases, they can fail to pass the capillaries in the microcirculation and be removed from the circulation (48). All these changes might thus contribute to enhanced eryptosis. Further studies are needed to quantify the permeation of glycerol in the RBCs in vivo and to extend our knowledge on its osmotic effects. In the present study, in all groups but M100 the RBC shape alterations were not observed in blood samples collected 12 hrs after the glucogenic mixture administration. Nevertheless, reported results suggest the need to lowering the dose of the glycerol administered both as replacer of corn grain and as feed supplement, especially in long-lasting protocols. sampling hour and their first-order interactions as fixed effects. As post-hoc test, for plasma glycerol, glucose, insulin, NEFA and urea, a one-tailed Dunnett’s test was used to highlight negative or positive differences, in respect of M100 group, as appropriate.

livestock production (1). These waste streams can source from different agro-industrial processes, and their use as feed ingredient or feed supplement can contribute to the environmental sustainability of livestock sectors. Recently, there has been an increasing global demand for biofuels, resulting in an increased demand for the feedstuffs (corn, wheat, and oilseeds) used for fuel extraction. The main by-product during the production of biodiesel is glycerol, which can be used as livestock feed ingredient or feed supplement (2). The EU legislation approved glycerol and another biodiesel co-product, namely propylene glycol, as animal feed additives with no restrictions on animal species or quantity that may be fed, while recommending from the European Food Safety Authority the collection of data on the presence of impurities and contaminants in crude glycerine from biodiesel production (3,4).
In small ruminants, glycerol, alone or in association with propylene glycol, has been used both as replacer of corn grain in the diet (5)(6)(7)(8)(9) and as a feed supplement to increase productive performance (10)(11)(12)(13)(14)(15)(16)(17). In these studies, glycerol was supplemented at a percentage of inclusion on dry matter (DM) basis ranging from 3% up to 45%. Overall the metabolic effects described include an increase in glucose and a decrease in NEFA circulating concentrations (6,11,18). Other studies in ewes also reported an increase in insulin (19)(20)(21), insulin-like growth factor-1 (20), and a decrease in urea circulating concentrations (20,21) after oral administration of 280 mL glycerol and 80 mL 1,2propylene glycol daily for 4 days. These changes have been associated with changes in follicular fluid composition during treatment period (20), and may contribute in creating a suitable systemic and local metabolic milieu for the promotion of ovarian function. Glucogenic mixtures based on glycerol, administered as short-term flushing can indeed increase ovulation rate (19) and oocyte quality (10) in sheep. The increase in ovulation rate following glycerol administration has been described also in other studies in which glycerol was administered in a single dose (14,15). Therefore, a short or very short (daily) flushing based on glycerol can represent a valid alternative to other short flushing techniques based on conventional feedstuff such lupin seeds (22) or soybean meal (23,24). Moreover, thanks to its boosting activity on ovarian functions, glycerol could possibly substitute for synthetic gonadotropins dosed to increase lambing and twinning rates in sheep. Among them, it is noteworthy the equine chorionic gonadotropin (eCG), a hormone commonly used in concert with progestogen to induce ovulation prior to natural mating or artificial insemination. The future use and availability of eCG appears to be strongly challenged by highly active animal-rights movement because the hormone is obtained from pregnant mares. Hence, there is a need for alternative synchronization protocols without eCG (25).
Although there is evidence that glycerol and glucogenic mixtures can improve production and reproductive performance in sheep, these additives may sometimes negatively impact on animal welfare, at least under some circumstances. In fact, after oral administration, a significant proportion of glycerol (44% according to (26) is rapidly absorbed through rumen wall by passive diffusion and transported via the bloodstream to serve directly as a substrate for glucose synthesis in the liver (27,28). In a previous study (21), we reported that the administration of a glycerol-based formulation (280 mL glycerol and 80 mL propylene glycol for 4 days; 23% DM) caused a 400-folds increase in circulating glycerol concentration in dairy ewes. Unfortunately, glycerol is a hyperosmotic agent and, consequently, plasma osmolality was greatly increased (21). This effect had been already observed in humans (29)(30)(31). Glycerol can readily permeate the red blood cells (RBCs) membrane at body temperature following a concentration gradient (32). Both the increase in plasma osmolality (33) and glycerol diffusion in the RBCs (34) can impact on cellular volume homeostasis which is critical for erythrocyte survival (35). Changes in RBCs volume can affect their membrane integrity, leading to impaired functionality and increased eryptosis (36,37). The critical importance of RBC volume regulation is demonstrated by several pathologies resulting from either overhydration or dehydration of RBCs (33).
Thus, the objective of this study was to investigate the metabolic and osmotic effects of different doses of glycerol and glycerol-propylene glycol mixtures in Sarda sheep in order to identify those able to positively modify ewe's metabolic status without significantly affecting RBC indices. The highest dose of the glucogenic mixture was taken from previous studies, being already evaluated for its effects on sheep metabolism and reproductive performance (8,16,18). Thereafter, to determine the ability of the selected glucogenic doses to promote ovarian function, we evaluated their effects on ewe's ovarian activity during an induced follicular phase, comparing them to the effects of the administration of a synthetic gonadotropin (eCG). Ovarian function was assessed by threedimensional ovarian ultrasound scanning and oestradiol production.
Before starting the titration feeding treatments, plasma osmolality, RBC indices, circulating concentrations of analysed metabolites and hormones were within the physiological ranges for the species (38) in all the ewes showing no differences between groups (D 0; Table 1).

Corresponding P-values are shown within parenthesis
The administration of the glucogenic formulations determined indeed a significant increase in plasma osmolality in the high and moderate dose groups (M100, G100, M75, G75; p<0.001), in which peak values were significantly higher than those recorded in the medium (M50 and G50) and low (M25 and G25;) dose groups (p<0.05, Table 3). RBC indices and their variations after treatment are shown in Table 4. The absence of microscopic platelet clumping or RBC agglutination or roulaux formation in whole blood smears allowed to consider blood samples processed at automatic cell counter as diagnostic for the purposes of this investigation. It is to point out that the absence of roulaux bodies may also have a diagnostic indication for particular conditions, including erythrocyte sedimentation rate (ESR) which may occur in the presence of an inflammatory condition, for instance, despite it is consolidated that the two parameters may not be strictly proportional (39). In view of the diagnostic aspects, all animals enrolled appeared clinically healthy throughout the trial.
Compared to pre-administration values, no changes in RBC indices were found in M25, G25 and M50 groups, while in G50 a significant decrease in mean cell haemoglobin concentration (MCHC) was found. In the high and moderate dose groups (M100, G100, M75, G75), the significant decrease in MCHC was accompanied by a significant increase in mean cell volume (MCV) and RBC distribution width (RDW-SD; Table 4). Representative pictures of the RBC shape alterations found in the high dose group 2 hrs after treatment administrations are shown in Figure 1. 00:00 h: immediately before the morning administration of glucogenic mixture. MCV: defines the size of the red blood cells and is expressed as femtoliters (10 −15 ; fl); MCH: quantifies the amount of haemoglobin per red blood cell and is expressed in picograms (pg); MCHC: indicates the amount of haemoglobin per unit volume and it is expressed as g/L of RBCs; RDW-SD: represents the coefficient of variation of the RBC volume distribution (size) and is expressed as femtoliters.
* Asterisks indicate significant differences within the same group between time 00.00 and time 02.00.
In line with these results, Pearson's correlation analyses showed that the moles of glycerol and propylene glycol administered daily were directly correlated with RDW-SD (r=0.539; p<0.001) and MCV (r=0.495; p<0.001) and indirectly correlated with MCHC (r=0.582; p<0.001; Table 5). It should be pointed out that in all the groups but M100 these alterations were not observed in blood samples collected on the D 3 before glucogenic mixture administration, i.e., after 12 hrs from the last administration (Table 2). Table 5. Pearson correlation coefficients between the moles of glucogenic compounds administered and RBC indices.

Corresponding P-values are shown within parenthesis
It should be pointed out that in all the groups but M100 these alterations were not observed in blood samples collected on the D 3 before glucogenic mixture administration, i.e., after 12 hrs from the last administration (Table 1). To identify the formulations and doses capable to induce a positive change in ewe's metabolic status we assessed their effects against those induced by M100. The latter dose has been indeed already evaluated for its effects on sheep metabolism and reproductive performance (8,16,18). Obtained results shows that the high and moderate dose groups (G100, M75 and G75) showed no differences in mean glycerol plasma concentrations compared to the M100 one ( Figure 2, panel a). Consequently, no differences were found between those groups and M100 in the circulating concentrations of glucose, insulin, urea and NEFA ( Figure 2, panels b,c,d,e). Pearson's correlation analyses showed indeed that the moles of glycerol and propylene glycol administered daily showed a fair direct correlation with glycerol, glucose, insulin and NEFA circulating concentration (r=0.592; p<0.001) and a fair indirect correlation with urea circulating concentration (r=0.574; p<0.001; Table   2).
In line with this results, the G50 and low dose groups (M25, G25) showed lower mean glycerol concentrations compared to M100, and consequently differed from M100 in the mean concentrations of glucose, insulin, and urea ( Figure 2, panels b,d). It is noteworthy that M50, while differing from M100 in glycerol, glucose and urea mean concentrations, showed no difference in mean insulin concentrations ( Figure 2, panel c). NEFA mean concentrations were higher than those found in M100 only in the low dose groups. Supplementary Figure 1 shows the concentrations-time data of the analysed metabolites and hormones.
Taken together, results show that the formulations M75 and G75 are able to induce metabolic changes comparable to those induced by M100. Hence, the dose of the glucogenic formulation can be lowered from 23% to 18% of diet on DM basis, without apparently inducing any significant change in ewe's metabolic response. However, with these moderate doses, plasma osmolality increased significantly, and RBC indices showed significant variations. Between the medium dose group, M50 showed lower glucose and higher urea mean concentrations than M100. However, insulin and NEFA mean concentrations did not differ from those found in the M100 group. Moreover, the M50 group did not show any significant variations in plasma osmolality and RBC indices thus suggesting that glycerol osmotic effects were less marked than in the high and moderate dose groups. Therefore, the glucogenic doses selected to determine the promotion of ovarian function (phase 2 of the present study) turned out to be M75, G75 and M50.

Phase 2
As in the previous phase, body weight did not change after the administration of the feeding treatments and no differences were observed between groups (M75 42.3 ± 1.7, G75 42.4 ± 1.8, M50 44.4 ± 1.8, GON 42.1 ± 3.5).
The number of follicles ranged between 1 to 3 mm in diameter, ovarian volume (mL), vascularity index, flow index and their relationship ( Table 6) did not show any differences between groups. The number of large follicles (4 to 7 mm in diameter), however, was higher in the M75 group compared to the GON group (P=<0.05; Table 6), while no differences were found between the other groups. were not statistically significant.

Discussion
The results of the present study showed, for the first time, that the administration of glycerol-based glucogenic formulations causes significant changes in RBC indices of ewes with a dose-related effect.
When the glucogenic formulations were administered at a percentage of inclusion on DM basis ranging from 18.2% to 27.4%, these changes indicated an increase in RBC volume (MCV), which was associated with a variation in the size of the erythrocytes (increased RDW-SD) and with a dilution in the amount of haemoglobin per volume unit (decreased MCHC). In the high and moderate dose groups the RBC shape alterations were found 2 hrs after the administration of glucogenic formulation administration when glycerol concentrations in the bloodstream were higher than ≈ 200 mg/dL. This rise was also associated with a significant increase in plasma osmolality, as already reported in lactating ewes when the glycogenic formulation was given at 23% of inclusion on DM basis (21). When the glycerol formulations were administered at lower doses (G50, 15.9% DM) a significant decrease in MCHC was the sole alteration observed. Below such value (≤ 12.9% DM; M50, M25 and G25), no significant alterations in RBC indices were found.
Glycerol is lipid soluble so it can diffuse by simple diffusion directly through the RBC membrane following a concentration gradient (32). It is known that glycerol enters human RBCs by two pathways: facilitated diffusion through aquaporins, and simple diffusion (40,41). By contrast, glycerol is transported in sheep RBC via simple diffusion only (41,42), and hence it has a lower permeability coefficient compared to humans (43). In humans, the effects of glycerol permeation on RBC functionality have been described in in vitro studies focusing on RBC cryopreservation, where glycerol is used as penetrating cryoprotectant (44). It has been reported that glycerol permeation in the RBC membrane resulted in shape alterations (increased MCV and RDW-SD) which appeared to be closely related to membrane alterations, increased permeability to ions and increased osmotic fragility (45).
These results are confirmed by another study in which the glycerol-dependent shape alterations (increased MCV and decreased MCHC) were also associated with an electrolyte shift (34). Glycerol is indeed a polar molecule which is known to alter the ionic strength and the dielectric constant of aqueous solutions (46). In sheep, incubation of RBCs in glycerol media causes haemolysis due to a gradual increase of the intracellular glycerol concentration (47). We can thus speculate that the RBC shape alterations observed in the high dose groups of this investigation may be caused by the glycerol permeating cell membrane following a concentration gradient (32)  RBCs (34,45), and might thus be associated with erythrocyte impaired functionality. In addition, when the average size of RBC's increases, they can fail to pass the capillaries in the microcirculation and be removed from the circulation (48). All these changes might thus contribute to enhanced eryptosis.
Further studies are needed to quantify the permeation of glycerol in the RBCs in vivo and to extend our knowledge on its osmotic effects. In the present study, in all groups but M100 the RBC shape alterations were not observed in blood samples collected 12 hrs after the glucogenic mixture administration. Nevertheless, reported results suggest the need to lowering the dose of the glycerol administered both as replacer of corn grain and as feed supplement, especially in long-lasting protocols.
By lowering the dose, however, the efficacy should be guaranteed. To verify whether the tested formulations were able to create a metabolic milieu suitable for the final follicular growth and thus for the conception period in the ewe, as described in previous studies (20), we compared the circulating concentrations of glycerol, glucose, insulin, urea and NEFA against those found in the M100 dose. This dose was taken from previous studies, being already evaluated for its effects on sheep metabolism and reproductive performance (8,16,18). In the low dose groups, showing no alterations in plasma osmolality and RBC indices, the metabolic milieu was markedly different compared to that found in M100 (lower glycerol, glucose and insulin circulating concentrations, higher NEFA and urea circulating concentrations). In the high and moderate dose groups (G100, G75 and M75) the metabolic milieu elicited did not differ significantly from M100. However, plasma osmolality and RBC indices showed significant alterations compared to pre-administration values.
The medium dose groups, despite providing the same amount of energy (0.6 Mcal/d of NE L ), differed both for the metabolic response and for the RBC indices alterations observed. When this amount of energy was supplied by an association of glycerol and propylene glycol (M50), the metabolic milieu elicited was more similar to that found in M100 and no alterations in plasma osmolarity and RBC indices were observed. If the same amount of energy was supplied by glycerol alone (G50), the total dose of glycerol administered daily increased from 70 to 110 mL. This higher dose was associated with RBC indices alterations (increased MCHC). In addition, insulin concentrations were lower than those found in M100. The association of glycerol with propylene glycol seems thus a better option than glycerol alone. The hyperinsulinemic effect of propylene glycol is related to its glucogenic effect consequent to the increased absorption of propanol, propionate, and propanal originating from ruminal metabolism of propylene glycol (49). This effect could be related to insulin resistance consequent to increased circulating levels of propylene glycol and propanol or to a decrease in the ratio of ketogenic to glucogenic metabolites in plasma (49).
However, if administered at high doses propylene glycol fermentation may produce sulphurcontaining gases (50) and provoke negative effects, including ataxia, salivation, hyperventilation and depression (51). At low doses, however, it allows reducing the amount of glycerol administered, and hence its osmotic effects, while acting synergistically in rising insulin concentrations.
The rise in glucose and insulin directly promotes follicular growth and development, mainly by altering FSH-induced effects on the synthesis of oestradiol by the granulosa cells, in accordance with the availability of glucose (52). The follicle has indeed a functional insulin-glucose-IGF-1 system which is affected by short-term nutritional treatments, and components of this metabolic system are nutritionally regulated in the follicle (53). Moreover, elevated circulating concentrations of NEFA, which are indicative of negative energy balance, affect follicular growth and fertility by acting directly on the different ovarian cells, including the oocyte within the growing follicle, and may also affect the development of the fertilized oocyte, morula and early blastocyst (52). Therefore, diets or metabolic states that favour high concentrations of NEFAs should be avoided during the cycle of conception and the early post-conception period (52). In the same way, high levels of urea in blood have been associated with lower fertility due to a changed uterine environment and poor embryo viability (54)(55)(56).
The administration of M50 was associated with an increase in the mean circulating concentrations of glucose, although less marked than that found in M100, and insulin and with a decrease in NEFA circulating concentrations. Taken together, these results suggest that the M50 may represent a safe glucogenic formulation in ewes, being able to elicit comparable metabolic conditions as found at higher doses, but without altering RBC indices and plasma osmolality. The effects at ovary levels were thus compared with those of moderate doses of glucogenic formulation (M75 and G75) and with those provoked by eCG, having an FSH/LH-like effect. The use of gonadotropins is routinely incorporated into synchronization systems used in ewes to induce ovulation. The most commonly used product is eCG because, besides inducing ovulation, such hormone reduces variability in the onset of oestrus among treated animals and increases ovulation rate as a consequence of a greater follicular development (57). In the present study, all glucogenic formulations were able to elicit an ovulatory response which did not differ significantly from what found in eCG treated ewes. It is remarkable that the high doses (G75 and M75) did not overcome the medium one (M50), which in contrast showed a numerically higher ovulation rate than the other treatments (n.s).
Glycerol is widely used as a feed supplement in ruminant nutrition. However, the results of the present study showed an alteration of RBC indices, and possibly of their functions, as a side effects of glycerol administration at moderate and high doses. Therefore, protocols foreseeing the administration of glycerol should be tested for their effects on RBC indices and functions. In general terms, medium doses of glucogenic mixtures should be preferred for flushing dairy ewes, as they proved to be effective at metabolic level without causing alterations in RBC indexes and possibly of their functionality.

Materials And Methods
The experiments were carried out at Bonassai research station of Agris Sardegna, located in northwestern Sardinia, Italy (40 °N, 8 °E, 32 m a.s.l.). This study was divided into two experimental phases.
Phase 1 aimed at determining the lowest doses of glycerol as a titration feeding trial able to modify ewe's metabolic status, as previously described using a high dose (20), without impacting on RBC indices. Glycerol was administered both alone and in combination with propylene glycol, as previously described (10,19,20). there is an 81% chance of correctly rejecting the null hypothesis of absence of differences between groups.

Animals and treatments
The experiment was run during March-April 2019, within the anoestrus season described for this breed at this latitude. Twenty adult (3-7 years) and multiparous non-lactating Sarda dairy ewes were used. Glycerol and propylene glycol had a purity grade of 99.5-100% and complied with EU Reg. 231/2012 for food additives (E422 and E1520 for glycerol and propylene glycol, respectively; Farmalabor srl, Assago, Milano, Italy). Glycerol was administered alone (G groups: 90% glycerol and 10% water; % v/v) or in combination with propylene glycol (M groups: 70% glycerol, 20% propylene glycol, 10% water; % v/v). Treatments were formulated in order to give 100%, 75%, 50% and 25% of the amount of energy supplied in previous experiments (10,20,21).   Figure 3 shows the exact timing at which treatments, samplings and measurements were performed.

Feedstuff composition
The compound feed was a commercial pelleted formulation based on cereals and soybean meal, with a net energy content (NE L ) of 1.43 Mcal/kg, 17% Crude protein (CP) and 32% starch, on DM basis. The home-grown hay was a late-cut Italian ryegrass with 6.5% CP and 75 % NDF and a NE L content of 0.86 Mcal/kg on DM basis. The hay was chopped before feeding.

Blood samplings and plasma osmolality determination
On day 0, in both periods, plasma concentrations of glycerol, glucose, insulin, NEFA, urea, plasma osmolality and, in whole blood, RBC indices were determined from samples drawn from jugular vein before morning feeding at 0800. On day 3, plasma concentrations of glycerol, glucose, insulin, NEFA, urea, and plasma osmolality were also determined from eight consecutive samples, while RBC indices from two samples (at 0800 and at 1000). Repeated sampling on day 3 were collected after jugular cannulation. At each sampling, from each ewe, two blood samples were collected: one using 2-mL vacuum collection tubes with glycolytic inhibitor (5.0 mg sodium fluoride, 4.0 mg pot. ox. -Vacutainer Systems Europe; Becton Dickinson, MeylanCedex, France) for glucose assay; one using 2.0 mL vacuum collection whole blood tube with spray-coated K 2 EDTA (Vacutainer Systems Europe; Becton Dickinson, MeylanCedex, France) for other metabolites, insulin and RBC indices determination.
Immediately after recovery, samples were cooled to 4°C. RBC indices were determined within 2 hrs of blood collection. The other blood samples were centrifuged at 1,500 g for 15 min at 4 °C degrees.
Individual plasma was removed and stored in vial at -20°C until assayed.

Complete blood count (CBC)
Haematological testing was integral to the interpretation of results to establish the safe amount of glycerol following the titration trial. Particular care was taken to prevent analytical biases responsible to affect blood sample diagnostic potentials. The control of non-experimental factors with osmotic effect started from the anticoagulant choice in our methodological approach. In light of the recommendations of the International Council for Standardization in Hematology (ICSH) the use of K 2 EDTA (ethylene diamine tetracetic acid) as anticoagulant was preferred, due to the less pronounced osmotic effect on blood cells than that exerted by K 2 EDTA (58,59). However, the K 2 EDTA requires the optimal blood to additive proportion and mixing to prevent clotting (58). That way, on sampling a same skilled veterinary surgeon collected blood samples in vacuum tubes until filled to the correct volume and gently mixed 5 to 10 times to allow complete interaction with the anticoagulants, before storage in the upright position and cooling to 4 °C. Time and temperature from blood sampling to haematological analyses were therefore crucial. All procedures adopted in this trial were carried out in respect of the recommendations to preserve diagnostic samples in the preanalytical phase (ICSH, 1993). All the samples were transported refrigerated to the laboratory where individual haematological profiles were assessed within 2 hours of blood sampling.
In brief, 15 μL of each whole blood sample were needed for the determination of haematological parameters through capillary analysis based on tri-angle scattering and chemical dying read through flow cytometry technology. The following parameters were determined using an automatic cell

Insulin
ELISA assays were performed using the Personal Lab Adaltis (Adaltissrl, Rome, Italy), which is a tool that performs automated ELISA protocols. Insulin concentration was measured in duplicate using a commercial Ovine Insulin ELISA Kit (Mercodia developing diagnostics, Germany) which is a solid-phase ELISA based on the direct sandwich technique. The kit is calibrated against an in-house reference preparation of ovine insulin, and it has been previously used for insulin determination in ovine plasma (61,62). The mean ovine insulin concentrations of the six reference solutions were 0, 0.05, 0.15, 0.5, 1.5, and 3 mg/L. The recovery on addition was 94%-114% (mean 103%). The analytical sensitivity was 0.025 mg/L and the intra-assay and interassay CV values were <7%.

Animals and treatments
The experiment was run during August-September, within the natural breeding season described for this breed at this latitude. Glucogenic treatments were administered following an oestrus synchronization protocol. In brief From D 0 to D 4, i.e. throughout the treatment periods, the groups were kept indoors in separate pens. Indoor daily feeding consisted of 200 g/head of a commercial pelleted feed and 150 g/head of maize grain divided in two equal meals (0800 and 1500) and individually fed in the milking parlour and 1500 g/head of hay fed in the afternoon. On the blood sampling day, the morning meal based on compound feeds was administered immediately after the first blood sampling. Water and mineral blocks were available ad libitum. Concentrates were completely consumed by the animals.

Feedstuff composition
Compound feeds and hay administered in the phase 2 were the same as those of phase 1. The whole maize grain had CP content of 7.80 %, starch level of 69.9 % and a NE L of 2.00 Mcal/kg on DM basis.

Determination of oestradiol plasma concentration
On day 5 of the experimental period (48 hrs after the second PGF2α administration), blood samples were drawn from jugular vein at fasting at 0800. Blood samples were collected using 10-mL vacuum collection tubes (BD Vacutainer® Rapid Serum Tube -Vacutainer Systems Europe; Becton Dickinson, MeylanCedex, France) for oestradiol assay. Immediately after recovery, blood samples were cooled at 4°C and centrifuged at 1,500 g for 15 minutes. Serum was removed and stored at -20°C until assayed.
Samples were measured in duplicate. Concentrations of oestradiol were measured after sample extraction by using a highly sensitive commercial enzimoimmunoassay kit for quantitative determination of estradiol-17β (Demeditec Diagnostics GmbH, Kiel-Wellsee, Germany). Sensitivity of the assay was 1.4 pg/mL and intra-assay variation coefficient was 5.7%.

Ovarian 3D ultrasound scanning
From day 1 to day 4, ovarian follicular population was studied daily using trans-vaginal 3D ultrasonography.
The images were recorded as digital files for later examination and the ovarian follicular population vascular index (VI), flow index (FI), blood vessels and blood flow index (VFI) by threedimensional Doppler histogram were calculated.
Thereafter, on day 14 ovulation rates were determined by counting the corpora lutea present in each ovary by transrectal ultrasonography with a real-time B-mode scanner (Aloka102 SSD 500; Aloka Co., Tokyo, Japan) fitted with a 7.5MHz linear-array probe.

Statistical analyses
Results are expressed as mean values (mean ± SE) or median values (median and range) and the differences were considered to be statistically significant at P < 0.05.
Differences in ewe's body weight at the beginning and at the end of the glucogenic-treatment period in phase 1 and 2 were analysed by a mono-factorial GLM. A SPSS GLM UNIANOVA model was also used to highlight changes among groups on plasma osmolality, RBC indices and circulating concentrations of analysed metabolites and hormones between day 0 and day 3.
Longitudinal data of plasma osmolality in the consecutive samples collected on day 3 (during treatment period) were analysed by a GLM UNIANOVA model in SPSS with treatment, sampling hour and their first-order interactions as fixed effects. A post-hoc test Tukey's test was used to highlight differences between and among groups.
A SPSS GLM UNIANOVA model was also used to evaluate within groups the effect of nutritional treatments on RBC indices before and 2 hrs after treatment administration on D 3.
Finally, the relationship between treatments administrated, concentration of metabolites and hormones and the relationship between treatments and RBC indices at two hours post treatment, were evaluated by Pearson's correlation analysis. The strength of the correlation was considered poor for r values ranging from ± 0.1 to ± 0.3, fair for r values ranging from ± 0.3 to ± 0.5, moderately strong for r values ranging from ± 0.6 to ± 0.8 and very strong for r values ranging from ± 0.8 to ± 1 (63).
In phase 2, longitudinal data on follicular population, ovarian vascularity and flow index and their relationship, as evaluated by ultrasound scanning, were analysed by a GLM UNIANOVA model in SPSS with treatment, day and their first-order interactions as fixed effects. As post-hoc test, Fisher LSD test was used to highlight differences between groups. A One-way ANOVA test was used to compare ovulation rate and oestradiol blood levels between the different experimental groups.

Consent for publication
Not applicable

Availability of data and materials
The datasets used and analyzed in the current study are available from the corresponding author on reasonable request.

Competing interests
The authors declare that they have no competing interests.

Supplementary Files
This is a list of supplementary files associated with this preprint. Click to download. Additionalfigure1.pdf NC3RsARRIVEGuidelinesChecklist.pdf