Comparison of models for describing the lactation curve of latxa sheep and an analysis of factors affecting milk yield

Feed efficiency and resilience in dairy ewes subjected to a nutritional challenge

It is currently uncertain how selection of more efficient animals might affect other traits such as resilience (which, in this context, is defined as the ability of an animal to sustain or revert quickly to its previous production level and health status after a disturbance), especially in small ruminants. However, improving, or at least maintaining, resilience is of utmost importance to ensure livestock production in the face of external perturbances, which are expected to become more prevalent in the near future due to climate change and global instability. This study was conducted to investigate whether a nutritional challenge consisting of animals receiving only 70% of their voluntary feed intake (DMI) for 26 d, might differentially affect the response of high and low feed efficiency (FE) in sheep. To meet this aim, residual feed intake (RFI), as a proxy of FE, was calculated in 40 lactating ewes fed a total mixed ration ad libitum. Using the RFI values, the ewes with the highest (H-FE; n = 8) and the lowest (L-FE; n = 8) FE were selected and subjected to the nutritional challenge. After the feed restriction, all sheep were again fed ad libitum for a further 25 d. Temporal patterns of variation in animal performance, ruminal fermentation, and blood and milk metabolites in response to, and recovery from, the underfeeding were examined in H-FE versus L-FE. More efficient ewes (H-FE) tended to eat less than less efficient ewes (L-FE) for the same level of production. Linear regressions for DMI and milk yield data showed that time-series response to the challenge was comparable in L-FE and H-FE. Despite temporal changes due to underfeeding, both FE divergent groups displayed a similar response to, and recovery from, the challenge with no significant differences in the temporal patterns of variation of the parameters analyzed, except for a tendency to higher plasma glucose in the L-FE ewes. These results challenge the general hypothesis that more efficient animals would be less resilient, but further research is still needed on the relationship and possible trade-offs between high efficiency and resilience.

Genetic and environmental determinants of immunoglobulin G in kid serum and adult colostrum of dairy goats

Immune competence is a key factor in the health and resilience of dairy ruminants. The aim of this study was to study IgG concentration in serum shortly after birth and in colostrum of lactating animals as proxy of passive immune transfer and colostrum quality, respectively, and thus provide new health and resilience indicators for selection. Data from 424 female kids and 203 primiparous goats, reared on an INRAE experimental farm, were analyzed. The IgG concentration was measured in kid serum and goat colostrum, respectively, using the radial immunodiffusion technique. Linear models were used to identify factors of variation in serum and colostrum IgG concentration, and health events were monitored for 5 wk after birth to measure the effect of passive immune transfer on kids' health. We found a significant effect of serum IgG level on hindquarters cleanliness score (as a proxy of digestive troubles) and nasal and ocular discharge. As expected, the absorption of IgG in the days following birth is essential to ensure the passive transfer of immunity from the mother to the youngster, and provide immune protection. Genetic parameters were estimated for serum IgG concentration in kids and for colostrum characteristics in primiparous goats. Heritability of serum IgG in kids was not significantly different from zero. This low value, with large SE, may be due to the limited sample size and large environmental variability likely not accounted for. Especially, we could not adjust for the amount of IgG that was ingested by the kid, which can have a significant effect on serum IgG and possibly on heritability estimates. Heritability of colostrum traits were moderate, with estimates of 0.25 (±0.16), 0.35 (±0.16), and 0.25 (±0.16) for IgG concentration, IgG mass, and colostrum volume, respectively. Genetic correlations with milk production showed high standard errors due to small numbers. However, estimates suggest a negative genetic association between production and IgG concentration in colostrum (genetic correlation = -0.44 ± 0.33), probably linked to a positive association with total colostrum volume. The significant heritabilities that were estimated for the first time demonstrate the existence of genetic control of colostrum quality in goats, and open up opportunities for selection in this species.

Evaluation of performance and estimation of genetic parameters for milk yield and reproductive traits in Awassi and Asaf sheep breeds and crosses in the Palestinian territory

This study aimed to estimate genetic parameters, including heritability and repeatability, for milk yield and prolificacy in sheep. It included 3682 records from 1837 ewes across various breeds: Awassi, Assaf, and Awassi x Assaf crosses, two Awassi lines: Improved Awassi, AFEC Awassi. The study measured total milk yield (TMY), yield up to 120 days (TMY120), and yield up to 150 days (TMY150), alongside reproductive traits: litter size (LS), number of lambs born alive (NLBA), and lambing interval (LI). The analysis utilized a mixed model and the REML procedure to estimate genetic parameters. Results indicated that litter size (LS) had no significant impact on milk traits, whereas breed, location, ewe parity, and lambing season showed significant effects. Lactation length also significantly influenced TMY. For reproductive traits, treatment was significant for NLBA, with location associated with breed, parity, and season affecting all traits. Heritability estimates for TMY ranged from 0.00 in Awassi x Assaf crosses to 0.11 in Awassi, and for TMY120 and TMY150, from 0.00 in Awassi x Assaf crosses to 0.16 in Awassi. Estimates for LS and NLBA varied similarly. The heritability for LI was 0.03 in Awassi and zero in other breeds and crosses. The findings suggest that in the northern West Bank regions of Nablus and Jenin, Assaf or Awassi x Assaf breeds are preferable. The study underscores the importance of comprehensive performance and pedigree recording for effective sheep farm management.

