Opportunities and limits of commercial farm data to study the genetic determinism of feed efficiency throughout lactation in dairy sheep.

The collective economic and environmental interest of the whole dairy sheep sector is to reduce feed costs and the negative impact of milk production on the environment. Thus, this study focused on the characterisation and genetic selection potential of feed efficiency in the Lacaune breed. Estimates for feed efficiency in dairy ewes are limited, mainly due to a lack of individual feed intake measurements in the sheepfold or in the pasture. We estimated the genetic parameters for two approximated (not entirely based on individual data) feed efficiency traits (lactation feed conversion ratio (LFCR) and residual energy intake (REI)) and daily milk yield (DMY) at different stages of lactation and throughout lactation. The accuracy of the efficiency traits was first evaluated on samples from Lacaune dairy ewes that were monitored individually, especially for their feed intake. Then, feed efficiency estimation methods were applied on eight commercial farms corresponding to 4 680 Lacaune dairy ewes over two milk lactations (30 854 records). Animals were collectively (for a large part of feed intake) or individually (for milk performance and dynamics of body fat reserves) monitored at different lactation stages. The heritabilities of LFCR and REI were estimated over lactations at 0.10 ± 0.01 and 0.11 ± 0.01, respectively. High genetic correlations were observed between the two efficiency traits and milk production traits, with a genetic correlation between LFCR and DMY of 0.74 ± 0.04 and between REI and DMY of -0.79 ± 0.04. A strong influence of environmental factors such as farm, year of milk production and lactation stage affected the genetic link between REI and milk production traits. Efficiency values observed in early lactation when animals were bred in the sheepfold were less genetically correlated with values obtained later in lactation when animals were grass-fed. However, individual characterisation of feed efficiency remains difficult due to the collective feeding context in dairy ewe farms.

Partitioning of the genetic trends of French dairy sheep in Mendelian samplings and long-term contributions.

The genetic trend of milk yield for 4 French dairy sheep breeds (Lacaune, Basco-Béarnaise, Manech Tête Noire, and Manech Tête Rousse) was partitioned in Mendelian sampling trends by categories of animals defined by sex and by selection pathways. Five categories were defined, as follows: (1) artificial insemination (AI) males (after progeny testing), (2) males discarded after progeny testing, (3) natural mating males, (4) dams of males, and (5) dams of females. Dams of males and AI males were the most important sources of genetic progress, as observed in the decomposition in Mendelian sampling trends. The yearly contributions were more erratic for AI males than for dams of males, as AI males are averaged across a smaller number of individuals. Natural mating males and discarded males did not contribute to the trend in terms of Mendelian sampling, as their estimated Mendelian sampling term is either null (natural mating males) or negative (discarded males). Overall, in terms of Mendelian sampling, females contributed more than males to the total genetic gain, and we interpret that this is because females constitute a larger pool of genetic diversity. In addition, we computed long-term contributions from each individual to the following pseudo-generations (one pseudo-generation spanning 4 years). With this information, we studied the selection decisions (selected or not selected) for females, and the contributions to the following generations. Mendelian sampling was more important than parent average to determine the selection of individuals and their long-term contributions. Long-term contributions were greater for AI males (with larger progeny sizes than females) and in Basco-Béarnaise than in Lacaune (with the latter being a larger population).

High genetic correlation for milk yield across Manech and Latxa dairy sheep from France and Spain.

