Convergent genomic signatures of domestication in sheep and goats

Large-scale genomic analysis reveals the genetic cost of chicken domestication

BACKGROUND

Species domestication is generally characterized by the exploitation of high-impact mutations through processes that involve complex shifting demographics of domesticated species. These include not only inbreeding and artificial selection that may lead to the emergence of evolutionary bottlenecks, but also post-divergence gene flow and introgression. Although domestication potentially affects the occurrence of both desired and undesired mutations, the way wild relatives of domesticated species evolve and how expensive the genetic cost underlying domestication is remain poorly understood. Here, we investigated the demographic history and genetic load of chicken domestication.

RESULTS

We analyzed a dataset comprising over 800 whole genomes from both indigenous chickens and wild jungle fowls. We show that despite having a higher genetic diversity than their wild counterparts (average π, 0.00326 vs. 0.00316), the red jungle fowls, the present-day domestic chickens experienced a dramatic population size decline during their early domestication. Our analyses suggest that the concomitant bottleneck induced 2.95% more deleterious mutations across chicken genomes compared with red jungle fowls, supporting the "cost of domestication" hypothesis. Particularly, we find that 62.4% of deleterious SNPs in domestic chickens are maintained in heterozygous states and masked as recessive alleles, challenging the power of modern breeding programs to effectively eliminate these genetic loads. Finally, we suggest that positive selection decreases the incidence but increases the frequency of deleterious SNPs in domestic chicken genomes.

CONCLUSION

This study reveals a new landscape of demographic history and genomic changes associated with chicken domestication and provides insight into the evolutionary genomic profiles of domesticated animals managed under modern human selection.

Herded and hunted goat genomes from the dawn of domestication in the Zagros Mountains

The Aceramic Neolithic (∼9600 to 7000 cal BC) period in the Zagros Mountains, western Iran, provides some of the earliest archaeological evidence of goat (Capra hircus) management and husbandry by circa 8200 cal BC, with detectable morphological change appearing ∼1,000 y later. To examine the genomic imprint of initial management and its implications for the goat domestication process, we analyzed 14 novel nuclear genomes (mean coverage 1.13X) and 32 mitochondrial (mtDNA) genomes (mean coverage 143X) from two such sites, Ganj Dareh and Tepe Abdul Hosein. These genomes show two distinct clusters: those with domestic affinity and a minority group with stronger wild affinity, indicating that managed goats were genetically distinct from wild goats at this early horizon. This genetic duality, the presence of long runs of homozygosity, shared ancestry with later Neolithic populations, a sex bias in archaeozoological remains, and demographic profiles from across all layers of Ganj Dareh support management of genetically domestic goat by circa 8200 cal BC, and represent the oldest to-this-date reported livestock genomes. In these sites a combination of high autosomal and mtDNA diversity, contrasting limited Y chromosomal lineage diversity, an absence of reported selection signatures for pigmentation, and the wild morphology of bone remains illustrates domestication as an extended process lacking a strong initial bottleneck, beginning with spatial control, demographic manipulation via biased male culling, captive breeding, and subsequently phenotypic and genomic selection.

The assembly of caprine Y chromosome sequence reveals a unique paternal phylogenetic pattern and improves our understanding of the origin of domestic goat

The mammalian Y chromosome offers a unique perspective on the male reproduction and paternal evolutionary histories. However, further understanding of the Y chromosome biology for most mammals is hindered by the lack of a Y chromosome assembly. This study presents an integrated in silico strategy for identifying and assembling the goat Y-linked scaffolds using existing data. A total of 11.5 Mb Y-linked sequences were clustered into 33 scaffolds, and 187 protein-coding genes were annotated. We also identified high abundance of repetitive elements. A 5.84 Mb subset was further ordered into an assembly with the evidence from the goat radiation hybrid map (RH map). The existing whole-genome resequencing data of 96 goats (worldwide distribution) were utilized to exploit the paternal relationships among bezoars and domestic goats. Goat paternal lineages were clearly divided into two clades (Y1 and Y2), predating the goat domestication. Demographic history analyses indicated that maternal lineages experienced a bottleneck effect around 2,000 YBP (years before present), after which goats belonging to the A haplogroup spread worldwide from the Near East. As opposed to this, paternal lineages experienced a population decline around the 10,000 YBP. The evidence from the Y chromosome suggests that male goats were not affected by the A haplogroup worldwide transmission, which implies sexually unbalanced contribution to the goat trade and population expansion in post-Neolithic period.

Genomics of Adaptations in Ungulates

Ungulates are a group of hoofed animals that have long interacted with humans as essential sources of food, labor, clothing, and transportation. These consist of domesticated, feral, and wild species raised in a wide range of habitats and biomes. Given the diverse and extreme environments inhabited by ungulates, unique adaptive traits are fundamental for fitness. The documentation of genes that underlie their genomic signatures of selection is crucial in this regard. The increasing availability of advanced sequencing technologies has seen the rapid growth of ungulate genomic resources, which offers an exceptional opportunity to understand their adaptive evolution. Here, we summarize the current knowledge on evolutionary genetic signatures underlying the adaptations of ungulates to different habitats.

Accessing Livestock Resources in Ensembl

Genome assembly is cheaper, more accurate and more automated than it has ever been. This is due to a combination of more cost-efficient chemistries, new sequencing technologies and better algorithms. The livestock community has been at the forefront of this new wave of genome assembly, generating some of the highest quality vertebrate genome sequences. Ensembl's goal is to add functional and comparative annotation to these genomes, through our gene annotation, genomic alignments, gene trees, regulatory, and variation data. We run computationally complex analyses in a high throughput and consistent manner to help accelerate downstream science. Our livestock resources are continuously growing in both breadth and depth. We annotate reference genome assemblies for newly sequenced species and regularly update annotation for existing genomes. We are the only major resource to support the annotation of breeds and other non-reference assemblies. We currently provide resources for 13 pig breeds, maternal and paternal haplotypes for hybrid cattle and various other non-reference or wild type assemblies for livestock species. Here, we describe the livestock data present in Ensembl and provide protocols for how to view data in our genome browser, download via it our FTP site, manipulate it via our tools and interact with it programmatically via our REST API.

A Whole-Genome Scan for Association with Invasion Success in the Fruit Fly Drosophila suzukii Using Contrasts of Allele Frequencies Corrected for Population Structure

Evidence is accumulating that evolutionary changes are not only common during biological invasions but may also contribute directly to invasion success. The genomic basis of such changes is still largely unexplored. Yet, understanding the genomic response to invasion may help to predict the conditions under which invasiveness can be enhanced or suppressed. Here, we characterized the genome response of the spotted wing drosophila Drosophila suzukii during the worldwide invasion of this pest insect species, by conducting a genome-wide association study to identify genes involved in adaptive processes during invasion. Genomic data from 22 population samples were analyzed to detect genetic variants associated with the status (invasive versus native) of the sampled populations based on a newly developed statistic, we called C2, that contrasts allele frequencies corrected for population structure. We evaluated this new statistical framework using simulated data sets and implemented it in an upgraded version of the program BayPass. We identified a relatively small set of single-nucleotide polymorphisms that show a highly significant association with the invasive status of D. suzukii populations. In particular, two genes, RhoGEF64C and cpo, contained single-nucleotide polymorphisms significantly associated with the invasive status in the two separate main invasion routes of D. suzukii. Our methodological approaches can be applied to any other invasive species, and more generally to any evolutionary model for species characterized by nonequilibrium demographic conditions for which binary covariables of interest can be defined at the population level.

British Sheep Breeds as a Part of World Sheep Gene Pool Landscape: Looking into Genomic Applications