The Impact of Age at First Lambing on Milk Yield and Lactation Length in a Population of Istrian Sheep under Semi-Intensive Management

This study aimed to examine the impact of ewe's age at first lambing (AFL) on days in milk (DIM), average daily milk yield (DMY), and total milk yield (TMY). Symmetrical bimodal distribution of AFL enabled classification of maidens in those mated in the first (47%) or second year of life (53%). After accounting for all available sources of phenotypic variability with the linear mixed model for repeated records, it was estimated that AFL had a statistically significant effect only on DIM (p < 0.001). The litter size had a significant effect only on TMY (p < 0.001), while the effect of the parity was significant for all the examined traits (p < 0.001). The results of the study suggest that prolongation of age at first mating to the second year of life is not justified in dairy-orientated sheep farms. However, more evidence on this issue is needed for generalization, especially considering some other traits that can impact profitability of dual-purpose sheep farms (reproduction traits, growth rate of lambs, etc.).

Lactation curve model with explicit representation of perturbations as a phenotyping tool for dairy livestock precision farming

In the context of dairy farming, ruminant females often face challenges inducing perturbations that affect their performance and welfare. A key issue is how to assess the effect of perturbations and provide metrics to quantify how animals cope with their environment. Milk production dynamics are good candidates to address this issue: i) they are easily accessible, ii) overall dynamics throughout lactation process are well described and iii) perturbations are visible through milk losses. In this study, a perturbed lactation model (PLM) with explicit representation of perturbations was developed. The model combines two components: i) the unperturbed lactation model that describes a theoretical lactation curve, assumed to reflect female production potential and ii) the perturbation model that describes all the deviations from the unperturbed lactation model with four parameters: starting date, intensity and shape (collapse and recovery). To illustrate the use of the PLM as a phenotyping tool, it was fitted on a data set of 319 complete lactations from 181 individual dairy goats. A total of 2 354 perturbations were detected, with an average of 7.40 perturbations per lactation. Loss of milk production for the whole lactation due to perturbations varied between 2 and 19% of the milk production predicted by the unperturbed lactation model. The number of perturbations was not the major factor explaining the loss of milk production, suggesting that there are different types of animal response to challenges. By incorporating explicit representation of perturbations in a lactation model, it was possible to determine for each female the potential milk production, characteristics of each perturbation and milk losses induced by perturbations. Further, it was possible to compare animals and analyze individual variability. The indicators produced by the PLM are likely to be useful to move from raw data to decision support tools in dairy production.

Milk production losses in Latxa dairy sheep associated with small ruminant lentivirus infection

Visna/Maedi is a disease of sheep caused by small ruminant lentivirus (SRLV) infection that is widespread throughout the world and that has been recognized to be present in the Basque Country (Spain) since the early 1980's. Nearly seven decades of studies have improved the knowledge on its clinical signs and epidemiology. However, its slow progressive nature, subclinical most of the time, makes difficult to assess its real impact on productive traits, a question of critical importance to balance out the economic costs it causes and the benefits of designing and deploying an eradication program. Development of a dairy breeding program since the 90 s in the local Latxa sheep population has provided data on milk productivity in several flocks where SRLV infection prevalence has been continuously monitored. This study analyses retrospectively the association between SRLV prevalence and production variables during ten yearly lactations in three Latxa dairy flocks with medium-high SRLV seroprevalence. Our results indicate that average standard lactation of seropositive sheep was 6.7 % lower than controls. The largest differences (p < 0.001) were observed at the ewe lifetime peak of production between second and fourth lactations. Lifelong milk and lamb production data indicated even a higher impact, with costs rising up to nearly 50 €/ewe/year. This substantial production decrease associated with subclinical SRLV infection in Latxa dairy sheep supports the benefit of establishing a SRLV control program. A rough cost-benefit analysis indicated that even in a medium-yielding breed, testing expenses would be largely covered by milk production improvement.

Bulk-Tank Milk Longitudinal Serosurvey Reveals Decreasing Schmallenberg Virus Circulation in a Confined Population of Portuguese Sheep