Spanish Latxa and French Manech are dairy sheep breeds that split into Blond (Latxa Cara Rubia, LCR; Manech Tête Rousse, MTR) and Black (Latxa Cara Negra of Navarre, LCN; Manech Tête Noire, MTN) strains. Exchange of genetic material (artificial insemination doses) is becoming more and more frequent across these breeds, within color, to boost both genomic precision using a larger reference population and genetic progress using a larger selection base. This exchange leads to some rams having descendance across both countries. However, additional gains can only be achieved if the selected traits are genetically similar across countries. The objective of this work was to estimate the genetic correlation across breeds for milk yield. We combine across-country, within-color records, pedigree, and marker information. The number of animals with records oscillates from 65,000 (LCN) to 544,000 (MTR), whereas the number of connecting artificial insemination rams (with more than 10 daughters in the other country) is 381 MTR rams in LCR and 58 MTN rams in LCN. Blond strains had a stronger and more extended-in-time connection. The number of genotyped rams goes from 328 (LCN) to 4,901 (MTR). The relatedness of populations was assessed by principal component analysis and Fst coefficients. The genetic correlation was estimated using 2 (one per color) 2-trait models (each country a trait), including all available data (records, pedigree and genotypes), by maximum profile likelihood while fixing other variance components to within-population estimates. Results showed a closer genetic relationship of Blond strains than of Black strains (Fst: 0.01 vs. 0.05, respectively). Genetic correlation estimates for milk yield were 0.70 in both cases. Based on Fst distances, we expected a lower correlation for Black strains than for Blond ones if dominance or epistasis are important. Thus, we attribute the value of this correlation not being close to 1 mostly to genotype-by-environment interaction, including on-farm management and trait modeling. Regardless, the correlation of 0.7 across populations is encouraging for future joint work of Latxa and Manech breeders, including joint genetic evaluations.

Removing data and using metafounders alleviates biases for all traits in Lacaune dairy sheep predictions.

Bias in dairy genetic evaluations, when it exists, has to be understood and properly addressed. The origin of biases is not always clear. We analyzed 40 yr of records from the Lacaune dairy sheep breeding program to evaluate the extent of bias, assess possible corrections, and emit hypotheses on its origin. The data set included 7 traits (milk yield, fat and protein contents, somatic cell score, teat angle, udder cleft, and udder depth) with records from 600,000 to 5 million depending on the trait, ∼1,900,000 animals, and ∼5,900 genotyped elite artificial insemination rams. For the ∼8% animals with missing sire, we fit 25 unknown parent groups. We used the linear regression method to compare "partial" and "whole" predictions of young rams before and after progeny testing, with 7 cut-off points, and we obtained estimates of their bias, (over)dispersion, and accuracy in early proofs. We tried (1) several scenarios as follows: multiple or single trait, the "official" (routine) evaluation, which is a mixture of both single and multiple trait, and "deletion" of data before 1990; and (2) several models as follows: BLUP and single-step genomic (SSG)BLUP with fixed unknown parent groups or metafounders, where, for metafounders, their relationship matrix gamma was estimated using either a model for inbreeding trend, or base allele frequencies estimated by peeling. The estimate of gamma obtained by modeling the inbreeding trend resulted in an estimated increase of inbreeding, based on markers, faster than the pedigree-based one. The estimated genetic trends were similar for most models and scenarios across all traits, but were shrunken when gamma was estimated by peeling. This was due to shrinking of the estimates of metafounders in the latter case. Across scenarios, all traits showed bias, generally as an overestimate of genetic trend for milk yield and an underestimate for the other traits. As for the slope, it showed overdispersion of estimated breeding values for all traits. Using multiple-trait models slightly reduced the overestimate of genetic trend and the overdispersion, as did including genomic information (i.e., SSGBLUP) when the gamma matrix was estimated by the model for inbreeding trend. However, only deletion of historical data before 1990 resulted in elimination of both kind of biases. The SSGBLUP resulted in more accurate early proofs than BLUP for all traits. We considered that a snowball effect of small errors in each genetic evaluation, combined with selection, may have resulted in biased evaluations. Improving statistical methods reduced some bias but not all, and a simple solution for this data set was to remove historical records.

Genomic and pedigree estimation of inbreeding depression for semen traits in the Basco-Béarnaise dairy sheep breed.