Sheep farming has been an important sector of the UK's economy and rural life for many centuries. It is the favored source of wool, meat and milk products. In the era of exponential progress in genomic technologies, we can now address the questions of what is special about UK sheep breed genotypes and how they differ genetically form one another and from other countries. We can reflect how their natural history has been determined at the level of their genetic code and what traces have been left in their genomes because of selection for phenotypic traits. These include adaptability to certain environmental conditions and management, as well as resistance to disease. Application of these advancements in genetics and genomics to study sheep breeds of British domestic selection has begun and will continue in order to facilitate conservation solutions and production improvement.

Genome-Wide Analyses Reveal Genetic Convergence of Prolificacy between Goats and Sheep

The litter size of domestic goats and sheep is an economically important trait that shows variation within breeds. Strenuous efforts have been made to understand the genetic mechanisms underlying prolificacy in goats and sheep. However, there has been a paucity of research on the genetic convergence of prolificacy between goats and sheep, which likely arose because of similar natural and artificial selection forces. Here, we performed comparative genomic and transcriptomic analyses to identify the genetic convergence of prolificacy between goats and sheep. By combining genomic and transcriptomic data for the first time, we identified this genetic convergence in (1) positively selected genes (CHST11 and SDCCAG8), (2) differentially expressed genes (SERPINA14, RSAD2, and PPIG at follicular phase, and IGF1, GPRIN3, LIPG, SLC7A11, and CHST15 at luteal phase), and (3) biological pathways (genomic level: osteoclast differentiation, ErbB signaling pathway, and relaxin signaling pathway; transcriptomic level: the regulation of viral genome replication at follicular phase, and protein kinase B signaling and antigen processing and presentation at luteal phase). These results indicated the potential physiological convergence and enhanced our understanding of the overlapping genetic makeup underlying litter size in goats and sheep.

GGVD: A goat genome variation database for tracking the dynamic evolutionary process of selective signatures and ancient introgressions

Understanding the evolutionary history and adaptive process depends on the knowledge that we can acquire from both ancient and modern genomic data. With the availability of a deluge of whole-genome sequencing data from ancient and modern goat samples, a user-friendly database making efficient reuse of these important resources is needed. Here, we use the genomes of 208 modern domestic goats, 24 bezoars, 46 wild ibexes, and 82 ancient goats to present a comprehensive goat genome variation database (GGVD). GGVD hosts a total of ∼41.44 million SNPs, ∼5.14 million indels, 6,193 selected loci, and 112 introgression regions. Users can freely visualize the frequency of genomic variations in geographical maps, selective sweeps in interactive tables, Manhattan plots, or line charts, as well as the heatmap patterns of the SNP genotype. Ancient data can be shown in haplotypes to track the state of genetic variants of selection and introgression events in the early, middle, and late stages. For facilitating access to sequence features, the UCSC Genome Browser, BLAT, BLAST, LiftOver, and pcadapt are also integrated into GGVD. GGVD will be a convenient tool for population genetic studies and molecular marker designing in goat breeding programs, and it is publicly available at http://animal.nwsuaf.edu.cn/GoatVar.

Historical Introgression from Wild Relatives Enhanced Climatic Adaptation and Resistance to Pneumonia in Sheep

How animals, particularly livestock, adapt to various climates and environments over short evolutionary time is of fundamental biological interest. Further, understanding the genetic mechanisms of adaptation in indigenous livestock populations is important for designing appropriate breeding programs to cope with the impacts of changing climate. Here, we conducted a comprehensive genomic analysis of diversity, interspecies introgression, and climate-mediated selective signatures in a global sample of sheep and their wild relatives. By examining 600K and 50K genome-wide single nucleotide polymorphism data from 3,447 samples representing 111 domestic sheep populations and 403 samples from all their seven wild relatives (argali, Asiatic mouflon, European mouflon, urial, snow sheep, bighorn, and thinhorn sheep), coupled with 88 whole-genome sequences, we detected clear signals of common introgression from wild relatives into sympatric domestic populations, thereby increasing their genomic diversities. The introgressions provided beneficial genetic variants in native populations, which were significantly associated with local climatic adaptation. We observed common introgression signals of alleles in olfactory-related genes (e.g., ADCY3 and TRPV1) and the PADI gene family including in particular PADI2, which is associated with antibacterial innate immunity. Further analyses of whole-genome sequences showed that the introgressed alleles in a specific region of PADI2 (chr2: 248,302,667-248,306,614) correlate with resistance to pneumonia. We conclude that wild introgression enhanced climatic adaptation and resistance to pneumonia in sheep. This has enabled them to adapt to varying climatic and environmental conditions after domestication.

Goats show higher behavioural flexibility than sheep in a spatial detour task

The ability to adapt to changing environments is crucial for survival and has evolved based on socio-ecological factors. Goats and sheep are closely related, with similar social structures, body sizes and domestication levels, but different feeding ecologies, i.e. goats are browsers and sheep are grazers. We investigated whether goats' reliance on more patchily distributed food sources predicted an increased behavioural flexibility compared to sheep. We tested 21 goats and 28 sheep in a spatial A-not-B detour task. Subjects had to navigate around a straight barrier through a gap at its edge. After one, two, three or four of these initial A trials, the gap was moved to the opposite end and subjects performed four B trials. Behaviourally more flexible individuals should move through the new gap faster, while those less behaviourally flexible should show greater perseveration. While both species showed an accuracy reduction following the change of the gap position, goats recovered from this perseveration error from the second B trial onwards, whereas sheep did so only in the fourth B trial, indicating differences in behavioural flexibility between the species. This higher degree of flexibility in goats compared to sheep could be linked to differences in their foraging strategies.

Documentation of extensive genetic diversity in the Ovar-DRB1 gene in native Turkish sheep

Indigenous breeds have a high level of genetic diversity that might contribute to develop animal breeds with desired traits such as disease resistance and high productivity. Major histocompatibility complex (MHC) is a key component of adaptive immune system and consists of highly polymorphic genes that take part in adaptive immune response and disease resistance. Exploring and understanding the effect of polymorphisms in MHC could be beneficial to future animal breeding strategies. In this study, we sequenced the highly polymorphic Exon2 of the ovine DRB1 gene using Sanger sequencing to explore the diversity of this gene in six indigenous Turkish sheep breeds and two crossbreeds. In total, 894 haplotypes from 447 sheep were investigated, and 69 different haplotypes including 27 novel ones were identified. Among the identified haplotypes there were common and breed specific haplotypes. There was a relatively high diversity of the alleles within indigenous breeds. Allelic diversity patterns were mostly associated with geographical differences. The results of this study highlight the genetic variation within indigenous breeds which has important implications for biodiversity and the adaptability of breeds to specific environments. There is value to further studies which include other genomic regions and traits, and these could guide breeding strategies.

Darwin's Fancy Revised: An Updated Understanding of the Genomic Constitution of Pigeon Breeds

Through its long history of artificial selection, the rock pigeon (Columba livia Gmelin 1789) was forged into a large number of domestic breeds. The incredible amount of phenotypic diversity exhibited in these breeds has long held the fascination of scholars, particularly those interested in biological inheritance and evolution. However, exploiting them as a model system is challenging, as unlike with many other domestic species, few reliable records exist about the origins of, and relationships between, each of the breeds. Therefore, in order to broaden our understanding of the complex evolutionary relationships among pigeon breeds, we generated genome-wide data by performing the genotyping-by-sequencing (GBS) method on close to 200 domestic individuals representing over 60 breeds. We analyzed these GBS data alongside previously published whole-genome sequencing data, and this combined analysis allowed us to conduct the most extensive phylogenetic analysis of the group, including two feral pigeons and one outgroup. We improve previous phylogenies, find considerable population structure across the different breeds, and identify unreported interbreed admixture events. Despite the reduced number of loci relative to whole-genome sequencing, we demonstrate that GBS data provide sufficient analytical power to investigate intertwined evolutionary relationships, such as those that are characteristic of animal domestic breeds. Thus, we argue that future studies should consider sequencing methods akin to the GBS approach as an optimal cost-effective approach for addressing complex phylogenies.

