Skip to main content
  • Research article
  • Open access
  • Published:

Genome-wide candidate regions for selective sweeps revealed through massive parallel sequencing of DNA across ten turkey populations

Abstract

Background

The domestic turkey (Meleagris gallopavo) is an important agricultural species that is largely used as a meat-type bird. Characterizing genetic variation in populations of domesticated species and associating these variation patterns with the evolution, domestication, and selective breeding is critical for understanding the dynamics of genomic change in these species. Intense selective breeding and population bottlenecks are expected to leave signatures in the genome of domesticated species, such as unusually low nucleotide diversity or the presence of exceptionally extended haplotype homozygosity. These patterns of variation in selected populations are highly useful to not only understand the consequences of selective breeding and population dynamics, but also to provide insights into biological mechanisms that may affect physiological processes important to bring changes in phenotype of interest.

Results

We observed 54 genomic regions in heritage and commercial turkey populations on 14 different chromosomes that showed statistically significant (P < 0.05) reduction in genomic variation indicating candidate selective sweeps. Areas with evidence of selective sweeps varied from 1.5 Mb to 13.8 Mb in length. Out of these 54 sweeps, 23 overlapped at least partially between two or more populations. Overlapping sweeps were found on 13 different chromosomes. The remaining 31 sweeps were population-specific and were observed on 12 different chromosomes, with 26 of these regions present only in commercial populations. Genes that are known to affect growth were enriched in the sweep regions.

Conclusion

The turkey genome showed large sweep regions. The relatively high number of sweep regions in commercial turkey populations compared to heritage varieties and the enrichment of genes important to growth in these regions, suggest that these sweeps are the result of intense selection in these commercial lines, moving specific haplotypes towards fixation.

Background

Characterizing genetic variation in populations of domesticated species and associating these variation patterns with the evolution, domestication, and selective breeding is critical for understanding the dynamics of genomic change in these species. Intense selective breeding and population bottlenecks are expected to leave signatures in the genome of domesticated species, such as unusually low nucleotide diversity or the presence of exceptionally extended haplotype homozygosity [1]-[3]. Genome-wide characterization of many different breeds and populations for these selective sweeps, along with the functional knowledge of the region, can reveal which genes are linked to traits or diseases with a complex genetic basis [4]. These patterns of variation in selected populations are highly useful to not only understand the consequences of selective breeding and population dynamics, but also to provide insights into biological mechanisms that may affect physiological processes important to bring changes in phenotype of interest [5],[6].

The turkey (Meleagris gallopavo) is an important agricultural species that is largely used as a meat-type bird. All domesticated turkeys descend from the wild turkeys indigenous to North and South America. There are seven subspecies of the wild form [7] distinguished by geographic range and plumage differences. They are: South Mexican (M. g. gallopavo), Rio Grande (M. g. intermedia), Merriam’s (M. g. merriami), Gould’s (M. g. mexicana), Eastern (M. g. silverstris), Moore’s (M. g. oneusta) and Florida (M. g. osceola). Three of the seven are believed to play an important role in domestication. It is generally accepted that domestication of turkey involved South Mexican turkey [8]. The earliest signs of turkey domestication dates to approximately 2000 years ago at Mayan sites in Southern Mexico such as Cobá [9]. Domestic turkey stocks were established by at least 180 AD within the Tehuacán valley [10]. A separate domestication event likely occurred in what is now the Southwest United States, where the first strong archaeological evidence for domestic stocks dates to similar time (ca. 200 BC-AD 500), although the wild progenitor has been long debated [11].

The modern domestic turkey has been recognized by the American Standard of Perfection since 1971 [12], and is registered as a single breed with eight varieties defined primarily by plumage color. Out of these eight heritage turkey varieties, five (Bronze, Narragansett, White Holland, Spanish Black and Blue Slate) were registered in 1874 [12], while the remaining three (Beltsville Small White, Bourbon Red, and Royal Palm) were registered in 1951, 1909, and 1971 respectively [12]. These domestic turkeys are presumed to be highly inbred [12], and to have undergone intensive selection for traits of economic importance such as body weight, meat quality and egg production [9],[11].

Recent census data show that turkey is the second largest contributor in worldwide poultry meat production [13]. Global production of turkeys has experienced a massive expansion over the past 40 years. In 2008, turkeys represented 6.65% of the world poultry meat production [13]. Global turkey stocks nearly tripled from 178 million in 1970 to over 482 million in 2008 [13]. Astonishingly, in those four decades, average meat production per bird doubled from 6.7 to 12.7 Kg, showing the result of intensive selection in turkeys.

An important genomic indicator of a selective sweep involves local reduction in genetic variation within the selected gene(s) and in nearby single nucleotide polymorphism (SNP) variants [14]. Selection affects all the genomic variability in the genome, including SNPs, microsatellites and several types of structural variations (SVs). The SV category includes large insertions and deletions, inversions, duplications and balanced or unbalanced inter-chromosomal translocations. Next generation sequencing (NGS) is an efficient approach for a large-scale, genome-wide SNP discovery and genotyping of individuals [15],[16]. Availability of a high quality reference genome sequence [17] and resequencing of individuals or groups with appropriate genome coverage are key prerequisites for whole-genome SNP discovery [15],[16]. Genomic sequences of individuals are aligned to a reference genome to detect nucleotide variations, i.e., differences in genotype of individuals at specific positions of the genome [18],[19].

Our search was aimed at finding genomic regions where selection or domestication has changed the frequency of favourable alleles towards fixation. Genomic regions where these changes are observed elucidate the effect from the selective pressure of domestication or selection that was applied to the domesticated turkey.

Methods

Populations