Inbreeding depression is associated with a decrease in performance and fitness of the animals. The goal of this study was to evaluate pedigree-based and genomic methods to estimate the level of inbreeding and inbreeding depression for 3 semen traits (volume, concentration, and motility score) in the Basco-Béarnaise sheep breed. Data comprised 16,196 (or 15,071) phenotypic records from 620 rams (of which 533 rams had genotypes of 36,464 SNPs). The pedigree included 8,266 animals, composed of the 620 rams and their ancestors. The number of equivalent complete generations for the 620 rams was 7.04. Inbreeding coefficients were estimated using genomic and pedigree-based information. Genomic inbreeding coefficients were estimated from individual SNP and using segments of homozygous SNP (runs of homozygosity, ROH). Short ROH are of old origin, whereas long ROH are due to recent inbreeding. Considering that the equivalent number of generations in Basco-Béarnaise was 6, inbreeding coefficients for ROH with a length >4 Mb refer to all (recent + old) inbreeding, those with a length >17 Mb correspond to recent inbreeding, and the difference between them indicates old inbreeding. Pedigree-based inbreeding coefficients were also estimated classically, or accounting for nonzero relationships for unknown parents, or including metafounder relationships (estimated using markers) to account for missing pedigree information. Finally, inbreeding coefficients combining genotyped and nongenotyped animal information were computed from matrix H of the single-step approach, also including metafounders. Inbreeding depression was estimated differently depending on the approach used to compute inbreeding coefficients. These 8 estimators of inbreeding coefficients were included as covariates in different animal models. No inbreeding depression was detected for sperm volume or sperm concentration. Inbreeding depression was significant for the motility of spermatozoa. The effect of old and recent inbreeding on motility was null and negative, respectively, demonstrating the existence of purging by selection of deleterious recessive alleles affecting motility. A 10% increase in inbreeding would result in a reduction in mean motility ranging between 0.09 and 0.22 points in the score (from 0 to 5). Motility is unfavorably affected by increasing recent inbreeding but the impact is very small. Runs of homozygosity and metafounders allow us to accurately estimate inbreeding depression and detect recent inbreeding.

Resistance to gastrointestinal nematodes in dairy sheep: Genetic variability and relevance of artificial infection of nucleus rams to select for resistant ewes on farms.

Breeding sheep for enhanced resistance to gastrointestinal parasites is a promising strategy to limit the use of anthelmintics due to the now widespread resistance of parasites to these molecules. This paper reports the genetic parameters estimated for parasite resistance and resilience traits in the Blond-faced Manech dairy sheep breed and the putative impacts of the selection for resistance to gastrointestinal nematodes (GIN) on farms. Two datasets were used. First, the rams of the selection scheme were artificially infected twice with L3 Haemonchus contortus larvae. Faecal egg counts (FEC) and packed cell volume (PCV) loss were measured 30 days after each infection. Secondly, the FEC, PCV and body condition score (BCS) (1-6 measures per ewe) of naturally infected ewes on farms were measured in the spring, summer and autumn over a two-year period. Genetic parameters were estimated for each dataset independently but also globally based on the pedigree connections between the two datasets. For the experimentally infected sires, the FEC following the second infection was moderately heritable (heritability: 0.35) and strongly correlated with FEC after the first infection (genetic correlation: 0.92). For the naturally infected ewes, FEC was also heritable (0.18). Using the two datasets together, a genetic correlation of 0.56-0.71 was estimated between the FEC values of the experimentally infected rams and naturally infected ewes. Consequently, the genetic variability of parasite resistance is similar whatever the physiological status (males or milking/pregnant ewes) and the infection conditions (experimental infection with one parasite or natural infection with several parasites). In practice, when the sire population is divided into two groups based on their genetic value, the FEC of the ewes born to the 50% most resistant sires is half that of the ewes born to the 50% most susceptible sires. Our study shows the feasibility and efficiency of genetic selection for parasitism resistance based on the sires' FEC records to improve parasite resistance in naturally grazing ewes. For breed improvement, and to increase the selection pressure on parasite resistance, it seems more appropriate to measure FEC values on rams after experimental infection rather than on ewes in natural infection conditions because this limits the number and standardizes the conditions of FEC measurements.

Optimal mating strategies to manage a heterozygous advantage major gene in sheep.