Genomic scans for selective sweeps through haplotype homozygosity and allelic fixation in 14 indigenous sheep breeds from Middle East and South Asia

The performance and productivity of livestock have consistently improved by natural and artificial selection over the centuries. Both these selections are expected to leave patterns on the genome and lead to changes in allele frequencies, but natural selection has played the major role among indigenous populations. Detecting selective sweeps in livestock may assist in understanding the processes involved in domestication, genome evolution and discovery of genomic regions associated with economically important traits. We investigated population genetic diversity and selection signals in this study using SNP genotype data of 14 indigenous sheep breeds from Middle East and South Asia, including six breeds from Iran, namely Iranian Balochi, Afshari, Moghani, Qezel, Zel, and Lori-Bakhtiari, three breeds from Afghanistan, namely Afghan Balochi, Arabi, and Gadik, three breeds from India, namely Indian Garole, Changthangi, and Deccani, and two breeds from Bangladesh, namely Bangladeshi Garole and Bangladesh East. The SNP genotype data were generated by the Illumina OvineSNP50 Genotyping BeadChip array. To detect genetic diversity and population structure, we used principal component analysis (PCA), admixture, phylogenetic analyses, and Runs of homozygosity. We applied four complementary statistical tests, FST (fixation index), xp-EHH (cross-population extended haplotype homozygosity), Rsb (extended haplotype homozygosity between-populations), and FLK (the extension of the Lewontin and Krakauer) to detect selective sweeps. Our results not only confirm the previous studies but also provide a suite of novel candidate genes involved in different traits in sheep. On average, FST, xp-EHH, Rsb, and FLK detected 128, 207, 222, and 252 genomic regions as candidates for selective sweeps, respectively. Furthermore, nine overlapping candidate genes were detected by these four tests, especially TNIK, DOCK1, USH2A, and TYW1B which associate with resistance to diseases and climate adaptation. Knowledge of candidate genomic regions in sheep populations may facilitate the identification and potential exploitation of the underlying genes in sheep breeding.

Genome-wide comparative analyses reveal selection signatures underlying adaptation and production in Tibetan and Poll Dorset sheep

The identification of genome-wide selection signatures can provide insights on the mechanisms of natural and/or artificial selection and uncover genes related to biological functions and/or phenotypes. Tibetan sheep are an important livestock in Tibet, providing meat and wool for Tibetans who are renown for breeding livestock that adapt well to high altitudes. Using whole-genome sequences with an effective sequencing depth of 5×, we investigated the genomic diversity and structure and, identified selection signatures of White Tibetan, Oula and Poll Dorset sheep. We obtained 30,163,679 Single Nucleotide Polymorphisms (SNPs) and 5,388,372 indels benchmarked against the ovine Oar_v4.0 genome assembly. Next, using F_ST, ZHp and XP-EHH approaches, we identified selection signatures spanning a set of candidate genes, including HIF1A, CAPN3, PRKAA1, RXFP2, TRHR and HOXA10 that are associated with pathways and GO categories putatively related to hypoxia responses, meat traits and disease resistance. Candidate genes and GO terms associated with coat color were also identified. Finally, quantification of blood physiological parameters, revealed higher levels of mean corpuscular hemoglobin measurement and mean corpuscular hemoglobin concentration in Tibetan sheep compared with Poll Dorset, suggesting a greater oxygen-carrying capacity in the Tibetan sheep and thus better adaptation to high-altitude hypoxia. In conclusion, this study provides a greater understanding of genome diversity and variations associated with adaptive and production traits in sheep.

An atlas of CNV maps in cattle, goat and sheep

Copy number variation (CNV) is the most prevalent type of genetic structural variation that has been recognized as an important source of phenotypic variation in humans, animals and plants. However, the mechanisms underlying the evolution of CNVs and their function in natural or artificial selection remain unknown. Here, we generated CNV region (CNVR) datasets which were diverged or shared among cattle, goat, and sheep, including 886 individuals from 171 diverse populations. Using 9 environmental factors for genome-wide association study (GWAS), we identified a series of candidate CNVRs, including genes relating to immunity, tick resistance, multi-drug resistance, and muscle development. The number of CNVRs shared between species is significantly higher than expected (P<0.00001), and these CNVRs may be more persist than the single nucleotide polymorphisms (SNPs) shared between species. We also identified genomic regions under long-term balancing selection and uncovered the potential diversity of the selected CNVRs close to the important functional genes. This study provides the evidence that balancing selection might be more common in mammals than previously considered, and might play an important role in the daily activities of these ruminant species.

A selection signatures study among Middle Eastern and European sheep breeds

Selection, both natural and artificial, leaves patterns on the genome during domestication of animals and leads to changes in allele frequencies among populations. Detecting genomic regions influenced by selection in livestock may assist in understanding the processes involved in genome evolution and discovering genomic regions related to traits of economic and ecological interests. In the current study, genetic diversity analyses were conducted on 34,206 quality-filtered SNP positions from 450 individuals in 15 sheep breeds, including six indigenous breeds from the Middle East, namely Iranian Balouchi, Afshari, Moghani, Qezel, Karakas and Norduz, and nine breeds from Europe, namely East Friesian Sheep, Ile de France, Mourerous, Romane, Swiss Mirror, Spaelsau, Suffolk, Comisana and Engadine Red Sheep. The SNP genotype data generated by the Illumina OvineSNP50 Genotyping BeadChip array were used in this analysis. We applied two complementary statistical analyses, F_ST (fixation index) and xp-EHH (cross-population extended haplotype homozygosity), to detect selection signatures in Middle Eastern and European sheep populations. F_ST and xp-EHH detected 629 and 256 genes indicating signatures of selection, respectively. Genomic regions identified using F_ST and xp-EHH contained the CIDEA, HHATL, MGST1, FADS1, RTL1 and DGKG genes, which were reported earlier to influence a number of economic traits. Both F_ST and xp-EHH approaches identified 60 shared genes as the signatures of selection, including four candidate genes (NT5E, ADA2, C8A and C8B) that were enriched for two significant Gene Ontology (GO) terms associated with the adenosine metabolic procedure. Knowledge about the candidate genomic regions under selective pressure in sheep breeds may facilitate identification of the underlying genes and enhance our understanding on these genes role in local adaptation.

Runs of homozygosity analysis of South African sheep breeds from various production systems investigated using OvineSNP50k data

Background

Population history, production system and within-breed selection pressure impacts the genome architecture resulting in reduced genetic diversity and increased frequency of runs of homozygosity islands. This study tested the hypothesis that production systems geared towards specific traits of importance or natural or artificial selection pressures influenced the occurrence and distribution of runs of homozygosity (ROH) in the South African sheep population. The Illumina OvineSNP50 BeadChip was used to genotype 400 sheep belonging to 13 breeds from South Africa representing mutton, pelt and mutton and wool dual-purpose breeds, including indigenous non-descript breeds that are reared by smallholder farmers. To get more insight into the autozygosity and distribution of ROH islands of South African breeds relative to global populations, 623 genotypes of sheep from worldwide populations were included in the analysis. Runs of homozygosity were computed at cut-offs of 1–6 Mb, 6–12 Mb, 12–24 Mb, 24–48 Mb and > 48 Mb, using the R package detectRUNS. The Golden Helix SVS program was used to investigate the ROH islands.