Ten turkey populations that included seven commercial lines and three heritage varieties were used for whole genome sequencing (WGS). The seven commercial lines, L1 through L7, were provided by two breeding companies. Commercial lines were selected for different objectives including higher adult body weight and rapid growth except L5 which is a female line that was selected for medium adult body weight, conformation and egg production. The heritage varieties were Beltsville Small White (BvSW), Royal Palm (RP) and Narragansett (Nset) [20]-[22]. In total, 29 individuals were selected for WGS, with three individuals per population except for RP, which was represented by two individuals.

Genomic DNA Extraction, Library Preparation and Sequencing

Genomic DNA was extracted from whole blood with the QIAamp DNA blood Midi Kit (Qiagen, Valencia, CA); the procedure included a proteinase K digestion followed by column purification. Integrity of high molecular weight DNA following the extraction was confirmed by agarose gel analysis. Genomic DNA was sheared using the Covaris S2 to yield an average fragment size of 450 bp, as determined with the Agilent Bioanalyzer 2100 (Agilent, Santa Clara, CA).

Genomic libraries were prepared with the Paired-end Sequencing Sample Preparation Kit (Illumina, San Diego, CA) with 5 μg of genomic DNA according to the manufacturer’s instructions. All genomic DNA libraries were validated with the Agilent Bioanalyzer (model 2100). The automated cBot Cluster Generation System (Illumina) was used to generate clusters on the flow cell. Each individual was sequenced (paired-end; read length 120 bp) in a single lane of a flow cell using the Illumina GAIIx.

Sequence mapping

Sequence reads of each turkey were filtered on base quality; reads were trimmed if three consecutive bases had an average Phred-like quality score of less than 13. Both paired-end sequences of a fragment were required to be at least 40 bp long after trimming to be retained for analyses. Retained reads were aligned against the turkey reference genome (UMD 2.01) using the MOSAIK aligner [23]. Mapping of reads from each individual to the reference genome sequence was performed with hash size (hs) of 15, maximum hash positions (mhp) of 100, an alignment candidate threshold (act) of 20, and a maximum mismatch percentage (mmp) of 5. Banded Smith-Waterman algorithm (bw = 41) was used to increase the speed of alignments. The algorithm implemented in MOSAIK calculates a mapping quality for each sequence that measures the probability that a sequence belongs to a specific target. The alignments were filtered for ambiguously mapped reads and sorted using MosaikSort. Finally, the file was converted to BAM format [16] using MosaikText. All BAM files have been uploaded to NCBI's Sequence Read Archive (SRA) database under the study accession number “sRP012021” [24].

Heterozygosity

Genome wide nucleotide diversity across the whole genome was assessed for each individual of the different turkey populations. The pileup function of SamTools version 0.1.12a [15] was used to perform SNP genotype calling, after which the nucleotide diversity was estimated across the whole genome for each individual separately. Nucleotide diversity was estimated by calculating the number of heterozygous SNP as well as the number of homozygous non-reference genotypes within each 300 Kb window. Windows of 300 Kb were necessary to avoid large random fluctuations in heterozygosity that were observed in a preliminary analysis with smaller windows. The random fluctuations with smaller windows were due to the low SNP detection rate. For calling SNPs, coverage per base was limited to 5-10 fold to avoid analysing repetitive regions of the genome as the average sequence depth per animal, at bases covered by at least one read, ranged from 2.07 to 6.72 [24]. In addition genotypes were only called when the genotype quality was at least 20. Observed number of heterozygous SNPs per nucleotide position was then averaged for each population within the window size of 300 Kb.

Estimation of threshold values for calling sweeps

Turkey chromosomes were divided into bins of 300 Kb, and these bins were used to estimate threshold values to determine significance levels of sweep regions in the genome. Patterns of heterozygosity among these bins were investigated separately for each turkey population. A sweep region was defined when heterozygosity was below the threshold for at least 5 consecutive bins. To obtain the genome wide significance threshold (P <0.05), heterozygosity values of the bins were randomly permuted across the genome. Subsequently the threshold that would lead to exactly one significant region of 5 consecutive bins was determined for each of 7000 replicates. The distribution of these 7000 thresholds was used to obtain the 5% genome wide threshold. With this 5% threshold heterozygosity value, each population had a 5% probability of finding 1 sweep region by chance. A threshold of five consecutive bins was selected because preliminary results showed large regions of homozygosity in the turkey genome, and also to obtain stable statistics for heterozygosity. Using these threshold values, each turkey population was investigated for regions of low heterozygosity indicative of the presence of a sweep. Subsequently, turkey populations were compared with each other for the overlap in putative sweep regions. Overlapping sweep regions were identified when a sweep was replicated in more than one population. The overlapping sweep regions were defined as the genomic region covered by the sweeps from at least two populations.

Heat map plots

Heat maps for the individual turkey chromosomes and for the whole turkey genome, including all turkey autosomes, were plotted separately to visualize overlapping sweeps in different turkey populations using the “heatmap.plus” package in R [25]. The color scale was based on the square root of heterozygosity values, for visualization and distinction of sweep areas in the genomic regions.

Functional annotation analysis

All genes lying within the overlapping sweep regions of turkey were used for functional annotation analysis. Functional annotation analysis was performed using DAVID (Database for Annotation, Visualization, and Integrated Discovery) with default parameters [26]. DAVID is a web-based bioinformatics application that systematically identifies enriched biology associated with large gene list(s) derived from high-throughput genomic experiments [26]. Correction for multiple comparison was done by the Benjamini-Hochberg method [27]. Annotation for turkey and chicken genes is very limited; therefore we used one to one orthologs of turkey to human to perform this functional annotation analysis.

Ethical approval for the use of animals in this study

