Whole-genome resequencing reveals domestication and signatures of selection in Ujimqin,Sunit, and Wu Ranke Mongolian sheep breeds

ObjectiveThe current study aimed to perform whole-genome resequencing of Chinese indigenous Mongolian sheep breeds including Ujimqin, Sunit, and Wu Ranke sheep breeds (UJMQ, SNT, WRK) and deeply analyze genetic variation, population structure, domestication, and selection for domestication traits among these Mongolian sheep breeds.MethodsBlood samples were collected from a total of 60 individuals comprising 20 WRK, 20 UJMQ, and 20 SNT. For genome sequencing, about 1.5 μg of genomic DNA was used for library construction with an insert size of about 350 bp. Pair-end sequencing were performed on Illumina NovaSeq platform, with the read length of 150 bp at each end. We then investigated the domestication and signatures of selection in these sheep breeds.ResultsAccording to the population and demographic analyses, WRK and SNT populations were very similar, which were different from UJMQ populations. Genome wide association study identified 468 and 779 significant loci from SNT vs UJMQ, and UJMQ vs WRK, respectively. However, only 3 loci were identified from SNT vs WRK. Genomic comparison and selective sweep analysis among these sheep breeds suggested that genes associated with regulation of secretion, metabolic pathways including estrogen metabolism and amino acid metabolism, and neuron development have undergone strong selection during domestication.ConclusionOur findings will facilitate the understanding of Chinese indigenous Mongolian sheep breeds domestication and selection for complex traits and provide a valuable genomic resource for future studies of sheep and other domestic animal breeding.     

Genomic Diversity in Pig (Sus scrofa) and its Comparison with Human and other Livestock

We have reviewed the current pig (Sus scrofa) genomic diversity within and between sites and compared them with human and other livestock. The current Porcine 60K single nucleotide polymorphism (SNP) panel has an average SNP distance in a range of 30 - 40 kb. Most of genetic variation was distributed within populations, and only a small proportion of them existed between populations. The average heterozygosity was lower in pig than in human and other livestock. Genetic inbreeding coefficient (F(IS)), population differentiation (F(ST)), and Nei's genetic distance between populations were much larger in pig than in human and other livestock. Higher average genetic distance existed between European and Asian populations than between European or between Asian populations. Asian breeds harboured much larger variability and higher average heterozygosity than European breeds. The samples of wild boar that have been analyzed displayed more extensive genetic variation than domestic breeds. The average linkage disequilibrium (LD) in improved pig breeds extended to 1 - 3 cM, much larger than that in human (~ 30 kb) and cattle (~ 100 kb), but smaller than that in sheep (~ 10 cM). European breeds showed greater LD that decayed more slowly than Asian breeds. We briefly discuss some processes for maintaining genomic diversity in pig, including migration, introgression, selection, and drift. We conclude that, due to the long time of domestication, the pig possesses lower heterozygosity, higher F(IS), and larger LD compared with human and cattle. This implies that a smaller effective population size and less informative markers are needed in pig for genome wide association studies.     

Conservation genetics of cattle, sheep, and goats

Cattle, sheep and goats were domesticated about 10,000 years ago, spread out of the domestication centers in Europe, Asia, and Africa during the next few thousands years, and gave many populations locally adapted. After a very long period of soft selection, the situation changed dramatically 200 years ago with the emergence of the breed concept. The selection pressure strongly increased, and the reproduction among breeds was seriously reduced, leading to the fragmentation of the initial gene pool. More recently, the selection pressure was increased again via the use of artificial insemination, leading to a few industrial breeds with very high performances, but with low effective population sizes. Beside this performance improvement of industrial breeds, genetic resources are being lost, because of the replacement of traditional breeds by high performance industrial breeds at the worldwide level, and because of the loss of genetic diversity in these industrial breeds. Many breeds are already extinct, and genetic resources in cattle, sheep, and goats are thus highly endangered, particularly in developed countries. The recent development of next generation sequencing technologies opens new avenues for properly characterizing the genetic resources, not only in the very diverse domestic breeds, but also in their wild relatives. Based on sound genetic characterization, urgent conservation measures must be taken to avoid an irremediable loss of farm animal genetic resources, integrating economical, sociological, and political parameters.     

