Boulbaba Rekik

Fitting lactation curves of dairy cattle in different types of herds in Tunisia

The incomplete gamma function was used to fit lactation curves of Holstein–Friesian cows in four types of herds in Tunisia. A total of 8640 records were used in the analysis. These included 1269, 637, 239, and 498 first lactation and 2986, 1441, 650, and 920 second and later lactation records in four herd groups namely investors, state, cooperative, and farmers’ herds, respectively. The effects of environmental variables, production sector, herd, parity, first test-day date, calving year, and calving season on the main lactation curve traits were analysed. The factors associated with milk yield at the beginning of lactation and the decreasing phase of the curve, persistency, and peak yield varied significantly (P<0.01) with all variables. The ascending phase of the lactation curve was not affected by parity and calving year, while days in milk until peak depended only on the rank of lactation. The state herds had the lowest peak and total yields. The summer season was unfavourable for milk production. In contrast to first lactation cows, third lactation cows had the highest peak and total yields. Milk yield was highly correlated with peak yield (r=0.79) and was not related to persistency measure.

Genetic parameters for Tunisian Holsteins using a test-day random regression model.

Genetic parameters of milk, fat, and protein yields were estimated in the first 3 lactations for registered Tunisian Holsteins. Data included 140,187; 97,404; and 62,221 test-day production records collected on 22,538; 15,257; and 9,722 first-, second-, and third-parity cows, respectively. Records were of cows calving from 1992 to 2004 in 96 herds. (Co)variance components were estimated by Bayesian methods and a 3-trait-3-lactation random regression model. Gibbs sampling was used to obtain posterior distributions. The model included herd x test date, age x season of calving x stage of lactation [classes of 25 days in milk (DIM)], production sector x stage of lactation (classes of 5 DIM) as fixed effects, and random regression coefficients for additive genetic, permanent environmental, and herd-year of calving effects, which were defined as modified constant, linear, and quadratic Legendre coefficients. Heritability estimates for 305-d milk, fat and protein yields were moderate (0.12 to 0.18) and in the same range of parameters estimated in management systems with low to medium production levels. Heritabilities of test-day milk and protein yields for selected DIM were higher in the middle than at the beginning or the end of lactation. Inversely, heritabilities of fat yield were high at the peripheries of lactation. Genetic correlations among 305-d yield traits ranged from 0.50 to 0.86. The largest genetic correlation was observed between the first and second lactation, potentially due to the limited expression of genetic potential of superior cows in later lactations. Results suggested a lack of adaptation under the local management and climatic conditions. Results should be useful to implement a BLUP evaluation for the Tunisian cow population; however, results also indicated that further research focused on data quality might be needed.

Genetic analysis of heat stress effects on yield traits, udder health, and fatty acids of Walloon Holstein cows

Genetic parameters that considered tolerance for heat stress were estimated for production, udder health, and milk composition traits. Data included 202,733 test-day records for milk, fat, and protein yields, fat and protein percentages, somatic cell score (SCS), 10 individual milk fatty acids (FA) predicted by mid-infrared spectrometry, and 7 FA groups. Data were from 34,468 first-lactation Holstein cows in 862 herds in the Walloon region of Belgium and were collected between 2007 and 2010. Test-day records were merged with daily temperature-humidity index (THI) values based on meteorological records from public weather stations. The maximum distance between each farm and its corresponding weather station was 21km. Linear reaction norm models were used to estimate the intercept and slope responses of 23 traits to increasing THI values. Most yield and FA traits had phenotypic and genetic declines as THI increased, whereas SCS, C18:0, C18:1 cis-9, and 4 FA groups (unsaturated FA, monounsaturated FA, polyunsaturated FA, and long-chain FA) increased with THI. Moreover, the latter traits had the largest slope-to-intercept genetic variance ratios, which indicate that they are more affected by heat stress at high THI levels. Estimates of genetic correlations within trait between cold and hot environments were generally high (>0.80). However, lower estimates (<=0.67) were found for SCS, fat yield, and C18:1 cis-9, indicating that animals with the highest genetic merit for those traits in cold environments do not necessarily have the highest genetic merit for the same traits in hot environments. Among all traits, C18:1 cis-9 was the most sensitive to heat stress. As this trait is known to reflect body reserve mobilization, using its variations under hot conditions could be a very affordable milk biomarker of heat stress for dairy cattle expressing the equilibrium between intake and mobilization under warm conditions.

Genotype x environment interaction for milk yield in Holsteins using Luxembourg and Tunisian populations.