Although animals were used in this work, no direct experiments were performed on them. Blood sample collection was carried out by highly skilled and experienced personnel from the breeding companies. No approval from the ethics committee was necessary according to local legislation.

Results

In order to identify candidate selective sweeps, threshold values were estimated for heterozygosity in each of the different turkey populations. These threshold heterozygosity values ranged from 1.0E-5 to 5.1E-5 (Table 1). The highest threshold value was obtained for the L3 commercial line while the lowest threshold value was obtained for BvSW.

Table 1 Estimates of heterozygosity threshold (P ≤0.05) of analyzed turkey populations

A whole genome view of the sweep regions in the different turkey populations is presented in 1. In total, we observed 54 genomic regions where heterozygosity was significantly reduced (P <0.05). These candidate selective sweeps were found on 14 different chromosomes across turkey populations (Additional file 1). Areas with evidence of candidate selective sweeps varied from 1.5 Mb to 11.1 Mb in length (Additional file 1). Out of these 54 sweep regions, 31 were population-specific (Additional file 1) and observed on 12 different turkey chromosomes, while 23 were overlapping sweep regions that were observed in multiple populations and distributed across 13 different chromosomes (Table 2 & Additional file 1).

Figure 1
figure 1

Turkey chromosomes showing sweep regions shared by different turkey populations. A) Turkey autosomes (1-30) showing variation in pattern of heterozygosity, colour pattern from white to red indicating a low to high level of heterozygosity. B) Turkey chromosome 5 with 2 candidate selective sweep regions from positions 7.8-8.4 Mb and 41.1-42.0 Mb shared by 5 different turkey populations (L1, L4, L6, BvSW, RP and L1, L3, L5, L6, L7 respectively). C) Turkey chromosome 7 with sweep region from positions 9.9-11.7 Mb shared by 4 different turkey populations (L1, L4, BvSW and RP). D) Turkey chromosome 9 with candidate selective sweep region from positions 17.4-18.6 Mb shared by 5 different turkey populations (L3, L5, L6, Nset and BvSW). E) Turkey chromosome 14 with sweep region from positions 3.3-4.5 Mb shared by 4 different turkey populations (L1, L3, L6 and BvSW).

Table 2 Turkey candidate selective sweeps showing overlap in multiple turkey populations

The majority of the population-specific regions, 26 in total, were observed in the commercial populations (L1-L7), on average nearly 4 per population; whereas heritage populations (BvSW, Nset and RP) showed 1.6 population-specific sweeps per population. Differences between commercial populations were considerable, with as many as 8 sweep regions observed in population L3 and only one population-specific sweep region observed in population L6. Five population-specific sweep regions were observed in heritage varieties with 1 (RP) or 2 (BvSW and Nset) sweeps per population.

Out of 23 sweep regions that showed overlap in multiple populations, one was observed only in the heritage varieties (Nset and RP) while 13 were observed only in the commercial lines (Table 2). Commercial line L1 had the largest sweep region, 11.1 Mb, (Additional file 1) as well as the highest number (10) of overlapping sweep regions. The lowest number (3) of overlapping sweep regions was observed in the heritage variety Nset (Table 2).

Differences were observed along the turkey genome, regarding the presence of sweeps at different chromosomes. Out of 54 observed sweep regions at different turkey chromosomes, chromosome 2 showed the highest number of significant regions, 8 in total, while chromosome 14 showed the lowest number, 2 in total. Chromosomes 5, 7, 9 and 14 had five candidate selective sweep regions that showed an overlap in at least 4 different turkey populations (Table 2; Figure 1). Chromosome 5 had two overlapping sweep regions that were each shared by at least five populations, and one of these two regions was presented by commercial lines only (Table 2). Chromosome 9 also had a sweep region that was shared by five populations (Table 2 and Figure 1).

