The easy road to genome‐wide medium density SNP screening in a non‐model species: development and application of a 10 K SNP‐chip for the house sparrow (Passer domesticus)

With the advent of next generation sequencing, new avenues have opened to study genomics in wild populations of non‐model species. Here, we describe a successful approach to a genome‐wide medium density Single Nucleotide Polymorphism (SNP) panel in a non‐model species, the house sparrow (Passer domesticus), through the development of a 10 K Illumina iSelect HD BeadChip. Genomic DNA and cDNA derived from six individuals were sequenced on a 454 GS FLX system and generated a total of 1.2 million sequences, in which SNPs were detected. As no reference genome exists for the house sparrow, we used the zebra finch (Taeniopygia guttata) reference genome to determine the most likely position of each SNP. The 10 000 SNPs on the SNP‐chip were selected to be distributed evenly across 31 chromosomes, giving on average one SNP per 100 000 bp. The SNP‐chip was screened across 1968 individual house sparrows from four island populations. Of the original 10 000 SNPs, 7413 were found to be variable, and 99% of these SNPs were successfully called in at least 93% of all individuals. We used the SNP‐chip to demonstrate the ability of such genome‐wide marker data to detect population sub‐division, and compared these results to similar analyses using microsatellites. The SNP‐chip will be used to map Quantitative Trait Loci (QTL) for fitness‐related phenotypic traits in natural populations.     

A high‐density linkage map enables a second‐generation collared flycatcher genome assembly and reveals the patterns of avian recombination rate variation and chromosomal evolution

Detailed linkage and recombination rate maps are necessary to use the full potential of genome sequencing and population genomic analyses. We used a custom collared flycatcher 50 K SNP array to develop a high‐density linkage map with 37 262 markers assigned to 34 linkage groups in 33 autosomes and the Z chromosome. The best‐order map contained 4215 markers, with a total distance of 3132 cM and a mean genetic distance between markers of 0.12 cM . Facilitated by the array being designed to include markers from most scaffolds, we obtained a second‐generation assembly of the flycatcher genome that approaches full chromosome sequences (N50 super‐scaffold size 20.2 Mb and with 1.042 Gb (of 1.116 Gb) anchored to and mostly ordered and oriented along chromosomes). We found that flycatcher and zebra finch chromosomes are entirely syntenic but that inversions at mean rates of 1.5–2.0 event (6.6–7.5 Mb) per My have changed the organization within chromosomes, rates high enough for inversions to potentially have been involved with many speciation events during avian evolution. The mean recombination rate was 3.1 cM /Mb and correlated closely with chromosome size, from 2 cM /Mb for chromosomes >100 Mb to >10 cM /Mb for chromosomes <10 Mb. This size dependence seemed entirely due to an obligate recombination event per chromosome; if 50 cM was subtracted from the genetic lengths of chromosomes, the rate per physical unit DNA was constant across chromosomes. Flycatcher recombination rate showed similar variation along chromosomes as chicken but lacked the large interior recombination deserts characteristic of zebra finch chromosomes.     

The effect of single nucleotide polymorphism identification strategies on estimates of linkage disequilibrium

At present there is tremendous interest in characterizing the magnitude and distribution of linkage disequilibrium (LD) throughout the human genome, which will provide the necessary foundation for genome-wide LD analyses and facilitate detailed evolutionary studies. To this end, a human high-density single-nucleotide polymorphism (SNP) marker map has been constructed. Many of the SNPs on this map, however, were identified by sampling a small number of chromosomes from a single population, and inferences drawn from studies using such SNPs may be influenced by ascertainment bias (AB). Through extensive simulations, we have found that AB is a potentially significant problem in estimating and comparing LD within and between populations. Specifically, the magnitude of AB is a function of the SNP discovery strategy, number of chromosomes used for SNP discovery, population genetic characteristics of the particular genomic region considered, amount of gene flow between populations, and demographic history of the populations. We demonstrate that a balanced SNP discovery strategy (where equal numbers of chromosomes are sampled from multiple subpopulations) is the optimal study design for generating broadly applicable SNP resources. Finally, we validate our theoretical predictions by comparing our results to publicly available data from ten genes sequenced in 24 African American and 23 European American individuals.     

The identification of SNPs with indeterminate positions using the Equine SNP50 BeadChip