Test-day (TD) milk yield records of first-lactation Holstein cows in Luxembourg and Tunisia were analyzed using within-and between-country random regression TD models. Edited data used for within-country analysis included 661,453 and 281,913 TD records in Luxembourg and Tunisia, respectively. The joint data included 730,810 TD records of 87,734 cows and 231 common sires. Both data sets covered calving years 1995 to 2006. Fourth-order Legendre polynomials for random effects and a Gibbs sampling method were used to estimate variance components of lactation curve parameters in separate and joint analyses. Genetic variances of the first 3 coefficients from Luxembourg data were 46 to 69% larger than corresponding estimates from the Tunisian data. Inversely, the Tunisian permanent environment variances for the same coefficients were 52 to 65% larger than the Luxembourg ones. Posterior mean heritabilities of 305-d milk yield and persistency, defined as estimated breeding values (EBV) at 280 days in milk-EBV at 80 days in milk, from between-country analysis were 0.42 and 0.12 and 0.19 and 0.08 in Luxembourg and Tunisia, respectively. Heritability estimates for the same traits from within-country analyses, mainly from the Tunisian data, were lower than those from the joint analysis. Genetic correlations for 305-d milk yield and persistency between countries were 0.60 and 0.36. Product moment and rank correlations between EBV of common sires for 305-d milk yield and persistency from within-country analyses were 0.38 and 0.41 and 0.27 and 0.26, respectively. Differences between genetic variances found in both countries reflect different milk production levels. Moreover, low genetic and rank correlations suggest different ranking of sires in the 2 environments, which implies the existence of a genotype x environment interaction for milk yield in Holsteins.

Accessing genotype by environment interaction using within- and across-country test-day random regression sire models.

First-lactation test-day (TD) milk records of Luxembourg and Tunisian Holsteins were analysed for evidence of genotype by environment interaction (G x E). The joint data included 730 810 TD records of 87 734 cows and 231 common sires. Random regression TD sire models with fourth-order Legendre polynomials were used to estimate genetic parameters via within- and across-country analyses. Daily heritability estimates of milk yield from within-country analysis were between 0.11 and 0.32, and 0.03 and 0.13 in Luxembourg and Tunisia, respectively. Heritability estimates for 305-day milk yield and persistency (defined as the breeding value for milk yield on DIM 280 minus the breeding value on DIM 80) were lower for Tunisian Holsteins compared with the Luxembourg population. Specifically, heritability for 305-day milk yield was 0.16 for within- and 0.11 for across-country analyses for Tunisian Holsteins and 0.38 for within- and 0.40 for across-country analyses for Luxembourg Holsteins. Heritability for apparent persistency was 0.02 for both within- and across-country analyses for Tunisian Holsteins and 0.08 for within- and 0.09 for across-country analyses for Luxembourg Holsteins. Genetic correlations between the two countries were 0.50 for 305-day milk yield and 0.43 for apparent persistency. Moreover, rank correlations between the estimated breeding values of common sires for 305-day milk yield and persistency, estimated separately in each country, were low. Low genetic correlations are evidence for G x E for milk yield production while low rank correlations suggest different rankings of sires in both environments. Results from this study indicate that milk production of daughters of the same sires depends greatly on the production environment and that importing high merit semen for limited input systems might not be an effective strategy to improve milk production.

Genetic analysis of milk production traits of Tunisian Holsteins using random regression test-day model with Legendre polynomials

Objective The objective of this study was to estimate genetic parameters of milk, fat, and protein yields within and across lactations in Tunisian Holsteins using a random regression test-day (TD) model. Methods A random regression multiple trait multiple lactation TD model was used to estimate genetic parameters in the Tunisian dairy cattle population. Data were TD yields of milk, fat, and protein from the first three lactations. Random regressions were modeled with third-order Legendre polynomials for the additive genetic, and permanent environment effects. Heritabilities, and genetic correlations were estimated by Bayesian techniques using the Gibbs sampler. Results All variance components tended to be high in the beginning and the end of lactations. Additive genetic variances for milk, fat, and protein yields were the lowest and were the least variable compared to permanent variances. Heritability values tended to increase with parity. Estimates of heritabilities for 305-d yield-traits were low to moderate, 0.14 to 0.2, 0.12 to 0.17, and 0.13 to 0.18 for milk, fat, and protein yields, respectively. Within-parity, genetic correlations among traits were up to 0.74. Genetic correlations among lactations for the yield traits were relatively high and ranged from 0.78±0.01 to 0.82±0.03, between the first and second parities, from 0.73±0.03 to 0.8±0.04 between the first and third parities, and from 0.82±0.02 to 0.84±0.04 between the second and third parities. Conclusion These results are comparable to previously reported estimates on the same population, indicating that the adoption of a random regression TD model as the official genetic evaluation for production traits in Tunisia, as developed by most Interbull countries, is possible in the Tunisian Holsteins.