Overlapping sweep regions covered 5,452 genes, 34.7% of the total number of genes that were identified in turkey genome sequence [17]. BioMart website version 0.7 (http://www.biomart.org) was used to identify human orthologs for turkey genes. Out of these turkey genes, 3,858 were one to one orthologs with human genes and 3,832 turkey genes had a corresponding HUGO Gene Nomenclature Committee (HGNC) symbol in human genebuild (GRC37.p7). Finally, 3,718 of these genes with HGNC symbol had annotation information available in DAVID and were used in the functional annotation analysis. Functional annotation analyses resulted in 514 gene ontology (GO) terms with an Expression Analysis Systematic Explorer (EASE) P-value [28] of less than 0.1 (Additional file 2) which is a rather liberal threshold because it does not correct for multiple testing. The EASE P-value is a modified Fisher Exact P-value. GO terms that passed the significant threshold of 0.05 after Benjamini Hochberg (B-H) correction [27] are shown in Table 3. The most enriched (B-H corrected P <0.0005) was embryonic morphogenesis, while the other terms in Table 3 are related to nucleic acid binding. The nominally significant GO terms (P <0.10, Additional file 2) included a few more terms related with morphogenesis or growth but were not significant after B-H correction.

Table 3 Gene ontology (GO) terms that passed significant threshold of 0.05 after Benjamini Hochberg correction

Discussion

We aimed at finding genomic regions with reduced heterozygosity, either resulting from strong selection in favor of specific alleles or other reasons such as genetic drift. For the discovery of these regions in different turkey populations (commercial lines and the heritage varieties), we used a modified whole genome heterozygosity distribution approach [2]. In a particular population, the occurrence of heterozygosity values equal or less than the threshold value (Table 1) within at least 5 consecutive bins (each with 300 Kb size) indicates a significant reduction in heterozygosity in that region. Use of large window size might have limited our access to highlight smaller significant sweep regions. This large window size was chosen due to the detection of a large number of consecutive sweep deserts in our preliminary analyses which might be due to species specific low heterozygosity and/or overall low sequence depth [24]. In general, heterozygosity in turkey is low with an estimated average heterozygosity of 1.07 SNPs Kb-1 [24], much lower than the observed heterozygosity in chicken, with 4.28 and 2.24 SNPs Kb-1 reported in two different studies [2],[29]. We estimated threshold values separately for each turkey population. The threshold values (Table 1) can also be regarded as a measure of the level of genetic diversity in a particular population. In our study, we found the highest threshold value for commercial population L3, which is concordant with the highest observed genetic diversity and the highest number of SNPs discovered in this population in our previous study [24]. Similarly, the lowest threshold value was observed for BvSW, also concordant with the previously observed lowest genetic diversity and the lowest number of SNPs discovered in this population [24].

In our study, 48 significant regions (population-specific and overlapping) were observed in the commercial populations, while only 6 significant regions (population-specific and overlapping) were observed in the heritage populations (Additional file 1 & Table 2). The small number of individuals (2-3) used per population could not reveal the complete variation of a particular population but each of these individuals still belonged to a specific population, therefore population specific terminology was used for the group of individuals that belong to a same population. The high number of candidate selective sweeps in commercial lines can be explained as a result of the high selection intensities applied to these populations [30]. A lower number of sweep regions in heritage varieties may be due to a number of reasons, such as the admixture of populations, relatively high effective populations size in heritage varieties, or relatively less intensive and less specific directional selection applied to these varieties in comparison to commercial turkeys. Specific information about population admixture or effective population size of heritage varieties is limited, but these varieties were likely pure lines given the anecdotal information from the turkey breeders.

In our previous study, among the heritage varieties, Nset showed the highest heterozygosity followed by RP and BvSW respectively [24]. A consistent pattern was observed with a lower number of sweep regions and a higher threshold heterozygosity value for Nset compared to BvSW and RP. These differences in number of sweeps and threshold heterozygosity values for the different populations may also be an indication of difference in level of admixture or effective population size. The heritage variety BvSW showed the lowest threshold heterozygosity value and also the highest number of sweeps of all heritage varieties, which is consistent with the severe bottleneck that this population went through in 2000 (Alexandra Scupham, Personal communication). Similarly, Nset population showed highest threshold heterozygosity value and the lowest number of sweeps of all domesticated turkey populations which could represent a higher level of admixture or comparatively larger effective population size for this population. However, no historical information is available to support this.

Regions of sweep with variable but large sizes (1.5-11.1 Mb) were observed. Reduction in genetic diversity/heterozygosity at different locations in the genome can persist for a long time, and indicate selection across a long genomic region [31]. The size of a sweep region may vary with history of domestication, the type of population (inbred or outbred), intensity of selection within a particular population, population dynamics such as bottlenecks and drift. SNP analyses of domestic dogs and cats show large stretches of alternating heterozygous and homozygous regions as a consequence of domestication and breed development [32],[33]. In most outbred species, a selected region would display local SNP homozygosity, compared to abundant polymorphism elsewhere in the genome [34].

Uneven distribution of homozygous regions can be expected across the genome due to selection pressure through natural or artificial means [1]-[3],[35]. Chromosome 5, 7, 9 and 14 are highly distinct with overlapping regions in at least four different turkey populations (Table 2). This suggests that genomic regions on these chromosomes contain gene(s) which affect the traits that are important for turkey production. Turkey populations that showed overlap in sweeps on these chromosomes may be selected for specific objectives that all populations had in common or, alternatively, may have been developed from parents that already were homozygous for these sweep regions. Two candidate selective sweep regions discovered on chromosome 5 and chromosome 22 show overlapping stretches only in commercial populations (Additional file 1). These regions may contain genes involved in commercially important traits. The regions, however, are too large to identify the individual genes that may have been selected.

We could not use museum samples (South Mexican turkeys) in our current data which were included in our previous study [24] due to their very low available sequence depth. Average sequence depth at bases covered by at least one read in museum samples ranged from 1.38 to 1.81 [24] which is less than half the depth (5 to 10 -fold) that was used as the criterion for calling SNPs in all individuals of the current study. However, even though coverage was low, in our previous study [24] we identified genomic regions at four chromosomes with increased homozygosity of non-reference alleles in the museum samples. The domesticated populations were found to be fixed for the reference alleles at those same loci [24]. These genomic regions with high non-reference allele homozygosity were aligned with the candidate selective sweep regions of current study to find any overlap. Besides the region at chromosome 3, the other regions at chromosomes 9, 14 and 22 showed overlap with the detected sweep regions (Additional file 1) of current study. These sweep regions of chromosome 9, 14 and 22 show overlap in 5, 4 and 3 populations respectively. This concordance of results supports our hypothesis that these candidate sweep regions are likely result of selection in commercial populations.

Chromosome studies have revealed that the karyotype is more conserved among avian species than in other taxa, such as mammals, with most avian species showing a diploid chromosome number between 76 and 80 (http://www.genomesize.com). This suggests that chromosomal evolution or large-scale rearrangements affecting chromosome number occur at a low rate in birds, and as a result many chromosomes have remained more or less intact during avian evolution [36]. Comparative cytogenetic and linkage maps between turkey and chicken showed conserved synteny and close ancestral relationships [37],[38] that support the hypothetical ancestral Galliform karyotype [39]. The strong structural and functional conservation between the turkey and the chicken genomes [40],[41], as well as the similarities in breeding objectives, suggest that overlap in selective sweep regions between the two species could be expected. To test whether selective sweep regions are conserved between chicken and turkey, the orthology to chicken for all significant overlapping sweep regions of turkey was determined. These genomic regions were then examined for the presence of sweeps, based on two different studies in the chicken [2],[42]. Selective sweep studies reported about 400 sweep regions [2],[42] which is about 0.38 sweep per Mb in chicken genome. Thirteen out of the 23 overlapping candidate selective sweep regions identified in turkey also harbored a selective sweep reported in chicken. Rubin et al. [2] reported 40 highly significant chicken sweep regions with very low Z transformed heterozygosity (ZHp < -6). Two of these highly significant chicken sweeps mapped within the syntenic regions of turkey sweeps on chromosomes 7 and 11 (Additional file 1). Overall, the concordance of chicken sweep regions with turkey sweep regions was low. Approximately 0.32 chicken sweeps were observed per Mb within the total overlapping sweep length of turkey. This result shows no enrichment of chicken sweeps within the overlapping sweep regions of turkey.

Selective sweep regions are expected to have been involved in producing phenotypic variation for the traits of interest and intensive selection leads these regions towards fixation. To investigate the variation explained by these regions, we looked for available turkey quantitative trait loci (QTL) information within these regions [41]. We did not find overlap between the QTL regions from our previous study [41] and the candidate sweep regions in the current analyses. This discordance could be explained if QTL regions were still too much variable to be identified in a search for selective sweeps. Due to the limited availability of information on turkey QTLs and the presence of structural and functional conservation between the turkey and chicken genomes [24],[38],[40], overlapping regions of candidate selective sweeps (Table 2) of turkey were aligned with chicken genome sequence (WASHUC2) to determine their positions in the chicken genome (Additional file 3). The orthologous chicken regions were subsequently examined for the presence of reported chicken QTL for growth [43]. Many QTL were found to be overlapped with these genomic regions (Additional file 3). The frequency of chicken growth QTL for which the confidence interval overlapped with the turkey sweep regions was found to be 11.33 growth QTL per Mb of sweep region. This high frequency of chicken growth QTL overlapping with the turkey candidate selective sweep regions was however a result of the high number of growth QTL discovered in chicken. The sweep regions did not show an enrichment of chicken QTL compared to other parts of the genome.

Production censuses of turkeys from the last four decades show that turkeys have doubled in size. We had therefore expected to see a sweep in the region of the somatomedin, insulin-like growth factor 1 (IGF-1), which is well known to play an important role in muscle growth and development in various domesticated species [44]-[46]. However, we did not find a candidate sweep near the IGF-1 region on turkey chromosome 1 (56348061 bp-56402610 bp). This observation suggests that the sequence variation at the IGF-1 locus itself is not involved in regulating the level in turkeys. Previously, two QTL were detected for IGF-1 level in blood plasma in chicken at chromosome 1 and 2 [47],[48]. These two chicken QTL regions, both are syntenic with turkey overlapping candidate sweep regions at chromosome 1 and 6, respectively (Additional file 3), showing that some genes are present within the candidate sweep regions that appear to affect the level of IGF-1 hormone in blood, which has been shown to regulate growth, reproduction, energy balance, cell proliferation and cell death [49].

Given the large increase in production per bird from 6.7 to 12.7 Kg in a 40 years period [13], intensive selection for growth must have taken place in turkeys. The likely candidate genes such as IGF2, Pit1, AFABP, PRKAG3, GDF8 etc. that have been previously reported to affect growth were not present within the candidate sweep regions. Gene ontology (GO) enrichment analysis was therefore performed to see if the complete set of genes within the candidate sweep regions has been enriched for association with growth. We performed gene functional annotation analysis using DAVID. Gene-based enrichment analysis showed some enrichment of genes for regulation of development and morphogenesis within the candidate sweep regions (Additional file 2). We found significantly (Benjamini P <0.05) enriched GO term with embryonic morphogenesis (Table 3) and other suggestive terms (1 < P >0.05) with embryonic organ morphogenesis, body development, maintenance of growth etc. (Additional file 2). This shows that the observed candidate selective sweep regions of turkey are enriched with genes that are important for some factors in growth and development.

Conclusion

The turkey genome showed large candidate sweep regions. The relatively high number of candidate selective sweep regions in commercial turkey populations compared to heritage varieties provided evidence of intense selection in these commercial lines. In addition, the enrichment of these candidate sweep regions with genes of importance to growth indicates that these regions may have been targets of selection for growth in these commercial lines, moving specific haplotypes towards fixation.

Availability of supporting data

BAM files of all the individuals used in current study are available online at NCBI's Sequence Read Archive (SRA) database under the study accession number “sRP012021” [http://www.ncbi.nlm.nih.gov/Traces/sra/?study=SRP012021].

Additional files

References

  1. Qanbari S, Pimentel ECG, Tetens J, Thaller G, Lichtner P, Sharifi AR, Simianer H: A genome-wide scan for signatures of recent selection in Holstein cattle. Anim Genet. 2010, 41: 377-389.

    PubMed  CAS  Google Scholar 

  2. Rubin C-J, Zody MC, Eriksson J, Meadows JRS, Sherwood E, Webster MT, Jiang L, Ingman M, Sharpe T, Ka S, Hallbook F, Besnier F, Carlborg O, Bed'hom B, Tixier-Boichard M, Jensen P, Siegel P, Lindblad-Toh K, Andersson L: Whole-genome resequencing reveals loci under selection during chicken domestication. Nature. 2010, 464 (7288): 587-591. 10.1038/nature08832.

    Article  PubMed  CAS  Google Scholar 

  3. Qanbari S, Gianola D, Hayes B, Schenkel F, Miller S, Moore S, Thaller G, Simianer H: Application of site and haplotype-frequency based approaches for detecting selection signatures in cattle. BMC Genomics. 2011, 12: 318-10.1186/1471-2164-12-318.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Nielsen R, Bustamante C: A scan for positively selected genes in the genomes of humans and chimpanzees. PloS Biology. 2005, 3 (6): e170-10.1371/journal.pbio.0030170.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Lettre G, Jackson AU, Gieger C, Schumacher FR, Berndt SI, Sanna S, Eyheramendy S, Voight BF, Butler JL, Guiducci C, Illig T, Hackett R, Heid IM, Jacobs KB, Lyssenko V, Uda M, Boehnke M, Chanock SJ, Groop LC, Hu FB, Isomaa B, Kraft P, Peltonen L, Salomaa V, Schlessinger D, Hunter DJ, Hayes RB, Abecasis GR, Wichmann HE, Mohlke KL, et al: Identification of ten loci associated with height highlights new biological pathways in human growth. Nat Genet. 2008, 40 (5): 584-591. 10.1038/ng.125.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  6. Vaysse A, Ratnakumar A, Derrien T, Axelsson E, Rosengren Pielberg G, Sigurdsson S, Fall T, Seppälä EH, Hansen MST, Lawley CT, Karlsson EK, Bannasch D, Vilà C, Lohi H, Galibert F, Fredholm M, Häggström J, Hedhammar A, André C, Lindblad-Toh K, Hitte C, Webster MT:: Identification of genomic regions associated with phenotypic variation between dog breeds using selection mapping. PLoS Genet. 2011, 7 (10): e1002316-10.1371/journal.pgen.1002316.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  7. Howard R, Moore A: A Complete Checklist of Birds of the World. 1984, Macmillan, London, Revised

    Google Scholar 

  8. Speller CF, Kemp BM, Wyatt SD, Monroe C, Lipe WD, Arndt UM, Yang DY: Ancient mitochondrial DNA analysis reveals complexity of indigenous North American turkey domestication. Proc Natl Acad Sci U S A. 2010, 107 (7): 2807-2812. 10.1073/pnas.0909724107.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  9. Hirst KK: Turkey ( Melagris gallapavo and spp.). In History of the Domestication of Turkeys. United States. About.com Education Archaeology. Accessed October 30, 2014. [http://archaeology.about.com/od/tterms/qt/turkey.htm]

  10. Flannery KV: Vertebrate fauna and hunting practices. Prehistory of the Tehuacan Valley. Environment and Subsistence. Edited by: Byers DS. 1967, Univ of Texas Press, Austin, 132-177.

    Google Scholar 

  11. Schorger AW: The Wild Turkey: Its History and Domestication. 1966, Univ of Oklahoma Press, Norman

    Google Scholar 

  12. The American Standard of Perfection. 2001, APA, Troy, NY

  13. Food and agriculture organization statistical division (FAOSTAT) of the United Nations. [http://faostat.fao.org/]

  14. Maynard Smith J, Haigh J: The hitchhiking effect of a favorable gene. Genet Res. 1974, 23: 23-35. 10.1017/S0016672300014634.

    Article  Google Scholar 

  15. Kijas JM, Townley D, Dalrymple BP, Heaton MP, Maddox JF, McGrath A, Wilson P, Ingersoll RG, McCulloch R, McWilliam S, Tang D, McEwan J, Cockett N, Oddy VH, Nicholas FW, Raadsma H: A genome wide survey of SNP variation reveals the genetic structure of sheep breeds. PLoS One. 2009, 4 (3): e4668-10.1371/journal.pone.0004668.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, Subgroup GPDP: The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009, 25 (16): 2078-2079. 10.1093/bioinformatics/btp352.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Dalloul RA, Long JA, Zimin AV, Aslam L, Beal K, Ann Blomberg L, Bouffard P, Burt DW, Crasta O, Crooijmans RPMA, Cooper K, Coulombe RA, De S, Delany ME, Dodgson JB, Dong JJ, Evans C, Frederickson KM, Flicek P, Florea L, Folkerts O, Groenen MAM, Harkins TT, Herrero J, Hoffmann S, Megens H-J, Jiang A, de Jong P, Kaiser P, Kim H, et al: Multi-platform next-generation sequencing of the domestic turkey (Meleagris gallopavo): genome assembly and analysis. PLoS Biol. 2010, 8 (9): e1000475-10.1371/journal.pbio.1000475.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Kerstens HHD, Crooijmans RPMA, Veenendaal A, Dibbits BW, Chin-A-Woeng TFC, Dunnen JTD, Groenen MAM: Large scale single nucleotide polymorphism discovery in unsequenced genomes using second generation high throughput sequencing technology: applied to turkey. BMC Genomics. 2009, 10: 479-10.1186/1471-2164-10-479.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Li G, Ma L, Song C, Yang Z, Wang X, Huang H, Li Y, Li R, Zhang X, Yang H, Wang J, Wang J: The YH database: the first Asian diploid genome database. Nucleic Acids Res. 2009, 37: 1025-1028. 10.1093/nar/gkn966.

    Article  Google Scholar 

  20. Marsden SJ: The Beltsville small white turkey. World’s Poult Sci J. 1967, 23 (1): 32-42. 10.1079/WPS19670008.

    Article  Google Scholar 

  21. The Livestock Conservancy: Narragansett Turkey. North Carolina, USA; Accessed October 30, 2014. [http://www.albc-usa.org/cpl/narragansett.html]

  22. The Livestock Conservancy: Royal Palm Turkey. North Carolina, USA; Accessed October 30, 2014. [http://albc-usa.org/cpl/royalpalm.html]

  23. Stromberg M: Mosaik Assembler. In 110014. Edited by Lee W-P. Boston College; 2010. [https://code.google.com/p/mosaik-aligner/downloads/list]

  24. Aslam ML, Bastiaansen JWM, Elferink M, Megens H-J, Crooijmans R, Long J, Groenen MAM: Whole genome SNP discovery and the genetic diversity of different turkey populations. BMC Genomics. 2012, 13: 391-10.1186/1471-2164-13-391.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  25. Day A: heatmap.plus: Heatmap with more sensible behavior. In R Package Version 1.3. Vienna; 2007. [http://cran.r-project.org/web/packages/heatmap.plus/index.html]

  26. Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2008, 4 (1): 44-57. 10.1038/nprot.2008.211.

    Article  Google Scholar 

  27. Benjamini Y, Hochberg Y: Controlling the false discovery rate - a practical and powerful approach to multiple testing. J R Stat Soc Series B Stat Methodol. 1995, 57: 289-300.

    Google Scholar 

  28. Hosack D, Dennis G, Sherman B, Lane H, Lempicki R: Identifying biological themes within lists of genes with EASE. Genome Biol. 2003, 4 (10): R70-10.1186/gb-2003-4-10-r70.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Wong GK-S, Liu B, Wang J, Zhang Y, Yang X, Zhang Z, Meng Q, Zhou J, Li D, Zhang J, Ni P, Li S, Ran L, Li H, Zhang J, Li R, Li S, Zheng H, Lin W, Li G, Wang X, Zhao W: A genetic variation map for chicken with 2.8 million single-nucleotide polymorphisms. Nature. 2004, 432 (7018): 717-722. 10.1038/nature03156.

    Article  PubMed  CAS  Google Scholar 

  30. McKay JC: The genetics of modern commercial poultry. In In proceedings of 23 rdWorld’s Poultry Congress: 2008; Brisbane, Australia. 2008. CD-ROM.,

  31. Initial sequence of the chimpanzee genome and comparison with the human genome. Nature. 2005, 437 (7055): 69-87. 10.1038/nature04072.

  32. Lindblad-Toh K, Wade CM, Mikkelsen TS, Karlsson EK, Jaffe DB, Kamal M, Clamp M, Chang JL, Kulbokas EJ, Zody MC, Mauceli E, Xie X, Breen M, Wayne RK, Ostrander EA, Ponting CP, Galibert F, Smith DR, de Jong PJ, Kirkness E, Alvarez P, Biagi T, Brockman W, Butler J, Chin C-W, Cook A, Cuff J, Daly MJ, DeCaprio D, Gnerre S, et al: Genome sequence, comparative analysis and haplotype structure of the domestic dog. Nature. 2005, 438 (7069): 803-819. 10.1038/nature04338.

    Article  PubMed  CAS  Google Scholar 

  33. Pontius JU, Mullikin JC, Smith DR, Lindblad-Toh K, Gnerre S, Clamp M, Chang J, Stephens R, Neelam B: Initial sequence and comparative analysis of the cat genome. Genome Res. 2007, 17: 1675-1689. 10.1101/gr.6380007.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  34. Oleksyk TK, Zhao K, De La Vega FM, Gilbert DA, O'Brien SJ, Smith MW: Identifying selected regions from heterozygosity and divergence using a light-coverage genomic dataset from two human populations. PLoS One. 2008, 3: e1712-10.1371/journal.pone.0001712.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Davidson S, Starkey A, MacKenzie A: Evidence of uneven selective pressure on different subsets of the conserved human genome; implications for the significance of intronic and intergenic DNA. BMC Genomics. 2009, 10: 614-10.1186/1471-2164-10-614.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Nie W, O’Brien PCM, Ng BL, Fu B, Volobouev V, Carter NP, Ferguson-Smith MA, Yang F: Avian comparative genomics: reciprocal chromosome painting between domestic chicken (Gallus gallus) and the stone curlew (Burhinus oedicnemus, Charadriiformes)-An atypical species with low diploid number. Chromosome Res. 2009, 17 (1): 99-113. 10.1007/s10577-009-9021-6.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  37. Chaves LD, Rowe JA, Reed KM: Survey of a cDNA library from the turkey (Meleagris gallopavo). Genome. 2005, 48: 12-17. 10.1139/g04-088.

    Article  PubMed  CAS  Google Scholar 

  38. Griffin D, Robertson L, Tempest H, Vignal A, Fillon V, Crooijmans R, Groenen M, Deryusheva S, Gaginskaya E, Carre W, Waddington D, Talbot R, Volker M, Masabanda J, Burt D: Whole genome comparative studies between chicken and turkey and their implications for avian genome evolution. BMC Genomics. 2008, 9 (1): 168-10.1186/1471-2164-9-168.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Griffin DK, Robertson LBW, Tempest HG, Skinner BM: The evolution of the avian genome as revealed by molecular cytogenetics. Cytogenet Genome Res. 2007, 117: 64-77. 10.1159/000103166.

    Article  PubMed  CAS  Google Scholar 

  40. Aslam ML, Bastiaansen JWM, Crooijmans RPMA, Vereijken A, Megens H-J, Groenen MAM: A SNP based linkage map of the turkey genome reveals multiple intrachromosomal rearrangements between the turkey and chicken genomes. BMC Genomics. 2010, 11: 647-10.1186/1471-2164-11-647.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Aslam ML, Bastiaansen JWM, Crooijmans RPMA, Vereijken A, Groenen MAM: Whole genome QTL mapping for growth, meat quality and breast meat yield traits in turkey. BMC Genet. 2011, 12: 61-10.1186/1471-2156-12-61.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Elferink MG, Megens H-J, Vereijken A, Hu X, Crooijmans RPMA, Groenen MAM: Signatures of Selection in the Genomes of Commercial and Non-Commercial Chicken Breeds. PLoS One. 2012, 7 (2): e32720-10.1371/journal.pone.0032720.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  43. Hu Z-L, Reecy JM: Animal QTLdb: beyond a repository - a public platform for QTL comparisons and integration with diverse types of structural genomic information. Mamm Genome. 2007, 18: 1-4. 10.1007/s00335-006-0105-8.

    Article  PubMed  Google Scholar 

  44. Klindt J, Yen JT, Buonomo FC, Roberts AJ, Wise T: Growth, body composition, and endocrine responses to chronic administration of insulin-like growth factor I and(or) porcine growth hormone in pigs. J Anim Sci. 1998, 76 (9): 2368-2381.

    PubMed  CAS  Google Scholar 

  45. Wood BJ, Archer JA, van der Werf JHJ: Response to selection in beef cattle using IGF-1 as a selection criterion for residual feed intake under different Australian breeding objectives. Liv Prod Sci. 2004, 91 (1-2): 69-81. 10.1016/j.livprodsci.2004.06.009.

    Article  Google Scholar 

  46. Sato S, Ohtake T, Uemoto Y, Okumura Y, Kobayashi E: Polymorphism of insulin-like growth factor 1 gene is associated with breast muscle yields in chickens. Anim Sci J. 2012, 83 (1): 1-6. 10.1111/j.1740-0929.2011.00917.x.

    Article  PubMed  CAS  Google Scholar 

  47. Park H-B, Jacobsson L, Wahlberg P, Siegel PB, Andersson L: QTL analysis of body composition and metabolic traits in an intercross between chicken lines divergently selected for growth. Physiol Genomics. 2006, 25 (2): 216-223. 10.1152/physiolgenomics.00113.2005.

    Article  PubMed  CAS  Google Scholar 

  48. Nadaf J, Pitel F, Gilbert H, Duclos MJ, Vignoles F, Beaumont C, Vignal A, Porter TE, Cogburn LA, Aggrey SE, Simon J, Le Bihan-Duval E: QTL for several metabolic traits map to loci controlling growth and body composition in an F2 intercross between high- and low-growth chicken lines. Physiol Genomics. 2009, 38 (3): 241-249. 10.1152/physiolgenomics.90384.2008.

    Article  PubMed  CAS  Google Scholar 

  49. McMurtry JP, Francis GL, Upton FZ: Insulinlike growth factors in poultry. Domest Anim Endocrinol. 1997, 14: 199-229. 10.1016/S0739-7240(97)00019-2.

    Article  PubMed  CAS  Google Scholar 

Download references

Acknowledgements

This project was supported by Agriculture and Food Research Initiative Competitive Grant no. 2010-65205-20428 from the USDA National Institute of Food and Agriculture. The authors thank Hybrid Turkeys, a division of Hendrix Genetics, and Aviagen Turkeys for supplying blood samples from pedigree turkey lines. The authors thank the USDA’s National Animal Disease Center for providing blood samples from Beltsville Small White turkeys, and the Department of Animal and Poultry Sciences at Virginia Polytechnic Institute for providing blood samples from Narragansett and Royal Palm heritage varieties. The authors thank Timothy L. Conn, and Lori Schreier for DNA extraction and library preparation, as well as Alicia Beavers for library sequencing.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Muhammad L Aslam.

Additional information

Competing interest

The authors declare that they have no competing interests.

Authors’ contributions

MLA, JWMB, and HJM analysed the data. RPMAC and FN helped in improving manuscript by their useful suggestions on data analyses and write-up. LAB assisted and trained laboratory personnel in genomic DNA isolation from blood and preparation of all libraries. CPVT assisted with project design and coordination. TSS helped develop sequencing strategy and supervised sequence production. SGS was responsible for sequencing processing. JAL organized blood sample collection and DNA sequencing, and was the Principal Investigator (Agriculture and Food Research Initiative Competitive Grant no. 2010-65205-20428). MLA wrote the paper and all other authors gave suggestions and comments for the improvement of paper. All authors read and approved the final manuscript. Overall coordination of the project was by JWMB, JAL and MAMG.

Electronic supplementary material

12863_2014_117_MOESM1_ESM.xlsx

Additional file 1: Position of turkey genomic regions with candidate selective sweeps in different turkey populations. Description: This file contains the start and the end positions of turkey sweep regions on different chromosomes for each of the populations.This file also shows the syntenic selective sweep regions in chicken that coincide with turkey overlapping sweeps. (XLSX 20 KB)

12863_2014_117_MOESM2_ESM.xlsx

Additional file 2: Gene ontology (GO) terms observed with functional annotation analysis performed using the online tool DAVID. Description: This file contains all GO terms included in the functional annotation analysis with their biological functions and enrichment P-values. This file also contains the names of genes in the sweep regions that are annotated with each of the GO terms. (XLSX 118 KB)

12863_2014_117_MOESM3_ESM.xlsx

Additional file 3: Syntenic positions of turkey overlapping sweep regions with chicken. Description: This file contains the start and the end positions of turkey overlapping sweep regions and their syntenic positions in chicken genome. This file also contains chicken growth QTL for which the confidence interval overlapped with the turkey sweep regions. (XLSX 40 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Rights and permissions

This article is published under license to BioMed Central Ltd. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Aslam, M.L., Bastiaansen, J.W., Megens, HJ. et al. Genome-wide candidate regions for selective sweeps revealed through massive parallel sequencing of DNA across ten turkey populations. BMC Genet 15, 117 (2014). https://doi.org/10.1186/s12863-014-0117-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12863-014-0117-4

Keywords