Detection of selection signatures in dairy and beef cattle using high-density genomic information
© Zhao et al. 2015
Received: 23 February 2014
Accepted: 19 May 2015
Published: 19 June 2015
Artificial selection for economically important traits in cattle is expected to have left distinctive selection signatures on the genome. Access to high-density genotypes facilitates the accurate identification of genomic regions that have undergone positive selection. These findings help to better elucidate the mechanisms of selection and to identify candidate genes of interest to breeding programs.
Information on 705 243 autosomal single nucleotide polymorphisms (SNPs) in 3122 dairy and beef male animals from seven cattle breeds (Angus, Belgian Blue, Charolais, Hereford, Holstein-Friesian, Limousin and Simmental) were used to detect selection signatures by applying two complementary methods, integrated haplotype score (iHS) and global fixation index (FST). To control for false positive results, we used false discovery rate (FDR) adjustment to calculate adjusted iHS within each breed and the genome-wide significance level was about 0.003. Using the iHS method, 83, 92, 91, 101, 85, 101 and 86 significant genomic regions were detected for Angus, Belgian Blue, Charolais, Hereford, Holstein-Friesian, Limousin and Simmental cattle, respectively. None of these regions was common to all seven breeds. Using the FST approach, 704 individual SNPs were detected across breeds. Annotation of the regions of the genome that showed selection signatures revealed several interesting candidate genes i.e. DGAT1, ABCG2, MSTN, CAPN3, FABP3, CHCHD7, PLAG1, JAZF1, PRKG2, ACTC1, TBC1D1, GHR, BMP2, TSG1, LYN, KIT and MC1R that play a role in milk production, reproduction, body size, muscle formation or coat color. Fifty-seven common candidate genes were found by both the iHS and global FST methods across the seven breeds. Moreover, many novel genomic regions and genes were detected within the regions that showed selection signatures; for some candidate genes, signatures of positive selection exist in the human genome. Multilevel bioinformatic analyses of the detected candidate genes suggested that the PPAR pathway may have been subjected to positive selection.
This study provides a high-resolution bovine genomic map of positive selection signatures that are either specific to one breed or common to a subset of the seven breeds analyzed. Our results will contribute to the detection of functional candidate genes that have undergone positive selection in future studies.
Artificial selection in cattle has resulted in divergent breeds that are specialized for either milk or meat production or raised as dual-purpose breeds. Such selection strategies are likely to have imposed selection pressures on particular regions of the genome that control these traits as well as other important animal characteristics such as disease resistance or general immune competence. Under positive selection pressure, the frequency of favorable alleles in the genome will rapidly increase. If intensive selection pressure occurred only over a few generations, it is unlikely that recombination had an impact on haplotype structure, and thus it resulted in (extended) linkage disequilibrium (LD) patterns between the mutation and neighboring loci . Analysis of these selection signatures can reveal genomic regions of interest for selection and provide insights into the mechanisms of evolution [2, 3].
Various statistical approaches have been proposed for the detection of selection signatures. Such tests include Tajima’s D-statistic , Fay and Wu’s H-statistic , extended haplotype homozygosity (EHH) , integrated haplotype score (iHS) , the Ka/Ks test , and the McDonald and Kreitman test . The EHH test is particularly useful to detect signatures of positive selection within a population using single nucleotide polymorphism (SNP) data [9–11]. This method that was first developed by Sabeti et al.  exploits knowledge on the relationship between the frequency of an allele and the measures of LD with neighboring alleles. An EHH is defined as the probability that two randomly chosen chromosomes that carry the core haplotype of interest are identical by descent for the entire interval between the core region and a certain locus . To overcome the influence of heterogeneous recombination rates across the genome, Voight et al.  developed the iHS approach, which is an extension of the EHH method and is based on the comparison of EHH between derived and ancestral alleles within a population. The iHS achieves maximal power when a selected allele segregates at intermediate frequencies in the population.
An alternative approach to the detection of selection signatures is based on the measure of population differentiation due to locus-specific allele frequencies between populations, which is quantified using the FST statistic . The fixation index, FST was first defined by Wright  to quantify the degree of genetic differentiation among populations based on differences in allele frequencies. FST provides information on the genomic variation at a locus among populations relative to that within populations. Thus, FST is also a test for evidence of selection i.e. high FST values indicate local positive adaptation while low FST values suggest negative or neutral selection .
Both iHS and FST statistics are useful to detect selection signatures . Previous analyses suggested that they are largely complementary; iHS has good power to detect selection signatures within breeds, while global FST is useful to detect selection signatures (i.e., loci that were differentially fixed in different breeds) across breeds . Global FST is also used to determine how divergent selection has impacted the genome of these breeds. The objective of our study was to detect signatures of selection using a large dataset of beef and dairy cattle with high-density SNP genotyping data. Potential biological functions of the genes that are present in the identified selection signatures were also examined using multi-level bioinformatic analyses.
Animal Care and Use Committee approval was not obtained for this study because all the data used were from the pre-existing database infrastructure operated by the Irish Cattle Breeding Federation (ICBF, Bandon, Co. Cork, Ireland).
Illumina (http://www.illumina.com) high-density genotypes (777 962 SNPs) were available on 3122 dairy and beef bulls; all animals had a genotype call rate of at least 95 %. The number of bulls per breed was 269, 196, 710, 234, 719, 730, and 264 for Angus, Belgian Blue, Charolais, Hereford, Holstein-Friesian, Limousin and Simmental, respectively. Mendelian inconsistencies among autosomal genotypes were used to validate animal identification through parentage assessment but also to discard 2816 reportedly autosomal SNPs that did not adhere to Mendelian inheritance patterns. An additional 11 654 autosomal SNPs with GenTrain scores less than 0.55 (i.e., a measure of genotype call quality) and a call rate less than greater than 90 % were also discarded as well as 29 939 SNPs that were monomorphic across all breeds or for which the position on the genome was unknown. The UMD3.1 genome build was used. Missing genotypes were imputed and genotypes were phased using Beagle Version 3.1.0 (http://faculty.washington.edu/browning/beagle/beagle.html) [17, 18]. After quality control, 705 243 SNPs were available with a mean distance of 3.56 kb between adjacent SNPs [See Additional file 1: Table S1].
Calculation of inbreeding coefficients
The pedigrees of all animals were traced back to the founder populations and mean inbreeding coefficients per breed were calculated using the algorithm in .
Detection of genomic regions with selection signatures
Integrated haplotype score (iHS) test
The iHS score is based on a ratio of extended haplotype homozygosities (EHH) associated with each allele. Thus, the iHS method requires information on the status of the ancestral and derived alleles for each SNP. Before computing iHS, the ancestral allele of all bovine SNPs was established from http://genome.jouy.inra.fr/downloads/Bovine_Ancestral_Allele/ . The iHS score was computed for each autosomal SNP using the R package “rehh” .
where iHH A and iHH D represent the integrated EHH score for ancestral and derived core alleles, respectively. Values of iHS were standardized so that they followed a standard normal distribution . To calculate the P value at the genomic level, iHS scores for each SNP were further transformed as p iHS = − log[1 − 2|Φ(iHS) − 0.5|], where Φ(x) represents the Gaussian cumulative distribution function (under neutrality) and p iHS is the two sided P-value associated with the neutral hypothesis (i.e., no selection) . In order to control for false positives, the R package “fdrtool”  was used with its default options for “statistic = p-value”, which uses the empirical data below the 75th percentile to determine the null distribution of the test statistics. After false discovery rate (FDR) adjustment within a breed, the genome-wide significance level was equal to approximately 0.003.
To better understand the genetic divergence among all breeds, FST was calculated using the HierFstat R package  with the unbiased estimator proposed by Weir and Cockerham . The negative FST values obtained for 24 800 SNPs were set to 0, since negative values have no biological interpretation . Raw global FST values were ranked and used to identify regions under positive selection. The empirical P-value was calculated for each SNP as a proportion of the total number of SNPs [26, 27]. As in , the genome-wide significance level was set to 0.001. i.e., only the top 0.1 % FST values were considered to represent a selection signature. Hence, no adjustment was made for multiple-testing for this statistic.
A gene was considered as being under selection if it overlapped with significant genomic windows based on iHS or if it contained an unexpectedly high proportion of highly differentiated SNPs based on FST values. Gene annotation was performed by exploiting the knowledge on UMD3.1 locations of genes from the NCBI (ftp://ftp.ncbi.nih.gov/genomes/Bos_taurus/mapview/seq_gene.md.gz). Because the annotation of the bovine genome is still incomplete, BioMart (www.ensembl.org/biomart) was used to determine the orthologous human gene ID for each gene detected. Enrichment analysis of these genes was performed using DAVID 6.7 by aligning the detected genes to human genes . Functional annotations (Gene Ontology (GO) Biological Process, GO Cellular Component, GO Molecular Function and Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway) were assigned to genes using the functional annotation tool.
Inbreeding coefficients per breed
Mean inbreeding coefficients of 0.0059, 0.0163, 0.0046, 0.0118, 0.0333, 0.0043 and 0.0106 were found for Angus, Belgian Blue, Charolais, Hereford, Holstein-Friesian, Limousin and Simmental breeds, respectively. In addition, genomic relationships were determined by calculating the Euclidean distances between alleles among all animals of the seven breeds analyzed (Figure S1 [See Additional file 2: Figure S1]). This figure shows that all the breeds can be clearly distinguished except Holstein and Friesian, for which animals were assigned to either of the two breeds based on their greatest breed proportion but many of the animals were actually crosses between Holstein and Friesian.
Genomic region and associated genes of the top 5 significant iHS for each of the seven bovine breeds
Genomic region (kb)
Mean |iHS| value
1.6 × 10−13
SPATA6, SLC5A9, LOC100295301, SLINT1, LOC100337055, LCO787081
3.6 × 10−12
CMPK1, STIL, TAL1, PDZK1IP1, LOC100847677, LOC784358, CYP4A11, LOC784417, LOC787638, CYP4A22, LOC787656, CYP4B1
1.7 × 10−11
LCO513399, ERRFI1, PARK7, TNFRSF9, UTS2, LOC100848366, PER3, VAMP3, CAMTA1
3.8 × 10−10
KIAA0494, ATPAF1-AS1, ATPAF1, MOB3C, MKNK1, KNCN, DMBX1, LOC513210, FAAH, NSUN4, UQCRH, LOC100847274, LRRC41, RAD54L, POMGNT1, C3H1orf190, TSPAN1
5.1 × 10−9
ECHDC2, ZYG11A, ZYG11B, SELRC1, FAM159A, GPX7, ZCCHC11, LOC100138140
1.26 × 10−9
AK5, PIGK, ST6GALNAC5
5.44 × 10−9
1.34 × 10−8
TJP1, TARSL2, TM2D3, LOC100335373
6.61 × 10−8
LOC781635, LOC100140130, LOC100848009, LOC784799
2 × 10−12
SGK3, C14H80rf45, LOC784087, LOC100847363, TCF24, PPP1R42, COPS5, CSPP1, APFGEF1
2.56 × 10−12
CD109, LOC100294729, LOC100336449
1.54 × 10−10
DOCK3, MAPKAPK3, CISH, HEMK1, C22H3orf18, CACNA2D2
2.95 × 10−9
8.15 × 10−10
SCFD2, FIP1L1, LNX1
3.79 × 10−9
LRRC66, SGCB, LOC100335977, SPATA18, LOC100847183, USP46
6.88 × 10−9
USP46, MIR2445, LOC100847282, RASL11B, SCFD2
8.43 × 10−9
OCIAD1, OCIAD2, CWH43, DCUN1D4
2.06 × 10−8
1.36 × 10−12
DIP2C, ZMYND11, PRNP, PRND, RASSF2
2.55 × 10−11
GPCPD1, LOC513580, LOC100140729, C13H20orf196, CHGB, TRMT6, MCM8
5.88 × 10−10
LOC100337129, LOC100299180, ABCC4, LOC530437
2.08 × 10−9
ADARB2, LOC100297660, WDR37, IDI1, GTPBP4, LARP4B, DIP2C
2.67 × 10−9
LOC507930, PROC, MAP3K2, ERCC3, CYP27C1, LOC784980, LOC524236, BIN1, MIR2350
HIBCH, C2H2orf88, LOC100335775, MSTN, OLC100335809, PMS1
1.18 × 10−11
2.31 × 10−11
2.98 × 10−11
LOC781635, LOC100140130, LOC100848009, LOC784799
1.05 × 10−12
LOC781635, LOC100140130, LOC100848009, LOC784799
7.55 × 10−12
8.18 × 10−12
DGCR8, TRMT2A, RANBP1, LOC526847, LOC100848428, RTN4R, LOC100336451, LOC100138815, LOC100301173, HSFY2, LOC786340, LOC100336511
1.07 × 10−11
ESYT1, ZC3H10, LOC100848780, PA2G4, ERBB3, RPS26, IKZF4, SUOX, RAB5B, CDK2, PMEL, DGKA, WIBG, LOC785991, MMP19, MGC142702, DNAJC14, ORMDL2, SARNP, GDF11, CD63, RDH5, BLOC1S1, ITGA7, METTL7B, LOC520938, OR10P1, OR6C4, LOC530539, LOC781363, LOC515967
3.34 × 10−9
Number of candidate genomic regions for each breed (on the diagonal) and number of overlapping candidate genomic regions between pair-wise breeds (below the diagonal) detected by iHS
Enriched pathway terms for genes in regions under selection with count, ratio and P- value for seven bovine breeds
Pathway terms description
hsa03320:PPAR signaling pathway
hsa00071:Fatty acid metabolism
hsa04514:Cell adhesion molecules (CAM)
hsa03320:PPAR signaling pathway
In total, 357 of the 704 significant FST values were located in 120 gene regions. The locus with the greatest FST value (0.914) was within an RNA region of the SCFD2 gene [See Additional file: 3 Table S6]. Seventeen other SNPs with significant FST values also resided within this genomic region. The genes that overlapped with the three other FST peaks were FANCA, a candidate for breast cancer susceptibility , PLAG1 that is associated with stature and body weight , and BIN1 that is associated with Alzheimer’s disease . Some of these 120 genes were overrepresented in the PPAR signaling pathway (Table 3). In addition, 13 GO terms that are associated with bone development, metabolic response and reproductive traits [See Additional file: 3 Table S5] were represented by these 120 genes.
Genes that were detected by both methods
Common genes identified by the complementary approaches iHS and FST
BPIFA4, BPIFB1, DCUN1D4, EPB41L1, LOC100337489, NECAB3, OCIAD1, PRDM5, SCP2, SNTA1
LOC530539, MGC142702, R3HDM1, TPRG1
DCUN1D4, FKBP5, FRYL, IMPACT, KCNN2, LNX1, LOC786242, LOC789547, LOC789558, MAML3, OCIAD1, PTK2, SCFD2, SGCB, SPATA18, TTC39C, ZNF521
ACSF3, ANKRD11, C18H16orf7, CPNE7, FANCA, LOC100299071, LOC100848941, LYN, SPG7, TGS1, TMEM68, UBE3A, XKR4, ZNF276
BIN1, COL5A2, HERC2, INPP1, MFSD6, OSGEPL1, PLAG1, PMS1, SDR16C5, SLC40A1, TPRG1
CNPY2, DCUN1D4, KIT, LOC530539, MGC142702, OCIAD1, SCFD2
In this work, our aim was to detect selection signatures in cattle using high-density genotypes (777 962 SNP) for seven breeds, with a mean distance between adjacent SNPs of 3.56 kb. Qanbari et al.  suggested that, in cattle, core regions may be more accurately detected by the relative EHH statistic if more than 50 000 SNPs were used. Barendse et al.  also suggested that more than 150 000 evenly spaced SNPs on the bovine genome would be an ideal number to accurately detect selection signatures using the FST approach.
Using simulations, Vilas et al.  recommended caution regarding the extent of false positive selection signatures which could be in fact false positive results. They detected apparent selection signatures on a large proportion of the simulated chromosomes for which actually no QTL had been simulated. In order to control the false positive rate in our study, we applied FDR adjustment within each breed. For FST, only the top 0.1 % FST values were considered as representing selection signatures as recommended in previous studies [26–28]. Nonetheless, it is likely that some of the apparent selection signatures detected in our study may represent false positive results. However, since many of the selection signatures that we detected are consistent with other reports on selection signatures in independent cattle populations [9, 34, 37–41] and since, overall, they confirm the results of genome-wide association studies for the same traits [42–44], we consider that the number of false positive results in our study is small. Applying a stricter FDR would inevitably reduce the number of true positive selection signatures detected.
Identifying recent positive selection signatures in domesticated animals can provide information on genomic regions that are under the influence of both artificial and natural selection, and thus, can help the identification of beneficial mutations and underlying biological pathways for economically important traits. Here, we used two different, yet complementary, statistical approaches, iHS and global FST, to detect selection signatures. The iHS approach appears to be the most powerful for detecting ongoing selection processes for which the target allele has a moderate to high frequency (0.4 < P < 0.8) within a population . Once an allele becomes fixed, the iHS statistic can still identify selection signatures in the regions of strong LD surrounding the selected site rather than the region itself because fixation eliminates variability at and near the selected site. If the iHS method detects a genomic region, this region can contain several loci that may actually be undergoing selection within the breed. Therefore, the iHS method can detect breed-specific candidate genes under positive selection. For example, the MC1R gene that plays a role in coat color types, including black and white coats and spotted phenotypes, was only detected by the iHS analysis in the Holstein-Friesian population. Previously, this method was successfully applied in human [6, 46] and pig  populations as well as in other international cattle populations [9, 22, 48].
Global FST values are useful to detect selection signatures across breeds (i.e., loci for which alleles are differentially fixed in different breeds) . Global FST analysis identifies selection signatures that are common to different breeds, and determines how divergent selection may have affected the genomic pattern of these breeds. In our study, global FST analysis highlighted genes that are associated with phenotypes that differ among breeds, in particular, coat color and body size. It has been shown that the KIT gene is associated with the level of white coat spotting in cattle [49, 50], which differentiates the breeds included here. For the PLAG1, LYN and TGS1 genes, several studies have documented their association with stature in both cattle and human populations [32, 51–54], thus, they may explain the differences in height between the breeds investigated in our study. One genomic region was detected by both iHS and global FST analyses, which indicates that both positive and divergent selection is acting on this region. Furthermore, complete concordance was found between the genes within regions of selection signatures identified by both methods (Table 4), which probably indicates true positive selection signatures. Integrating these two complementary approaches provides a valuable tool for positioning genomic regions that have undergone positive selection with more confidence.
Genotyping SNPs, which were discovered in another study, can result in ascertainment bias  and if the protocol used to identify the SNPs is not known, it will not be possible to directly correct for ascertainment bias. Based on simulations, Voight et al.  generated SNPs that had the same allele frequencies as in the real dataset in order to control, at least partially, for the effects of ascertainment bias. Nevertheless, in most studies, little or no attempt has been made to correct for ascertainment bias, and its effect is currently unknown. Previous genome-wide studies to detect positive selection in cattle have used the Bovine SNP50 BeadChip, which can also suffer from ascertainment bias due to the protocol used to discover the SNPs as well as to limited resolution. High-density SNP panels such as the Illumina BovineHD SNP chip have been designed to be less sensitive to ascertainment bias . Furthermore, iHS analysis exploits information on allele frequencies of both selected and neighboring SNPs, which increases its power to detect selection signatures . This method is more suited to genotyping data generated from SNP chips than to whole-genome sequence data, which minimizes the problems of ascertainment bias [34, 57]. To completely overcome ascertainment bias, a large-scale whole-genome sequencing project across breeds is necessary.
Although ascertainment bias may occur with the genotyping data used in our study, its effect is probably the same along the whole genome, unlike selection pressure that acts at certain genomic regions and will impact FST values for those regions only. The mean (± standard deviation) minor allele frequency (MAF) per SNP was similar in all breeds and ranged from 0.224 ± 0.162 (Angus) to 0.245 ± 0.158 (Hereford). In general, methods for the detection of selection signatures are designed to analyze non-related animals. Our animal sample consisted of 3122 dairy and beef animals from seven breeds. All animals were selected for genotyping with the high-density SNP panel to maximize imputation accuracy of their descendants from lower-density genotypes. Thus, although these animals were related, they were chosen to have as many progeny as possible and to be as genetically diverse as possible. Furthermore, they were highly selected animals and although it is likely that some false positive selection signatures may result from random genetic drift (and other factors), many of the detected signals probably reflect true selection signatures. Inbreeding will increase the extent of LD, which may result in false positives or type I errors. However the level of inbreeding in the populations analyzed here was relatively low.
Apart from inbreeding, the demographic history of a population can also influence the variome, i.e. the whole set of genetic variations found for a population of a given species, which complicates the interpretation of selection signatures. Expansion of a population increases the frequency of alleles that originally have a low frequency compared to expectations under a neutral model. Similarly, recent positive selection for an allele may have begun from a set of beneficial alleles with a higher initial frequency . Such alleles may have been introgressed into a population through historical crossbreeding and, thus, be included in various haplotypes, which prevents LD-based estimators to detect the selection signature. Furthermore, crossbreeding can also generate false selection signatures, if for example a large conserved region of the genome from another breed is mixed with many smaller segments from the genome of the original breed .
Genome-wide analyses of selection signatures were reported for several international Holstein populations including Chinese , German [9, 34] and Israeli  Holsteins. Several of the genes that we identified here confirm previously documented selection signatures in Holstein cattle populations, such as ACTC1 , FABP3 , RORA , GHR  and LACTB . Of particular interest is the region on BTA20 that was detected by the iHS method as having a strong selection signature in Holstein-Friesian cattle (Fig. 1) and [See Additional file: 2 Figure S2]. This result confirms selection signatures reported in a population of Israeli Holstein cows  and to a lesser extent in German Holstein cows . This region on BTA20 between 20 and 40 Mb harbors many genes including GDNF, WDR70, NUP155, GHR, ITGA2, LOC100847619, ITGA1, PELO, NDUFS4, FST, LOC100847646, LOC782165, MOCS2, ITGA2, NIM1, ZNF131, LOC100848437, LOC100336494, LOC785615, LOC785744, LOC100139184, LOC100848479, LOC783463, LOC527137, SEPP1, CCDC152, LOC100848533, PARP8, EMB and LOC785429.
For the beef cattle populations studied here, several of the candidate genes found for body size were previously reported in horse , human [15, 53], dog  and/or cattle  populations. These genes included CHCHD7, PLAG1 and SMAD2 for the Limousin breed, GDF5 for the Angus and Simmental breeds, CDK6 for Simmental, JAZF1 and PRKG2 for Belgian Blue but selection signatures that overlapped with these genes were not found for the Holstein-Friesian population. Furthermore, some of the genes that overlapped with selection signatures in our study were consistent with those identified in other beef cattle populations [9, 38, 39, 41], such as ACTC1 in the Charolais and Holstein-Friesian populations, a gene that is related to muscle formation .
In addition, some of the candidate genes that we detected were previously found by GWAS on cattle populations. Several of the genes that were identified here by the iHS analysis were previously suggested to be associated with milk production, fertility, body size or body conformation [42–44]. Apart from the aforementioned genes associated with body size, body weight and feed intake, we also identified DGAT1, ABCG2, MSTN, GHR, CAPN3, PDGFRA, GAS1, ZNF521 and TMEM130. Biological justifications of why many of these genes reside within selection signatures were discussed in detail elsewhere . Detection of selection signatures and GWAS are two different approaches to identify candidate genes of interest . GWAS evaluates the relationship between genotype and phenotype, while detection of selection signatures relies on population genetic and evolutionary parameters that are obtained only from genomic information.
In our study, it should be noted that DGAT1 and ABCG2 were not found within a selection signature in the Holstein-Friesian population, whereas they were detected in the Limousin and Charolais populations. Both genes exhibited selection signatures only in the beef breeds and not in the dairy breed. Previously, DGAT1 and ABCG2 were detected in selection sweeps [39, 60] and by GWAS for performance traits  in cattle. The results of the present study are nonetheless consistent with the results of Kemper et al.  who detected selection signatures in Limousin and Charolais populations for DGAT1 and ABCG2, but not in a Holstein population. The reason for not finding these genes in the Holstein-Friesian population used in our study may be that their alleles are no longer segregating in the population and therefore could not be detected by the iHS statistic. Possibly, the alleles that are still segregating, even after the intensive artificial selection during domestication, may have unfavorable pleiotropic effects that prevent their frequency from increasing in the Holstein-Friesian population. In addition, selection is likely to have affected standing variation. If the selected mutations were segregating on multiple different haplotypes before selection began, the iHS statistic may have too little power to detect the selection signature.
Some of the genomic regions that we identified here were previously documented to be under selection in human and other livestock populations. One selection signature region that was observed in the global FST analysis was on BTA 2 between 61881578 and 62129511 bp and contained the R3HDM1 and LCT genes . These two genes are associated with energy homeostasis; R3HDM1 has a role in efficient food conversion and intramuscular fat content in some breeds [62, 63], while LCT is involved in the digestion of lactose in human adults . These two genes have also been shown to be under positive selection in human populations . In addition, some of the candidate genes (TBC1D1, WIF1, LEMD3, KIT, and BMP2) that we detected here were previously found within selection signatures in pig [65, 66], sheep [28, 67] and horse  populations.
It should also be noted that we detected several poorly annotated genomic regions that appear to have undergone strong selection. For example, genomic regions that had the greatest P-value estimated by the iHS method were on BTA20 in the Holstein-Friesian population [See Additional file: 2 Figure S2] but no genes in this region have been documented (Table 2). Similar patterns were also observed for FST signatures [See Additional file: 3 Table S6]. This observation is consistent with other genome-wide analyses of selection signatures in cattle , thoroughbred horses  and humans . Thus, these results suggest that regions that do not appear to contain genes may also have an important role in adaptive evolution. Another reason, particularly in cattle, may be due to the relatively poor annotation of the bovine genome. Priority should be given to an improved annotation of the genomic regions that are suspected to be within positive selection signatures.
To better understand the molecular functions of these genes, we examined their GO classifications. Many of the genes detected in our study are consistent with expectations since they are involved in fatty acid metabolism, reproductive traits, and both meat and milk production. An intriguing candidate pathway that we identified is the PPAR signaling pathway which is known to be associated with meat quality and production traits in pigs  and cattle . These observations need to be explored and verified in an independent population.
We used two complementary methods (iHS and global FST) to detect selection signatures across the whole bovine genome and across seven diverse cattle breeds using high-density genotypes. Our analyses revealed multiple genes under positive selection, which are related to milk production, reproduction, body size, muscle formation and coat color. Moreover, we identified the PPAR signaling pathway, which is an intriguing candidate pathway. Our results can contribute to the identification of the variants that underlie the detected selection signatures. In most cases, further studies are required to distinguish between selection signatures that are due to breed-specific characteristics or traits of practical interest for agriculture. However, both types of selection signatures are relevant to better understand the mechanisms and identify the targets of natural and artificial selection in domesticated cattle.
We greatly thank two anonymous reviewers and the editor for their useful comments and suggestions on an earlier version of the manuscript. This study was financially supported by a grant from the Irish Department of Agriculture, Food and Marine Research Stimulus Fund (11/S/112), the Agricultural Science and Technology Innovation Program (No. ASTIP-IAS-TS-6) and the Natural Science Foundation of China (No. 31200927).
- Sabeti PC, Reich DE, Higgins JM, Levine HZ, Richter DJ, Schaffner SF, et al. Detecting recent positive selection in the human genome from haplotype structure. Nature. 2002;419:832–7.PubMedView ArticleGoogle Scholar
- Akey JM, Zhang G, Zhang K, Jin L, Shriver MD. Interrogating a high-density SNP map for signatures of natural selection. Genome Res. 2002;12:1805–14.PubMed CentralPubMedView ArticleGoogle Scholar
- Andersson L, Georges M. Domestic-animal genomics: deciphering the genetics of complex traits. Nat Rev Genet. 2004;5:202–12.PubMedView ArticleGoogle Scholar
- Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123:585–95.PubMed CentralPubMedGoogle Scholar
- Fay JC, Wu CI. Hitchhiking under positive Darwinian selection. Genetics. 2000;155:1405–13.PubMed CentralPubMedGoogle Scholar
- Voight BF, Kudaravalli S, Wen X, Pritchard JK. A map of recent positive selection in the human genome. PLoS Biol. 2006;4:e72.PubMed CentralPubMedView ArticleGoogle Scholar
- Hughes AL, Nei M. Pattern of nucleotide substitution at major histocompatibility complex class I loci reveals overdominant selection. Nature. 1988;335:167–70.PubMedView ArticleGoogle Scholar
- McDonald JH, Kreitman M. Adaptive protein evolution at the Adh locus in Drosophila. Nature. 1991;351:652–4.PubMedView ArticleGoogle Scholar
- Qanbari S, Gianola D, Hayes B, Schenkel F, Miller S, Moore S, et al. Application of site and haplotype-frequency based approaches for detecting selection signatures in cattle. BMC Genomics. 2011;12:318.PubMed CentralPubMedView ArticleGoogle Scholar
- Walsh EC, Sabeti P, Hutcheson HB, Fry B, Schaffner SF, de Bakker PI, et al. Searching for signals of evolutionary selection in 168 genes related to immune function. Hum Genet. 2006;119:92–102.PubMedView ArticleGoogle Scholar
- Zhang C, Bailey DK, Awad T, Liu G, Xing G, Cao M, et al. A whole genome long-range haplotype (WGLRH) test for detecting imprints of positive selection in human populations. Bioinformatics. 2006;22:2122–8.PubMedView ArticleGoogle Scholar
- Lewontin RC, Krakauer J. Distribution of gene frequency as a test of the theory of the selective neutrality of polymorphisms. Genetics. 1973;74:175–95.PubMed CentralPubMedGoogle Scholar
- Wright S. The genetical structure of populations. Ann Eugen. 1951;15:323–54.PubMedView ArticleGoogle Scholar
- Kullo IJ, Ding K. Patterns of population differentiation of candidate genes for cardiovascular disease. BMC Genet. 2007;8:48.PubMed CentralPubMedView ArticleGoogle Scholar
- Pybus M, Dall’Olio GM, Luisi P, Uzkudun M, Carreño-Torres A, Pavlidis P, et al. 1000 Genomes Selection Browser 1.0: a genome browser dedicated to signatures of natural selection in modern humans. Nucleic Acids Res. 2013;42:D903–9.PubMed CentralPubMedView ArticleGoogle Scholar
- Biswas S, Akey JM. Genomic insights into positive selection. Trends Genet. 2006;22:437–46.PubMedView ArticleGoogle Scholar
- Browning BL, Browning SR. A unified approach to genotype imputation and haplotype-phase inference for large data sets of trios and unrelated individuals. Am J Hum Genet. 2009;84:210–23.PubMed CentralPubMedView ArticleGoogle Scholar
- Browning SR. Browning BL Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. Am J Hum Genet. 2007;81:1084–97.PubMed CentralPubMedView ArticleGoogle Scholar
- Meuwissen THE, Luo Z. Computing inbreeding coefficients in large populations. Genet Sel Evol. 1992;24:305–13.PubMed CentralView ArticleGoogle Scholar
- Rocha D, Billerey C, Samson F, Boichard D, Boussaha M. Identification of the putative ancestral allele of bovine single-nucleotide polymorphisms. J Anim Breed Genet. 2014;131:483–6.PubMedView ArticleGoogle Scholar
- Gautier M, Vitalis R. rehh: an R package to detect footprints of selection in genome-wide SNP data from haplotype structure. Bioinformatics. 2012;28:1176–7.PubMedView ArticleGoogle Scholar
- Gautier M, Naves M. Footprints of selection in the ancestral admixture of a New World Creole cattle breed. Mol Ecol. 2011;20:3128–43.PubMedView ArticleGoogle Scholar
- Strimmer K. fdrtool: a versatile R package for estimating local and tail area-based false discovery rates. Bioinformatics. 2008;24:1461–2.PubMedView ArticleGoogle Scholar
- Goudet J. Hierfstat, a package for R to compute and test hierarchical F‐statistics. Molec Ecol Notes. 2005;5:184–6.View ArticleGoogle Scholar
- Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution. 1984;38:1358–70.View ArticleGoogle Scholar
- Moradi MH, Nejati-Javaremi A, Moradi-Shahrbabak M, Dodds KG, McEwan JC. Genomic scan of selective sweeps in thin and fat tail sheep breeds for identifying of candidate regions associated with fat deposition. BMC Genet. 2012;13:10.PubMed CentralPubMedView ArticleGoogle Scholar
- Dong K, Yao N, Pu Y, He X, Zhao Q, Luan Y, et al. Genomic scan reveals loci under altitude adaptation in Tibetan and Dahe pigs. PLoS One. 2014;9:e110520.PubMed CentralPubMedView ArticleGoogle Scholar
- Kijas JW, Lenstra JA, Hayes B, Boitard S, Porto-Neto LR, San Cristobal M, et al. Genome-wide analysis of the world’s sheep breeds reveals high levels of historic mixture and strong recent selection. PLoS Biol. 2012;10:e1001258.PubMed CentralPubMedView ArticleGoogle Scholar
- da Huang W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37:1–13.PubMed CentralView ArticleGoogle Scholar
- Wright S. Evolution and the genetics of population. In Variability within and among natural populations. Volume 4. Chicago: University of Chicago Press; 1978.Google Scholar
- Litim N, Labrie Y, Desjardins S, Ouellette G, Plourde K, Belleau P, et al. Polymorphic variations in the FANCA gene in high-risk non-BRCA1/2 breast cancer individuals from the French Canadian population. Mol Oncol. 2013;7:85–100.PubMedView ArticleGoogle Scholar
- Karim L, Takeda H, Lin L, Druet T, Arias JA, Baurain D, et al. Variants modulating the expression of a chromosome domain encompassing PLAG1 influence bovine stature. Nat Genet. 2011;43:405–13.PubMedView ArticleGoogle Scholar
- Kingwell K. Alzheimer disease: BIN1 variant increases risk of Alzheimer disease through tau. Nat Rev Neurol. 2013;9:184.PubMedView ArticleGoogle Scholar
- Qanbari S, Pimentel E, Tetens J, Thaller G, Lichtner P, Sharifi A, et al. A genome‐wide scan for signatures of recent selection in Holstein cattle. Anim Genet. 2010;41:377–89.PubMedGoogle Scholar
- Barendse W, Harrison BE, Bunch RJ, Thomas MB, Turner LB. Genome wide signatures of positive selection: The comparison of independent samples and the identification of regions associated to traits. BMC Genomics. 2009;10:178.PubMed CentralPubMedView ArticleGoogle Scholar
- Vilas A, Pérez‐Figueroa A, Caballero A. A simulation study on the performance of differentiation‐based methods to detect selected loci using linked neutral markers. J Evol Biol. 2012;25:1364–76.PubMedView ArticleGoogle Scholar
- Pan D, Zhang S, Jiang J, Jiang L, Zhang Q, Liu J. Genome-wide detection of selective signature in Chinese Holstein. PLoS One. 2013;8:e60440.PubMed CentralPubMedView ArticleGoogle Scholar
- Xu L, Bickhart DM, Cole JB, Schroeder SG, Song J, Van Tassell CP, et al. Genomic signatures reveal new evidences for selection of important traits in domestic cattle. Mol Biol Evol. 2015;32:711–25.PubMedView ArticleGoogle Scholar
- Kemper KE, Saxton SJ, Bolormaa S, Hayes BJ, Goddard ME. Selection for complex traits leaves little or no classic signatures of selection. BMC Genomics. 2014;15:246.PubMed CentralPubMedView ArticleGoogle Scholar
- Glick G, Shirak A, Uliel S, Zeron Y, Ezra E, Seroussi E, et al. Signatures of contemporary selection in the Israeli Holstein dairy cattle. Anim Genet. 2012;43 Suppl 1:45–55.PubMedView ArticleGoogle Scholar
- Ramey HR, Decker JE, McKay SD, Rolf MM, Schnabel RD, Taylor JF. Detection of selective sweeps in cattle using genome-wide SNP data. BMC Genomics. 2013;14:382.PubMed CentralPubMedView ArticleGoogle Scholar
- Wu XP, Fang M, Liu L, Wang S, Liu JF, Ding XD, et al. Genome wide association studies for body conformation traits in the Chinese Holstein cattle population. BMC Genomics. 2013;14:897.PubMed CentralPubMedView ArticleGoogle Scholar
- Jiang L, Liu JF, Sun DX, Ma PP, Ding XD, Yu Y, et al. Genome wide association studies for milk production traits in Chinese Holstein population. PLoS One. 2010;5:e13661.PubMed CentralPubMedView ArticleGoogle Scholar
- Doran AG, Berry DB, Creevey CJ. Whole genome association study identifieds regions of the bovine genome and biological pathways involved in carcass trait performance in Holstein-Friesian cattle. BMC Genomics. 2014;15:837.PubMed CentralPubMedView ArticleGoogle Scholar
- Simianer H. Statistical problems in livestock population genomics. In Proceedings of the 10th World Congress on Genetics Applied to Livestock Production: 17–22 August 2014; Vancouver. 2014. https://asas.org/docs/default-source/wcgalp-proceedings-oral/202_paper_10373_manuscript_1346_0.pdf?sfvrsn=2.
- Grossman SR, Shylakhter I, Karlsson EK, Byrne EH, Morales S, Frieden G, et al. A composite of multiple signals distinguishes causal variants in regions of positive selection. Science. 2010;327:883–6.PubMedView ArticleGoogle Scholar
- Ma Y, Zhang H, Zhang Q, Ding X. Identification of selection footprints on the X chromosome in pig. PLoS One. 2014;9:e94911.PubMed CentralPubMedView ArticleGoogle Scholar
- Qanbari S, Pausch H, Jansen S, Somel M, Strom TM, Fries R, et al. Classic selective sweeps revealed by massive sequencing in cattle. PLoS Genet. 2014;10:e1004148.PubMed CentralPubMedView ArticleGoogle Scholar
- Hayes BJ, Pryce J, Chamberlain AJ, Bowman PJ, Goddard ME. Genetic architecture of complex traits and accuracy of genomic prediction: coat colour, milk-fat percentage, and type in Holstein cattle as contrasting model traits. PLoS Genet. 2010;6:e1001139.PubMed CentralPubMedView ArticleGoogle Scholar
- Fontanesi L, Tazzoli M, Russo V, Beever J. Genetic heterogeneity at the bovine KIT gene in cattle breeds carrying different putative alleles at the spotting locus. Anim Genet. 2010;41:295–303.PubMedView ArticleGoogle Scholar
- Pryce JE, Hayes BJ, Bolormaa S, Goddard ME. Polymorphic regions affecting human height also control stature in cattle. Genetics. 2011;187:981–4.PubMed CentralPubMedView ArticleGoogle Scholar
- Gudbjartsson DF, Walters GB, Thorleifsson G, Stefansson H, Halldorsson BV, Zusmanovich P, et al. Many sequence variants affecting diversity of adult human height. Nat Genet. 2008;40:609–15.PubMedView ArticleGoogle Scholar
- Lettre G, Jackson AU, Gieger C, Schumacher FR, Berndt SI, Sanna S, et al. Identification of ten loci associated with height highlights new biological pathways in human growth. Nat Genet. 2008;40:584–91.PubMed CentralPubMedView ArticleGoogle Scholar
- Littlejohn M, Grala T, Sanders K, Walker C, Waghorn G, Macdonald K, et al. Genetic variation in PLAG1 associates with early life body weight and peripubertal weight and growth in Bos taurus. Anim Genet. 2012;43:591–4.PubMedView ArticleGoogle Scholar
- Nielsen R, Hellmann I, Hubisz M, Bustamante C, Clark AG. Recent and ongoing selection in the human genome. Nat Rev Genet. 2007;8:857–68.PubMed CentralPubMedView ArticleGoogle Scholar
- Porto-Neto LR, Sonstegard TS, Liu GE, Bickhart DM, Da Silva MVB, Machado MA, et al. Genomic divergence of zebu and taurine cattle identified through high-density SNP genotyping. BMC Genomics. 2013;14:876.PubMed CentralPubMedView ArticleGoogle Scholar
- Tang K, Thornton KR, Stoneking M. A new approach for using genome scans to detect recent positive selection in the human genome. PLoS Biol. 2007;5:e171.PubMed CentralPubMedView ArticleGoogle Scholar
- Petersen JL, Mickelson JR, Rendahl AK, Valberg SJ, Andersson LS, Axelsson J, et al. Genome-wide analysis reveals selection for important traits in domestic horse breeds. PLoS Genet. 2013;9:e1003211.PubMed CentralPubMedView ArticleGoogle Scholar
- Akey JM. Constructing genomic maps of positive selection in humans: Where do we go from here? Genome Res. 2009;19:711–22.PubMed CentralPubMedView ArticleGoogle Scholar
- Flori L, Fritz S, Jaffrézic F, Boussaha M, Gut I, Heath S, et al. The genome response to artificial selection: a case study in dairy cattle. PLoS One. 2009;4:e6595.PubMed CentralPubMedView ArticleGoogle Scholar
- Bouwman AC, Bovenhuis H, Visker MHPW, van Arendonk JAM. Genome-wide association of milk fatty acids in Dutch dairy cattle. BMC Genet. 2011;12:43.PubMed CentralPubMedView ArticleGoogle Scholar
- Bovine HapMap C, Gibbs RA, Taylor JF, Van Tassell CP, Barendse W, Eversole KA, et al. Genome-wide survey of SNP variation uncovers the genetic structure of cattle breeds. Science. 2009;324:528–32.View ArticleGoogle Scholar
- Barendse W, Reverter A, Bunch RJ, Harrison BE, Barris W, Thomas MB. A validated whole-genome association study of efficient food conversion in cattle. Genetics. 2007;176:1893–905.PubMed CentralPubMedView ArticleGoogle Scholar
- Sabeti PC, Varilly P, Fry B, Lohmueller J, Hostetter E, Cotsapas C, et al. Genome-wide detection and characterization of positive selection in human populations. Nature. 2007;449:913–8.PubMed CentralPubMedView ArticleGoogle Scholar
- Wilkinson S, Lu ZH, Megens H-J, Archibald AL, Haley C, Jackson IJ, et al. Signatures of diversifying selection in European pig breeds. PLoS Genet. 2013;9:e1003453.PubMed CentralPubMedView ArticleGoogle Scholar
- Ai H, Huang L, Ren J. Genetic diversity, linkage disequilibrium and selection signatures in Chinese and Western pigs revealed by genome-wide SNP markers. PLoS One. 2013;8:e56001.PubMed CentralPubMedView ArticleGoogle Scholar
- Fariello MI, Servin B, Tosser-Klopp G, Rupp R, Moreno C, San Cristobal M, et al. Selection signatures in worldwide sheep populations. PLoS One. 2014;9:e103813.PubMed CentralPubMedView ArticleGoogle Scholar
- Gu J, Orr N, Park SD, Katz LM, Sulimova G, MacHugh DE, et al. A genome scan for positive selection in thoroughbred horses. PLoS One. 2009;4:e5767.PubMed CentralPubMedView ArticleGoogle Scholar
- He K, Wang Q, Wang Z, Pan Y. Association study between gene polymorphisms in PPAR signaling pathway and porcine meat quality traits. Mamm Genome. 2013;24:322–31.PubMedView ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.