We have used linkage disequilibrium (LD) to identify single nucleotide polymorphisms (SNPs) on the Illumina Equine SNP50 BeadChip, which may be incorrectly positioned on the genome map. A total of 1201 Thoroughbred horses were genotyped using the Illumina Equine SNP50 BeadChip. LD was evaluated in a pairwise fashion between all autosomal SNPs, both within and across chromosomes. Filters were then applied to the data, firstly to identify SNPs that may have been mapped to the wrong chromosome and secondly to identify SNPs that may have been incorrectly positioned within chromosomes. We identified a single SNP on ECA28, which showed low LD with neighbouring SNPs but considerable LD with a group of SNPs on ECA10. Furthermore, a cluster of SNPs on ECA5 showed unusually low LD with surrounding SNPs. A total of 39 SNPs met the criteria for unusual within-chromosome LD. The results of this study indicate that some SNPs may be misplaced. This finding is significant, as misplaced SNPs may lead to difficulties in the application of genomic methods, such as homozygosity mapping, for which SNP order is important.     

Replicated high-density genetic maps of two great tit populations reveal fine-scale genomic departures from sex-equal recombination rates

Linking variation in quantitative traits to variation in the genome is an important, but challenging task in the study of life-history evolution. Linkage maps provide a valuable tool for the unravelling of such trait−gene associations. Moreover, they give insight into recombination landscapes and between-species karyotype evolution. Here we used genotype data, generated from a 10k single-nucleotide polymorphism (SNP) chip, of over 2000 individuals to produce high-density linkage maps of the great tit (Parus major), a passerine bird that serves as a model species for ecological and evolutionary questions. We created independent maps from two distinct populations: a captive F2-cross from The Netherlands (NL) and a wild population from the United Kingdom (UK). The two maps contained 6554 SNPs in 32 linkage groups, spanning 2010 cM and 1917 cM for the NL and UK populations, respectively, and were similar in size and marker order. Subtle levels of heterochiasmy within and between chromosomes were remarkably consistent between the populations, suggesting that the local departures from sex-equal recombination rates have evolved. This key and surprising result would have been impossible to detect if only one population was mapped. A comparison with zebra finch Taeniopygia guttata, chicken Gallus gallus and the green anole lizard Anolis carolinensis genomes provided further insight into the evolution of avian karyotypes.     

Construction of a dense SNP map for bovine chromosome 6 to assist the assembly of the bovine genome sequence

A linkage map was constructed for bovine chromosome 6 (BTA6), using 399 single nucleotide polymorphisms (SNPs) detected primarily from PCR-resequencing. The efficiency of SNP detection was highly dependent on the source of sequence information chosen for primer design (BAC-end sequences, introns or promoters). The SNPs were used to build a linkage map comprising 104 cM on BTA6. The SNP order in the linkage map corresponded very well with radiation hybrid (RH) maps available for BTA6 as well as with expected positions in the human comparative map, but diverged significantly from the current assembly of the bovine genome (Btau_3.1). When performing linkage analysis with the marker order suggested from the Btau_3.1 we observed an expansion of the genetic map from 104 cM to 137 cM, strongly suggesting a reordering of scaffolds in the current version of the bovine genome assembly. The extent of LD on BTA6 was evaluated by calculating the average r ² for SNP pairs separated by given distances. The decline of LD was rapid with distance, such that r ² was 0.1 at 100 kb. Our results indicate that linkage mapping will be a valuable source of information for correcting errors in the current bovine assembly. These errors were sufficiently frequent to be of concern for the accuracy of mapping QTL with panels of SNPs whose positions are based on the current assembly.     

A first‐generation microsatellite linkage map of the ruff

A linkage map of the ruff (Philomachus pugnax) genome was constructed based on segregation analysis of 58 microsatellite loci from 381 captive‐bred individuals spanning fourteen breeding years and comprising 64 families. Twenty‐eight of the markers were resolved into seven linkage groups and five single marker loci, homologous to known chicken (Gallus gallus) and zebra finch (Taeniopygia guttata) chromosomes. Linkage groups range from 10.1 to 488.7 cM in length and covered a total map distance of 641.6 cM, corresponding to an estimated 30–35% coverage of the ruff genome, with a mean spacing of 22.9 cM between loci. Through comparative mapping, we are able to assign linkage groups Ppu1, Ppu2, Ppu6, Ppu7, Ppu10, Ppu13, and PpuZ to chromosomes and identify several intrachromosomal rearrangements between the homologs of chicken, zebra finch, and ruff microsatellite loci. This is the first linkage map created in the ruff and is a major step toward providing genomic resources for this enigmatic species. It will provide an essential framework for mapping of phenotypically and behaviorally important loci in the ruff.     

