Gene-based mapping and pathway analysis of metabolic traits in dairy cows

Dissecting loci that underpin the genetic correlations between production, fertility, and urea traits in Australian Holstein cattle

Unfavorable genetic correlations between milk production, fertility, and urea traits have been reported. However, knowledge of the genomic regions associated with these unfavorable correlations is limited. Here, we used the correlation scan method to identify and investigate the regions driving or antagonizing the genetic correlations between production vs. fertility, urea vs. fertility, and urea vs. production traits. Driving regions produce an estimate of correlation that is in the same direction as the global correlation. Antagonizing regions produce an estimate in the opposite direction of the global estimates. Our dataset comprised 6567, 4700, and 12,658 Holstein cattle with records of production traits (milk yield, fat yield, and protein yield), fertility (calving interval) and urea traits (milk urea nitrogen and blood urea nitrogen predicted using milk-mid-infrared spectroscopy), respectively. Several regions across the genome drive the correlations between production, fertility, and urea traits. Antagonizing regions were confined to certain parts of the genome and the genes within these regions were mostly involved in preventing metabolic dysregulation, liver reprogramming, metabolism remodeling, and lipid homeostasis. The driving regions were enriched for QTL related to puberty, milk, and health-related traits. Antagonizing regions were mostly related to muscle development, metabolic body weight, and milk traits. In conclusion, we have identified genomic regions of potential importance for dairy cattle breeding. Future studies could investigate the antagonizing regions as potential genomic regions to break the unfavorable correlations and improve milk production as well as fertility and urea traits.

Nutrigenomic Interventions to Address Metabolic Stress and Related Disorders in Transition Cows

For dairy cattle, the period involving a shift from late pregnancy to early lactation termed transition or periparturient is an excruciating phase. Health-related disorders are likely to happen in this time frame. Timely postpartum and metabolic adjustments to this new physical state demands correct management strategies to fulfill the cow's needs for a successful transition to this phase. Among the management strategies, one of the most researched methods for managing transition-related stress is nutritional supplementation. Dietary components directly or indirectly affect the expression of various genes that are believed to be involved in various stress-related responses during this phase. Nutrigenomics, an interdisciplinary approach that combines nutritional science with omics technologies, opens new avenues for studying the genome's complicated interactions with food. This revolutionary technique emphasizes the importance of food-gene interactions on various physiological and metabolic mechanisms. In animal sciences, nutrigenomics aims to promote the welfare of livestock animals and enhance their commercially important qualities through nutritional interventions. To this end, an increasing volume of research shows that nutritional supplementation can be effectively used to manage the metabolic stress dairy cows undergo during the transition period. These nutritional supplements, including polyunsaturated fatty acids, vitamins, dietary amino acids, and phytochemicals, have been shown to modulate energy homeostasis through different pathways, leading to addressing metabolic issues in transition cows.

Structural equation modeling for unraveling the multivariate genomic architecture of milk proteins in dairy cattle

The aims of this study were to investigate potential functional relationships among milk protein fractions in dairy cattle and to carry out a structural equation model (SEM) GWAS to provide a decomposition of total SNP effects into direct effects and effects mediated by traits that are upstream in a phenotypic network. To achieve these aims, we first fitted a mixed Bayesian multitrait genomic model to infer the genomic correlations among 6 milk nitrogen fractions [4 caseins (CN), namely κ-, β-, α_S1-, and α_S2-CN, and 2 whey proteins, namely β-lactoglobulin (β-LG) and α-lactalbumin (α-LA)], in a population of 989 Italian Brown Swiss cows. Animals were genotyped with the Illumina BovineSNP50 Bead Chip v.2 (Illumina Inc.). A Bayesian network approach using the max-min hill-climbing (MMHC) algorithm was implemented to model the dependencies or independence among traits. Strong and negative genomic correlations were found between β-CN and α_S1-CN (-0.706) and between β-CN and κ-CN (-0.735). The application of the MMHC algorithm revealed that κ-CN and β-CN seemed to directly or indirectly influence all other milk protein fractions. By integrating multitrait model GWAS and SEM-GWAS, we identified a total of 127 significant SNP for κ-CN, 89 SNP for β-CN, 30 SNP for α_S1-CN, and 14 SNP for α_S2-CN (mostly shared among CN and located on Bos taurus autosome 6) and 15 SNP for β-LG (mostly located on Bos taurus autosome 11), whereas no SNP passed the significance threshold for α-LA. For the significant SNP, we assessed and quantified the contribution of direct and indirect paths to total marker effect. Pathway analyses confirmed that common regulatory mechanisms (e.g., energy metabolism and hormonal and neural signals) are involved in the control of milk protein synthesis and metabolism. The information acquired might be leveraged for setting up optimal management and selection strategies aimed at improving milk quality and technological characteristics in dairy cattle.