Results

A total of 121,399 ROH with mean number of ROH per animal per breed ranging from 800 (African White Dorper) to 15,097 (Australian Poll Dorset) were obtained. Analysis of the distribution of ROH according to their size showed that, for all breeds, the majority of the detected ROH were in the short (1–6 Mb) category (88.2%). Most animals had no ROH > 48 Mb. Of the South African breeds, the Nguni and the Blackhead Persian displayed high ROH based inbreeding (F_ROH) of 0.31 ± 0.05 and 0.31 ± 0.04, respectively. Highest incidence of common runs per SNP across breeds was observed on chromosome 10 with over 250 incidences of common ROHs. Mean proportion of SNPs per breed per ROH island ranged from 0.02 ± 0.15 (island ROH224 on chromosome 23) to 0.13 ± 0.29 (island ROH175 on chromosome 15). Seventeen (17) of the islands had SNPs observed in single populations (unique ROH islands). The MacArthur Merino (MCM) population had five unique ROH islands followed by Blackhead Persian and Nguni with three each whilst the South African Mutton Merino, SA Merino, White Vital Swakara, Karakul, Dorset Horn and Chinese Merino each had one unique ROH island. Genes within ROH islands were associated with predominantly metabolic and immune response traits and predomestic selection for traits such as presence or absence of horns.

Conclusions

Overall, the frequency and patterns of distribution of ROH observed in this study corresponds to the breed history and implied selection pressures exposed to the sheep populations under study.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07314-2.

Association of TMEM8B and SPAG8 with Mature Weight in Sheep

Simple Summary

Domestication and the subsequent selection of animals for either economic or morphological features can impact the legacy genome of a population in myriad ways. In sheep, the rs426272889 single nucleotide polymorphism (SNP) was identified as the peak of a signature of selection. We examined phenotypic data and identified associations for the Transmembrane protein 8B (TMEM8B) rs426272889 SNP and its genetically linked Sperm-associated antigen 8 (SPAG8) rs160159557 SNP with ewe mature weight in four sheep breeds. These data provided the first production-relevant phenotypes, as well as the first organism-level (as opposed to cellular or tumor-derived) phenotypes, associated with TMEM8B, and in so doing, improved the annotation of this gene and genomic region by adding body weight implications. Once validated, these data can be applied in genetic or genomic selection aiming to achieve desired mature body weight.

Abstract

Signature of selection studies have identified many genomic regions with known functional importance and some without verified functional roles. Multiple studies have identified Transmembrane protein 8B (TMEM8B)rs426272889 as having been recently under extreme selection pressure in domesticated sheep, but no study has provided sheep phenotypic data clarifying a reason for extreme selection. We tested rs426272889 for production trait association in 770 U.S. Rambouillet, Targhee, Polypay, and Suffolk sheep. TMEM8Brs426272889 was associated with mature weight at 3 and 4 years (p < 0.05). This suggested selection for sheep growth and body size might explain the historical extreme selection pressure in this genomic region. We also tested Sperm-associated antigen 8 (SPAG8) rs160159557 encoding a G493C substitution. While this variant was associated with mature weights at ages 3 and 4, it was not as strongly associated as TMEM8Brs426272889. Transmembrane protein 8B has little functional information except as an inhibitor of cancer cell proliferation. To our knowledge, this is the first study linking TMEM8B to whole organism growth and body size under standard conditions. Additional work will be necessary to identify the underlying functional variant(s). Once identified, such variants could be used to improve sheep production through selective breeding.

Donkey genomes provide new insights into domestication and selection for coat color

Current knowledge about the evolutionary history of donkeys is still incomplete due to the lack of archeological and whole-genome diversity data. To fill this gap, we have de novo assembled a chromosome-level reference genome of one male Dezhou donkey and analyzed the genomes of 126 domestic donkeys and seven wild asses. Population genomics analyses indicate that donkeys were domesticated in Africa and conclusively show reduced levels of Y chromosome variability and discordant paternal and maternal histories, possibly reflecting the consequences of reproductive management. We also investigate the genetic basis of coat color. While wild asses show diluted gray pigmentation (Dun phenotype), domestic donkeys display non-diluted black or chestnut coat colors (non-Dun) that were probably established during domestication. Here, we show that the non-Dun phenotype is caused by a 1 bp deletion downstream of the TBX3 gene, which decreases the expression of this gene and its inhibitory effect on pigment deposition.

A new donkey reference genome and comparisons with wild asses yields insights into the evolutionary history of donkey domestication and identifies a genetic variant that results in the non-Dun coat colours of domestic donkeys.

Gaze Following in Ungulates: Domesticated and Non-domesticated Species Follow the Gaze of Both Humans and Conspecifics in an Experimental Context

Gaze following is the ability to use others' gaze to obtain information about the environment (e.g., food location, predators, and social interactions). As such, it may be highly adaptive in a variety of socio-ecological contexts, and thus be widespread across animal taxa. To date, gaze following has been mostly studied in primates, and partially in birds, but little is known on the gaze following abilities of other taxa and, especially, on the evolutionary pressures that led to their emergence. In this study, we used an experimental approach to test gaze following skills in a still understudied taxon, ungulates. Across four species (i.e., domestic goats and lamas, and non-domestic guanacos and mouflons), we assessed the individual ability to spontaneously follow the gaze of both conspecifics and human experimenters in different conditions. In line with our predictions, species followed the model's gaze both with human and conspecific models, but more likely with the latter. Except for guanacos, all species showed gaze following significantly more in the experimental conditions (than in the control ones). Despite the relative low number of study subjects, our study provides the first experimental evidence of gaze following skills in non-domesticated ungulates, and contributes to understanding how gaze following skills are distributed in another taxon-an essential endeavor to identify the evolutionary pressures leading to the emergence of gaze following skills across taxa.

Evolutionary history and current distribution of the West Mediterranean lineage of Brucella melitensis in Italy

Ovine and caprine brucellosis, caused by Brucella melitensis, is one of the world’s most widespread zoonoses and is a major cause of economic losses in domestic ruminant production. In Italy, the disease remains endemic in several southern provinces, despite an ongoing brucellosis eradication programme. In this study, we used whole-genome sequencing to detail the genetic diversity of circulating strains, and to examine the origins of the predominant sub-lineages of B. melitensis in Italy. We reconstructed a global phylogeny of B. melitensis, strengthened by 339 new whole-genome sequences, from Italian isolates collected from 2011 to 2018 as part of a national livestock surveillance programme. All Italian strains belonged to the West Mediterranean lineage, which further divided into two major clades that diverged roughly between the 5th and 7th centuries. We observed that Sicily serves as a brucellosis burden hotspot, giving rise to several distinct sub-lineages. More than 20 putative outbreak clusters of ovine and caprine brucellosis were identified, several of which persisted over the 8 year survey period despite an aggressive brucellosis eradication campaign. While the outbreaks in Central and Northern Italy were generally associated with introductions of single clones of B. melitensis and their subsequent dissemination within neighbouring territories, we observed weak geographical segregation of genotypes in the southern regions. Biovar determination, recommended in routine analysis of all Brucella strains by the World Organisation for Animal Health (OIE), could not discriminate among the four main global clades. This demonstrates a need for updating the guidelines used for monitoring B. melitensis transmission and spread, both at the national and international level, and to include whole-genome-based typing as the principal method for identification and tracing of brucellosis outbreaks.

A Hu sheep genome with the first ovine Y chromosome reveal introgression history after sheep domestication

