Understanding the genetic changes underlying phenotypic variation in sheep (Ovis aries) may facilitate our efforts towards further improvement. Here, we report the deep resequencing of 248 sheep including the wild ancestor (O. orientalis), landraces, and improved breeds. We explored the sheep variome and selection signatures. We detected genomic regions harboring genes associated with distinct morphological and agronomic traits, which may be past and potential future targets of domestication, breeding, and selection. Furthermore, we found non-synonymous mutations in a set of plausible candidate genes and significant differences in their allele frequency distributions across breeds. We identified PDGFD as a likely causal gene for fat deposition in the tails of sheep through transcriptome, RT-PCR, qPCR, and Western blot analyses. Our results provide insights into the demographic history of sheep and a valuable genomic resource for future genetic studies and improved genome-assisted breeding of sheep and other domestic animals.

Domestic sheep and their wild relatives harbor substantial genetic variants that can form the backbone of molecular breeding, but their genome landscapes remain understudied. Here, we present a comprehensive genome resource for wild ovine species, landraces and improved breeds of domestic sheep, comprising high-coverage (∼16.10×) whole genomes of 810 samples from 7 wild species and 158 diverse domestic populations. We detected, in total, ∼121.2 million single nucleotide polymorphisms, ∼61 million of which are novel. Some display significant (P < 0.001) differences in frequency between wild and domestic species, or are private to continent-wide or individual sheep populations. Retained or introgressed wild gene variants in domestic populations have contributed to local adaptation, such as the variation in the HBB associated with plateau adaptation. We identified novel and previously reported targets of selection on morphological and agronomic traits such as stature, horn, tail configuration, and wool fineness. We explored the genetic basis of wool fineness and unveiled a novel mutation (chr25: T7,068,586C) in the 3′-UTR of IRF2BP2 as plausible causal variant for fleece fiber diameter. We reconstructed prehistorical migrations from the Near Eastern domestication center to South-and-Southeast Asia and found two main waves of migrations across the Eurasian Steppe and the Iranian Plateau in the Early and Late Bronze Ages. Our findings refine our understanding of genome variation as shaped by continental migrations, introgression, adaptation, and selection of sheep.

Goat domestication was critical for agriculture and civilization, but its underlying genetic changes and selection regimes remain unclear. Here, we analyze the genomes of worldwide domestic goats, wild caprid species, and historical remains, providing evidence of an ancient introgression event from a West Caucasian tur-like species to the ancestor of domestic goats. One introgressed locus with a strong signature of selection harbors the MUC6 gene, which encodes a gastrointestinally secreted mucin. Experiments revealed that the nearly fixed introgressed haplotype confers enhanced immune resistance to gastrointestinal pathogens. Another locus with a strong signal of selection may be related to behavior. The selected alleles at these two loci emerged in domestic goats at least 7200 and 8100 years ago, respectively, and increased to high frequencies concurrent with the expansion of the ubiquitous modern mitochondrial haplogroup A. Tracking these archaeologically cryptic evolutionary transformations provides new insights into the mechanisms of animal domestication.

Pistachio (Pistacia vera), one of the most important commercial nut crops worldwide, is highly adaptable to abiotic stresses and is tolerant to drought and salt stresses. Here, we provide a draft de novo genome of pistachio as well as large-scale genome resequencing. Comparative genomic analyses reveal stress adaptation of pistachio is likely attributable to the expanded cytochrome P450 and chitinase gene families. Particularly, a comparative transcriptomic analysis shows that the jasmonic acid (JA) biosynthetic pathway plays an important role in salt tolerance in pistachio. Moreover, we resequence 93 cultivars and 14 wild P. vera genomes and 35 closely related wild Pistacia genomes, to provide insights into population structure, genetic diversity, and domestication. We find that frequent genetic admixture occurred among the different wild Pistacia species. Comparative population genomic analyses reveal that pistachio was domesticated about 8000 years ago and suggest that key genes for domestication related to tree and seed size experienced artificial selection. Our study provides insight into genetic underpinning of local adaptation and domestication of pistachio. The Pistacia genome sequences should facilitate future studies to understand the genetic basis of agronomically and environmentally related traits of desert crops.

This study investigated the effects of a SNP in the myostatin gene (MSTN or growth differentiation factor 8, GDF8) on birth, growth, carcass, and beef quality traits in Australia (Aust.) and New Zealand (NZ). The SNP is a cytosine to adenine transversion in exon 1, causing an amino acid substitution of leucine for phenylalanine⁹⁴ (F94L). The experiment used crosses between the Jersey and Limousin breeds, with the design being a backcross using first-cross bulls of Jersey × Limousin or Limousin × Jersey breeding, mated to Jersey and Limousin cows. Progeny were genotyped for the myostatin SNP and phenotyped in Aust., with finishing on feedlot (366 calves, over 3 birth years) and in NZ with finishing on pasture (416 calves, over 2 birth years). The effect of the F94L allele (A allele) on birth and growth traits was not significant. The F94L allele in Limousin backcross calves was associated with an increase in meat weight (7.3 and 5.9% of the trait mean in Aust. and NZ, respectively, P < 0.001), and a reduction in fat depth (−13.9 and −18.7% of the trait means on live calves (600 d) and carcasses, respectively, Aust. only, P < 0.001), intramuscular fat content (−8.2% of the trait mean in Aust., P < 0.05; −7.1% in NZ, not significant), total carcass fat weight (−16.5 and −8.1% of the trait mean, Aust. and NZ; P < 0.001 and P < 0.05, respectively). Meat tenderness, pH, and cooking loss of the M. longissimus dorsi were not affected by the F94L variant. In the Jersey backcross calves, additive and dominance effects were confounded because the F94L allele was not segregating in the Jersey dams. The combined effects, however, were significant on LM area (4.4% in both Aust., P < 0.05, and NZ, P < 0.01), channel fat (−11.7%, NZ only, P < 0.01), rib fat depth (−11.2%, NZ only, P < 0.05), and carcass fat weight (−7.1%, NZ only, P < 0.05). The results provide strong evidence that this myostatin F94L variant provides an intermediate and more useful phenotype than the more severe double-muscling phenotype caused by knockout mutations in the myostatin gene. Copyright © 2008. . Copyright 2008 Journal of Animal Science