Recombination rates across porcine autosomes inferred from high-density linkage maps

Studies of the variation in recombination rate across the genome provide a better understanding of evolutionary genomics and are also an important step towards mapping and dissecting complex traits in domestic animals. With the recent completion of the porcine genome sequence and the availability of a high‐density porcine single nucleotide polymorphism (SNP) array, it is now possible to construct a high‐density porcine linkage map and estimate recombination rate across the genome. A total of 416 animals were genotyped with the Porcine SNP60BeadChip, and high‐density chromosome linkage maps were constructed using CRI‐MAP, assuming the physical order of the Sscrofa10 assembly. The total linkage map length was 2018.79 cM, using 658 meioses and 14 503 SNPs. The estimated average recombination rate across the porcine autosomes was 0.86 cM/Mb. However, a large variation in recombination rate was observed among chromosomes. The estimated average recombination rates (cM/Mb) per chromosome ranged from 0.48 in SSC1 to 1.48 in SSC10, displaying a significant negative correlation with the chromosome sizes. In addition, the analysis of the variation in the recombination rates taking 1‐Mb sliding windows has allowed us to demonstrate the variation in recombination rates within chromosomes. In general, a larger recombination rate was observed in the extremes than in the centre of the chromosome. Finally, the ratio between female and male recombination rates was also inferred, obtaining a value of 1.38, with the heterogametic sex having the least recombination.     

High density SNP mapping and QTL analysis for fruit quality characteristics in peach (Prunus persica L.)

Single nucleotide polymorphisms (SNPs) were used to construct an integrated SNP linkage map of peach (Prunus persica (L.) Batsch). A set of 1,536 SNPs were evaluated with the GoldenGate® Genotyping assay in two mapping populations, Pop-DF, and Pop-DG. After genotyping and filtering, a final set of 1,400 high quality SNPs in Pop-DF and 962 in Pop-DG with full map coverage were selected and used to construct two linkage maps with JoinMap®4.0. The Pop-DF map covered 422 cM of the peach genome and included 1,037 SNP markers, and Pop-DG map covered 369 cM and included 738 SNPs. A consensus map was constructed with 588 SNP markers placed in eight linkage groups (n?=?8 for peach), with map coverage of 454 cM and an average distance of 0.81 cM/marker site. Placements of SNPs on the “peach v1.0” physical map were compared to placement on the linkage maps and several differences were observed. Using the SNP linkage map of Pop-DG and phenotypic data collected for three harvest seasons, a QTL analysis for fruit quality traits and chilling injury symptoms was carried out with the mapped SNPs. Significant QTL effects were detected for mealiness (M) and flesh bleeding (FBL) QTLs on linkage group 4 and flesh browning (FBr) on linkage group 5. This study represents one of the first examples of QTL detection for quality traits and chilling injury symptoms using a high-density SNP map in a single peach F1 family.     

Genetic linkage maps of Populus nigra L. including AFLPs, SSRs, SNPs, and sex trait

Black poplar (Populus nigra L.) is a tree of ecological and economic interest. A better knowledge of P. nigra genome is needed for an effective protection and use of its genetic resources. The main objective of this study is the construction of a highly informative genetic map of P. nigra species including genes of adaptive and economic interest. Two genotypes originated from contrasted natural Italian populations were crossed to generate a F₁ mapping pedigree of 165 individuals. Amplification fragment length polymorphism (AFLP), simple sequence repeat (SSR), and single nucleotide polymorphism (SNP) markers were used to genotype 92 F₁ individuals, and the pseudo-test-cross strategy was applied for linkage analysis. The female parent map included 368 markers (274 AFLPs, 91 SSRs, and 3 SNPs) and spanned 2,104 cM with 20 linkage groups, and the male parent map, including 317 markers (205 AFLPs, 106 SSRs, 5 SNPs, and sex trait), spanned 2,453 cM with 23 main linkage groups. The sex, as morphological trait, was mapped on the linkage group XIX of the male parent map. The generated maps are among the most informative in SSRs when compared to the Populus maps published so far and allow a complete alignment with the 19 haploid chromosomes of Populus sequence genome. These genetic maps provide informative tools for a better understanding of P. nigra genome structure and genetic improvement of this ecologically and economically important European tree species. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Construction of a high-density,high-resolution genetic map and its integration with BAC-based physical map in channel catfish