The Y chromosome plays key roles in male fertility and reflects the evolutionary history of paternal lineages. Here, we present a de novo genome assembly of the Hu sheep with the first draft assembly of ovine Y chromosome (oMSY), using nanopore sequencing and Hi-C technologies. The oMSY that we generated spans 10.6 Mb from which 775 Y-SNPs were identified by applying a large panel of whole genome sequences from worldwide sheep and wild Iranian mouflons. Three major paternal lineages (HY1a, HY1b and HY2) were defined across domestic sheep, of which HY2 was newly detected. Surprisingly, HY2 forms a monophyletic clade with the Iranian mouflons and is highly divergent from both HY1a and HY1b. Demographic analysis of Y chromosomes, mitochondrial and nuclear genomes confirmed that HY2 and the maternal counterpart of lineage C represented a distinct wild mouflon population in Iran that diverge from the direct ancestor of domestic sheep, the wild mouflons in Southeastern Anatolia. Our results suggest that wild Iranian mouflons had introgressed into domestic sheep and thereby introduced this Iranian mouflon specific lineage carrying HY2 to both East Asian and Africa sheep populations.

Whole genome analysis of water buffalo and global cattle breeds highlights convergent signatures of domestication

More people globally depend on the water buffalo than any other domesticated species, and as the most closely related domesticated species to cattle they can provide important insights into the shared evolutionary basis of domestication. Here, we sequence the genomes of 79 water buffalo across seven breeds and compare patterns of between breed selective sweeps with those seen for 294 cattle genomes representing 13 global breeds. The genomic regions under selection between cattle breeds significantly overlap regions linked to stature in human genetic studies, with a disproportionate number of these loci also shown to be under selection between water buffalo breeds. Investigation of potential functional variants in the water buffalo genome identifies a rare example of convergent domestication down to the same mutation having independently occurred and been selected for across domesticated species. Cross-species comparisons of recent selective sweeps can consequently help identify and refine important loci linked to domestication.

The comparative genomics of domesticated lineages can yield insights into the signatures of artificial selection. This study sequences 79 water buffalo genomes from 7 breeds and reveals examples of convergent domestication at the genetic level between water buffalo and cattle.

The growth of domestic goats and sheep: A meta study with Bertalanffy-Pütter models

•

We studied 122 mass-at-age data for goats, sheep, and wildlife.

•

We used Bertalanffy-Pütter models to obtain better fitting growth curves.

•

The Brody model had an acceptable fit to 70% of the data.

•

For 39% of data the best-fitting BP-model had a discernable inflection-point.

•

For these models, maximal weight gain per day was 55% higher than natal weight gain.

Growth literature often uses the Brody, Gompertz, Verhulst, and von Bertalanffy models. Is there a rationale for the preference of these classical named models? The versatile five-parameter Bertalanffy-Pütter (BP) model generalizes these models. We revisited peer-reviewed publications from the years 1970–2019 that fitted growth models to together 122 mass-at-age data of sheep and goats from 19 countries and studied the best-fit BP-models using the least-squares method. None of the named models was ever best-fitting. However, for 70% of the data a single non-sigmoidal model had an acceptable fit (normalized root mean squared error 〈 5% and F-ratio test 〉 5% in comparison to the best-fit): the Brody model. The inherently non-sigmoidal character was further underlined, as there were only 39% of the data, where the best-fitting BP-model had a discernible inflection point. For these data, conclusions of biological interest could be drawn from the sigmoidal best-fit BP-models: the maximal weight gain per day was about 55% higher than the natal weight gain per day.

863 genomes reveal the origin and domestication of chicken

Despite the substantial role that chickens have played in human societies across the world, both the geographic and temporal origins of their domestication remain controversial. To address this issue, we analyzed 863 genomes from a worldwide sampling of chickens and representatives of all four species of wild jungle fowl and each of the five subspecies of red jungle fowl (RJF). Our study suggests that domestic chickens were initially derived from the RJF subspecies Gallus gallus spadiceus whose present-day distribution is predominantly in southwestern China, northern Thailand and Myanmar. Following their domestication, chickens were translocated across Southeast and South Asia where they interbred locally with both RJF subspecies and other jungle fowl species. In addition, our results show that the White Leghorn chicken breed possesses a mosaic of divergent ancestries inherited from other subspecies of RJF. Despite the strong episodic gene flow from geographically divergent lineages of jungle fowls, our analyses show that domestic chickens undergo genetic adaptations that underlie their unique behavioral, morphological and reproductive traits. Our study provides novel insights into the evolutionary history of domestic chickens and a valuable resource to facilitate ongoing genetic and functional investigations of the world's most numerous domestic animal.

Genetic diversity in Iranian indigenous sheep vis-à-vis selected exogenous sheep breeds and wild mouflon

The heterogeneity of climate and different agro-ecological conditions in Iran have resulted in development of 27 indigenous sheep breeds. Wild Asiatic mouflon (Ovis orientalis) is believed to be the ancestor of Iranian sheep. Evaluation of genetic diversity and population structure within and among domestic breeds has important implications for animal breeding programs and genetic resources management. Based on 50K SNP genotype data, we studied the genetic diversity of five indigenous Iranian sheep breeds: Afshari (n = 37), Moghani (n = 34), Qezel (n = 35), Zel (n = 46) and Lori-Bakhtiari (n = 46), and Asiatic mouflon (n = 8) sampled from Iran. Furthermore, genetic diversity and the breed admixture of Iranian sheep were assessed on a larger geographic scale using a reference panel comprising: three indigenous Afghan breeds - Arabi (n = 15), Balouchi (n = 15) and Gadik (n = 15); three indigenous breeds from Turkey and Cyprus - Cyprus Fat Tail (n = 30), Karakas (n = 18) and Norduz (n = 20); and three commercial European breeds - Suffolk (n = 19), Comisana (n = 24) and Engadine Red Sheep (n = 24). The results revealed that the investigated breeds are divided into five genetically distinct clusters according to their geographic origin. Afshari was closest to the local mouflon population and showed signs of mouflon admixture. Qezel was identified as a hybrid sheep breed. Much evidence supported the Afghan breeds being identical. Inbreeding values, which were estimated based on ROHs, were highest for Suffolk (F_ROH  = 0.0544) and lowest for Balouchi (F_ROH  = 0.0078). In conclusion, analysis of selected breeds from neighboring countries along with Asiatic mouflon gave a deeper insight into the evolutionary history and origin of Iranian sheep with important implications for future breed management.

Genomic, Transcriptomic and Epigenomic Tools to Study the Domestication of Plants and Animals: A Field Guide for Beginners

In the last decade, genomics and the related fields of transcriptomics and epigenomics have revolutionized the study of the domestication process in plants and animals, leading to new discoveries and new unresolved questions. Given that some domesticated taxa have been more studied than others, the extent of genomic data can range from vast to nonexistent, depending on the domesticated taxon of interest. This review is meant as a rough guide for students and academics that want to start a domestication research project using modern genomic tools, as well as for researchers already conducting domestication studies that are interested in following a genomic approach and looking for alternate strategies (cheaper or more efficient) and future directions. We summarize the theoretical and technical background needed to carry out domestication genomics, starting from the acquisition of a reference genome and genome assembly, to the sampling design for population genomics, paleogenomics, transcriptomics, epigenomics and experimental validation of domestication-related genes. We also describe some examples of the aforementioned approaches and the relevant discoveries they made to understand the domestication of the studied taxa.

Ancient Genomes Reveal the Evolutionary History and Origin of Cashmere-Producing Goats in China

