Geographical contrasts of Y-chromosomal haplogroups from wild and domestic goats reveal ancient migrations and recent introgressions.

By their paternal transmission, Y-chromosomal haplotypes are sensitive markers of population history and male-mediated introgression. Previous studies identified biallelic single-nucleotide variants in the SRY, ZFY and DDX3Y genes, which in domestic goats identified four major Y-chromosomal haplotypes, Y1A, Y1B, Y2A and Y2B, with a marked geographical partitioning. Here, we extracted goat Y-chromosomal variants from whole-genome sequences of 386 domestic goats (75 breeds) and seven wild goat species, which were generated by the VarGoats goat genome project. Phylogenetic analyses indicated domestic haplogroups corresponding to Y1B, Y2A and Y2B, respectively, whereas Y1A is split into Y1AA and Y1AB. All five haplogroups were detected in 26 ancient DNA samples from southeast Europe or Asia. Haplotypes from present-day bezoars are not shared with domestic goats and are attached to deep nodes of the trees and networks. Haplogroup distributions for 186 domestic breeds indicate ancient paternal population bottlenecks and expansions during migrations into northern Europe, eastern and southern Asia, and Africa south of the Sahara. In addition, sharing of haplogroups indicates male-mediated introgressions, most notably an early gene flow from Asian goats into Madagascar and the crossbreeding that in the 19th century resulted in the popular Boer and Anglo-Nubian breeds. More recent introgressions are those from European goats into the native Korean goat population and from Boer goat into Uganda, Kenya, Tanzania, Malawi and Zimbabwe. This study illustrates the power of the Y-chromosomal variants for reconstructing the history of domestic species with a wide geographical range.

Across-country genomic predictions in Norwegian and New Zealand Composite sheep populations with similar development history.

The goal of this study was to assess the feasibility of across-country genomic predictions in Norwegian White Sheep (NWS) and New Zealand Composite (NZC) sheep populations with similar development history. Different training populations were evaluated (i.e., including only NWS or NZC, or combining both populations). Predictions were performed using the actual phenotypes (normalized) and the single-step GBLUP via Bayesian inference. Genotyped NWS animals born in 2016 (N = 267) were used to assess the accuracy and bias of genomic estimated breeding values (GEBVs) predicted for birth weight (BW), weaning weight (WW), carcass weight (CW), EUROP carcass classification (EUC), and EUROP fat grading (EUF). The accuracy and bias of GEBVs differed across traits and training population used. For instance, the GEBV accuracies ranged from 0.13 (BW) to 0.44 (EUC) for GEBVs predicted including only NWS, from 0.06 (BW) to 0.15 (CW) when including only NZC, and from 0.10 (BW) to 0.41 (EUC) when including both NWS and NZC animals in the training population. The regression coefficients used to assess the spread of GEBVs (bias) ranged from 0.26 (BW) to 0.64 (EUF) for only NWS, 0.10 (EUC) to 0.52 (CW) for only NZC, and from 0.42 (WW) to 2.23 (EUC) for both NWS and NZC in the training population. Our findings suggest that across-country genomic predictions based on ssGBLUP might be possible for NWS and NZC, especially for novel traits.

Genetic Connectedness Between Norwegian White Sheep and New Zealand Composite Sheep Populations With Similar Development History.

The Norwegian White sheep (NWS) and New Zealand Terminal Sire Composite (NZC) sheep breeds have been developed based on crossing of multiple breeds, mainly of Northern European origin. A close genetic relationship between these populations could enable across-country genomic evaluations. The main objectives of this study were to assess the genetic connectedness between Norwegian and New Zealand sheep populations and estimate numerous genetic diversity metrics for these two populations. A total of 792 NWS and 16,912 NZC animals were genotyped using a high-density Illumina SNP chip panel (∼606K SNPs). The NZC animals were grouped based on their breed composition as: Finn, Lamb Supreme, Primera, Texel, "Other Dual Purpose", and "Other Terminal Sire". The average level of linkage disequilibrium ranged from 0.156 (for Primera) to 0.231 (for Finn). The lowest consistency of gametic phase was estimated between NWS and Finn (0.397), and between NWS and Texel (0.443), respectively. Similar consistency of gametic phase was estimated between NWS and the other NZC populations (∼ 0.52). For all composite sheep populations analyzed in this study, the majority of runs of homozygosity (ROH) segments identified had short length (<2,500 kb), indicating ancient (instead of recent) inbreeding. The variation in the number of ROH segments observed in the NWS was similar to the variation observed in Primera and Lamb Supreme. There was no clear discrimination between NWS and NZC based on the first few principal components. In addition, based on admixture analyses, there seems to be a significant overlap of the ancestral populations that contributed to the development of both NWS and NZC. There were no evident signatures of selection in these populations, which might be due to recent crossbreeding. In conclusion, the NWS composite breed was shown to be moderately related to NZC populations, especially Primera and Lamb Supreme. The findings reported here indicate a promising opportunity for collaborative genomic analyses involving NWS and NZC sheep populations.

Genomic selection strategies to improve maternal traits in Norwegian White Sheep.

This study tested and compared different implementation strategies for genomic selection for Norwegian White Sheep, aiming to increase genetic gain for maternal traits. These strategies were evaluated for their genetic gain ingrowth, carcass and maternal traits, total genetic gain, a weighted sum of the gain in each trait and rates of inbreeding through a full-scale stochastic simulation. Results showed genomic selection schemes to increase genetic gain for maternal traits but reduced genetic gain for other traits. This could also be obtained by selecting rams for artificial selection at a higher age. Implementation of genomic selection in the current breeding structure increased genetic gain for maternal traits up to 57%, outcompeted by reducing the generation interval for artificial insemination rams from current 3 to 2 years. Then, total genetic gain for maternal traits increased by 65%-77% and total genetic gain by18%-20%, but at increased rates of inbreeding.

Impact of two myostatin (MSTN) mutations on weight gain and lamb carcass classification in Norwegian White Sheep (Ovis aries).

BACKGROUND: Our aim was to estimate the effect of two myostatin (MSTN) mutations in Norwegian White Sheep, one of which is close to fixation in the Texel breed. METHODS: The impact of two known MSTN mutations was examined in a field experiment with Norwegian White Sheep. The joint effect of the two MSTN mutations on live weight gain and weaning weight was studied on 644 lambs. Carcass weight gain from birth to slaughter, carcass weight, carcass conformation and carcass fat classes were calculated in a subset of 508 lambs. All analyses were carried out with a univariate linear animal model. RESULTS: The most significant impact of both mutations was on conformation and fat classes. The largest difference between the genotype groups was between the wild type for both mutations and the homozygotes for the c.960delG mutation. Compared to the wild types, these mutants obtained a conformation score 5.1 classes higher and a fat score 3.0 classes lower, both on a 15-point scale. CONCLUSIONS: Both mutations reduced fatness and increased muscle mass, although the effect of the frameshift mutation (c.960delG) was more important as compared to the 3'-UTR mutation (c.2360G>A). Lambs homozygous for the c.960delG mutation grew more slowly than those with other MSTN genotypes, but had the least fat and the largest muscle mass. Only c.960delG showed dominance effects.