Benoit Auvray

Genomic prediction and genome-wide association study for dagginess and host internal parasite resistance in New Zealand sheep

BackgroundDagginess (faecal soiling of the perineum region) and host nematode parasite resistance are important animal welfare traits in New Zealand sheep. Genomic prediction (GP) estimates the genetic merit, as a molecular breeding value (mBV), for each trait based on many SNPs. The additional information the mBV provides (as determined by its accuracy) has led to its incorporation into breeding schemes. Some GP methods give SNP effects, which provide additional information to identify genome-wide associations (GWAS) for a trait of interest. Here we report results from a GP and GWAS study for dagginess and host nematode parasite resistance in a New Zealand sheep industry resource.ResultsGenomic prediction analysis was performed using 50K SNP chip data and parent average-removed, de-regressed BVs for five traits, from a resource of 8705 pedigree recorded animals. The five traits were dag score at three and eight months (DAG3, DAG8) and nematode faecal egg count in summer (FEC1), autumn (FEC2) and as an adult (AFEC). The resource consisted of Romney, Coopworth, Perendale, Texel and various breed crosses (designated: CompRCP, CompRCPT and CompCRP). The pure breeds, apart from Texel, plus CompRCP were used to develop the GP. The resulting SNP effects were used to identify genetic regions associated with dagginess and parasite resistance. Accuracies of the weighted correlation between mBV and true BV ranged between −0.07 (Texel) and 0.56 (Coopworth) for DAG3 and DAG8. For FEC1, FEC2 and AFEC accuracies ranged between −0.22 (CompRCPT) and 0.69 (Coopworth). The weighted average individual accuracy (calculated from theory) ranges were 0.13 (Texel) to 0.52 (Coopworth) and 0.11 (Texel) to 0.55 (Coopworth) respectively, for dagginess and parasite traits. There was one SNP for DAG8 that reached Bonferroni significance threshold (P < 1 × 10−6) on OAR15, the same two SNPs for each of the parasite traits (OAR26) and none for DAG3. A notable peak was also observed on OAR7 for all the parasite traits, however, it did not reach the Bonferroni significance threshold.ConclusionsThis study presents the first results of a GWAS on dagginess and faecal egg count traits in New Zealand sheep. The results suggest that there are quantitative trait loci on OAR 15 for dagginess and on OAR26 and seven for faecal egg count.

Estimation of linkage disequilibrium and effective population size in New Zealand sheep using three different methods to create genetic maps

BackgroundInvestments in genetic selection have played a major role in the New Zealand sheep industry competitiveness. Selection may erode genetic diversity, which is a crucial factor for the success of breeding programs. Better understanding of linkage disequilibrium (LD) and ancestral effective population size (Ne) through quantifying this diversity and comparison between populations allows for more informed decisions with regards to selective breeding taking population genetic diversity into account. The estimation of Ne can be determined via genetic markers and requires knowledge of genetic distances between these markers. Single nucleotide polymorphisms (SNP) data from a sample of 12,597 New Zealand crossbred and purebred sheep genotyped with the Illumina Ovine SNP50 BeadChip was used to perform a genome-wide scan of LD and Ne. Three methods to estimate genetic distances were investigated: 1) M1: a ratio fixed across the whole genome of one Megabase per centiMorgan; 2) M2: the ratios of genetic distance (using M3, below) over physical distance fixed for each chromosome; and, 3) M3: a genetic map of inter-SNP distances estimated using CRIMAP software (v2.503).ResultsThe estimates obtained with M2 and M3 showed much less variability between autosomes than those with M1, which tended to give lower Ne results and higher LD decay. The results suggest that Ne has decreased since the development of sheep breeds in Europe and this reduction in Ne has been accelerated in the last three decades. The Ne estimated for five generations in the past ranged from 71 to 237 for Texel and Romney breeds, respectively. A low level of genetic kinship and inbreeding was estimated in those breeds suggesting avoidance of mating close relatives.ConclusionsM3 was considered the most accurate method to create genetic maps for the estimation of LD and Ne. The findings of this study highlight the history of genetic selection in New Zealand crossbred and purebred sheep and these results will be very useful to understand genetic diversity of the population with respect to genetic selection. In addition, it will help geneticists to identify genomic regions which have been preferentially selected within a variety of breeds and populations.

Derivation of a new lamb survival trait for the New Zealand sheep industry

Previous research identified that a review of the current industry New Zealand lamb survival trait was necessary as its recording accuracy was reliant on farmers notifying their Sheep Improvement Limited bureau of lamb deaths. This paper reports the decision rules and genetic parameters for a new lamb survival trait for the New Zealand sheep industry. These rules define the new lamb survival trait (NEWSUR) using lamb birth fate (BFATE) codes and the presence/absence of lamb weight measurements. Six univariate animal models were tested and used to estimate variance or covariance components and the resulting direct and maternal heritabilities for NEWSUR. The models differed in the way they adjust for the effect of day of birth, the exclusion or inclusion of a litter (dam/year of birth) random effect, and the application or not of a logit transformation of the phenotypes. For both the linear and logistic methods, models including the random effect of litter provided the best fit for NEWSUR according to log-likelihood values. Log-likelihoods for the linear and logistic models cannot be compared; therefore, a cross-validation method was used to assess whether the logit transformation was appropriate by analyzing the predictive ability of the models. The mean square errors were slightly lower for the linear compared with the logistic model, and therefore, the linear model was recommended for industry use. The heritability attributed to direct effects ranged from 2 to 5.5%. A direct heritability of 5.5% resulted from a linear model without litter effect and omitting the effect of day of birth on survival, whereas a direct heritability of 2% resulted from the logistic model fitting a random litter effect. The heritability attributed to maternal genetic effects ranged from 1.9 to 7.7%. A maternal genetic heritability of 7.7% resulted from the logistic model omitting the litter effect, whereas a maternal genetic heritability of 1.9% resulted from the linear model fitting a random litter effect. The addition of the litter random effect substantially decreased the maternal heritabilities in all cases and was recommended for industry use to avoid overestimation of the maternal genetic variance. Sheep Improvement Limited has implemented NEWSUR and the associated genetic evaluation model based on information described in this paper. Industry-wide implementation will enable sheep breeders to produce more accurate genetic evaluations to their commercial clients.