Goats are one of the most widespread farmed animals across the world; however, their migration route to East Asia and local evolutionary history remain poorly understood. Here, we sequenced 27 ancient Chinese goat genomes dating from the Late Neolithic period to the Iron Age. We found close genetic affinities between ancient and modern Chinese goats, demonstrating their genetic continuity. We found that Chinese goats originated from the eastern regions around the Fertile Crescent, and we estimated that the ancestors of Chinese goats diverged from this population in the Chalcolithic period. Modern Chinese goats were divided into a northern and a southern group, coinciding with the most prominent climatic division in China, and two genes related to hair follicle development, FGF5 and EDA2R, were highly divergent between these populations. We identified a likely causal de novo deletion near FGF5 in northern Chinese goats that increased to high frequency over time, whereas EDA2R harbored standing variation dating to the Neolithic. Our findings add to our understanding of the genetic composition and local evolutionary process of Chinese goats.

Microevolutionary Dynamics of Chicken Genomes under Divergent Selection for Adiposity

Decades of artificial selection have significantly improved performance and efficiency of animal production systems. However, little is known about the microevolution of genomes due to intensive breeding. Using whole-genome sequencing, we document dynamic changes of chicken genomes under divergent selection on adiposity over 19 generations. Directional selection reduced within-line but increased between-line genomic differences. We observed that artificial selection tended to result in recruitment of preexisting variations of genes related to adipose tissue growth. In addition, novel mutations contributed to divergence of phenotypes under selection but contributed significantly less than preexisting genomic variants. Integration of 15 generations genome sequencing, genome-wide association study, and multi-omics data further identified that genes involved in signaling pathways important to adipogenesis, such as autophagy and lysosome (URI1, MBL2), neural system (CHAT), and endocrine (PCSK1) pathways, were under strong selection. Our study provides insights into the microevolutionary dynamics of domestic animal genomes under artificial selection.

Assessment of genetic diversity and differentiation among four indigenous Turkish sheep breeds using microsatellites

Conservation and breeding programmes of livestock species depend on determination of genetic diversity. Today in livestock species, microsatellite markers are commonly used to reveal population structure and genetic diversity in both breeds and varieties. In this study, population structure, genetic diversity, and differentiation among four native Turkish sheep breeds including Güney Karaman, Kangal, Norduz, and Karakas were assessed by using 21 microsatellite loci. By genotyping 120 individuals belonging to four sheep breeds, a total of 275 different alleles, 37 of which were private alleles, were observed across all loci. The mean number of alleles per breed ranged from 7.28 (Güney Karaman) to 8.09 (Karakas), while allelic richness ranged from 7.22 (Güney Karaman) to 7.87 (Karakas). Mean observed heterozygosity varied from 0.60 (Kangal) to 0.66 (Norduz and Karakas). The lowest pairwiseF ST value (0.084) was between Kangal and Karakas populations, while the highest pairwiseF ST value (0.142) was between Norduz and Karakas populations. Polymorphic information content (PIC) values, ranging from 0.71 (ETH10) to 0.91 (OarFCB304), were highly polymorphic (PIC  >  0.5) and informative in studied populations. In the present study, the results of phylogenetic analysis were of importance, since all studied populations have been accepted as Akkaraman varieties till today. However, factorial correspondence and structure analysis, pairwiseF ST values, and an unweighted pair group method with arithmetic mean analysis (UPGMA) dendrogram revealed that Güney Karaman and Norduz populations have became genetically different from the Akkaraman breed due being raised in different parts of Turkey under different climatic conditions together with their breeding practices. Therefore, we recommend that more comprehensive molecular studies should be conducted to clarify genetic differentiation of Akkaraman sheep varieties.

Whole-genome resequencing of wild and domestic sheep identifies genes associated with morphological and agronomic traits

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.

Insight into the Possible Formation Mechanism of the Intersex Phenotype of Lanzhou Fat-Tailed Sheep Using Whole-Genome Resequencing

Simple Summary

Individuals with hermaphroditism are a serious hazard to animal husbandry and production due to their abnormal fertility. The molecular mechanism of sheep intersex formation was unclear. This study was the first to locate the homologous sequence of the goat polled intersex syndrome (PIS) region in sheep and found that the intersex traits of Lanzhou fat-tailed sheep were not caused by the lack of this region. By detecting the selective sweep regions, vital genes associated with androgen biosynthesis and the follicle stimulating hormone response entry were found, including steroid 5 alpha-reductase 2 (SRD5A2), and pro-apoptotic WT1 regulator (PAWR). Additionally, the copy number variations of the four regions on chr9, chr1, chr4, and chr16 may affect the expression of the gonadal development genes, zinc finger protein, FOG family member 2 (ZFPM2), LIM homeobox 8 (LHX8), inner mitochondrial membrane peptidase subunit 2 (IMMP2L) and slit guidance ligand 3 (SLIT3), respectively, and further affect the formation of intersex traits.

Abstract

Intersex, also known as hermaphroditism, is a serious hazard to animal husbandry and production. The mechanism of ovine intersex formation is not clear. Therefore, genome-wide resequencing on the only two intersex and two normal Lanzhou fat-tailed (LFT) sheep, an excellent but endangered Chinese indigenous sheep breed, was performed. Herein, the deletion of homologous sequences of the goat polled intersex syndrome (PIS) region (8787 bp, 247747059–247755846) on chromosome 1 of the LFT sheep was not the cause of the ovine intersex trait. By detecting the selective sweep regions, we found that the genes related to androgen biosynthesis and follicle stimulating hormone response items, such as steroid 5 alpha-reductase 2 (SRD5A2), steroid 5 alpha-reductase 3 (SRD5A3), and pro-apoptotic WT1 regulator (PAWR), may be involved in the formation of intersex traits. Furthermore, the copy number variations of the four regions, chr9: 71660801–71662800, chr1: 50776001–50778000, chr4: 58119201–58121600, and chr16: 778801–780800, may affect the expression of the zinc finger protein, FOG family member 2 (ZFPM2), LIM homeobox 8 (LHX8), inner mitochondrial membrane peptidase subunit 2 (IMMP2L) and slit guidance ligand 3 (SLIT3) genes, respectively, which contribute to the appearance of intersex traits. These results may supply a theoretical basis for the timely detection and elimination of intersex individuals in sheep, which could accelerate the healthy development of animal husbandry.

Common Gene Modules Identified for Chicken Adiposity by Network Construction and Comparison

Excessive fat deposition can cause chicken health problem, and affect production efficiency by causing great economic losses to the industry. However, the molecular underpinnings of the complex adiposity trait remain elusive. In the current study, we constructed and compared the gene co-expression networks on four transcriptome profiling datasets, from two chicken lines under divergent selection for abdominal fat contents, in an attempt to dissect network compositions underlying adipose tissue growth and development. After functional enrichment analysis, nine network modules important to adipogenesis were discovered to be involved in lipid metabolism, PPAR and insulin signaling pathways, and contained hub genes related to adipogenesis, cell cycle, inflammation, and protein synthesis. Moreover, after additional functional annotation and network module comparisons, common sub-modules of similar functionality for chicken fat deposition were identified for different chicken lines, apart from modules specific to each chicken line. We further validated the lysosome pathway, and found TFEB and its downstream target genes showed similar expression patterns along with chicken preadipocyte differentiation. Our findings could provide novel insights into the genetic basis of complex adiposity traits, as well as human obesity and related metabolic diseases.

On the origin of European sheep as revealed by the diversity of the Balkan breeds and by optimizing population-genetic analysis tools

Background

In the Neolithic, domestic sheep migrated into Europe and subsequently spread in westerly and northwesterly directions. Reconstruction of these migrations and subsequent genetic events requires a more detailed characterization of the current phylogeographic differentiation.