Some mutations (or 'major genes') have a desirable effect in heterozygous carriers but an undesirable effect in homozygous carriers. When these mutations affect a trait of significant economic importance, their eradication, depending on their effect and frequency, may be counterproductive. This is especially the case of major genes affecting the ovulation rate and thus the prolificacy in meat sheep populations. To manage such situations, a mating design based on the major genotypes of reproducers has to be optimized. Both the effect of the major gene and the cost of genotyping candidates at this locus influence the expected genetic progress and profitability of the breeding plan. The aim of this study was to determine the optimal combination of matings that maximizes profitability at the level of the whole population (nucleus + commercial flocks). A deterministic model was developed and, using sequential quadratic programming methodology, the optimal strategy (optimal combination of matings) that maximized the economic gain achieved by the population across a range of genotype effects and genotyping costs was determined. The optimal strategy was compared with simpler and more practical strategies based on a limited number of parental genotype mating types. Depending on the genotype effect and genotyping costs, the optimal strategy varied, such that either the heterozygous frequency and/or polygenic gain was maximized with a large number of animals genotyped, or when genotyping costs were higher, the optimization led to lower heterozygous frequency and/or polygenic gain with fewer animals genotyped. Comparisons showed that some simpler strategies were close to the optimal strategy. An overlapping model was then derived as an application of the real case of the French Lacaune meat sheep OVI-TEST breeding program. Results showed that a practical strategy based on mating non-carriers to heterozygous carriers was only slightly less effective than the optimal strategy, with a reduction in efficiency from 3% to 8%, depending on the genotyping costs. Based on only two different parental genotype mating types, this strategy would be easy to implement.

Genetic and economic effects of the increase in female paternal filiations by parentage assignment in sheep and goat breeding programs.

In sheep and goat breeding programs, the proportion of females for which the sire is known (known paternity rate [KPR]) can be very low. In this context, paternity assignment using SNP is an attractive tool. The annual genetic gain (AGG) is impacted by the accuracy of the EBV. In populations with a low KPR, the number of known relatives for a given individual is low, and the EBV that are based on this information are imprecise. However, the impact of partially known paternal filiations, in terms of potential genetic and economic losses, has never been quantitatively evaluated in situations where natural mating is the main reproductive mode. A deterministic model was developed to assess, for a panel of real breeding programs, the influence of the female KPR on the AGG and economic benefit. First, males were divided into categories according to their status (natural mating or AI sire) and breeding cycle and females according to parity, sire status (including unknown sire), and breeding cycle of the sire. Second, a demographic model described, for each category, the accumulation of known records for individuals and their close relatives. The output from this model was used to compute the average accuracy of the EBV per category. Then, a genetic model based on the gene flow between categories over time was described. Using the average accuracy of EBV per category, it provided the asymptotic AGG of the nucleus given its KPR. In the economic studies, changes to the mean genetic values in the nucleus and the commercial population after an increase in KPR and various gain:cost ratios (monetary gain due to an extra genetic SD of the selected trait divided by the cost of 1 assignment) were considered. Relative profit and payback periods were computed. We showed that SNP-based parentage assignment aimed at increasing the female KPR was not always profitable and that the type of breeding program and the size of the commercial population should be taken into consideration. Notably, achieving a profit was largely dependent on obtaining a favorable gain:cost ratio. The maximum supplementary AGG (16.9%) was obtained for breeding programs using only natural mating. In such programs without AI, a gain:cost ratio of 5 was needed to make assignment profitable at the nucleus level whereas a gain:cost ratio of 2 was sufficient if the nucleus represented a third of the total population.

Economic evaluation of genomic selection in small ruminants: a sheep meat breeding program.

Recent genomic evaluation studies using real data and predicting genetic gain by modeling breeding programs have reported moderate expected benefits from the replacement of classic selection schemes by genomic selection (GS) in small ruminants. The objectives of this study were to compare the cost, monetary genetic gain and economic efficiency of classic selection and GS schemes in the meat sheep industry. Deterministic methods were used to model selection based on multi-trait indices from a sheep meat breeding program. Decisional variables related to male selection candidates and progeny testing were optimized to maximize the annual monetary genetic gain (AMGG), that is, a weighted sum of meat and maternal traits annual genetic gains. For GS, a reference population of 2000 individuals was assumed and genomic information was available for evaluation of male candidates only. In the classic selection scheme, males breeding values were estimated from own and offspring phenotypes. In GS, different scenarios were considered, differing by the information used to select males (genomic only, genomic+own performance, genomic+offspring phenotypes). The results showed that all GS scenarios were associated with higher total variable costs than classic selection (if the cost of genotyping was 123 euros/animal). In terms of AMGG and economic returns, GS scenarios were found to be superior to classic selection only if genomic information was combined with their own meat phenotypes (GS-Pheno) or with their progeny test information. The predicted economic efficiency, defined as returns (proportional to number of expressions of AMGG in the nucleus and commercial flocks) minus total variable costs, showed that the best GS scenario (GS-Pheno) was up to 15% more efficient than classic selection. For all selection scenarios, optimization increased the overall AMGG, returns and economic efficiency. As a conclusion, our study shows that some forms of GS strategies are more advantageous than classic selection, provided that GS is already initiated (i.e. the initial reference population is available). Optimizing decisional variables of the classic selection scheme could be of greater benefit than including genomic information in optimized designs.

