Genome wide association analysis on semen volume and milk yield using different strategies of imputation to whole genome sequence in French dairy goats.

BACKGROUND: Goats were domesticated 10,500 years ago to supply humans with useful resources. Since then, specialized breeds that are adapted to their local environment have been developed and display specific genetic profiles. The VarGoats project is a 1000 genomes resequencing program designed to cover the genetic diversity of the Capra genus. In this study, our main objective was to assess the use of sequence data to detect genomic regions associated with traits of interest in French Alpine and Saanen breeds. RESULTS: Direct imputation from the GoatSNP50 BeadChip genotypes to sequence level was investigated in these breeds using FImpute and different reference panels: within-breed, all Capra hircus sequenced individuals, European goats and French mainland goats. The best results were obtained with the French goat panel with allele and genotype concordance rates reaching 0.86 and 0.75 in the Alpine and 0.86 and 0.73 in the Saanen breed respectively. Mean correlations tended to be low in both breeds due to the high proportion of variants with low frequencies. For association analysis, imputation was performed using FImpute for 1129 French Alpine and Saanen males using within-breed and French panels on 23,338,436 filtered variants. The association results of both imputation scenarios were then compared. In Saanen goats, a large region on chromosome 19 was significantly linked to semen volume and milk yield in both scenarios. Significant variants for milk yield were annotated for 91 genes on chromosome 19 in Saanen goats. For semen volume, the annotated genes include YBOX2 which is related to azoospermia or oligospermia in other species. New signals for milk yield were detected on chromosome 2 in Alpine goats and on chromosome 5 in Saanen goats when using a multi-breed panel. CONCLUSION: Even with very small reference populations, an acceptable imputation quality can be achieved in French dairy goats. GWAS on imputed sequences confirmed the existence of QTLs and identified new regions of interest in dairy goats. Adding identified candidates to a genotyping array and sequencing more individuals might corroborate the involvement of identified regions while removing potential imputation errors.

Alternative methods improve the accuracy of genomic prediction using information from a causal point mutation in a dairy sheep model.

BACKGROUND: Genomic evaluation is usually based on a set of markers assumed to be linked with causal mutations. Selection and precise management of major genes and the remaining polygenic component might be improved by including causal polymorphisms in the evaluation models. In this study, various methods involving a known mutation were used to estimate prediction accuracy. The SOCS2 gene, which influences body growth, milk production and somatic cell scores, a proxy for mastitis, was studied as an example in dairy sheep. METHODS: The data comprised 1,503,148 phenotypes and 9844 54K SNPs genotypes. The SOCS2 SNP was genotyped for 4297 animals and imputed in the above 9844 animals. Breeding values and their accuracies were estimated for each of nine traits by using single-step approaches. Pedigree-based BLUP, single-step genomic BLUP (ssGBLUP) involving the 54K ovine SNPs chip, and four weighted ssGBLUP (WssGBLUP) methods were compared. In WssGBLUP methods, weights are assigned to SNPs depending on their effect on the trait. The ssGBLUP and WssGBLUP methods were again tested after including the SOCS2 causal mutation as a SNP. Finally, the Gene Content approach was tested, which uses a multiple-trait model that considers the SOCS2 genotype as a trait. RESULTS: EBV accuracies were increased by 14.03% between the pedigree-based BLUP and ssGBLUP methods and by 3.99% between ssGBLUP and WssGBLUP. Adding the SOCS2 SNP to ssGBLUP methods led to an average gain of 0.26%. Construction of the kinship matrix and estimation of breeding values was generally improved by placing emphasis on SNPs in regions with a strong effect on traits. In the absence of chip data, the Gene Content method, compared to pedigree-based BLUP, efficiently accounted for partial genotyping information on SOCS2 as accuracy was increased by 6.25%. This method also allowed dissociation of the genetic component due to the major gene from the remaining polygenic component. CONCLUSIONS: Causal mutations with a moderate to strong effect can be captured with conventional SNP chips by applying appropriate genomic evaluation methods. The Gene Content method provides an efficient way to account for causal mutations in populations lacking genome-wide genotyping.

A validation study of loci associated with mastitis resistance in two French dairy sheep breeds.

BACKGROUND: The identification of loci associated with resistance to mastitis or of the causative mutations may be helpful in breeding programs for dairy sheep as it is for cattle worldwide. Seven genomic regions that control milk somatic cell counts, an indirect indicator of udder infection, have already been identified in sheep (Spanish Churra, French Lacaune and Italian Sardinian-Lacaune backcross populations). In this study, we used a 960 custom-designed ovine single nucleotide polymorphism (SNP) chip in Lacaune and Manech Tête Rousse dairy sheep to validate these seven genomic regions associated with mastitis. RESULTS: The most significant SNP (rs868996547) on Ovis aries chromosome (OAR) 3 was a previously described mutation in the suppressor of cytokine signalling 2 (SOCS2) gene. An antagonist effect of this causal candidate between health and growth in Lacaune sheep was confirmed. Effects of the mutation on the infectious status of the udder, i.e. increases in milk somatic cell counts and bacteria shedding, were also identified. This SNP was not present in the data available on Manech Tête Rousse. Three other regions associated with mastitis were also confirmed on OAR16 (Manech Tête Rousse), 19 (Lacaune) and 2 (both breeds). For the OAR2 region, we validated previously detected SNPs in several other breeds (Sarda, Churra, and Chios). For significant SNPs in the four mastitis regions, the effect varied from 0.24 to 0.67 phenotypic standard deviation of the traits. Two of the mastitis quantitative trait loci (QTL) regions (OAR2 and 16) that we validated here were also associated in opposite ways with milk production traits in both populations. CONCLUSIONS: These results indicate, at least in part, a genomic basis for the trade-off between milk production and mastitis resistance. Four of the seven mastitis QTL regions that were previously identified in independent populations, were confirmed in this study, which demonstrates partial sharing of mastitis-related genetic mechanisms between different distant dairy sheep populations.

Signatures of selection and environmental adaptation across the goat genome post-domestication.

BACKGROUND: Since goat was domesticated 10,000 years ago, many factors have contributed to the differentiation of goat breeds and these are classified mainly into two types: (i) adaptation to different breeding systems and/or purposes and (ii) adaptation to different environments. As a result, approximately 600 goat breeds have developed worldwide; they differ considerably from one another in terms of phenotypic characteristics and are adapted to a wide range of climatic conditions. In this work, we analyzed the AdaptMap goat dataset, which is composed of data from more than 3000 animals collected worldwide and genotyped with the CaprineSNP50 BeadChip. These animals were partitioned into groups based on geographical area, production uses, available records on solid coat color and environmental variables including the sampling geographical coordinates, to investigate the role of natural and/or artificial selection in shaping the genome of goat breeds. RESULTS: Several signatures of selection on different chromosomal regions were detected across the different breeds, sub-geographical clusters, phenotypic and climatic groups. These regions contain genes that are involved in important biological processes, such as milk-, meat- or fiber-related production, coat color, glucose pathway, oxidative stress response, size, and circadian clock differences. Our results confirm previous findings in other species on adaptation to extreme environments and human purposes and provide new genes that could explain some of the differences between goat breeds according to their geographical distribution and adaptation to different environments. CONCLUSIONS: These analyses of signatures of selection provide a comprehensive first picture of the global domestication process and adaptation of goat breeds and highlight possible genes that may have contributed to the differentiation of this species worldwide.