A Gene-Set Enrichment and Protein-Protein Interaction Network-Based GWAS with Regulatory SNPs Identifies Candidate Genes and Pathways Associated with Carcass Traits in Hanwoo Cattle

Non-synonymous SNPs and protein coding SNPs within the promoter region of genes (regulatory SNPs) might have a significant effect on carcass traits. Imputed sequence level data of 10,215 Hanwoo bulls, annotated and filtered to include only regulatory SNPs (450,062 SNPs), were used in a genome-wide association study (GWAS) to identify loci associated with backfat thickness (BFT), carcass weight (CWT), eye muscle area (EMA), and marbling score (MS). A total of 15, 176, and 1 SNPs were found to be significantly associated (p < 1.11 × 10-7) with BFT, CWT, and EMA, respectively. The significant loci were BTA4 (CWT), BTA6 (CWT), BTA14 (CWT and EMA), and BTA19 (BFT). BayesR estimated that 1.1%~1.9% of the SNPs contributed to more than 0.01% of the phenotypic variance. So, the GWAS was complemented by a gene-set enrichment (GSEA) and protein-protein interaction network (PPIN) analysis in identifying the pathways affecting carcass traits. At p < 0.005 (~2,261 SNPs), 25 GO and 18 KEGG categories, including calcium signaling, cell proliferation, and folate biosynthesis, were found to be enriched through GSEA. The PPIN analysis showed enrichment for 81 candidate genes involved in various pathways, including the PI3K-AKT, calcium, and FoxO signaling pathways. Our finding provides insight into the effects of regulatory SNPs on carcass traits.

Longitudinal Phenotypes Improve Genotype Association for Hyperketonemia in Dairy Cattle

Simple Summary

Dairy cows have differing success in supporting their physiological functions while in energy deficit right after calving. Identification of genomic regions associated with different concentrations of non–esterified fatty acids and β–hydroxybutyrate in early postpartum Holstein cows provide insight into an animal’s genetic susceptibility to these conditions. Longitudinal phenotypes may provide a different perspective than cross-sectional phenotype variation and their association with genotypes in the study of complex metabolic diseases in dairy cows. This might allow us to reinforce preventative measures that decrease the incidence of hyperketonemia and improve genetic selection criteria.

Abstract

The objective of our study was to identify genomic regions associated with varying concentrations of non-esterified fatty acid (NEFA), β-hydroxybutyrate (BHB), and the development of hyperketonemia (HYK) in longitudinally sampled Holstein dairy cows. Our study population consisted of 147 multiparous cows intensively characterized by serial NEFA and BHB concentrations. To identify individuals with contrasting combinations in longitudinal BHB and NEFA concentrations, phenotypes were established using incremental area under the curve (AUC) and categorized as follows: Group (1) high NEFA and high BHB, group (2) low NEFA and high BHB), group (3) low NEFA and low BHB, and group (4) high NEFA and low BHB. Cows were genotyped on the Illumina Bovine High-density (777 K) beadchip. Genome-wide association studies using mixed linear models with the least-related animals were performed to establish a genetic association with HYK, BHB-AUC, NEFA-AUC, and the comparisons of the 4 AUC phenotypic groups using Golden Helix software. Nine single-nucleotide polymorphisms were associated with high longitudinal concentrations of BHB and further investigated. Five candidate genes related to energy metabolism and homeostasis were identified. These results provide biological insight and help identify susceptible animals thus improving genetic selection criteria thereby decreasing the incidence of HYK.

On the genomic regions associated with milk lactose in Fleckvieh cattle