Genetic diversity of <Emphasis Type="Italic">Ovis aries</Emphasis> populations near domestication centers and in the New World

Domestic sheep in Kazakhstan may provide an interesting source of genetic variability due to their proximity to the center of domestication and the Silk Route. Additionally, those breeds have never been compared to New World sheep populations. This report compares genetic diversity among five Kazakhstan (KZ) and 13 United States (US) sheep breeds (N = 442) using 25 microsatellite markers from the FAO panel. The KZ breeds had observed and expected measures of heterozygosity greater than 0.60 and an average number of alleles per locus of 7.8. In contrast, US sheep breeds had observed heterozygosity ranged from 0.37 to 0.62 and had an average number of alleles of 5.7. A Bayesian analysis indicated there were two primary populations (K = 2). Surprisingly, the US breeds were near evenly split between the two clusters, while all of the KZ breeds were placed in one of the two clusters. Pooling breeds within country of sample origin showed KZ and US populations to have similar levels of expected heterozygosity and the average number of alleles per locus. The results of breeds pooled within country suggest that there was no difference between countries for these diversity measures using this set of neutral markers. This finding suggests that populations’ geographically isolated from centers of domestication can be more diverse than previously thought, and as a result, conservation strategies can be adjusted accordingly. Furthermore, these results suggest there may be limited need for countries to alter the protocols for trade and exchange of animal genetic resources that are in place today, since no one population has a unique set of private alleles.     

Linkage disequilibrium over short physical distances measured in sheep using a high‐density SNP chip

The extent of linkage disequilibrium (LD) between genetic loci has implications for both association studies and the accuracy of genomic prediction. To characterise the persistence of LD in diverse sheep breeds, two SNP genotyping platforms were used. First, existing SNP genotypes from 63 breeds obtained using the ovine SNP50 BeadChip (49 034 loci) were used to estimate LD decay in populations with contrasting levels of genetic diversity. Given the paucity of marker pairs separated by short physical distances on the SNP50 BeadChip, genotyping was subsequently performed for four breeds using the recently developed ovine HD BeadChip that assays approximately 600 000 SNPs with an average genomic spacing of 5 kb. This facilitated a highly accurate estimate of LD over short genomic distances (<30 kb) and revealed LD varies considerably between sheep breeds. Further, sheep appear to contain generally lower levels of LD than do other domestic species, likely a reflection of aspects of their past population history.     

Origin, genetic diversity, and population structure of Chinese domestic sheep

To characterize the origin, genetic diversity, and phylogeographic structure of Chinese domestic sheep, we here analyzed a 531-bp fragment of mtDNA control region of 449 Chinese autochthonous sheep from 19 breeds/populations from 13 geographic regions, together with previously reported 44 sequences from Chinese indigenous sheep. Phylogenetic analysis showed that all three previously defined lineages A, B, and C were found in all sampled Chinese sheep populations, except for the absence of lineage C in four populations. Network profiles revealed that the lineages B and C displayed a star-like phylogeny with the founder haplotype in the centre, and that two star-like subclades with two founder haplotypes were identified in lineage A. The pattern of genetic variation in lineage A, together with the divergence time between the two central founder haplotypes suggested that two independent domestication events have occurred in sheep lineage A. Considerable mitochondrial diversity was observed in Chinese sheep. Weak structuring was observed either among Chinese indigenous sheep populations or between Asian and European sheep and this can be attributable to long-term strong gene flow induced by historical human movements. The high levels of intra-population diversity in Chinese sheep and the weak phylogeographic structuring indicated three geographically independent domestication events have occurred and the domestication place was not only confined to the Near East, but also occurred in other regions.     

Genetic structure of European sheep breeds