Introduction: In 2011, Schmallenberg virus (SBV) was first detected in dairy cattle herds in The Netherlands and Germany having since then spread across Europe. Today studies are starting to show a decrease in new SBV infections, a circumstance that raises alerts for possible re-emergence if ideal conditions for vector development occur. To assess the potential decrease in SBV circulation, we performed a 2-year longitudinal serological investigation for SBV infection at the herd level by using bulk-tank milk of a specific sheep breed from central Portugal. Materials and Methods: Bulk-tank milk samples from 68 flocks were collected in both 2015 and 2016, and lactosera were tested for IgG anti-SBV by EIA. Results and Discussion: Results show that in 2015, 92.6% (95% confidence interval [CI]: 83.9-96.8) of the bulk-tank milk samples were positive, whereas in 2016 only 77.9% (95% CI: 66.7-86.1 of the samples from the same flocks were positive. Differences in the 2015/2016 seroprevalences showed to be statistically significant (p = 0.027). This significant decrease at the herd level seems to be in agreement with reported data from other European countries and raise alerts, since increasingly favorable conditions (higher number of susceptible animals) are now present, potentially favoring SBV epidemics if improved conditions for midge replication occur in the future.

Evaluation of efficacy of a biofilm-embedded bacteria-based vaccine against staphylococcal mastitis in sheep-A randomized, placebo-controlled field study

Our objective was to evaluate the efficacy of a vaccine against staphylococcal mastitis in 5 dairy sheep farms, with 316 ewes in the vaccinated (V) group and 307 in the control (C) group studied throughout a lactation period. Two administrations of the vaccine were performed during the last stage of gestation of ewes. Starting 15 d after lambing and at monthly intervals thereafter, up to 9 milk samplings were performed for bacteriological and cytological examinations. Staphylococcal isolates recovered were examined for biofilm formation. Blood samples were collected for measurement of IgG poly-N-acetylglucosamine-specific antibodies. The most frequently isolated bacteria were staphylococci: 56.4 and 76.1%, respectively, of total isolates recovered from ewes of group V and C, respectively; staphylococci as causal agents of mastitis were isolated less frequently from V (5.3%) than in ewes in C (10.3%). Among mastitis-associated staphylococcal isolates recovered from V ewes, a smaller proportion was biofilm-forming than among ones from C: 53.2% versus 74.9% of isolates; biofilm-forming staphylococci as causal agents of mastitis were isolated less frequently from ewes in group V (2.3%) than in ewes in group C (6.0%). Anti-poly-N-acetylglucosamine-specific antibody values increased in V ewes and were higher than in C; a greater proportion of ewes with low antibody titers developed staphylococcal mastitis (41.4%) than of V ewes with high antibody titers (17.0%). Incidence risk of mastitis, staphylococcal mastitis, and biofilm-associated staphylococcal mastitis was smaller in V than in C: 36.7, 17.1, and 8.0% versus 44.3, 30.9, and 18.9%, respectively. The first case of staphylococcal mastitis occurred later in V than in C: third versus second sampling point. Overall, efficacy of the vaccine was 44.6% for staphylococcal mastitis, 57.7% for biofilm-associated staphylococcal mastitis, 33.1% for staphylococcal intramammary infection, and 51.5% for biofilm-associated staphylococcal intramammary infection. Nevertheless, vaccination should not be the only means for controlling mastitis; other udder health management measures should be included therein to improve control of the infection.

Factors affecting the milk yield and composition over lactation of prolific D'man ewes in Tunisian oases

This work aimed to investigate milk potentiality and the influence of some non-genetic factors both on milk yield production and composition in prolific D'man sheep. Observations were performed in eighty lactating D'man ewes maintained under an intensive oases system in Southern Tunisia. Daily milk yield was determined at 7-day intervals using oxytocin + hand milking method. Chemical milk composition was determined by Lactoscan. Total productions over lactation (milk, fat, and protein) were estimated using the Fleischmann method. Daily milk yield averaged 1.64 ± 0.02 l. It increased to a plateau from the third to the fourth week of lactation, and decreased gradually to the tenth week. D'man ewes produced 128.91 ± 3.14 l of milk, 9.23 ± 0.26 kg of fat, and 5.26 ± 0.12 kg of total proteins. Overall, milk contained 16.17% of total solids, 7.08% of fat, 4.04% of proteins, 4.32% of lactose, 9.09% of solids-not-fat, and 0.73% of ash. Lambing season affected the daily milk yield and milk composition, where fat and total proteins were highest in autumn. Total solids, lactose, solids-not-fat, and ash were higher in winter than in autumn and summer. Ewes reared multiple lambs produced more milk than those reared singles. Dams nursing triplets or more had the highest milk total proteins, and the lowest total solids and fat levels. Adult ewes produced more milk compared to younger and older ewes, and their milk contained higher total solids, fat, and ash levels than milk from young ewes. Ewes with asymmetric udder produced less and more concentrated milk than ewes with symmetrical udder. Daily milk yield was negatively correlated to total solids, fat, and total proteins amounts. In contrast, litter survival at weaning was positively correlated to the amounts of total milk, fat, and total proteins.

Causal effect of prolificacy on milk yield in dairy sheep using propensity score

In animal production, it is often important to investigate causal relationships among variables. The gold standard tool for such investigation is randomized experiments. However, randomized experiments may not always be feasible, possible, or cost effective or reflect real-world farm conditions. Sometimes it is necessary to infer effects from farm-recorded data. Inferring causal effects between variables from field data is challenging because the association between them may arise not only from the effect of one on another but also from confounding background factors. Propensity score (PS) methods address this issue by correcting for confounding in different levels of the causal variable, which allows unbiased inference of causal effects. Here the objective was to estimate the causal effect of prolificacy on milk yield (MY) in dairy sheep using PS based on matched samples. Data consisted of 4,319 records from 1,534 crossbred ewes. Confounders were lactation number (first, second, and third through sixth) and dairy breed composition (<0.5, 0.5-0.75, and >0.75 of East Friesian or Lacaune). The causal variable prolificacy was considered as 2 levels (single or multiple lambs at birth). The outcome MY represented the volume of milk produced in the whole lactation. Pairs of single- and multiple-birth ewes (1,166) with similar PS were formed. The matching process diminished major discrepancies in the distribution of prolificacy for each confounder variable indicating bias reduction (cutoff standardized bias = 20%). The causal effect was estimated as the average difference within pairs. The effect of prolificacy on MY per lactation was 20.52 L of milk with a simple matching estimator and 12.62 L after correcting for remaining biases. A core advantage of causal over probabilistic approaches is that they allow inference of how variables would react as a result of external interventions (e.g., changes in the production system). Therefore, results imply that management and decision-making practices increasing prolificacy would positively affect MY, which is important knowledge at the farm level. Farm-recorded data can be a valuable source of information given its low cost, and it reflects real-world herd conditions. In this context, PS methods can be extremely useful as an inference tool for investigating causal effects. In addition, PS analysis can be implemented as a preliminary evaluation or a hypothesis generator for future randomized trials (if the trait analyzed allows randomization).

Equid milk production: evaluation of Martina Franca jennies and IHDH mares by Wood’s model application

The present work aims to evaluate milk production aptitude of two equid species: horse and donkey. The trial was conducted on the Italian Heavy Draught breed horse (IHDH) and on the Martina Franca breed donkey. Data were collected from 23 multiparous nursing IHDH mares aged from 6 to 10 years, with an average bodyweight (BW) (mean ± s.d.) of 822.3 ± 61.8 kg and from 28 multiparous nursing donkeys of Martina Franca breed aged from 7 to 12 years, with an average BW (mean ± s.d.) of 389.5 ± 21.2 kg. Milk yield was recorded on 15, 30, 45, 60, 75, 90, 105, 120, 135, 150, 165 and 180 days in milk (DIM). At each milk yield evaluation day, a milk sample was collected and its composition was analysed. All data were normalised according metabolic weight (MBW) of animals, so to make comparable results of the two considered species. Wood’s lactation model was used to describe milk, fat, protein and lactose yield. Peak lactation was reached on DIM 69 in IHDH horse and on DIM 48 in Martina Franca donkey with a milk production of 10.693 kg/100 kg MBW (R2 = 0.9139) and 5.264 kg/100 kg MBW (R2 = 0.9428), respectively. Persistency index of milk lactation curve was higher in Martina Franca donkey then in IHDH mares (7.0338 vs 6.2672). IHDH mares showed higher fat and protein production. Lactose curves showed a similar trend in both the species although IHDH’s curve was characterised by higher values. Data reported in the present work underline that both breeds considered are characterised by good milk production aptitude. Moreover, high R2 values obtained applying Wood’s model to milk, fat, protein and lactose production, calculated on the metabolic bodyweight basis, are indices of goodness of model fit.

Invited review: Cessation of lactation: Effects on animal welfare

The forced cessation of milk production, or dry-off, is a routine management practice in dairy cattle, sheep, and goats. This practice initiates a dry period, during which the animal is not milked. Milking begins again after parturition. Most of the literature on the dry period has focused on how various drying-off strategies affect milk production and disease; little work to date has addressed how dry-off affects the overall welfare of the dairy animal. The first aim of this review was to present an overview of the importance of dry-off and how it is commonly achieved. Our review shows that much scientific progress has been made in improving health status between lactations. The second aim was to identify important gaps in the literature, of which 2 key research disparities have been identified. We find that much of the work to date has focused on cattle and very little research has examined dry-off in dairy sheep and goats. We also find a lack of research addressing how common dry-off methodologies affect animal welfare on more than just a biological level, regardless of species.

Effects of crude protein level in the concentrate and time allotment on pasture on milk yield, urinary nitrogen, and purine derivative excretion in lactating Latxa ewes

This study assessed the influence of reducing the crude protein (CP) content (from 190 to 130 g/kg CP/kg) in the supplementary concentrate and time allotment on pasture (TAP) on the milk yield, bodyweight, and the urinary nitrogen (N) and purine derivatives (PD) in lactating Latxa ewes. Animals were reared in a production system that restricted the amount of time that ewes spent on pasture. In mid-April, at the start of the 42 days experiment, 40 dairy Latxa ewes were assigned to one of four groups on the basis of their initial milk yield, days in lactation (DIL), bodyweight (BW), and condition score; thereafter, each group grazed in its own fenced paddock. The experiment was based on a 2 × 2 factorial design that included two CP levels and two pasture grazing regimes: 4 h continuous grazing (CG) in the morning or 2 h grazing in the morning and 2 h in the afternoon (DG). Individual milk yield was recorded three times a week, and time spent grazing and BW was recorded weekly. In the middle (day 15–17; P1) and at the end (day 36–38; P2) of the experiment, urinary spot samples were collected using a catheter. At the end of the experiment, ewes were confined to metabolic cages and urine was collected. CP level of the concentrate was not correlated with time spent grazing; however, the ewes that were permitted access to pasture twice per day spent more time grazing (223 min/day vs 207 min/day, P < 0.01) and were more efficient with their time (56 min/h grazing vs 52 min/h, P < 0.05) than the ewes that were permitted access to pasture once per day. Concentrate CP levels were not correlated with milk yield or composition, although ewes that received the high protein (HP) lost more weight than did those that received the low protein (LP) concentrate. Concentrate CP level and TAP were not correlated with creatinine (CR) excretion rate (mean = 315 µmol/kg LW0.75 s.e. 0.0161). Urea-N was the largest component of urinary-N (68.7 s.e. 2.33%; P > 0.05), and urea-N waste was higher in HP ewes (202.7 mmol/day) than it was in LP ewes (159.5 mmol/day) (s.d. 27.83; P < 0.01). Ewes subjected to the DG regime had significantly (P < 0.01) higher urinary PD excretion (23.6 vs 21.4 mmol/day s.d. 4.01; P < 0.01), and tended (P < 0.1) to excrete less urea-N (175.1 vs 188.3 mmol, s.d. 27.83) than did the CG ewes. A reduction in the CP in the supplementary concentrate led to a reduction in N waste without having a detrimental effect on performance or milk production.

Application of the Wood model to analyse lactation curves of organic dairy sheep farming

The aim of the present study was to evaluate the performance of the Wood model to describe the characteristics of lactation curves of dairy ewes under organic management in Mexico. In total, 4861 weekly test-day milk yield records from 194 lactations of crossbred dairy ewes were analysed to assess the performance of an empirical model to fit their lactation curve. We used the mathematical model proposed by Wood. The evaluation criteria were the correlation coefficient (r) between the values of total milk yield observed and estimated, the coefficient of determination (R2), and the mean square prediction error (MSPE). In addition, the peak yield (PYest) and time at peak yield (PTest) were calculated. The Wood model showed adequate goodness of fit (r = 0.95, R2 = 0.92 and MSPE = 0.024). The Wood model detected that 52.06% of lactation curves had a continuously decreasing shape (atypical curve), probably as a consequence of the characteristic management of the organic system, mainly due to the genotype used and the nutritional management. Residuals were greater for atypical curves than for typical ones, indicating differences in the ability of the Wood model to fit the two types of shapes. In typical curves, the Wood model showed adequate estimates of total milk yield and time at peak yield. The peak yield was underestimated both in typical and atypical curves. The Wood model in atypical curves underestimated the time at peak yield and milk yields in late lactation. The Wood model showed a reasonable fit of lactation curve in dairy sheep in organic systems but presented deficiencies of fit in atypical curves; therefore, estimates should be interpreted carefully.

Effects of crossbreeding on milk production and composition in dairy sheep under organic management

The crossbreeding of local sheep breeds with dairy breeds is an option to improve dairy production parameters in organic sheep dairy systems. Weekly milk yield (WMY) was recorded and individual samples of milk for chemical analysis were taken during 17 weeks from 45 dairy ewes of the following three genotypes: 15 East Friesian (EF), 15 EF × Suffolk (EF × SF) and 15 EF × Pelibuey (EF × PL) under organic management. For analysis of the lactation curve the Wood gamma model was used. The effect of genotype on the WMY was analysed using repeated-measures. The comparison of the least square means among genotypes for total milk yield (TMY), daily milk yield, protein content, protein yield, fat content, fat yield, non-fat solids concentration, non-fat solids yield, total solids yield and acidity was analysed using a general linear model. The genetic group influenced only in the ascent phase of the lactation curve, with values of the Parameter b of model Wood higher in EF (P = 0.01). There were no differences (P > 0.05) between genotypes in relation to the WMY, TMY, protein content and acidity; however, the effects of week of lactation trial and the interaction of genotype and week of lactation trial on WMY were significant (P < 0.05). The values of daily milk yield, fat yield, protein yield and total solids yield were higher (P < 0.005) in EF and EF × SF than EF × PL. Fat content was higher in EF × PL. EF × SF had similar values of TMY than EF and better chemical composition, which places this genotype as an option of crossbreeding in dairy sheep systems under organic management with similar agro climatic characteristics to the present study.

Factors affecting the lactation curves of intensively managed sheep based on a clustering approach

This study investigated factors affecting milk production and lactation curves from complete lactations of Lacaune dairy sheep. Animals were part of a single flock under intensive management and were milked twice daily starting at lambing. The results of the analyses of 7788 complete lactations showed an average total milk yield of 434±183 l from lactations 234±63 d long, with an average lambing interval of 302±65 d. A Pollott additive mathematical model was used to estimate complete lactation curves. Clustering analysis identified four lactation types among Lacaune dairy sheep differing mainly in productivity i.e. milk yield per lactation (MY) and length of lactation (DIM). The so-called SL type involved short, less productive lactations (n=2137; 27·4%; MY=222±75·5 l and DIM=182±52·9 d). The SN type involved short lactations of normal productivity (n=2039; 26·2%; MY=396±73·7 l and DIM=205±33·1 d). The LP type involved long and productive lactations (n=2169; 27·9%; MY=487±70·5 l and DIM=265±40·7 d), while the LVP type included long and extremely productive lactations (n=1443; 18·5%; MY=694±114·0 l and DIM=295±54·7 d). Sheep showing the best lactation curves were usually younger than other sheep, and they had higher yield during the previous lactation, a shorter previous dry period (55±50·4 for LP and 61±55·0 d for LVP types) and longer lambing intervals. In addition, they tended to be born in September and to lamb in March, October and December. Sheep were remarkably stable in their lactation curve behaviour: the curve type observed for the first lactation was highly likely to persist in subsequent lactations (P<0·0001). These results suggest that farmers can use the shape of the first lactation curve to guide their selection of ewes for breeding and retention on the farm, thereby improving flock productivity.

Comparison of mathematical models applied to f1 dairy sheep lactations in organic farm and environmental factors affecting lactation curve parameter

THE OBJECTIVE OF THIS STUDY WAS TO COMPARE THE GOODNESS OF FIT OF FOUR LACTATION CURVE MODELS: Wood's Gamma model (WD), Wilmink (WL), and Pollott's multiplicative two (POL2) and three parameters (POL3) and to determine the environmental factors affecting the complete lactation curve of F1 dairy sheep under organic management. A total of 5,382 weekly milk yields records from 150 ewes, under organic management were used. Residual mean square (RMS), determination coefficients (R(2)), and correlation (r) analysis were used as an indicator of goodness of fit for each model. WL model best fitted the lactation curves as indicated by the lower RMS values (0.019), followed by WD (0.023), POL2 (0.025) and POL3 (0.029). The four models provided total milk yield (TMY) estimations that were highly correlated (0.93 to 0.97) with observed TMY (89.9 kg). The four models under estimated peak yield (PY), whereas POL2 and POL3 gave nearer peak time lactation estimations. Ewes lambing in autumn had higher TMY and showed a typical curve shape. Higher TMY were recorded in second and third lambing. Season of lambing, number of lambing and type of lambing had a great influenced over TMY shaping the complete lactation curve of F1 dairy sheep. In general terms WL model showed the best fit to the F1 dairy sheep lactation curve under organic management.

Influence of dry period length on reproductive performance and productivity of Lacaune dairy sheep under an intensive management system

Intensive management is almost the only way to ensure dairy farm profitability. The dry period length (DPL) is a key factor in the productivity and health of dairy cows, but whether the same is true of dairy sheep is unclear. This study investigated the effects of DPL on the performance of Lacaune sheep under intensive management. We recorded 8136 lactations from 4220 ewes on one farm for the period 2005-2010, and data from a total of 6762 complete lactations 1-4 were included in the study. The length of the dry period following the current lactation was studied. The larger the total milk yield (MY) and daily milk yield (DMY), the shorter was the DPL before the next lactation. DPL correlated with MY (r=-0·384), DMY (r=-0·277) and the lambing-to-conception interval (LC; r=0·201, P<0·0001) in the global analysis of all lactations (lactations 1-4). The influence of previous-DPL (P-DPL), or the length of the period prior to the start of the next lactation, was studied for 4318 lactations. P-DPL was classified into five intervals: very short (P-DPL-XS), 1-30 d; short (P-DPL-S), 31-60 d; medium (P-DPL-M), 61-90 d; long (P-DPL-L), 91-120 d; and very long (P-DPL-XL), >120 d. P-DPL positively correlated with lambing-to-next conception interval (LNC; r=0·095, P<0·0001) for lactations 1-4. LNC was significantly shorter for P-DPLs that were very short, short, or long (P-PDL-XS, 144·2±67·8 d; P-PDL-S, 149·1±57·2 d; P-PDL-L, 152·0±53·7 d) than for groups with very long or medium P-PDLs (P-DPL-XL, 161·5±62·9 d; P-DPL-M, 169·0±74·8 d; P<0·0001). Moreover, P-DPLs that were very short, long, or very long were associated with the lowest milk yields (P-PDL-XS, 377±215 l; P-PDL-l, 370±168 l; P-PDL-XL, 396±196 l). These yields were significantly lower than the yields for short and medium P-DPLs (P-DPL-S, 432±187 l; P-DPL-M, 436±191 l; P<0·0001) when averages of lactations 1-4 were analysed. These results indicate that lactations with larger MY are followed by a shorter dry period, and that a dry period of 30-90 d leads to larger yields in the next lactation. The best LNC was associated with the shortest Previous-DPL. Hence, 30-60 d should be the optimal dry period length for Lacaune sheep under intensive conditions.

Application of Wood's model to lactation curve of Italian Heavy Draft horse mares

The aim of this work was to study, using Wood's model, the lactation curve in Italian Heavy Draft horse (IHDH) mares, both for milk yield and composition. Interest in mare's milk has grown for use in human nutrition, as a substitute for human and cow milk for premature newborns and allergic children. For this study, 7 IHDH mares were used. Milk yield was evaluated from d 15 to 180 postpartum (every 15 d). Wood's lactation model was used to describe milk yield, fat, protein, and lactose lactation curves in the first 180 d in milk (DIM), with a prediction to 210 DIM. Mean values recorded across the considered lactation period for milk yield and constituents were 16.11 kg of milk/d, 11.38 g of fat/kg, 15.54 g of protein/kg, and 69.72 g of lactose/kg. The lactation curve showed a peak at 69 DIM, with a persistence index of 6.26. Fat (g/kg) and protein (g/kg) lactation curves decreased during the lactation period; instead, the lactose (g/kg) curve increased from the first to the last DIM. The coefficient of determination values indicated a good model fit by Wood's model application to IHDH mare lactation. The research showed good aptitude of a heavy horse breed, such as IHDH, for milk production.

Mapping quantitative trait loci (QTL) in sheep. II. Meta-assembly and identification of novel QTL for milk production traits in sheep

An (Awassi × Merino) × Merino backcross family of 172 ewes was used to map quantitative trait loci (QTL) for different milk production traits on a framework map of 200 loci across all autosomes. From five previously proposed mathematical models describing lactation curves, the Wood model was considered the most appropriate due to its simplicity and its ability to determine ovine lactation curve characteristics. Derived milk traits for milk, fat, protein and lactose yield, as well as percentage composition and somatic cell score were used for single and two-QTL approaches using maximum likelihood estimation and regression analysis. A total of 15 significant (P < 0.01) and additional 25 suggestive (P < 0.05) QTL were detected across both single QTL methods and all traits. In preparation of a meta-analysis, all QTL results were compared with a meta-assembly of QTL for milk production traits in dairy ewes from various public domain sources and can be found on the ReproGen ovine gbrowser http://crcidp.vetsci.usyd.edu.au/cgi-bin/gbrowse/oaries_genome/. Many of the QTL for milk production traits have been reported on chromosomes 1, 3, 6, 16 and 20. Those on chromosomes 3 and 20 are in strong agreement with the results reported here. In addition, novel QTL were found on chromosomes 7, 8, 9, 14, 22 and 24. In a cross-species comparison, we extended the meta-assembly by comparing QTL regions of sheep and cattle, which provided strong evidence for synteny conservation of QTL regions for milk, fat, protein and somatic cell score data between cattle and sheep.

Winter/spring changes in fatty acid composition of farmhouse Idiazabal cheese due to different flock management systems

Typically, two different flock managements are employed by basque sheepherders in winter and spring. Thus, seasonal changes in the fatty acid (FA) composition of Idiazabal PDO farmhouse cheeses were studied. Ewe's raw milk cheeses elaborated in winter and spring were collected after 120 days of ripening from 10 Idiazabal PDO farmhouses. In winter, concentrate and conserved forages were fed, whereas a part-time grazing system was adopted from spring onward. Spring cheeses had less (P <or= 0.05) saturated FA and higher (P <or= 0.05) content of unsaturated FA, including trans-FA (mainly trans-vaccenic acid) and conjugated linoleic acid (CLA), branched-chain FA (BCFA), and n-3 FA. Principal component analysis (PCA) separated winter and spring cheeses into two groups by the combination of two principal components (84.2% of variance). Fresh pasture in the diet enhanced desirable FA and lowered atherogenicity index in cheeses, supporting the benefits of using a part-time grazing system for the consumer.

Modeling extended lactations of Holsteins

Modeling extended lactations for the US Holsteins is useful because a majority (>55%) of the cows in the present population produce lactations longer than 305 d. In this study, 9 empirical and mechanistic models were compared for their suitability for modeling 305-d and 999-d lactations of US Holsteins. A pooled data set of 4,266,597 test-day yields from 427,657 (305-d complete) lactation records from the AIPL-USDA database was used for model fitting. The empirical models included Wood (WD), Wilmink (WIL), Rook (RK), monophasic (MONO), diphasic (DIPH), and lactation persistency (LPM) functions; Dijkstra (DJ), Pollott (POL), and new-multiphasic (MULT) models comprised the mechanistic counterparts. Each model was separately tested on 305-d (>280 days in milk) and 999-d (>800 days in milk) lactations for cows in first parity and those in third and greater parities. All models were found to produce a significant fit for all 4 scenarios (2 parity groups and 2 lactation lengths). However, the resulting parameter estimates for the 4 scenarios were different. All models except MONO, DIPH, and LPM yielded residuals with absolute values smaller than 2 kg for the entire period of the 305-d lactations. For the extended lactations, the prediction errors were larger. However, the RK, DJ, POL, and MULT models were able to predict daily yield within a +/- 3 kg range for the entire 999-d period. The POL and MULT models (having 6 and 12 parameters, respectively) produced the lowest mean square error and Bayesian information criteria values, although the differences from the other models were small. Conversely, POL and MULT were often associated with poor convergence and highly correlated, unreliable, or biologically atypical parameter estimates. Considering the computational problems of large mechanistic models and the relative predictive ability of the other models, smaller models such as RK, DJ, and WD were recommended as sufficient for modeling extended lactations unless mechanistic details on the extended curves are needed. The recommended models were also satisfactory in describing fat and protein yields of 305-d and 999-d lactations of all parities.

Modeling extended lactation curves of dairy cattle: a biological basis for the multiphasic approach

Objectives of this study are to describe the biological basis for multiphasic milk production and to propose a new empirical model for the lactation curve. To illustrate this model, we used data on 3573 first-lactation Holsteins having lactations of various lengths (285, 345, 405, 465, and 525 d) and with various days open (45, 105, 165, 225, and 285 d). The model describes an increasing first phase of milk yield and a series of decreasing phases of yield. The increasing phase, described by an increasing logistic function of time, is associated with increase in number of active mammary gland cells and increase in yield per cell. The decreasing phases, described by three decreasing logistic functions of time, are associated with decreases in cell number due to apoptosis and in yield per cell due to pregnancy. The new model is [equation: see text] where yDIM is milk yield at each day in milk (DIM), a1 is upper level for the increasing first phase, and p2, p3, and p4 = (1 - p2 - p3) are proportions of a1 for the decreasing second, third, and fourth phases; b's are proportional to duration of each phase; and c's are time of maximum increase or decrease. Nonlinear regression was used to fit average milk yield for each of nine datasets, four with 180 d carried calf and five with 240 d carried calf. Average results indicated that for the first phase, upper level of milk yield was about 22 kg. Duration was about 120 d, centered on time of maximum increase, which was about 11 d before calving. For the second phase (first phase of apoptosis), decrease in yield was relatively large (about 20%) and duration was relatively long (about 375 d). Time of maximum decrease was about 107 d after calving. For the third phase (pregnancy), decrease in yield was relatively small (about 6%) and duration was relatively short (about 200 d). Time of maximum decrease was about 300 d after calving. For each additional day open, time of maximum decrease increased about 1 d. For the fourth phase (second phase of apoptosis), decrease in yield was relatively large (about 74%) and duration was relatively long (about 765 d). Duration for the lactation length of 525 d was exceptionally long. Time of maximum decrease was about 382 d after calving. For each additional day of lactation, duration increased about 5.4 d and day of maximum decrease increased about 0.82 d. We believe that it is possible to model empirically standard and extended lactation curves of dairy cows, based on biological theory and predicated on the multiphasic approach. Further research to understand better the biology of extended lactations, using the proposed multiphasic model, should use planned extended lactations that are at least 525 d in milk and have at least 240 d carried calf.