Lactose is a sugar uniquely found in mammals' milk and it is the major milk solid in bovines. Lactose yield (LY, kg/d) is responsible for milk volume, whereas lactose percentage (LP) is thought to be more related to epithelial integrity and thus to udder health. There is a paucity of studies that have investigated lactose at the genomic level in dairy cows. This paper aimed to improve our knowledge on LP and LY, providing new insights into the significant genomic regions affecting these traits. A genome-wide association study for LP and LY was carried out in Fleckvieh cattle by using bulls' deregressed estimated breeding values of first lactation as pseudo-phenotypes. Heritabilities of first-lactation test-day LP and LY estimated using linear animal models were 0.38 and 0.25, respectively. A total of 2,854 bulls genotyped with a 54K SNP chip were available for the genome-wide association study; a linear mixed model approach was adopted for the analysis. The significant SNP of LP were scattered across the whole genome, with signals on chromosomes 1, 2, 3, 7, 12, 16, 18, 19, 20, 28, and 29; the top 4 significant SNP explained 4.90% of the LP genetic variance. The signals were mostly in regions or genes with involvement in molecular intra- or extracellular transport; for example, CDH5, RASGEF1C, ABCA6, and SLC35F3. A significant region within chromosome 20 was previously shown to affect mastitis or somatic cell score in cattle. As regards LY, the significant SNP were concentrated in fewer regions (chromosomes 6 and 14), related to mastitis/somatic cell score, immune response, and transport mechanisms. The 5 most significant SNP for LY explained 8.45% of genetic variance and more than one-quarter of this value has to be attributed to the variant within ADGRB1. Significant peaks in target regions remained even after adjustment for the 2 most significant variants previously detected on BTA6 and BTA14. The present study is a prelude for deeper investigations into the biological role of lactose for milk secretion and volume determination, stressing the connection with genes regulating intra- or extracellular trafficking and immune and inflammatory responses in dairy cows. Also, these results improve the knowledge on the relationship between lactose and udder health; they support the idea that LP and its derived traits are potential candidates as indicators of udder health in breeding programs aimed to enhance cows' resistance to mastitis.

Genome-wide association analysis for β-hydroxybutyrate concentration in Milk in Holstein dairy cattle

BACKGROUND

Ketosis in dairy cattle has been shown to cause a high morbidity in the farm and substantial financial losses to dairy farmers. Ketosis symptoms, however, are difficult to identify, therefore, the amount of ketone bodies (mainly β-hydroxybutyric acid, BHB) is used as an indicator of subclinical ketosis in cows. It has also been shown that milk BHB concentrations have a strong correlation with ketosis in dairy cattle. Mid-infrared spectroscopy (MIR) has recently became a fast, cheap and high-throughput method for analyzing milk components. The aim of this study was to perform a genome-wide association study (GWAS) on the MIR-predicted milk BHB to identify genomic regions, genes and pathways potentially affecting subclinical ketosis in North American Holstein dairy cattle.

RESULTS

Several significant regions were identified associated with MIR-predicted milk BHB concentrations (indicator of subclinical ketosis) in the first lactation (SCK1) and second and later lactations (SCK2) in Holstein dairy cows. The strongest association was located on BTA6 for SCK1 and BTA14 on SCK2. Several SNPs on BTA6 were identified in regions and variants reported previously to be associated with susceptibility to ketosis and clinical mastitis in Jersey and Holstein dairy cattle, respectively. One highly significant SNP on BTA14 was found within the DGAT1 gene with known functions on fat metabolism and inflammatory response in dairy cattle. A region on BTA6 and three SNPs on BTA20 were found to overlap between SCK1 and SCK2. However, a novel region on BTA20 (55-63 Mb) for SCK2 was also identified, which was not reported in previous association studies. Enrichment analysis of the list of candidate genes within the identified regions for MIR-predicted milk BHB concentrations yielded molecular functions and biological processes that may be involved in the inflammatory response and lipid metabolism in dairy cattle.

CONCLUSIONS

The results of this study confirmed several SNPs and genes identified in previous studies as associated with ketosis susceptibility and immune response, and also found a novel region that can be used for further analysis to identify causal variations and key regulatory genes that affect clinical/ subclinical ketosis.

Genome-wide association analysis and gene set enrichment analysis with SNP data identify genes associated with 305-day milk yield in Holstein dairy cows