Construction of genetic linkage map is essential for genetic and genomic studies. Recent advances in sequencing and genotyping technologies made it possible to generate high-density and high-resolution genetic linkage maps, especially for the organisms lacking extensive genomic resources. In the present work, we constructed a high-density and high-resolution genetic map for channel catfish with three large resource families genotyped using the catfish 250K single-nucleotide polymorphism (SNP) array. A total of 54,342 SNPs were placed on the linkage map, which to our knowledge had the highest marker density among aquaculture species. The estimated genetic size was 3,505.4 cM with a resolution of 0.22 cM for sex-averaged genetic map. The sex-specific linkage maps spanned a total of 4,495.1 cM in females and 2,593.7 cM in males, presenting a ratio of 1.7 : 1 between female and male in recombination fraction. After integration with the previously established physical map, over 87% of physical map contigs were anchored to the linkage groups that covered a physical length of 867 Mb, accounting for ∼90% of the catfish genome. The integrated map provides a valuable tool for validating and improving the catfish whole-genome assembly and facilitates fine-scale QTL mapping and positional cloning of genes responsible for economically important traits.     

A genetic linkage map and comparative genome analysis of common carp (Cyprinus carpio L.) using microsatellites and SNPs

A genetic linkage map is a powerful research tool for mapping traits of interest and is essential to understanding genome evolution. The aim of this study is to provide an expanded genetic linkage map of common carp to effectively carry out quantitative trait loci analysis and conduct comparative mapping analysis between lineages. Here, we constructed a genetic linkage map of common carp (Cyprinus carpio L.) using microsatellite and single-nucleotide polymorphism (SNP) markers in a 159 sibling family. A total of 246 microsatellites and 306 SNP polymorphic markers were genotyped in this family. Linkage analysis using JoinMap 4.0 organized 427 markers (186 microsatellites and 241 SNPs) to 50 linkage groups, ranging in size from 1.4 to 130.1 cM. Each group contained 2-30 markers. The linkage map covered a genetic distance of 2,039.2 cM and the average interval for markers within the linkage groups was approximately 6.4 cM. In addition, comparative genome analysis within five model teleost fish revealed a high percentage (74.7%) of conserved loci corresponding to zebrafish chromosomes. In most cases, each zebrafish chromosome comprised two common carp linkage groups. The comparative analysis also revealed independent chromosome rearrangements in common carp and zebrafish. The linkage map will be of great assistance in mapping genes of interest and serve as a reference to approach comparative mapping and enable further insights into the comprehensive investigations of genome evolution of common carp.     

A gene-based genetic linkage map of the collared flycatcher (Ficedula albicollis) reveals extensive synteny and gene-order conservation during 100 million years of avian evolution

By taking advantage of a recently developed reference marker set for avian genome analysis we have constructed a gene-based genetic map of the collared flycatcher, an important "ecological model" for studies of life-history evolution, sexual selection, speciation, and quantitative genetics. A pedigree of 322 birds from a natural population was genotyped for 384 single nucleotide polymorphisms (SNPs) from 170 protein-coding genes and 71 microsatellites. Altogether, 147 gene markers and 64 microsatellites form 33 linkage groups with a total genetic distance of 1787 cM. Male recombination rates are, on average, 22% higher than female rates (total distance 1982 vs. 1627 cM). The ability to anchor the collared flycatcher map with the chicken genome via the gene-based SNPs revealed an extraordinary degree of both synteny and gene-order conservation during avian evolution. The great majority of chicken chromosomes correspond to a single linkage group in collared flycatchers, with only a few cases of inter- and intrachromosomal rearrangements. The rate of chromosomal diversification, fissions/fusions, and inversions combined is thus considerably lower in birds (0.05/MY) than in mammals (0.6-2.0/MY). A dearth of repeat elements, known to promote chromosomal breakage, in avian genomes may contribute to their stability. The degree of genome stability is likely to have important consequences for general evolutionary patterns and may explain, for example, the comparatively slow rate by which genetic incompatibility among lineages of birds evolves.     