Large-scale evaluations of genetic diversity in domestic livestock populations are necessary so that region-specific conservation measures can be implemented. We performed the first such survey in European sheep by analysing 820 individuals from 29 geographically and phenotypically diverse breeds and a closely related wild species at 23 microsatellite loci. In contrast to most other domestic species, we found evidence of widespread heterozygote deficit within breeds, even after removing loci with potentially high frequency of null alleles. This is most likely due to subdivision among flocks (Wahlund effect) and use of a small number of rams for breeding. Levels of heterozygosity were slightly higher in southern than in northern breeds, consistent with declining diversity with distance from the Near Eastern centre of domestication. Our results highlight the importance of isolation in terms of both geography and management in augmenting genetic differentiation through genetic drift, with isolated northern European breeds showing the greatest divergence and hence being obvious targets for conservation. Finally, using a Bayesian cluster analysis, we uncovered evidence of admixture between breeds, which has important implications for breed management.     

Evidence of three maternal lineages in Near Eastern sheep supporting multiple domestication events

The variability of mtDNA was analysed in local sheep breeds reared throughout Turkey, for which a fragment of the D-loop region and the complete cytochrome b were sequenced. Phylogenetic analyses performed independently for the D-loop and the Cyt b gene revealed three clearly separated clusters indicating three major maternal lineages, two of which had been previously described as types B and A. The new type, C, was present in all the breeds analysed and showed considerable mtDNA variability. Divergence time was obtained on the basis of Cyt b gene and was estimated to be around 160,000-170,000 years ago for lineages B and A, whereas the divergence of lineage C proved to have occurred earlier (between 450,000 and 750,000 years ago). These times greatly predate domestication and suggest that the origin of modern sheep breeds was more complex than previously thought and that at least three independent sheep domestication events occurred. Our results, together with archaeological information and the current wild sheep populations in the Near East region support the high importance of this area in the sheep domestication process. Finally, the evidence of a third maternal lineage has important implications regarding the history of modern sheep.     

Genetic diversity of three Chinese native sheep breeds

The sheep (Ovis aries L.) has been an important farm animal species since its domestication. A wide array of indigenous sheep breeds with abundant phenotypic diversity exists for domestication and selection. Therefore, assessing the genetic diversity of a local sheep resource using a multi-molecular system is helpful for maintaining and conserving those breeds. This study aimed to investigate the genetic diversity of three native Chinese sheep breeds (Tibetan sheep, Sishui Fur sheep, and Small-tailed Han sheep) using 15 microsatellite markers and the second exon of the DRA gene. In regards to the microsatellites, on average, 19 alleles per loci were observed among all individuals. Across loci, the HO within the population was 0.652 ± 0.022 in Tibetan sheep, 0.603 ± 0.023 in Small-tailed Han sheep and 0.635 ± 0.022 in SFS, and for most populations, the H _E and H _O were inconsistent. In addition, affluent private alleles within the breed indicated that the breeds have different domestication histories or sites. In regards to the 2 exon of the DRA gene, three haplotypes were constructed by seven single-nucleotide polymorphisms (SNPs), which were identified in the second DRA exon and inferred the potential for phenotypic variety in these Chinese native sheep. In summary, the current study reveals the importance of implementing effective conservation strategies for these three native Chinese sheep.     

Conservation genetics in Chinese sheep: diversity of fourteen indigenous sheep (Ovis aries) using microsatellite markers

The domestic sheep (Ovis aries) has been an economically and culturally important farm animal species since its domestication around the world. A wide array of sheep breeds with abundant phenotypic diversity exists including domestication and selection as well as the indigenous breeds may harbor specific features as a result of adaptation to their environment. The objective of this study was to investigate the population structure of indigenous sheep in a large geographic location of the Chinese mainland. Six microsatellites were genotyped for 611 individuals from 14 populations. The mean number of alleles (±SD) ranged from 7.00 ± 3.69 in Gangba sheep to 10.50 ± 4.23 in Tibetan sheep. The observed heterozygote frequency (±SD) within a population ranged from 0.58 ± 0.03 in Gangba sheep to 0.71 ± 0.03 in Zazakh sheep and Minxian black fur sheep. In addition, there was a low pairwise difference among the Minxian black fur sheep, Mongolian sheep, Gansu alpine merino, and Lanzhou fat‐tailed sheep. Bayesian analysis with the program STRUCTURE showed support for 3 clusters, revealing a vague genetic clustering pattern with geographic location. The results of the current study inferred high genetic diversity within these native sheep in the Chinese mainland.     