Milk production traits, such as 305-day milk yield (305MY), have been under direct selection to improve production in dairy cows. Over the past 50 years, the average milk yield has nearly doubled, and over 56% of the increase is attributable to genetic improvement. As such, additional improvements in milk yield are still possible as new loci are identified. The objectives of this study were to detect SNPs and gene sets associated with 305MY in order to identify new candidate genes contributing to variation in milk production. A population of 781 primiparous Holstein cows from six central Washington dairies with records of 305MY and energy corrected milk were used to perform a genome-wide association analysis (GWAA) using the Illumina BovineHD BeadChip (777 962 SNPs) to identify QTL associated with 305MY (P < 1.0 × 10^-5 ). A gene set enrichment analysis with SNP data (GSEA-SNP) was performed to identify gene sets (normalized enrichment score > 3.0) and leading edge genes (LEGs) influencing 305MY. The GWAA identified three QTL comprising 34 SNPs and 30 positional candidate genes. In the GSEA-SNP, five gene sets with 58 unique and 24 shared LEGs contributed to 305MY. Identification of QTL and LEGs associated with 305MY can provide additional targets for genomic selection to continue to improve 305MY in dairy cattle.

Invited review: Metabolic challenges and adaptation during different functional stages of the mammary gland in dairy cows: Perspectives for sustainable milk production

Milk production of dairy cows has increased markedly during recent decades and continues to increase further. The evolutionarily conserved direction of nutrients to the mammary gland immediately after calving provided the basis for successful selective breeding toward higher performance. Considerable variation in adaptive responses toward energy and nutrient shortages exists; however, this variation in adaptability recently gained interest for identifying more metabolically robust dairy cows. Metabolic challenges during periods of high milk production considerably affect the immune system, reproductive performance, and product quality as well as animal welfare. Moreover, growing consumer concerns need to be taken into consideration because the public perception of industrialized dairy cow farming, the high dependency on feed sources suitable for human nutrition, and the apparently abundant use of antibiotics may affect the sales of dairy products. Breeding for high yield continues, but the metabolic challenges increasingly come close to the adaptational limits of meeting the mammary gland's requirements. The aim of the present review is to elucidate metabolic challenges and adaptational limitations at different functional stages of the mammary gland in dairy cows. From the challenges and adaptational limitations, we derive perspectives for sustainable milk production. Based on previous research, we highlight the importance of metabolic plasticity in adaptation mechanisms at different functional stages of the mammary gland. Metabolic adaptation and plasticity change among developing, nonlactating, remodeling, and lactational stages of the mammary gland. A higher metabolic plasticity in early-lactating dairy cows could be indicative of resilience, and a high performance level without an extraordinary occurrence of health disorders can be achieved.

Genome-wide association studies and genomic prediction of breeding values for calving performance and body conformation traits in Holstein cattle

Background

Our aim was to identify genomic regions via genome-wide association studies (GWAS) to improve the predictability of genetic merit in Holsteins for 10 calving and 28 body conformation traits. Animals were genotyped using the Illumina Bovine 50 K BeadChip and imputed to the Illumina BovineHD BeadChip (HD). GWAS were performed on 601,717 real and imputed single nucleotide polymorphism (SNP) genotypes using a single-SNP mixed linear model on 4841 Holstein bulls with breeding value predictions and followed by gene identification and in silico functional analyses. The association results were further validated using five scenarios with different numbers of SNPs.

Results

Seven hundred and eighty-two SNPs were significantly associated with calving performance at a genome-wise false discovery rate (FDR) of 5%. Most of these significant SNPs were on chromosomes 18 (71.9%), 17 (7.4%), 5 (6.8%) and 7 (2.4%) and mapped to 675 genes, among which 142 included at least one significant SNP and 532 were nearby one (100 kbp). For body conformation traits, 607 SNPs were significant at a genome-wise FDR of 5% and most of them were located on chromosomes 5 (30%), 18 (27%), 20 (13%), 6 (6%), 7 (5%), 14 (5%) and 13 (3%). SNP enrichment functional analyses for calving traits at a FDR of 1% suggested potential biological processes including musculoskeletal movement, meiotic cell cycle, oocyte maturation and skeletal muscle contraction. Furthermore, pathway analyses suggested potential pathways associated with calving performance traits including tight junction, oxytocin signaling, and MAPK signaling (P < 0.10). The prediction ability of the 1206 significant SNPs was between 78 and 83% of the prediction ability of the BovineSNP50 SNPs for calving performance traits and between 35 and 79% for body conformation traits.