Identification of the sex‐determining region in flathead grey mullet (Mugil cephalus)

Elucidation of the sex‐determination mechanism in flathead grey mullet (Mugil cephalus) is required to exploit its economic potential by production of genetically determined monosex populations and application of hormonal treatment to parents rather than to the marketed progeny. Our objective was to construct a first‐generation linkage map of the M. cephalus in order to identify the sex‐determining region and sex‐determination system. Deep‐sequencing data of a single male was assembled and aligned to the genome of Nile tilapia (Oreochromis niloticus). A total 245 M. cephalus microsatellite markers were designed, spanning the syntenic tilapia genome assembly at intervals of 10 Mb. In the mapping family of full‐sib progeny, 156 segregating markers were used to construct a first‐generation linkage map of 24 linkage groups (LGs), corresponding to the number of chromosomes. The linkage map spanned approximately 1200 cM with an average inter‐marker distance of 10.6 cM. Markers segregating on LG9 in two independent mapping families showed nearly complete concordance with gender (R² = 0.95). The sex determining locus was fine mapped within an interval of 8.6 cM on LG9. The sex of offspring was determined only by the alleles transmitted from the father, thus indicating an XY sex‐determination system.     

Association mapping of morphological traits in wild and captive zebra finches: reliable within,but not between populations

Identifying causal genetic variants underlying heritable phenotypic variation is a long‐standing goal in evolutionary genetics. We previously identified several quantitative trait loci (QTL) for five morphological traits in a captive population of zebra finches (Taeniopygia guttata) by whole‐genome linkage mapping. We here follow up on these studies with the aim to narrow down on the quantitative trait variants (QTN) in one wild and three captive populations. First, we performed an association study using 672 single nucleotide polymorphisms (SNPs) within candidate genes located in the previously identified QTL regions in a sample of 939 wild‐caught zebra finches. Then, we validated the most promising SNP–phenotype associations (n = 25 SNPs) in 5228 birds from four populations. Genotype–phenotype associations were generally weak in the wild population, where linkage disequilibrium (LD) spans only short genomic distances. In contrast, in captive populations, where LD blocks are large, apparent SNP effects on morphological traits (i.e. associations) were highly repeatable with independent data from the same population. Most of those SNPs also showed significant associations with the same trait in other captive populations, but the direction and magnitude of these effects varied among populations. This suggests that the tested SNPs are not the causal QTN but rather physically linked to them, and that LD between SNPs and causal variants differs between populations due to founder effects. While the identification of QTN remains challenging in nonmodel organisms, we illustrate that it is indeed possible to confirm the location and magnitude of QTL in a population with stable linkage between markers and causal variants.     

The development and characterization of a 57K single nucleotide polymorphism array for rainbow trout

In this study, we describe the development and characterization of the first high‐density single nucleotide polymorphism (SNP) genotyping array for rainbow trout. The SNP array is publically available from a commercial vendor (Affymetrix). The SNP genotyping quality was high, and validation rate was close to 90%. This is comparable to other farm animals and is much higher than previous smaller scale SNP validation studies in rainbow trout. High quality and integrity of the genotypes are evident from sample reproducibility and from nearly 100% agreement in genotyping results from other methods. The array is very useful for rainbow trout aquaculture populations with more than 40 900 polymorphic markers per population. For wild populations that were confounded by a smaller sample size, the number of polymorphic markers was between 10 577 and 24 330. Comparison between genotypes from individual populations suggests good potential for identifying candidate markers for populations' traceability. Linkage analysis and mapping of the SNPs to the reference genome assembly provide strong evidence for a wide distribution throughout the genome with good representation in all 29 chromosomes. A total of 68% of the genome scaffolds and contigs were anchored through linkage analysis using the SNP array genotypes, including ~20% of the genome assembly that has not been previously anchored to chromosomes.     

Genome wide linkage disequilibrium in Chinese asparagus bean (Vigna. unguiculata ssp. sesquipedialis) germplasm: implications for domestication history and genome wide association studies