Genetic variation of 5 SNPs of MC1R gene in Chinese indigenous sheep breeds

The purpose of this study was to assess genetic diversity, genetic differentiation relationship and population structure among 10 Chinese sheep populations using 5 single nucleotide polymorphisms (SNPs) in MC1R gene. The genetic diversity indices suggested that the intra-population variation levels of Chinese Merino and Large-tailed Han breeds were lowest than Kazakh Fat-Rumped. Chinese sheep breeds have maintained a high intra-population variation levels (95.23%). The genetic differentiation patterns and genetic relationships among Chinese sheep breeds displayed a high consistency with the traditional classification. The cluster trees were constructed by UPMGA method. The results showed that Chinese indigenous sheep populations have distinct genetic differentiation. The inter-population variation levels in Chinese sheep populations indicated three geographically independent domestication events have occurred. The Bayesian cluster analyses also showed a reliable clustering pattern, which revealed three major clusters in Chinese indigenous sheep populations (Mongolian group, Kazakh group and Tibetan group), except for Duolang and Minxian Black-fur. There were probably caused by different breeding history, geography isolation and different levels of inbreeding. The findings supported the related records in literature, ten sheep populations originated on different time stage from the primogenitor population and communicated genetically with each other in the process of natural and artificial selection, and in different ecological environment. It is concluded that Chinese indigenous sheep have higher genetic variation and diversity, genetic differentiation exist between Chinese sheep populations.     

Evolutionary genomics of dog domestication

We review the underlying principles and tools used in genomic studies of domestic dogs aimed at understanding the genetic changes that have occurred during domestication. We show that there are two principle modes of evolution within dogs. One primary mode that accounts for much of the remarkable diversity of dog breeds is the fixation of discrete mutations of large effect in individual lineages that are then crossed to various breed groupings. This transfer of mutations across the dog evolutionary tree leads to the appearance of high phenotypic diversity that in actuality reflects a small number of major genes. A second mechanism causing diversification involves the selective breeding of dogs within distinct phenotypic or functional groups, which enhances specific group attributes such as heading or tracking. Such progressive selection leads to a distinct genetic structure in evolutionary trees such that functional and phenotypic groups cluster genetically. We trace the origin of the nuclear genome in dogs based on haplotype-sharing analyses between dogs and gray wolves and show that contrary to previous mtDNA analyses, the nuclear genome of dogs derives primarily from Middle Eastern or European wolves, a result more consistent with the archeological record. Sequencing analysis of the IGF1 gene, which has been the target of size selection in small breeds, further supports this conclusion. Finally, we discuss how a black coat color mutation that evolved in dogs has transformed North American gray wolf populations, providing a first example of a mutation that appeared under domestication and selectively swept through a wild relative.     

Ear type in sheep is associated with the MSRB3 locus

Ear morphology is an important determinant of sheep breeds. It includes different variable traits such as ear size and erectness, suggesting a polygenic architecture. Here, we performed a comprehensive genome-wide analysis to identify regions under selection for ear morphology in 515 sheep from 17 breeds fixed for diverse ear phenotypes using 34k SNP genotyping data. GWASs for two ear type traits, size and erectness, revealed a single genome-wide significant association on ovine chromosome 3. The derived marker alleles were enriched in sheep with large and/or floppy ears. The GWAS signal harboured the MSRB3 gene encoding methionine sulphoxide reductase B3, which has already been found to be associated with different ear types in other species. We attempted whole-genome resequencing to identify causal variant(s) within a 1 Mb interval around MSRB3. This experiment excluded major copy number variants in the interval, but failed to identify a compelling candidate causal variant. Fine-mapping suggested that the causal variant for large floppy ears most likely resides in a 175 kb interval downstream of the MSRB3 coding region.     

The legacy of domestication: accumulation of deleterious mutations in the dog genome