Results

We collected 50 K single nucleotide polymorphism (SNP) profiles of Balkan sheep that are currently found near the major Neolithic point of entry into Europe, and combined these data with published genotypes from southwest-Asian, Mediterranean, central-European and north-European sheep and from Asian and European mouflons. We detected clines, ancestral components and admixture by using variants of common analysis tools: geography-informative supervised principal component analysis (PCA), breed-specific admixture analysis, across-breed \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{4}$$\end{document}f4 profiles and phylogenetic analysis of regional pools of breeds. The regional Balkan sheep populations exhibit considerable genetic overlap, but are clearly distinct from the breeds in surrounding regions. The Asian mouflon did not influence the differentiation of the European domestic sheep and is only distantly related to present-day sheep, including those from Iran where the mouflons were sampled. We demonstrate the occurrence, from southeast to northwest Europe, of a continuously increasing ancestral component of up to 20% contributed by the European mouflon, which is assumed to descend from the original Neolithic domesticates. The overall patterns indicate that the Balkan region and Italy served as post-domestication migration hubs, from which wool sheep reached Spain and north Italy with subsequent migrations northwards. The documented dispersal of Tarentine wool sheep during the Roman period may have been part of this process. Our results also reproduce the documented 18th century admixture of Spanish Merino sheep into several central-European breeds.

Conclusions

Our results contribute to a better understanding of the events that have created the present diversity pattern, which is relevant for the management of the genetic resources represented by the European sheep population.

Ensembl 2019