Within- and across-breed genomic predictions and genomic relationships for Western Pyrenees dairy sheep breeds Latxa, Manech, and Basco-Béarnaise.

Genotypes, phenotypes and pedigrees of 6 breeds of dairy sheep (including subdivisions of Latxa, Manech, and Basco-Béarnaise) from the Spain and France Western Pyrenees were used to estimate genetic relationships across breeds (together with genotypes from the Lacaune dairy sheep) and to verify by forward cross-validation single-breed or multiple-breed genetic evaluations. The number of rams genotyped fluctuated between 100 and 1,300 but generally represented the 10 last cohorts of progeny-tested rams within each breed. Genetic relationships were assessed by principal components analysis of the genomic relationship matrices and also by the conservation of linkage disequilibrium patterns at given physical distances in the genome. Genomic and pedigree-based evaluations used daughter yield performances of all rams, although some of them were not genotyped. A pseudo-single step method was used in this case for genomic predictions. Results showed a clear structure in blond and black breeds for Manech and Latxa, reflecting historical exchanges, and isolation of Basco-Béarnaise and Lacaune. Relatedness between any 2 breeds was, however, lower than expected. Single-breed genomic predictions had accuracies comparable with other breeds of dairy sheep or small breeds of dairy cattle. They were more accurate than pedigree predictions for 5 out of 6 breeds, with absolute increases in accuracy ranging from 0.05 to 0.30 points. They were significantly better, as assessed by bootstrapping of candidates, for 2 of the breeds. Predictions using multiple populations only marginally increased the accuracy for a couple of breeds. Pooling populations does not increase the accuracy of genomic evaluations in dairy sheep; however, single-breed genomic predictions are more accurate, even for small breeds, and make the consideration of genomic schemes in dairy sheep interesting.

Prediction of fatty acid profiles in cow, ewe, and goat milk by mid-infrared spectrometry.

Mid-infrared (MIR) spectrometry was used to estimate the fatty acid (FA) composition in cow, ewe, and goat milk. The objectives were to compare different statistical approaches with wavelength selection to predict the milk FA composition from MIR spectra, and to develop equations for FA in cow, goat, and ewe milk. In total, a set of 349 cow milk samples, 200 ewe milk samples, and 332 goat milk samples were both analyzed by MIR and by gas chromatography, the reference method. A broad FA variability was ensured by using milk from different breeds and feeding systems. The methods studied were partial least squares regression (PLS), first-derivative pretreatment + PLS, genetic algorithm + PLS, wavelets + PLS, least absolute shrinkage and selection operator method (LASSO), and elastic net. The best results were obtained with PLS, genetic algorithm + PLS and first derivative + PLS. The residual standard deviation and the coefficient of determination in external validation were used to characterize the equations and to retain the best for each FA in each species. In all cases, the predictions were of better quality for FA found at medium to high concentrations (i.e., for saturated FA and some monounsaturated FA with a coefficient of determination in external validation >0.90). The conversion of the FA expressed in grams per 100mL of milk to grams per 100g of FA was possible with a small loss of accuracy for some FA.

Assessment of accuracy of genomic prediction for French Lacaune dairy sheep.

