The success of Genome Wide Association Studies in the discovery of sequence variation linked to complex traits in humans has increased interest in high throughput SNP genotyping assays in livestock species. Primary goals are QTL detection and genomic selection. The purpose here was design of a 50–60,000 SNP chip for goats. The success of a moderate density SNP assay depends on reliable bioinformatic SNP detection procedures, the technological success rate of the SNP design, even spacing of SNPs on the genome and selection of Minor Allele Frequencies (MAF) suitable to use in diverse breeds. Through the federation of three SNP discovery projects consolidated as the International Goat Genome Consortium, we have identified approximately twelve million high quality SNP variants in the goat genome stored in a database together with their biological and technical characteristics. These SNPs were identified within and between six breeds (meat, milk and mixed): Alpine, Boer, Creole, Katjang, Saanen and Savanna, comprising a total of 97 animals. Whole genome and Reduced Representation Library sequences were aligned on >10 kb scaffolds of the de novo goat genome assembly. The 60,000 selected SNPs, evenly spaced on the goat genome, were submitted for oligo manufacturing (Illumina, Inc) and published in dbSNP along with flanking sequences and map position on goat assemblies (i.e. scaffolds and pseudo-chromosomes), sheep genome V2 and cattle UMD3.1 assembly. Ten breeds were then used to validate the SNP content and 52,295 loci could be successfully genotyped and used to generate a final cluster file. The combined strategy of using mainly whole genome Next Generation Sequencing and mapping on a contig genome assembly, complemented with Illumina design tools proved to be efficient in producing this GoatSNP50 chip. Advances in use of molecular markers are expected to accelerate goat genomic studies in coming years.

The diversity of populations in domestic species offers great opportunities to study genome response to selection. The recently published Sheep HapMap dataset is a great example of characterization of the world wide genetic diversity in sheep. In this study, we re-analyzed the Sheep HapMap dataset to identify selection signatures in worldwide sheep populations. Compared to previous analyses, we made use of statistical methods that (i) take account of the hierarchical structure of sheep populations, (ii) make use of linkage disequilibrium information and (iii) focus specifically on either recent or older selection signatures. We show that this allows pinpointing several new selection signatures in the sheep genome and distinguishing those related to modern breeding objectives and to earlier post-domestication constraints. The newly identified regions, together with the ones previously identified, reveal the extensive genome response to selection on morphology, color and adaptation to new environments.

Gastrointestinal nematodes are one of the main health issues in sheep breeding. To identify loci affecting the resistance to Haemonchus contortus, a genome scan was carried out using 1,275 Romane × Martinik Black Belly backcross lambs. The entire population was challenged with Haemonchus contortus in 2 consecutive experimental infections, and fecal egg counts (FEC) and packed cell volumes were measured. A subgroup of 332 lambs with extreme FEC was necropsied to determine the total worm burden, length of female worms, sex ratio in the worm population, abomasal pH, and serum and mucosal G immunoglobulins (IgG) responses. Pepsinogen concentration was measured in another subset of 229 lambs. For QTL detection, 160 microsatellite markers were used as well as the Illumina OvineSNP50 BeadChip that provided 42,469 SNP markers after quality control. Linkage, association, and joint linkage and association analyses were performed with the QTLMAP software. Linkage disequilibrium (LD) was estimated within each pure breed, and association analyses were carried out either considering or not the breed origin of the haplotypes. Four QTL regions on sheep chromosomes (OAR)5, 12, 13, and 21 were identified as key players among many other QTL with small to moderate effects. A QTL on OAR21 affecting pepsinogen concentration exactly matched the pepsinogen (PGA5) locus. A 10-Mbp region affecting FEC after the 1st and 2nd infections was found on OAR12. The SNP markers outperformed microsatellites in the linkage analysis. Taking advantage of the LD helped to refine the locations of the QTL mapped on OAR5 and 13.