The Ensembl project (https://www.ensembl.org) makes key genomic data sets available to the entire scientific community without restrictions. Ensembl seeks to be a fundamental resource driving scientific progress by creating, maintaining and updating reference genome annotation and comparative genomics resources. This year we describe our new and expanded gene, variant and comparative annotation capabilities, which led to a 50% increase in the number of vertebrate genomes we support. We have also doubled the number of available human variants and added regulatory regions for many mouse cell types and developmental stages. Our data sets and tools are available via the Ensembl website as well as a through a RESTful webservice, Perl application programming interface and as data files for download.

Understanding the Genetic Domestication History of the Jianchang Duck by Genotyping and Sequencing of Genomic Genes Under Selection

The Jianchang duck is mainly distributed in Southwest China, and has the characteristics of fast growth rate and strong abilities in lipid deposition in the liver. In order to investigate the effects of domestication process on formation of the unique characteristics of Jianchang duck, the whole genome of sixteen individuals and three pooling of Jianchang duck were re-sequenced, and genome data of 70 mallards and 83 domestic ducks from thirteen different places in China were obtained from NCBI. The population stratification and evolution analysis showed gene exchanges existed between the Jianchang and other domestic duck populations, as well as Jianchang ducks and mallards. Genomic comparison between mallards and Jianchang ducks showed genes, including CNTN1, CHRNA9, and SHANK2, which is involved in brain and nerve development, experienced strong positive selection in the process of Jianchang duck domestication. The genomic comparison between Jianchang and domestic duck populations showed that HSD17B12 and ESM1, which affect lipid metabolism, experienced strong positive selection during the domestication process. FST analysis among populations of Jianchang duck with different plumage colors indicated that MITF was related to the phenotype of a white feather, while MC1R was related to the phenotype of hemp feather. Our results provided a base for the domestication process of Jianchang duck and the genomic genes for unique traits.

The origin of domestication genes in goats

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.

Molecular cloning of the goose GnRH gene and identification of GnRH polymorphisms associated with laying traits

1. Egg-laying traits are important economic characteristics in goose production (Anser cygnoides). The gene GnRH, which encodes gonadotropin-releasing hormone, is a strong candidate gene for egg-laying traits in avian species. 2. In this study, a 3520 bp genomic sequence and a 279 bp mRNA sequence for GnRH, which encoded 92 amino acids, were determined. The GnRH DNA sequence contains four exons and three introns, and the DNA and deduced amino acid sequences were highly conserved across mammals (human, macaque, cow, and sheep) and avians (chicken, fulmar and quail). 3. Using a direct sequencing method, 46 single nucleotide polymorphisms (SNPs) were identified in the GnRH genomic sequence that were shared between two Sichuan White goose populations (217 and 208 individuals). Furthermore, 44 haplotypes were constructed using a sliding window approach. Association analysis between the SNPs and haplotypes and egg-laying traits showed that 10 SNPs affected the first egg weight, average egg weight, egg number at 48 weeks and egg number at 64 weeks. 4. These results lay the foundation for further studies of the function of GnRH in geese and provide a theoretical basis for marker-assisted selection of egg-laying traits in the Sichuan white goose population.

Comparing the diversity of the casein genes in the Asian mouflon and domestic sheep

We aimed to determine whether casein variants that are currently segregating in ovine populations existed before the domestication of sheep or, to the contrary, if their emergence is much more recent. To this end, we have retrieved whole-genome sequences from Iranian and domestic sheep from Africa, Europe, South and East Asia and West Asia. Population structure analysis based on 55,352,935 SNPs revealed a clear separation between Iranian mouflons and domestic sheep. Moreover, we also observed a strong genetic differentiation between Iranian mouflons sampled in geographic areas close to Tehran and Tabriz. Based on sequence data, hundreds of SNPs mapping to the casein α_S1 (CSN1S1, 248 SNPs), casein α_S2 (CSN1S2, 268 SNPs), casein ß (CSN2, 146 SNPs) and casein κ (CSN3, 112 SNPs) genes were identified. Approximately 25-63.02% of the casein variation was shared between Iranian mouflons and domestic sheep, and the four domestic sheep populations also shared 44.2-57.4% of the casein polymorphic sites. These findings suggest that an important fraction of the casein variation present in domestic sheep was already segregating in the mouflon prior to its domestication. Genomic studies performed in horses and dogs are consistent with this view, suggesting that much of the diversity that we currently detect in domestic animals comes from standing variation already segregating in their wild ancestors.

Changed Patterns of Genomic Variation Following Recent Domestication: Selection Sweeps in Farmed Atlantic Salmon

The introduction of wild Atlantic salmon into captivity, and their subsequent artificial selection for production traits, has caused phenotypic differences between domesticated fish and their wild counterparts. Identification of regions of the genome underling these changes offers the promise of characterizing the early biological consequences of domestication. In the current study, we sequenced a population of farmed European Atlantic salmon and compared the observed patterns of SNP variation to those found in conspecific wild populations. This identified 139 genomic regions that contained significantly elevated SNP homozygosity in farmed fish when compared to their wild counterparts. The most extreme was adjacent to versican, a gene involved in control of neural crest cell migration. To control for false positive signals, a second and independent dataset of farmed and wild European Atlantic salmon was assessed using the same methodology. A total of 81 outlier regions detected in the first dataset showed significantly reduced homozygosity within the second one, strongly suggesting the genomic regions identified are enriched for true selection sweeps. Examination of the associated genes identified a number previously characterized as targets of selection in other domestic species and that have roles in development, behavior and olfactory system. These include arcvf, sema6, errb4, id2-like, and 6n1-like genes. Finally, we searched for evidence of parallel sweeps using a farmed population of North American origin. This failed to detect a convincing overlap to the putative sweeps present in European populations, suggesting the factors that drive patterns of variation under domestication and early artificial selection were largely independent. This is the first analysis on domestication of aquaculture species exploiting whole-genome sequence data and resulted in the identification of sweeps common to multiple independent populations of farmed European Atlantic salmon.

How human behavior can impact the evolution of genetically-mediated behavior in wild non-human species

Humans intensely modify the ecosystems we inhabit. Many of the impacts that this behavior can have on other species also sharing these spaces are obvious. A prime example is the devastating current extinction crisis. Yet some populations of non-human, non-domesticated species survive or even appear to thrive in heavily disturbed or human-built habitats. Theoretically, this apparent paradox could be facilitated partly by the evolution of genetically-mediated trait adaptations to the impacts of human behavior. At the least, persistence in strongly modified habitats would provide requisite selection pressures for this process to potentially occur in the future. In fact, we have a growing number of well-characterized examples of morphological trait adaptations to human behavior. However, our knowledge of genetically-mediated behavioral adaptations in similar contexts is less well developed. In this review I set up and discuss several evolutionary scenarios by which human behavior might have impacted the evolution of genetically mediated behavior in non-human, non-domestic species and highlight several approaches that could be used in future studies of this process.

Phylogenetic analysis of small ruminant lentiviruses in Germany and Iran suggests their expansion with domestic sheep

Small ruminant lentiviruses (SRLVs) are found in sheep in Germany and Iran. SRLVs have been classified into four genotypes: A-C and E. Genotype A has been subdivided into 20 subtypes. Previous studies suggested that, first, the ancestors of genotype A are those SRLVs found in Turkey, second, the evolution of SRLVs is related to the domestication process, and, third, SRLV infection was first observed in sheep in Iceland and the source of that infection was a flock imported from Germany. This study generated, for the first time, partial SRLV sequence data from German and Iranian sheep, enhancing our knowledge of the genetic and evolutionary relationships of SRLVs, and their associations with the domestication process. Based on 54 SRLV sequences from German and Iranian sheep, our results reveal: (1) SRLV subtypes A4, A5, A11, A16 and A21 (new) are found in German sheep and A22 (new) in Iranian sheep. (2) Genotype A has potentially an additional ancestor (A22), found in Iran, Lebanon and Jordan. (3) Subtype A22 is likely an old version of SRLVs. (4) The transmission routes of some SRLVs are compatible with domestication pathways. (5) This study found no evidence of Icelandic subtype A1 in German sheep.

Adaptive introgression from indicine cattle into white cattle breeds from Central Italy

Cattle domestication occurred at least twice independently and gave rise to the modern taurine and indicine cattle breeds. European cattle diversity is generally dominated by taurine cattle, although elevated levels of indicine ancestry have been recorded in several breeds from southern Europe. Here we use genome-wide high-density SNP genotyping data to investigate the taurine and indicine ancestry in southern European cattle, based on a dataset comprising 508 individuals from 23 cattle breeds of taurine, indicine and mixed ancestry, including three breeds from Central Italy known to exhibit the highest levels of indicine introgression among southern European breeds. Based on local genomic ancestry analyses, we reconstruct taurine and indicine ancestry genome-wide and along chromosomes. We scrutinise local genomic introgression signals and identify genomic regions that have introgressed from indicine into taurine cattle under positive selection, harbouring genes with functions related to body size and feed efficiency. These findings suggest that indicine-derived traits helped enhance Central Italian cattle through adaptive introgression. The identified genes could provide genomic targets for selection for improved cattle performance. Our findings elucidate the key role of adaptive introgression in shaping the phenotypic features of modern cattle, aided by cultural and livestock exchange among historic human societies.

Selection signatures in goats reveal copy number variants underlying breed-defining coat color phenotypes

Domestication and human selection have formed diverse goat breeds with characteristic phenotypes. This process correlated with the fixation of causative genetic variants controlling breed-specific traits within regions of reduced genetic diversity, so called selection signatures or selective sweeps. Using whole genome sequencing of DNA pools (pool-seq) from 20 genetically diverse modern goat breeds and bezoars, we identified 2,239 putative selection signatures. In two Pakistani goat breeds, Pak Angora and Barbari, we found selection signatures in a region harboring KIT, a gene involved in melanoblast development, migration, and survival. The search for candidate causative variants responsible for these selective sweeps revealed two different copy number variants (CNVs) downstream of KIT that were exclusively present in white Pak Angora and white-spotted Barbari goats. Several Swiss goat breeds selected for specific coat colors showed selection signatures at the ASIP locus encoding the agouti signaling protein. Analysis of these selective sweeps revealed four different CNVs associated with the white or tan (A^Wt), Swiss markings (A^sm), badgerface (A^b), and the newly proposed peacock (A^pc) allele. RNA-seq analyses on skin samples from goats with the different CNV alleles suggest that the identified structural variants lead to an altered expression of ASIP between eumelanistic and pheomelanistic body areas. Our study yields novel insights into the genetic control of pigmentation by identifying six functionally relevant CNVs. It illustrates how structural changes of the genome have contributed to phenotypic evolution in domestic goats.

Domestic animals have been selected for hundreds or sometimes even thousands of years for traits that were appreciated by their human owners. This process correlated with the fixation of causative genetic variants controlling breed-specific traits within regions of reduced genetic diversity, so called selection signatures or selective sweeps. We conducted a comprehensive screen for selection signatures in 20 phenotypically and genetically diverse modern goat breeds and identified a total of 2,239 putative selection signatures in our dataset. Follow-up experiments on selection signatures harboring known candidate genes for coat color revealed six different copy number variants (CNVs). Two of these CNVs were located in the 3’-flanking region of KIT and associated with a completely white coat color phenotype in Pak Angora goats and a white-spotted coat color phenotype in Barbari goats, respectively. The other four CNVs were located at the ASIP locus. They were associated with four different types of coat color patterning in seven Swiss goat breeds. Their functional effect is mediated by region-specific quantitative changes in ASIP mRNA expression. Our study illustrates how structural changes of the genome have contributed to phenotypic evolution in domestic goats.

Comparative genome analyses reveal the unique genetic composition and selection signals underlying the phenotypic characteristics of three Chinese domestic goat breeds

BACKGROUND

As one of the important livestock species around the world, goats provide abundant meat, milk, and fiber to fulfill basic human needs. However, the genetic loci that underlie phenotypic variations in domestic goats are largely unknown, particularly for economically important traits. In this study, we sequenced the whole genome of 38 goats from three Chinese breeds (Chengdu Brown, Jintang Black, and Tibetan Cashmere) and downloaded the genome sequence data of 30 goats from five other breeds (four non-Chinese and one Chinese breed) and 21 Bezoar ibexes to investigate the genetic composition and selection signatures of the Chinese goat breeds after domestication.

RESULTS

Based on population structure analysis and FST values (average FST = 0.22), the genetic composition of Chengdu Brown goats differs considerably from that of Bezoar ibexes as a result of geographic isolation. Strikingly, the genes under selection that we identified in Tibetan Cashmere goats were significantly enriched in the categories hair growth and bone and nervous system development, possibly because they are involved in adaptation to high-altitude. In particular, we found a large difference in allele frequency of one novel SNP (c.-253G>A) in the 5'-UTR of FGF5 between Cashmere goats and goat breeds with short hair. The mutation at this site introduces a start codon that results in the occurrence of a premature FGF5 protein and is likely a natural causal variant that is involved in the long hair phenotype of cashmere goats. The haplotype tagged with the AGG-allele in exon 12 of DSG3, which encodes a cell adhesion molecule that is expressed mainly in the skin, was almost fixed in Tibetan Cashmere goats, whereas this locus still segregates in the lowland goat breeds. The pigmentation gene KITLG showed a strong signature of selection in Tibetan Cashmere goats. The genes ASIP and LCORL were identified as being under positive selection in Jintang Black goats.

CONCLUSIONS

After domestication, geographic isolation of some goat breeds has resulted in distinct genetic structures. Furthermore, our work highlights several positively selected genes that likely contributed to breed-related traits in domestic goats.