Genomic selection in Lacaune dairy sheep was investigated based on genotypes from the OvineSNP50 BeadChip (Illumina Inc., San Diego, CA). Historical artificial insemination progeny-tested rams formed a population of 2,892 genotyped rams. Additional ungenotyped rams and females were included by single-step genomic BLUP (ssGBLUP). Three prediction strategies were tried: pseudo-BLUP (using all rams and daughter yield deviations), pseudo-ssGBLUP (using all rams and daughter yield deviations), and regular ssGBLUP (using all phenotypes and pedigree in an animal model). The population linkage disequilibrium was determined, with an average squared correlation coefficient of 0.11 for markers closer than 0.1cM (lower than in dairy cattle). The estimated effective population is 370 individuals. Gain in accuracy of genomic selection over parent averages ranged from 0.10 to 0.20. Highest accuracies and lowest bias were found using regular ssGBLUP. Transition to a genomic breeding scheme is possible but costs need to be carefully evaluated.

Potential benefits of genomic selection on genetic gain of small ruminant breeding programs.

In conventional small ruminant breeding programs, only pedigree and phenotype records are used to make selection decisions but prospects of including genomic information are now under consideration. The objective of this study was to assess the potential benefits of genomic selection on the genetic gain in French sheep and goat breeding designs of today. Traditional and genomic scenarios were modeled with deterministic methods for 3 breeding programs. The models included decisional variables related to male selection candidates, progeny testing capacity, and economic weights that were optimized to maximize annual genetic gain (AGG) of i) a meat sheep breeding program that improved a meat trait of heritability (h(2)) = 0.30 and a maternal trait of h(2) = 0.09 and ii) dairy sheep and goat breeding programs that improved a milk trait of h(2) = 0.30. Values of ±0.20 of genetic correlation between meat and maternal traits were considered to study their effects on AGG. The Bulmer effect was accounted for and the results presented here are the averages of AGG after 10 generations of selection. Results showed that current traditional breeding programs provide an AGG of 0.095 genetic standard deviation (σa) for meat and 0.061 σa for maternal trait in meat breed and 0.147 σa and 0.120 σa in sheep and goat dairy breeds, respectively. By optimizing decisional variables, the AGG with traditional selection methods increased to 0.139 σa for meat and 0.096 σa for maternal traits in meat breeding programs and to 0.174 σa and 0.183 σa in dairy sheep and goat breeding programs, respectively. With a medium-sized reference population (nref) of 2,000 individuals, the best genomic scenarios gave an AGG that was 17.9% greater than with traditional selection methods with optimized values of decisional variables for combined meat and maternal traits in meat sheep, 51.7% in dairy sheep, and 26.2% in dairy goats. The superiority of genomic schemes increased with the size of the reference population and genomic selection gave the best results when nref > 1,000 individuals for dairy breeds and nref > 2,000 individuals for meat breed. Genetic correlation between meat and maternal traits had a large impact on the genetic gain of both traits. Changes in AGG due to correlation were greatest for low heritable maternal traits. As a general rule, AGG was increased both by optimizing selection designs and including genomic information.

Genomic selection in the French Lacaune dairy sheep breed.

Genomic selection aims to increase accuracy and to decrease generation intervals, thus increasing genetic gains in animal breeding. Using real data of the French Lacaune dairy sheep breed, the purpose of this study was to compare the observed accuracies of genomic estimated breeding values using different models (infinitesimal only, markers only, and joint estimation of infinitesimal and marker effects) and methods [BLUP, Bayes Cπ, partial least squares (PLS), and sparse PLS]. The training data set included results of progeny tests of 1,886 rams born from 1998 to 2006, whereas the validation set had results of 681 rams born in 2007 and 2008. The 3 lactation traits studied (milk yield, fat content, and somatic cell scores) had heritabilities varying from 0.14 to 0.41. The inclusion of molecular information, as compared with traditional schemes, increased accuracies of estimated breeding values of young males at birth from 18 up to 25%, according to the trait. Accuracies of genomic methods varied from 0.4 to 0.6, according to the traits, with minor differences among genomic approaches. In Bayes Cπ, the joint estimation of marker and infinitesimal effects had a slightly favorable effect on the accuracies of genomic estimated breeding values, and were especially beneficial for somatic cell counts, the less heritable trait. Inclusion of infinitesimal effects also improved slopes of predictive regression equations. Methods that select markers implicitly (Bayes Cπ and sparse PLS) were advantageous for some models and traits, and are of interest for further quantitative trait loci studies.