Association mapping of important traits of crop plants relies on first understanding the extent and patterns of linkage disequilibrium (LD) in the particular germplasm being investigated. We characterize here the genetic diversity, population structure and genome wide LD patterns in a set of asparagus bean (Vigna. unguiculata ssp. sesquipedialis) germplasm from China. A diverse collection of 99 asparagus bean and normal cowpea accessions were genotyped with 1127 expressed sequence tag-derived single nucleotide polymorphism markers (SNPs). The proportion of polymorphic SNPs across the collection was relatively low (39%), with an average number of SNPs per locus of 1.33. Bayesian population structure analysis indicated two subdivisions within the collection sampled that generally represented the 'standard vegetable' type (subgroup SV) and the 'non-standard vegetable' type (subgroup NSV), respectively. Level of LD (r(2)) was higher and extent of LD persisted longer in subgroup SV than in subgroup NSV, whereas LD decayed rapidly (0-2 cM) in both subgroups. LD decay distance varied among chromosomes, with the longest (≈ 5 cM) five times longer than the shortest (≈ 1 cM). Partitioning of LD variance into within- and between-subgroup components coupled with comparative LD decay analysis suggested that linkage group 5, 7 and 10 may have undergone the most intensive epistatic selection toward traits favorable for vegetable use. This work provides a first population genetic insight into domestication history of asparagus bean and demonstrates the feasibility of mapping complex traits by genome wide association study in asparagus bean using a currently available cowpea SNPs marker platform.     

SNP Discovery and Chromosome Anchoring Provide the First Physically-Anchored Hexaploid Oat Map and Reveal Synteny with Model Species

A physically anchored consensus map is foundational to modern genomics research; however, construction of such a map in oat (Avena sativa L., 2n = 6x = 42) has been hindered by the size and complexity of the genome, the scarcity of robust molecular markers, and the lack of aneuploid stocks. Resources developed in this study include a modified SNP discovery method for complex genomes, a diverse set of oat SNP markers, and a novel chromosome-deficient SNP anchoring strategy. These resources were applied to build the first complete, physically-anchored consensus map of hexaploid oat. Approximately 11,000 high-confidence in silico SNPs were discovered based on nine million inter-varietal sequence reads of genomic and cDNA origin. GoldenGate genotyping of 3,072 SNP assays yielded 1,311 robust markers, of which 985 were mapped in 390 recombinant-inbred lines from six bi-parental mapping populations ranging in size from 49 to 97 progeny. The consensus map included 985 SNPs and 68 previously-published markers, resolving 21 linkage groups with a total map distance of 1,838.8 cM. Consensus linkage groups were assigned to 21 chromosomes using SNP deletion analysis of chromosome-deficient monosomic hybrid stocks. Alignments with sequenced genomes of rice and Brachypodium provide evidence for extensive conservation of genomic regions, and renewed encouragement for orthology-based genomic discovery in this important hexaploid species. These results also provide a framework for high-resolution genetic analysis in oat, and a model for marker development and map construction in other species with complex genomes and limited resources.     

Characterizing linkage disequilibrium in pig populations

Knowledge of the extent and range of linkage disequilibrium (LD), defined as non-random association of alleles at two or more loci, in animal populations is extremely valuable in localizing genes affecting quantitative traits, identifying chromosomal regions under selection, studying population history, and characterizing/managing genetic resources and diversity. Two commonly used LD measures, r(2) and D', and their permutation based adjustments, were evaluated using genotypes of more than 6,000 pigs from six commercial lines (two terminal sire lines and four maternal lines) at ~4,500 autosomal SNPs (single nucleotide polymorphisms). The results indicated that permutation only partially removed the dependency of D' on allele frequency and that r(2) is a considerably more robust LD measure. The maximum r(2) was derived as a function of allele frequency. Using the same genotype dataset, the extent of LD in these pig populations was estimated for all possible syntenic SNP pairs using r(2) and the ratio of r(2) over its theoretical maximum. As expected, the extent of LD highest for SNP pairs was found in tightest linkage and decreased as their map distance increased. The level of LD found in these pig populations appears to be lower than previously implied in several other studies using microsatellite genotype data. For all pairs of SNPs approximately 3 centiMorgan (cM) apart, the average r(2) was equal to 0.1. Based on the average population-wise LD found in these six commercial pig lines, we recommend a spacing of 0.1 to 1 cM for a whole genome association study in pig populations.