Conclusions

Various SNPs that are significantly associated with calving performance are located within or nearby genes with potential roles in tight junction, oxytocin signaling, and MAPK signaling. Combining the significant SNPs or SNPs within or nearby gene(s) from the HD panel with the BovineSNP50 panel yielded a marginal increase in the accuracy of prediction of genomic estimated breeding values for all traits compared to the use of the BovineSNP50 panel alone.

Electronic supplementary material

The online version of this article (10.1186/s12711-017-0356-8) contains supplementary material, which is available to authorized users.

A Pathway-Centered Analysis of Pig Domestication and Breeding in Eurasia

Ascertaining the molecular and physiological basis of domestication and breeding is an active area of research. Due to the current wide distribution of its wild ancestor, the wild boar, the pig (Sus scrofa) is an excellent model to study these processes, which occurred independently in East Asia and Europe ca. 9000 yr ago. Analyzing genome variability patterns in terms of metabolic pathways is attractive since it considers the impact of interrelated functions of genes, in contrast to genome-wide scans that treat genes or genome windows in isolation. To that end, we studied 40 wild boars and 123 domestic pig genomes from Asia and Europe when metabolic pathway was the unit of analysis. We computed statistical significance for differentiation (Fst) and linkage disequilibrium (nSL) statistics at the pathway level. In terms of Fst, we found 21 and 12 pathways significantly differentiated at a q-value < 0.05 in Asia and Europe, respectively; five were shared across continents. In Asia, we found six significant pathways related to behavior, which involved essential neurotransmitters like dopamine and serotonin. Several significant pathways were interrelated and shared a variable percentage of genes. There were 12 genes present in >10 significant pathways (in terms of Fst), comprising genes involved in the transduction of a large number of signals, like phospholipase PCLB1, which is expressed in the brain, or ITPR3, which has an important role in taste transduction. In terms of nSL, significant pathways were mainly related to reproductive performance (ovarian steroidogenesis), a similarly important target trait during domestication and modern animal breeding. Different levels of recombination cannot explain these results, since we found no correlation between Fst and recombination rate. However, we did find an increased ratio of deleterious mutations in domestic vs. wild populations, suggesting a relaxed functional constraint associated with the domestication and breeding processes. Purifying selection was, nevertheless, stronger in significantly differentiated pathways than in random pathways, mainly in Europe. We conclude that pathway analysis facilitates the biological interpretation of genome-wide studies. Notably, in the case of pig, behavior played an important role, among other physiological and developmental processes.

Analysis of hepatic transcript profile and plasma lipid profile in early lactating dairy cows fed grape seed and grape marc meal extract

Background

It was recently reported that dairy cows fed a polyphenol-rich grape seed and grape marc meal extract (GSGME) during the transition period had an increased milk yield, but the underlying reasons remained unclear. As polyphenols exert a broad spectrum of metabolic effects, we hypothesized that feeding of GSGME influences metabolic pathways in the liver which could account for the positive effects of GSGME in dairy cows. In order to identify these pathways, we performed genome-wide transcript profiling in the liver and lipid profiling in plasma of dairy cows fed GSGME during the transition period at 1 week postpartum.

Results

Transcriptomic analysis of the liver revealed 207 differentially expressed transcripts, from which 156 were up- and 51 were down-regulated, between cows fed GSGME and control cows. Gene set enrichment analysis of the 155 up-regulated mRNAs showed that the most enriched gene ontology (GO) biological process terms were dealing with cell cycle regulation and the most enriched Kyoto Encyclopedia of Genes and Genomes pathways were p53 signaling and cell cycle. Functional analysis of the 43 down-regulated mRNAs revealed that a great part of these genes are involved in endoplasmic reticulum (ER) stress-induced unfolded protein response (UPR) and inflammatory processes. Accordingly, protein folding, response to unfolded protein, unfolded protein binding, chemokine activity and heat shock protein binding were identified as one of the most enriched GO biological process and molecular function terms assigned to the down-regulated genes. In line with the transcriptomics data the plasma concentrations of the acute phase proteins serum amyloid A (SAA) and haptoglobin were reduced in cows fed GSGME compared to control cows. Lipidomic analysis of plasma revealed no differences in the concentrations of individual species of major and minor lipid classes between cows fed GSGME and control cows.

Conclusions