Recombination is a complex biological process that results from a cascade of multiple events during meiosis. Understanding the genetic determinism of recombination can help to understand if and how these events are interacting. To tackle this question, we studied the patterns of recombination in sheep, using multiple approaches and data sets. We constructed male recombination maps in a dairy breed from the south of France (the Lacaune breed) at a fine scale by combining meiotic recombination rates from a large pedigree genotyped with a 50K SNP array and historical recombination rates from a sample of unrelated individuals genotyped with a 600K SNP array. This analysis revealed recombination patterns in sheep similar to other mammals but also genome regions that have likely been affected by directional and diversifying selection. We estimated the average recombination rate of Lacaune sheep at 1.5 cM/Mb, identified ∼50,000 crossover hotspots on the genome, and found a high correlation between historical and meiotic recombination rate estimates. A genome-wide association study revealed two major loci affecting interindividual variation in recombination rate in Lacaune, including the RNF212 and HEI10 genes and possibly two other loci of smaller effects including the KCNJ15 and FSHR genes. The comparison of these new results to those obtained previously in a distantly related population of domestic sheep (the Soay) revealed that Soay and Lacaune males have a very similar distribution of recombination along the genome. The two data sets were thus combined to create more precise male meiotic recombination maps in Sheep. However, despite their similar recombination maps, Soay and Lacaune males were found to exhibit different heritabilities and QTL effects for interindividual variation in genome-wide recombination rates. This highlights the robustness of recombination patterns to underlying variation in their genetic determinism.

Agroecology uses natural processes and local resources rather than chemical inputs to ensure production while limiting the environmental footprint of livestock and crop production systems. Selecting to achieve a maximization of target production criteria has long proved detrimental to fitness traits. However, since the 1990s, developments in animal breeding have also focussed on animal robustness by balancing production and functional traits within overall breeding goals. We discuss here how an agroecological perspective should further shift breeding goals towards functional traits rather than production traits. Breeding for robustness aims to promote individual adaptive capacities by considering diverse selection criteria which include reproduction, animal health and welfare, and adaptation to rough feed resources, a warm climate or fluctuating environmental conditions. It requires the consideration of genotype × environment interactions in the prediction of breeding values. Animal performance must be evaluated in low-input systems in order to select those animals that are adapted to limiting conditions, including feed and water availability, climate variations and diseases. Finally, we argue that there is no single agroecological animal type, but animals with a variety of profiles that can meet the expectations of agroecology. The standardization of both animals and breeding conditions indeed appears contradictory to the agroecological paradigm that calls for an adaptation of animals to local opportunities and constraints in weakly artificialized systems tied to their physical environment.

A French sheep case control study has been organised to estimate the effects of the PrP haplotypes on resistance to atypical scrapie. The ALHQ and AFRQ haplotypes are significantly more susceptible than the others.

The diversity of populations in domestic species offers great opportunities to study genome response to selection. The recently published Sheep HapMap dataset is a great example of characterization of the world wide genetic diversity in sheep. In this study, we re-analyzed the Sheep HapMap dataset to identify selection signatures in worldwide sheep populations. Compared to previous analyses, we made use of statistical methods that (i) take account of the hierarchical structure of sheep populations, (ii) make use of linkage disequilibrium information and (iii) focus specifically on either recent or older selection signatures. We show that this allows pinpointing several new selection signatures in the sheep genome and distinguishing those related to modern breeding objectives and to earlier post-domestication constraints. The newly identified regions, together with the ones previously identified, reveal the extensive genome response to selection on morphology, color and adaptation to new environments.

Gastrointestinal nematode infections are one of the main health/economic issues in sheep industries, worldwide. Indicator traits for resistance such as faecal egg count (FEC) are commonly used in genomic studies; however, published results are inconsistent among breeds. Meta (or joint)-analysis is a tool for aggregating information from multiple independent studies. The aim of this study was to identify loci underlying variation in FEC, as an indicator of nematode resistance, in a joint analysis using data from three populations (Scottish Blackface, Sarda × Lacaune and Martinik Black-Belly × Romane), genotyped with the ovine 50k SNP chip. The trait analysed was the average animal effect for Strongyles and Nematodirus FEC data. Analyses were performed with regional heritability mapping (RHM), fitting polygenic effects with either the whole genomic relationship matrix or matrices excluding the chromosome being interrogated. Across-population genomic covariances were set to zero. After quality control, 4123 animals and 38 991 SNPs were available for the analysis. RHM identified genome-wide significant regions on OAR4, 12, 14, 19 and 20, with the latter being the most significant. The OAR20 region is close to the major histocompatibility complex, which has often been proposed as a functional candidate for nematode resistance. This region was significant only in the Sarda × Lacaune population. Several other regions, on OAR1, 3, 4, 5, 7, 12, 19, 20 and 24, were significant at the suggestive level.