Dogs exhibit more phenotypic variation than any other mammal and are affected by a wide variety of genetic diseases. However, the origin and genetic basis of this variation is still poorly understood. We examined the effect of domestication on the dog genome by comparison with its wild ancestor, the gray wolf. We compared variation in dog and wolf genes using whole-genome single nucleotide polymorphism (SNP) data. The d(N)/d(S) ratio (omega) was around 50% greater for SNPs found in dogs than in wolves, indicating that a higher proportion of nonsynonymous alleles segregate in dogs compared with nonfunctional genetic variation. We suggest that the majority of these alleles are slightly deleterious and that two main factors may have contributed to their increase. The first is a relaxation of selective constraint due to a population bottleneck and altered breeding patterns accompanying domestication. The second is a reduction of effective population size at loci linked to those under positive selection due to Hill-Robertson interference. An increase in slightly deleterious genetic variation could contribute to the prevalence of disease in modern dog breeds.     

Runs of homozygosity and signatures of selection: a comparison among eight local Swiss sheep breeds

A dataset consisting of 787 animals with high‐density SNP chip genotypes (346 774 SNPs) and 939 animals with medium‐density SNP chip genotypes (33 828 SNPs) from eight indigenous Swiss sheep breeds was analyzed to characterize population structure, quantify genomic inbreeding based on runs of homozygosity and identify selection signatures. In concordance with the recent known history of these breeds, the highest genetic diversity was observed in Engadine Red sheep and the lowest in Valais Blacknose sheep. Correlation between F_PED and F_ROH was around 0.50 and thereby lower than that found in similar studies in cattle. Mean F_ROH estimates from medium‐density data and HD data were highly correlated (0.95). Signatures of selection and candidate gene analysis revealed that the most prominent signatures of selection were found in the proximity of genes associated with body size (NCAPG, LCORL, LAP3, SPP1, PLAG1, ALOX12, TP53), litter size (SPP1), milk production (ABCG2, SPP1), coat color (KIT, ASIP, TBX3) and horn status (RXFP2). For the Valais Blacknose sheep, the private signatures in proximity of genes/QTL influencing body size, coat color and fatty acid composition were confirmed based on runs of homozygosity analysis. These private signatures underline the genetic uniqueness of the Valais Blacknose sheep breed. In conclusion, we identified differences in the genetic make‐up of Swiss sheep breeds, and we present relevant candidate genes responsible for breed differentiation in locally adapted breeds.     

Insights into the Genetic History of French Cattle from Dense SNP Data on 47 Worldwide Breeds

Background

Modern cattle originate from populations of the wild extinct aurochs through a few domestication events which occurred about 8,000 years ago. Newly domesticated populations subsequently spread worldwide following breeder migration routes. The resulting complex historical origins associated with both natural and artificial selection have led to the differentiation of numerous different cattle breeds displaying a broad phenotypic variety over a short period of time.

Methodology/Principal Findings

This study gives a detailed assessment of cattle genetic diversity based on 1,121 individuals sampled in 47 populations from different parts of the world (with a special focus on French cattle) genotyped for 44,706 autosomal SNPs. The analyzed data set consisted of new genotypes for 296 individuals representing 14 French cattle breeds which were combined to those available from three previously published studies. After characterizing SNP polymorphism in the different populations, we performed a detailed analysis of genetic structure at both the individual and population levels. We further searched for spatial patterns of genetic diversity among 23 European populations, most of them being of French origin, under the recently developed spatial Principal Component analysis framework.

Conclusions/Significance

Overall, such high throughput genotyping data confirmed a clear partitioning of the cattle genetic diversity into distinct breeds. In addition, patterns of differentiation among the three main groups of populations—the African taurine, the European taurine and zebus—may provide some additional support for three distinct domestication centres. Finally, among the European cattle breeds investigated, spatial patterns of genetic diversity were found in good agreement with the two main migration routes towards France, initially postulated based on archeological evidence.     

Signatures of selection in the genomes of commercial and non-commercial chicken breeds