Analysis of hepatic transcript profile in cows fed GSGME during the transition period at 1 week postpartum indicates that polyphenol-rich feed components are able to inhibit ER stress-induced UPR and inflammatory processes, both of which are considered to contribute to liver-associated diseases and to impair milk performance in dairy cows, in the liver of dairy cows during early lactation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3638-1) contains supplementary material, which is available to authorized users.

Pathway-based genome-wide association analysis of milk coagulation properties, curd firmness, cheese yield, and curd nutrient recovery in dairy cattle

It is becoming common to complement genome-wide association studies (GWAS) with gene-set enrichment analysis to deepen the understanding of the biological pathways affecting quantitative traits. Our objective was to conduct a gene ontology and pathway-based analysis to identify possible biological mechanisms involved in the regulation of bovine milk technological traits: coagulation properties, curd firmness modeling, individual cheese yield (CY), and milk nutrient recovery into the curd (REC) or whey loss traits. Results from 2 previous GWAS studies using 1,011 cows genotyped for 50k single nucleotide polymorphisms were used. Overall, the phenotypes analyzed consisted of 3 traditional milk coagulation property measures [RCT: rennet coagulation time defined as the time (min) from addition of enzyme to the beginning of coagulation; k₂₀: the interval (min) from RCT to the time at which a curd firmness of 20 mm is attained; a₃₀: a measure of the extent of curd firmness (mm) 30 min after coagulant addition], 6 curd firmness modeling traits [RCT_eq: RCT estimated through the CF equation (min); CF_P: potential asymptotic curd firmness (mm); k_CF: curd-firming rate constant (% × min^-1); k_SR: syneresis rate constant (% × min^-1); CF_max: maximum curd firmness (mm); and t_max: time to CF_max (min)], 3 individual CY-related traits expressing the weight of fresh curd (%CY_CURD), curd solids (%CY_SOLIDS), and curd moisture (%CY_WATER) as a percentage of weight of milk processed and 4 milk nutrient and energy recoveries in the curd (REC_FAT, REC_PROTEIN, REC_SOLIDS, and REC_ENERGY calculated as the % ratio between the nutrient in curd and the corresponding nutrient in processed milk), milk pH, and protein percentage. Each trait was analyzed separately. In total, 13,269 annotated genes were used in the analysis. The Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway databases were queried for enrichment analyses. Overall, 21 Gene Ontology and 17 Kyoto Encyclopedia of Genes and Genomes categories were significantly associated (false discovery rate at 5%) with 7 traits (RCT, RCT_eq, k_CF, %CY_SOLIDS, REC_FAT, REC_SOLIDS, and REC_ENERGY), with some being in common between traits. The significantly enriched categories included calcium signaling pathway, salivary secretion, metabolic pathways, carbohydrate digestion and absorption, the tight junction and the phosphatidylinositol pathways, as well as pathways related to the bovine mammary gland health status, and contained a total of 150 genes spanning all chromosomes but 9, 20, and 27. This study provided new insights into the regulation of bovine milk coagulation and cheese ability that were not captured by the GWAS.

The 'heritability' of domestication and its functional partitioning in the pig

We propose to estimate the proportion of variance explained by regression on genome-wide markers (or genomic heritability) when wild/domestic status is considered the phenotype of interest. This approach differs from the standard Fst in that it can accommodate genetic similarity between individuals in a general form. We apply this strategy to complete genome data from 47 wild and domestic pigs from Asia and Europe. When we partitioned the total genomic variance into components associated to subsets of single nucleotide polymorphisms (SNPs) defined in terms of their annotation, we found that potentially deleterious non-synonymous mutations (9566 SNPs) explained as much genetic variance as the whole set of 25 million SNPs. This suggests that domestication may have affected protein sequence to a larger extent than regulatory or other kinds of mutations. A pathway-guided analysis revealed ovarian steroidogenesis and leptin signaling as highly relevant in domestication. The genomic regression approach proposed in this study revealed molecular processes not apparent through typical differentiation statistics. We propose that at least some of these processes are likely new discoveries because domestication is a dynamic process of genetic selection, which may not be completely characterized by a static metric like Fst. Nevertheless, and despite some particularly influential mutation types or pathways, our analyses tend to rule out a simplistic genetic basis for the domestication process: neither a single pathway nor a unique set of SNPs can explain the process as a whole.