Genetic correlation between female fertility and milk yield in Lacaune sheep.

A total of 416,670 lactations for 189,101 ewes from 3,603 sires and distributed across 1,978 herd-year groups were used to estimate genetic and environmental parameters of standardized milk yield (SMY(T)), fertility in ewe lambs (PR(1)), and fertility in adult ewes (PR(A)). Parameters were estimated with a multiple-trait sire linear model. Heritabilities for SMY(T), PR(1), and PR(A) were 0.27 (0.009), 0.04 (0.004), and 0.05 (0.004), respectively. These results were in accordance with the literature. The genetic correlation between PR(1) and PR(A) was 0.55, indicating that fertility is not the same trait in ewe lambs and adult ewes. The genetic correlation between milk yield and lamb fertility was not significantly different from zero. The genetic correlation between milk yield and fertility in adult ewe (-0.23) was in the range of antagonistic correlations reported in dairy cattle. Consequently, these results show that selection for milk yield can induce an indirect decrease in fertility. Nevertheless, no phenotypic decrease in fertility in artificial insemination matings has been observed in this population. This is the first time that correlation between milk yield and fertility is reported in sheep and further investigations are needed to confirm this result.

Genetic parameters for milk somatic cell scores and relationships with production traits in French Lacaune dairy sheep.

Records from 94,445 and 45,499 French Lacaune dairy ewes in first and second lactations, respectively, were used to estimate genetic parameters for somatic cell scores. Somatic cell count data came from an extensive recording scheme and sample testing that began in 1999 using the flocks enrolled in the official milk recording system. Somatic cell count data were from 2 to 4 test days per lactation. Lactation average and single test-day somatic cell scores were considered in multitrait sire models. The heritability estimate of lactation somatic cell score was close to 0.13 and similar for first and second parity. Heritabilities of somatic cell scores increased from first to fourth test day (from 0.07 to 0.11 in first lactation and from 0.05 to 0.13 in second lactation). Genetic correlations between somatic cell scores were high, usually more than 0.91, but lower between first test day and later test days in first lactation (0.64 to 0.88). The genetic correlations between lactation somatic cell score and milk yield, between lactation somatic cell score and fat content, and between lactation somatic cell score and protein content were 0.18, 0.04, and 0.03 in first lactation, respectively. The genetic antagonism between test day somatic cell score and milk yield measured in first lactation increased from beginning to the end of the lactation (0.05 to 0.23). This antagonism was slightly lower for somatic cell score in second lactation (from 0.09 to 0.14, and 0.08 for lactation mean). Environmental correlations in first lactation between lactation somatic cell score and milk yield, between lactation somatic cell score and fat content, and between lactation somatic cell score and protein content were -0.18, 0.13, and 0.30, respectively.

Genetic analysis for mastitis resistance and milk somatic cell score in French Lacaune dairy sheep.

Genetic analysis for mastitis resistance was studied from two data sets. Firstly, risk factors for different mastitis traits, i.e. culling due to clinical or chronic mastitis and subclinical mastitis predicted from somatic cell count (SCC), were explored using data from 957 first lactation Lacaune ewes of an experimental INRA flock composed of two divergent lines for milk yield. Secondly, genetic parameters for SCC were estimated from 5 272 first lactation Lacaune ewes recorded among 38 flocks, using an animal model. In the experimental flock, the frequency of culling due to clinical mastitis (5% ) was lower than that of subclinical mastitis (10% ) predicted from SCC. Predicted subclinical mastitis was unfavourably associated with the milk yield level. Such an antagonism was not detected for clinical mastitis, which could result, to some extent, from its low frequency or from the limited amount of data. In practice, however, selection for mastitis resistance could be limited in a first approach to selection against subclinical mastitis using SCC. The heritability estimate of SCC was 0.15 for the lactation mean trait and varied from 0.04 to 0.12 from the first to the fifth test-day. The genetic correlation between lactation SCC and milk yield was slightly positive (0.15) but showed a strong evolution during lactation, i.e. from favourable (-0.48) to antagonistic (0.27). On a lactation basis, our results suggest that selection for mastitis resistance based on SCC is feasible. Patterns for genetic parameters within first lactation, however, require further confirmation and investigation.