Identifying genomics regions that are affected by selection is important to understand the domestication and selection history of the domesticated chicken, as well as understanding molecular pathways underlying phenotypic traits and breeding goals. While whole-genome approaches, either high-density SNP chips or massively parallel sequencing, have been successfully applied to identify evidence for selective sweeps in chicken, it has been difficult to distinguish patterns of selection and stochastic and breed specific effects. Here we present a study to identify selective sweeps in a large number of chicken breeds (67 in total) using a high-density (58 K) SNP chip. We analyzed commercial chickens representing all major breeding goals. In addition, we analyzed non-commercial chicken diversity for almost all recognized traditional Dutch breeds and a selection of representative breeds from China. Based on their shared history or breeding goal we in silico grouped the breeds into 14 breed groups. We identified 396 chromosomal regions that show suggestive evidence of selection in at least one breed group with 26 of these regions showing strong evidence of selection. Of these 26 regions, 13 were previously described and 13 yield new candidate genes for performance traits in chicken. Our approach demonstrates the strength of including many different populations with similar, and breed groups with different selection histories to reduce stochastic effects based on single populations.     

Population structure in Indian sheep ascertained using microsatellite information

This study attempts to provide a comprehensive insight into the prevailing genetic status of Indian sheep breeds using microsatellite markers. Seventeen Indian sheep breeds from 3 agroecological zones were analysed using a panel of 25 microsatellite markers. All of the sheep breeds investigated were genetically diverse, as evident from the high allele (>6) and gene (>0.6) diversity values. The gene diversity values for all breeds ranged from 0.621 to 0.780. The within-population heterozygote deficit (F(IS)) varied from -0.098 to 0.234, reflecting significant levels for 12 of the 17 breeds investigated. The average genetic differentiation between all breeds (F(ST)) was 11.1%, revealing moderate discrimination between the indigenous sheep breeds. The genetic distance and principal component analysis revealed a separation of sheep breeds based on geographical propinquity. The Bayesian clustering approach suggested poor breed differentiation in the north-western arid and semi-arid region when compared to the breeds from the eastern and southern peninsular regions. The observed results mirror the divergent management strategies in the different agroecological regions, lack of specific selection policies, and intermixing of breeds in close proximity. Immediate steps to curb the intermixing and erosion of breed purity for some of these breeds need to be implemented, for example, by introducing measures like making proven rams available and ensuring their frequent exchange between flocks. The data generated here provides valuable information about the genetic structure of the 17 Indian sheep breeds and this can be used for designating priorities for their conservation.     

Genetic relationships among Vietnamese local pigs investigated using genome‐wide SNP markers

Vietnam is one of the most important countries for pig domestication, and a total of 26 local breeds have been reported. In the present study, genetic relationships among the various pig breeds were investigated using 90 samples collected from local pigs (15 breeds) in 15 distantly separated, distinct areas of the country and six samples from Landrace pigs in Hanoi as an out‐group of a common Western breed. All samples were genotyped using the Illumina Porcine SNP60 v2 Genotyping BeadChip. We used 15 160–15 217 SNPs that showed a high degree of polymorphism in the Vietnamese breeds for identifying genetic relationships among the Vietnamese breeds. Principal components analysis showed that most pigs indigenous to Vietnam formed clusters correlated with their original geographic locations. Some Vietnamese breeds formed a cluster that was genetically related to the Western breed Landrace, suggesting the possibility of crossbreeding. These findings will be useful for the conservation and management of Vietnamese local pig breeds.     

Signatures of Diversifying Selection in European Pig Breeds

Following domestication, livestock breeds have experienced intense selection pressures for the development of desirable traits. This has resulted in a large diversity of breeds that display variation in many phenotypic traits, such as coat colour, muscle composition, early maturity, growth rate, body size, reproduction, and behaviour. To better understand the relationship between genomic composition and phenotypic diversity arising from breed development, the genomes of 13 traditional and commercial European pig breeds were scanned for signatures of diversifying selection using the Porcine60K SNP chip, applying a between-population (differentiation) approach. Signatures of diversifying selection between breeds were found in genomic regions associated with traits related to breed standard criteria, such as coat colour and ear morphology. Amino acid differences in the EDNRB gene appear to be associated with one of these signatures, and variation in the KITLG gene may be associated with another. Other selection signals were found in genomic regions including QTLs and genes associated with production traits such as reproduction, growth, and fat deposition. Some selection signatures were associated with regions showing evidence of introgression from Asian breeds. When the European breeds were compared with wild boar, genomic regions with high levels of differentiation harboured genes related to bone formation, growth, and fat deposition.