- Research article
- Open Access
Extent of genome-wide linkage disequilibrium in Australian Holstein-Friesian cattle based on a high-density SNP panel
BMC Genomics volume 9, Article number: 187 (2008)
The extent of linkage disequilibrium (LD) within a population determines the number of markers that will be required for successful association mapping and marker-assisted selection. Most studies on LD in cattle reported to date are based on microsatellite markers or small numbers of single nucleotide polymorphisms (SNPs) covering one or only a few chromosomes. This is the first comprehensive study on the extent of LD in cattle by analyzing data on 1,546 Holstein-Friesian bulls genotyped for 15,036 SNP markers covering all regions of all autosomes. Furthermore, most studies in cattle have used relatively small sample sizes and, consequently, may have had biased estimates of measures commonly used to describe LD. We examine minimum sample sizes required to estimate LD without bias and loss in accuracy. Finally, relatively little information is available on comparative LD structures including other mammalian species such as human and mouse, and we compare LD structure in cattle with public-domain data from both human and mouse.
We computed three LD estimates, D', Dvol and r2, for 1,566,890 syntenic SNP pairs and a sample of 365,400 non-syntenic pairs. Mean D' is 0.189 among syntenic SNPs, and 0.105 among non-syntenic SNPs; mean r2 is 0.024 among syntenic SNPs and 0.0032 among non-syntenic SNPs. All three measures of LD for syntenic pairs decline with distance; the decline is much steeper for r2 than for D' and Dvol. The value of D' and Dvol are quite similar. Significant LD in cattle extends to 40 kb (when estimated as r2) and 8.2 Mb (when estimated as D'). The mean values for LD at large physical distances are close to those for non-syntenic SNPs. Minor allelic frequency threshold affects the distribution and extent of LD. For unbiased and accurate estimates of LD across marker intervals spanning < 1 kb to > 50 Mb, minimum sample sizes of 400 (for D') and 75 (for r2) are required. The bias due to small samples sizes increases with inter-marker interval. LD in cattle is much less extensive than in a mouse population created from crossing inbred lines, and more extensive than in humans.
For association mapping in Holstein-Friesian cattle, for a given design, at least one SNP is required for each 40 kb, giving a total requirement of at least 75,000 SNPs for a low power whole-genome scan (median r2 > 0.19) and up to 300,000 markers at 10 kb intervals for a high power genome scan (median r2 > 0.62). For estimation of LD by D' and Dvol with sufficient precision, a sample size of at least 400 is required, whereas for r2 a minimum sample of 75 is adequate.
Recent developments in high-throughput genotyping technology [1–3] and the discovery of large numbers of SNPs through sequencing of the cattle genome [4, 5] have generated enthusiasm and interest in genome-wide association mapping and marker-assisted selection in dairy cattle. Mapping by association requires a population-based sample rather than specific families. Association studies rely on the fact that alleles at loci that surround a quantitative trait nucleotide (QTN) tend to co-segregate. If the marker is sufficiently close to the QTN, the association remains intact in the majority of the individuals in the population, even after many generations. This non-random association between alleles at different loci is called linkage disequilibrium (LD). Such allelic associations are mostly due to physical proximity but are also affected by population history and evolutionary forces . The main limitation of association analyses is the requirement for sufficient markers to provide a high chance of identifying markers in LD with all important QTN. The marker density required for a successful association analysis, and for subsequent marker assisted selection (MAS), depends on the extent of LD across the genome.
Several measures of LD have been devised [7–12], and two measures, D' and r2, each with different statistical properties , are commonly used. Both range from 0 (no disequilibrium) to 1 (complete disequilibrium), but their interpretation is different. For biallelic markers, D' is equal to 1 if one or more of the four possible haplotypes is absent, and is < 1 if all four possible haplotypes are present. Most studies in livestock have reported the extent of LD based on D'. D' is most useful for representing historical recombination patterns, which are central to determining the extent and pattern of LD over a genome. D' is especially helpful in understanding long-range LD. Recently Chen et al.  suggested a volume measure of LD Dvol equivalent to D' and reported that this measure is more robust when estimated from small samples.
The measure r represents the statistical correlation between two sites, and takes the value of 1 (perfect LD) for a pair of biallelic markers if only two haplotypes are present within a population. Hence in order for r2 to be 1, allelic frequencies of the two SNPs in question need to be identical in addition to both being in LD. The power of association mapping is inversely proportional to r2, and to achieve the same power by typing a SNP in LD with a QTN versus typing the QTN directly, the sample size must be increased by a factor of 1/r2 . However, if more than one SNP is used to predict the effect of a QTN, then the association between the QTN and haplotypes is more informative.
The average extent of LD in the human genome has been extensively studied: it extends up to 50 kb but is highly variable, depending on the population and threshold used to measure LD [13, 6, 15, 16]. Studies on LD in cattle [17–22], sheep , pig [24, 25], horse , dog [27, 28] and chicken  have shown that LD in livestock populations is much more extensive than in humans. This is mainly due to smaller effective population size and, in some circumstances, stronger selection that has recently occurred in livestock populations .
The estimates of LD reported in cattle are mostly based on microsatellite markers. LD in humans estimated from SNP markers has been reported to be smaller than estimates from microsatellites . There is limited information available on the extent of LD between SNP markers in cattle and it is mainly limited to regions used in fine mapping [30, 31] or for a single chromosome [32, 22]. Furthermore, studies in cattle have usually involved relatively small sample sizes which are both subject to bias and loss of accuracy when estimating D' and r2; and such bias may vary with inter-marker distance. Finally, effects of minor allele frequency (MAF) of SNP loci on estimates of LD have not been reported in cattle.
In an ongoing gene discovery program in dairy cattle, a panel of 1546 Holstein-Friesian bulls was genotyped for 15,036 SNPs using a high-throughput genotyping service . This data set was recently described and used for the construction of a bovine HapMap . The present paper presents additional results of LD analyses of this genotype data. This study explores the effect of MAF and sample size on LD parameters and suggests the minimum sample size required for useful estimates. We also compare a recently suggested volume measure of LD  to the more commonly used statistics. Public-domain SNP data from mice and humans are also analysed to compare LD from different mammalian species using the same statistics and methods.
Of a total of 15,036 SNPs genotyped in this study, 9,195 survived filtering on MAF (> = 0.05) and HWE (P > 0.0001). The distribution of SNPs on different chromosomes is summarised in Table 1, ranging from 528 on BTA1 to 158 on BTA27. The overall mean MAF was 0.286 ± 0.0013. The distribution of MAF over the genome is approximately uniform (see Additional file 1, Figure S1), consistent with the ascertainment of SNP discovery identifying primarily the common SNPs within the bovine genome. The genomic coordinates of each SNP, together with its MAF, are available in Additional file 2.
Linkage disequilibrium among syntenic SNP pairs
Two of the most commonly used measures of LD, D' and r2, were estimated for each pair-wise combination of SNPs on each chromosome (syntenic SNPs): a total of 1,566,890 pairs were analyzed for all autosomes. The mean values of D' and r2 for individual autosomes are presented in Table 1. The mean D' and r2, pooled over autosomes (1–29) in different categories of map distances, are summarized in Table 2. Similar tables for individual chromosomes are provided in Additional file 3. The distribution of D' and r2 with respect to the physical distance separating loci is presented in Figures 1A and 1B, respectively. As expected, there is a gradual decline in D' with increasing physical distance between SNPs: for SNPs up to 1 kb apart, the mean D' is 0.99; for SNPs separated by 200–500 kb the mean D' is 0.46, and for SNPs separated by more than 50 Mb, the mean D' is 0.11. The distribution of expected D' obtained from fitting the Malécot model [34, 15] is also shown in Figure 1A. From this distribution, the estimated swept radius (the distance over which LD declines to ~37% of its initial value) is 8.2 Mb.
Compared to D', there is (as expected) a steeper decline in r2 with increasing distance between syntenic SNPs (Table 2, Figure 1B). The mean r2 across autosomes is 0.016 ± .00005. A mean r2 above 0.3 was observed only for SNPs less than 40 kb apart, and r2 declines rapidly with increasing distances. As can be seen in Table 2, starting from a maximum of 0.77 for SNPs less than 1 kb apart,r2 reduces to less than 10% of this value for SNPs greater than 500 kb apart. The value of r2 at a distance of more than 50 Mb (0.0034) is close to average r2 among the non-syntenic pairs (0.0032).
Similar trends for D' and r2 are evident for all the individual autosomes, although there is variation in the trend of decline (see Additional files 3, 4 and 5 for individual chromosomes). Heat maps for D' and r2 for individual chromosomes (see Additional files 6 and 7) clearly display variation in the pattern of LD across the genome. Figure 2A and Figure 2B compare the average LD in different distance bins of individual chromosomes with combined over genome. There is more variation in r2 across chromosomes at closer distance bins. The variation in LD across chromosomes suggests that inference of genome-wide LD based on single or few chromosomes can be biased. A different analysis of this variation in LD has recently been presented . The swept radius for D' of individual chromosomes shows that LD is higher for the longer chromosomes than the shorter chromosomes (ρ = 0.7). This can be explained by fewer recombinations per unit of physical distance in long chromosomes than in short chromosomes, as observed in the human genome.
A more robust measure of LD termed the volume measure of D' (Dvol) has recently been proposed . We computed Dvol for all the syntenic pairs. Since the value of D' and Dvol were quite similar (ρ = 0.99) in this data, the estimates of Dvol are not presented here, can be obtained from corresponding author.
LD among non-syntenic SNP pairs
The LD estimates for a sample of 365,400 non-syntenic pairs (two SNPs present on different chromosome) were computed. The mean D' for non-syntenic SNPs is 0.105 ± 0.0002 (Table 2). The distribution of D' for non-syntenic SNP pairs presented in Figure 3A indicates a larger proportion of non-syntenic pairs have small values of D'. The mean value r2 for non-syntenic SNPs is 0.0032 ± 0.000008. The distribution of r2 for non-syntenic pairs is shown in Figure 3B. The comparison of D' and Dvol estimated for non-syntenic SNP pairs shows that both estimates are quite similar (ρ = 0.99) for this sample (see Additional files 1, Figure S2).
The relationship between LD and average minor allelic frequency for non-syntenic SNPs
The estimate of D' seems to inflate with lower average minor allelic frequency (MAF) for non-syntenic SNP pairs (see Additional files 1, Figure S3). There are 1.27 percent of non-syntenic SNP pairs with D' values above 0.5. Most of these D' values were observed for SNP pairs with low average MAF. However, there is no relationship between MAF and LD expressed as r2 (see Additional files 1, Figure S4). It should be noted that all r2 values for non-syntenic pairs are very small. The maximum value for r2 observed for non-syntenic pairs is only 0.27, and only 9 pairs (0.002 %) have an r2 value of more than 0.1. Similarly, we also analysed chi-squared for non-syntenic SNP pairs and only small values of chi-squared were observed for the non-syntenic SNP pairs (data not shown). The proportion of non-syntenic pairs with a high D' value is small (1.27 %). However, of the pairs with a mean MAF less than 0.1, 18.9 percent have D' above 0.5 and there are no pairs with r2 > 0.1. These results suggest that D' is comparatively more sensitive to changes in MAF due to the increased likelihood of a missing haplotype when MAF is low. In contrast, r2 appears to be largely unaffected by low MAF. The relationship between Dvol and mean MAF shows a similar trend to that for D' (see Additional files 1, Figure S5), indicating that Dvol is also inflated for loci with rare alleles.
Effect of MAF on the extent of syntenic LD
The effect of MAF on the extent of LD amongst syntenic SNPs was studied using three different minimum MAF thresholds of 0.05, 0.1 and 0.2 (Figure 4A, B). Minimum MAF has a strong effect on mean r2 especially at short distances (up to 40 kb) (Figure 4B). For example, for SNPs within 1–10 kb, the mean r2 was 0.59 for MAF ≥ 0.05 but was as high as 0.70 for loci with MAF ≥ 0.2. A likely explanation is that as the MAF threshold increases, there is an increase in the number of SNP pairs with similar allele frequencies, and, consequently, an increase in r2. Compared to r2, there is an opposite effect of MAF threshold on D' (Figure 4A). There is almost no effect at close distances, especially up to 40 kb. However, at larger distances, higher MAF thresholds are associated with lower values of mean D', which is consistent with earlier studies .
Effect of sample size on accuracy and bias of LD estimators
Table 3, Figures 5A, B and 6 show the effect of sample size on accuracy and bias in estimation of the extent of LD (also see Additional files 1, Figure S6 and S7). In all cases of small samples sizes, estimates were biased upwards, and this trend was more noticeable for LD measured across marker intervals greater than 40 kb, and became exponentially worse for long-range LD estimated from inter-marker distances of > 10 Mb. Both estimates of LD based on the full sample of 1546 bulls were only slightly higher (%) as compared to estimates based on the reference sample of 1000 bulls. This may be because of higher average kinship between the 1546 bulls as compared with the subset of 1000 bulls. In Figures 5A and 5B, bias is expressed as the ratio of the LD estimate from the test sample of size n, to the LD estimate from the reference sample (n = 1000).
Figure 6 and Table 3 show the relative accuracy of LD estimates from different sample sizes, expressed as correlations between the LD parameters obtained from the test sample and the LD estimates obtained from the reference sample of 1000, for all chromosomes combined. For samples sizes of 75 or less, the accuracy of r2 appears to be significantly compromised. At samples sizes of 200, both the regression, correlation and residual variance of prediction shows goodness of fit consistent with estimates of the full sample of 1000. For D' a substantially larger sample size of 400 or more is required to obtain accurate estimates, compared with the estimates obtained for 1000 animals (Table 3). The accuracy of prediction, as indicated by correlations, is shown in Table 3. We compared the volume measure of D' (Dvol)  with D' estimated from different sample sizes for two chromosomes. As mentioned earlier, the estimates of D' and Dvol are quite similar especially for large samples. The correlation of Dvol estimated from different sample sizes was similar to that for D' (see Additional files 8, Table S2) indicating that accurate estimation of Dvol requires a sample size comparable to D'. Overall, we can conclude that a sample size of at least 75 is required for a reasonable estimate of r2, and a sample of at least 400 is required for D' and Dvol. In contrast to LD, MAF could be accurately estimated from samples sizes as small as 50 (Table 3). However, this comparison for MAF is based on the common SNPs with MAF more than 0.05 used in this study. Additional comparison of 1446 SNPs with less than 0.05 MAF (not included in the LD analysis) indicated that only large samples were able to detect the polymorphism in the majority of the SNPs (data not shown).
Comparison of LD in the cattle, mouse and human genomes
In order to compare the extent of LD in the bovine genome with mouse and human, we analysed some mouse and human SNP genotyping data available in the public domain. Although the extent of LD in cattle is much greater than in human (Figure 7A, B), the pattern of decay in LD with distance is similar. In contrast, LD in a mouse population derived from the intercrossing of 8 inbred lines, maintained for 50 generations by pseudorandom matting , shows extensive extended LD for both D' and r2. The extent of LD in this population extends up to 10 Mb for D', and for r2 useful LD of 0.3 extended more than 2 Mb. In the human data set, the distance to r2 = 0.3 extends only to 10 kb.
This is a comprehensive study of LD with a high density SNP panel in HF dairy cattle and currently reflects the best estimates of genome-wide LD in this breed based on number of animals screened and number of SNPs genotyped. The pairwise measures of LD decline over increasing distance between SNPs. The LD estimated by D' appears to be quite extensive and is much higher in cattle than in humans. This may be due to random drift caused by relatively small effective population in dairy cattle . Comparable estimates of extensive LD based on D' estimated from microsatellites have been reported in cattle [17–21], sheep , pig , horse , dog  and chicken . The LD between SNP markers reported in the present study is slightly smaller than that estimated between microsatellite markers in earlier studies [17, 21] and this may be explained in part by the differences in the mutation rate between these two types of markers and reflect the more recent origin of microsatellite polymorphisms. Secondly this may be due to the higher power for detecting LD when using markers with many alleles (e.g., microsatellites) as compared to biallelic SNP markers . The difference in LD detected using SNP and microsatellite loci was more pronounced in humans .
The D' metric has been suggested as a good measure to explain the extent of LD in population and variation of LD over the genome [22, 15]. However, individual values of D' are more influenced by variation in allele frequencies than for the r2 metric. This is reflected in the inflated D' values at low MAF . For non-syntenic SNP pairs we observed higher values of D' for the pairs with low mean MAF. It is known that SNPs of different MAFs have different LD properties . Higher D' between loci with rare alleles is mainly due to population genetic effects (rare alleles are, in general, younger than common alleles, and hence may still be in LD) and to effects of sampling. Smaller samples may fail to sample rare fourth gametes and, therefore, can inflate the D' estimate . On the other hand SNPs with rare alleles tend to have lower r2 values. The inflated D' estimates between non-syntenic SNPs with rare alleles are probably due to the effect of sampling caused by random drift or the loss of rare haplotypes in sampling in the present study.
The decline in LD is much steeper for r2 than for D'. Such differences in the pattern of D' and r2 over physical distances have previously been observed in humans [38, 39]. These two measures, which have been widely used in practice, have different statistical properties : D' focuses on historical recombination, which is central to defining the extent and pattern of LD over a genome. However, r2 is more useful for predicting the power of association mapping. To obtain the same power as obtained when testing the causative quantitative trait nucleotide (QTN)/mutation, the sample size required for association mapping is inversely proportional to r2 [13, 40, 41]. From the pattern of decline of r2 the average useful LD for single-point association mapping in this population extends only up to 40 kb which suggests that at least 75,000 SNPs are required for a whole genome association scan. At this distance r2 values between adjacent SNPs decrease to an average of ~0.3 (median r2 = 0.19). At this spacing a QTN would be at a maximum of 20 kb (located at the mid point of the interval) from an adjacent SNP which would give an average r2 of 0.4 (median r2 = 0.28) between SNP and QTN. However, if more stringent criteria are considered for higher power genome scans, then the number of SNPs required would be 150,158 (one SNP every 20 kb) and 300,631 (one SNP every 10 kb) to obtain average r2 values between adjacent SNPs of 0.4 (median r2 = 0.28), and 0.6 (median r2 = 0.62), respectively. In addition there is a lot of variation in r2, as indicated by the large standard deviation of r2 (Table 2), within an interval of extent of LD considered. The 25th percentile of r2 (Table 2) indicates that only 75 % of the pairs of SNPs in the 1–10 kb distance bin have r2 of more than 0.2. Similar low values of the 25th percentile are noted for other distance bins in Table 2. The variation in r2 in a distance bin is partly because of the variation in LD in different genomic regions. In addition r2 is dependent upon the matching allelic frequencies and are known to have very low values between markers even at very short distances . This variation has been ignored in most studies while suggesting the number of SNPs required for genome scans based on average r2 . To accommodate this variation, more SNPs will need to be genotyped in each interval which will increase the estimate of the total number of SNPs required for association mapping. However, these partial correlations can be exploited using multi-markers haplotype analysis which provides more discriminatory power as compared to individual SNPs to detect putative causal mutation. Recently McKay et al.  suggested one SNP every 100 kb to obtain an average r2 of 0.15-0.2 between adjacent SNPs. Gautier et al.  reported the LD analysis of 526 SNPs mostly located on one chromosome in 14 breeds and suggested a common panel of 300,000 SNPs (one SNP every 10 kb) for association mapping in different breeds similar to the requirement for a high power panel within-breed study as shown here.
The number of SNPs required for association mapping can be reduced by excluding some of the redundant SNPs, by optimally using the LD information present in the population. However, there are differences in the pattern of LD across the genome. This can be addressed more precisely by identifying tag SNPs based on haplotype block structure information, as was done in the human HapMap project . An attempt to construct the haplotype block map of the bovine genome and the concomitant identification of tag SNPs from this dataset was presented by Khatkar et al.  but at present only covers a relatively small proportion of the bovine genome based on 15,036 markers. Genome-wide identification of tags SNPs for the whole genome would be possible from a saturated haplotype block map of the bovine genome. Until such maps become available, the extent of LD as expressed by r2 can provide an interim guide for number and spacing of the SNPs over the genome.
The extent of LD within a genome can be affected by a number of factors. Our results confirm that MAF has direct effects on the estimation of extent of LD. The r2 between common SNPs is higher especially for SNPs at close physical distances. On the other hand D' between SNPs with low MAF is higher at longer distances. Similarly, sample size also affects the estimation of LD. The results in this study clearly indicate that the estimate of D' is affected most by sample size. It seems that for reliable estimates of D' a sample of 400 or more is required. Similar observations were also made for Dvol. The requirement of sample size would be even higher in human because of the larger effective population size. Hence, it may be suggested that analyses utilizing D' matrix (like construction of LD maps and HapMap) should be based on a large sample size and preferably from within-breed group samples.
However, the identification of tag SNPs, which is generally based on r2 values, can be accomplished using smaller samples. Similar estimates of correlations between the estimates of r2 were obtained from different samples in a simulation study by Visscher . Similarly, a small sample size of 50 and above can provide precise estimate of MAF for common SNPs in the population (Table 3).
We have conducted an extensive LD analysis using 9,195 SNPs genotyped on 1546 Holstein Frisian bulls. Overall the extent of LD over different chromosomes was similar but varied in cattle. We compared three different measures of LD. All are affected by sample size although to different extents. We suggest that the sample size required to compute reliable estimates of D' or for any analysis (HapMap and LD maps) based on the D' matrix should be at least 400 samples. A similar sample size is required for Dvol. However, it seems that r2 can be reliably estimated with smaller sample sizes of 75 individuals. Based on the extent of decline of r2 we suggest more than 75,000 SNPs would be required for low power association mapping in the Holstein Friesian population and up to 300,000 SNPs for a high power study. The extent of LD in cattle is higher than in human, but much less than in a mouse population.
A total of 1,546 Holstein Friesian bulls were genotyped for 15,036 SNPs. A subset of 1,000 bulls (called the reference sample) was selected for LD analysis based on minimizing their pedigree relationship with each other. However, when studying the effect of sample size on LD estimators, estimates of LD from the entire sample of 1,546 bulls were also used. Of the total SNPs, 10,000 (MegAllele Genotyping Bovine 10,000 SNP Panel, ParAllele) were generated as part of the community project of the International Bovine Genome Sequencing Consortium (IBGSC) . The remaining 5,036 custom SNPs were selected from the Interactive Bovine In Silico SNP (IBISS) database [47, 48], from in-house sequencing, and from publications. A high-throughput SNP assay service provided by Affymetrix, Inc. was used for the genotyping. The locations of the SNPs were determined on the bovine sequence assembly Btau 3.1  The details of the genotyping and mapping of SNPs were given in Khatkar et al. .
Data on 2202 mice genotyped for 13,459 SNPs were downloaded from the Wellcome Trust Centre for Human Genetics's web page [50, 35] and analysed for D' and r2 with Haploview  for all 19 autosomes. Pedigree structure was ignored for these computations and the results were pooled over all autosomes. This data set is not representative of natural populations: it resulted from the intercrossing of 8 inbred lines to create a population maintained for 50 generations of pseudo-random mating and referred to as heterogeneous stock mice. It is of interest, however, because of its central use in QTL mapping .
Data on 537 participants genotyped for 408,000 SNPs described by Fung et al.  were downloaded from NINDS Human Genetics DNA and Cell Line Repository . A sample of 1500 SNPs on one chromosome (HSA15) was selected for estimation of LD parameters and comparison to mouse and cattle.
Estimation of linkage disequilibrium
Estimates of the two standard descriptive linkage disequilibrium (LD) parameters, D'  and r2 , were obtained via Haploview software  for all pair-wise combinations of SNPs on each chromosome. SNPs showing significant deviations from HWE (P < 0.0001) were excluded from analysis, as were SNPs with minor allelic frequency (MAF) < 0.05, and those genotyped in less than 50% of the bulls. The Malecot model for the relationship between LD and distance between loci ;  was fitted to the D' data on individual autosomes and combined over all autosomes. As shown in the results section, D' above 50 Mb were close to non-syntenic (background) LD, hence SNPs pairs with up to 50 Mb distance between them were included for fitting Malecot model. This also allowed similar window size for comparison of LD over the autosomes of different length. Fitting this model enabled estimation of the swept radius, the extent of LD that is useful for mapping  defined as the distance in kb over which LD declines to approximately 37% of its initial value. Estimates of LD statistics among a sample of non-syntenic SNP pairs were computed using the genetics package  in the R statistical software suite . A random sample of 30 SNPs was chosen from each chromosome, and all non-syntenic pairs were used to assess the extent of non-syntenic LD.
We also computed a volume measure of D' (Dvol) as per the formulae given and implemented in the C codes of Chen et al. . The program computes Dvol from the four haplotype counts. Four haplotype counts were computed with GOLD .
Estimation of samples size on accuracy and bias of LD estimates
To examine the effect of sample size on the accuracy of estimation of LD parameters, samples of 25, 50, 75, 100, 150, 200, 300, 400, and 500 were randomly drawn from the total sample of 1,000 analysed in this study (Table 3, Figure 6). The effect of sample size on the accuracy and bias in estimation of the extent of LD obtained from the 1000 animal reference sample. The bias was expressed as a ratio of statistics from test sample size to the reference sample. The data consisting of all 1546 animals (which included some closely related animals) were also included for comparison here. The SNPs in each sample were evaluated for MAF and HWE. The LD between SNPs pairs for all autosomes was estimated for each sample, after excluding SNPs with MAF < 0.05 and HWE (P < 0.0001). Pair-wise estimates of LD and MAF were compared between different samples.
Gunderson KL, Steemers FJ, Lee G, Mendoza LG, Chee MS: A genome-wide scalable SNP genotyping assay using microarray technology. Nat Genet. 2005, 37 (5): 549-554. 10.1038/ng1547.
Hirschhorn JN, Daly MJ: Genome-wide association studies for common diseases and complex traits. Nat Rev Genet. 2005, 6 (2): 95-108. 10.1038/nrg1521.
Craig DW, Stephan DA: Applications of whole-genome high-density SNP genotyping. Expert Rev Mol Diagn. 2005, 5 (2): 159-170. 10.1586/14737188.8.131.52.
Snelling WM, Casas E, Stone RT, Keele JW, Harhay GP, Bennett GL, Smith TP: Linkage mapping bovine EST-based SNP. BMC Genomics. 2005, 6 (1): 74-10.1186/1471-2164-6-74.
Lee SH, Park EW, Cho YM, Lee JW, Kim HY, Lee JH, Oh SJ, Cheong IC, Yoon DH: Confirming single nucleotide polymorphisms from expressed sequence tag datasets derived from three cattle cDNA libraries. J Biochem Mol Biol. 2006, 39 (2): 183-188.
Ardlie KG, Kruglyak L, Seielstad M: Patterns of linkage disequilibrium in the human genome. Nat Rev Genet. 2002, 3 (4): 299-309. 10.1038/nrg777.
Lewontin RC: The interaction of selection and linkage. I. General considerations: heterotic models. Genetics. 1964, 49: 49-67.
Hedrick PW: Gametic disequilibrium measures: proceed with caution. Genetics. 1987, 117 (2): 331-341.
Devlin B, Risch N: A comparison of linkage disequilibrium measures for fine-scale mapping. Genomics. 1995, 29 (2): 311-322. 10.1006/geno.1995.9003.
Morton NE, Zhang W, Taillon-Miller P, Ennis S, Kwok PY, Collins A: The optimal measure of allelic association. Proc Natl Acad Sci USA. 2001, 98 (9): 5217-5221. 10.1073/pnas.091062198.
Jorde LB: Linkage disequilibrium and the search for complex disease genes. Genome Res. 2000, 10 (10): 1435-1444. 10.1101/gr.144500.
Zhao H, Nettleton D, Soller M, Dekkers JC: Evaluation of linkage disequilibrium measures between multi-allelic markers as predictors of linkage disequilibrium between markers and QTL. Genet Res. 2005, 86 (1): 77-87. 10.1017/S001667230500769X.
Pritchard JK, Przeworski M: Linkage disequilibrium in humans: models and data. Am J Hum Genet. 2001, 69 (1): 1-14. 10.1086/321275.
Chen Y, Lin CH, Sabatti C: Volume measures for linkage disequilibrium. BMC Genet. 2006, 7: 54-10.1186/1471-2156-7-54.
Maniatis N, Collins A, Xu CF, McCarthy LC, Hewett DR, Tapper W, Ennis S, Ke X, Morton NE: The first linkage disequilibrium (LD) maps: delineation of hot and cold blocks by diplotype analysis. Proc Natl Acad Sci USA. 2002, 99 (4): 2228-2233. 10.1073/pnas.042680999.
Weiss KM, Clark AG: Linkage disequilibrium and the mapping of complex human traits. Trends Genet. 2002, 18 (1): 19-24. 10.1016/S0168-9525(01)02550-1.
Farnir F, Coppieters W, Arranz JJ, Berzi P, Cambisano N, Grisart B, Karim L, Marcq F, Moreau L, Mni M: Extensive genome-wide linkage disequilibrium in cattle. Genome Res. 2000, 10 (2): 220-227. 10.1101/gr.10.2.220.
Vallejo RL, Li YL, Rogers GW, Ashwell MS: Genetic diversity and background linkage disequilibrium in the North American Holstein cattle population. J Dairy Sci. 2003, 86 (12): 4137-4147.
Tenesa A, Knott SA, Ward D, Smith D, Williams JL, Visscher PM: Estimation of linkage disequilibrium in a sample of the United Kingdom dairy cattle population using unphased genotypes. J Anim Sci. 2003, 81 (3): 617-623.
Odani M, Narita A, Watanabe T, Yokouchi K, Sugimoto Y, Fujita T, Oguni T, Matsumoto M, Sasaki Y: Genome-wide linkage disequilibrium in two Japanese beef cattle breeds. Anim Genet. 2006, 37 (2): 139-144. 10.1111/j.1365-2052.2005.01400.x.
Khatkar MS, Thomson PC, Tammen I, Cavanagh JA, Nicholas FW, Raadsma HW: Linkage disequilibrium on chromosome 6 in Australian Holstein-Friesian cattle. Genet Sel Evol. 2006, 38 (5): 463-477. 10.1051/gse:2006015.
Khatkar MS, Collins A, Cavanagh JA, Hawken RJ, Hobbs M, Zenger KR, Barris W, McClintock AE, Thomson PC, Nicholas FW: A first-generation metric linkage disequilibrium map of bovine chromosome 6. Genetics. 2006, 174 (1): 79-85. 10.1534/genetics.106.060418.
McRae AF, McEwan JC, Dodds KG, Wilson T, Crawford AM, Slate J: Linkage disequilibrium in domestic sheep. Genetics. 2002, 160 (3): 1113-1122.
Nsengimana J, Baret P, Haley CS, Visscher PM: Linkage disequilibrium in the domesticated pig. Genetics. 2004, 166 (3): 1395-1404. 10.1534/genetics.166.3.1395.
Jungerius BJ, Gu J, Crooijmans RP, Poel van der JJ, Groenen MA, van Oost BA, te Pas MF: Estimation of the extent of linkage disequilibrium in seven regions of the porcine genome. Anim Biotechnol. 2005, 16 (1): 41-54. 10.1081/ABIO-200053402.
Tozaki T, Hirota KI, Hasegawa T, Tomita M, Kurosawa M: Prospects for whole genome linkage disequilibrium mapping in thoroughbreds. Gene. 2005, 346: 127-132. 10.1016/j.gene.2004.10.011.
Lou XY, Todhunter RJ, Lin M, Lu Q, Liu T, Wang Z, Bliss SP, Casella G, Acland GM, Lust G: The extent and distribution of linkage disequilibrium in a multi-hierarchic outbred canine pedigree. Mamm Genome. 2003, 14 (8): 555-564. 10.1007/s00335-003-2272-1.
Sutter NB, Eberle MA, Parker HG, Pullar BJ, Kirkness EF, Kruglyak L, Ostrander EA: Extensive and breed-specific linkage disequilibrium in Canis familiaris. Genome Res. 2004, 14 (12): 2388-2396. 10.1101/gr.3147604.
Heifetz EM, Fulton JE, O'Sullivan N, Zhao H, Dekkers JC, Soller M: Extent and consistency across generations of linkage disequilibrium in commercial layer chicken breeding populations. Genetics. 2005, 171 (3): 1173-1181. 10.1534/genetics.105.040782. Epub 2005 Aug 1122.
Grisart B, Coppieters W, Farnir F, Karim L, Ford C, Berzi P, Cambisano N, Mni M, Reid S, Simon P: Positional candidate cloning of a QTL in dairy cattle: Identification of a missense mutation in the bovine DGAT1 gene with major effect on milk yield and composition. Genome Res. 2002, 12 (2): 222-231. 10.1101/gr.224202.
de Koning DJ: Conflicting candidates for cattle QTLs. Trends Genet. 2006, 22: 301-305. 10.1016/j.tig.2006.04.006.
Spelman RJ, Coppieters W: Linkage disequilibrium in the New Zealand Jersey population. Presented in 8th World Congress on Genetics Applied to Livestock Production, August 13–18, Belo Horizonte, MG, Brasil. 2006
Khatkar MS, Zenger KR, Hobbs M, Hawken RJ, Cavanagh JA, Barris W, McClintock AE, McClintock S, Thomson PC, Tier B: A primary assembly of a bovine haplotype block map based on a 15,036-single-nucleotide polymorphism panel genotyped in Holstein-Friesian cattle. Genetics. 2007, 176 (2): 763-772. 10.1534/genetics.106.069369. Epub 2007 Apr 2015.
Malecot G: Les Mathématiques de l'Hérédité. Maison et Cie, Paris. 1948
Valdar W, Solberg LC, Gauguier D, Burnett S, Klenerman P, Cookson WO, Taylor MS, Rawlins JN, Mott R, Flint J: Genome-wide genetic association of complex traits in heterogeneous stock mice. Nat Genet. 2006, 38 (8): 879-887. 10.1038/ng1840. Epub 2006 Jul 2009.
Zenger KR, Khatkar MS, Cavanagh JA, Hawken RJ, Raadsma HW: Genome-wide genetic diversity of Holstein Friesian cattle reveals new insights into Australian and global population variability, including impact of selection. Anim Genet. 2007, 38 (1): 7-14. 10.1111/j.1365-2052.2006.01543.x.
Slatkin M: Linkage disequilibrium in growing and stable populations. Genetics. 1994, 137 (1): 331-336.
Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, Higgins J, DeFelice M, Lochner A, Faggart M: The structure of haplotype blocks in the human genome. Science. 2002, 296 (5576): 2225-2229. 10.1126/science.1069424.
International HapMap Consortium: A haplotype map of the human genome. Nature. 2005, 437 (7063): 1299-1320. 10.1038/nature04226.
Wall JD, Pritchard JK: Haplotype blocks and linkage disequilibrium in the human genome. Nat Rev Genet. 2003, 4 (8): 587-597. 10.1038/nrg1123.
Zondervan KT, Cardon LR: The complex interplay among factors that influence allelic association. Nat Rev Genet. 2004, 5 (2): 89-100. 10.1038/nrg1270.
McKay SD, Schnabel RD, Murdoch BM, Matukumalli LK, Aerts J, Coppieters W, Crews D, Dias Neto E, Gill CA, Gao C: Whole genome linkage disequilibrium maps in cattle. BMC Genet. 2007, 8 (1): 74-10.1186/1471-2156-8-74.
Gautier M, Faraut T, Moazami-Goudarzi K, Navratil V, Foglio M, Grohs C, Boland A, Garnier JG, Boichard D, Lathrop GM: Genetic and haplotypic structure in 14 European and african cattle breeds. Genetics. 2007, 177 (2): 1059-1070. 10.1534/genetics.107.075804. Epub 2007 Aug 1024.
International HapMap Project. [http://www.hapmap.org/]
Visscher PM: Variation of estimates of SNP and haplotype diversity and linkage disequilibrium in samples from the same population due to experimental and evolutionary sample size. Ann Hum Genet. 2007, 71 (Pt 1): 119-126. 10.1111/j.1469-1809.2006.00305.x.
International Bovine Genome Sequencing Consortium (IBGSC). [http://www.hgsc.bcm.tmc.edu/projects/bovine/]
The Interactive Bovine In Silico SNP (IBISS) database. [http://www.livestockgenomics.csiro.au/ibiss/]
Hawken RJ, Barris WC, McWilliam SM, Dalrymple BP: An interactive bovine in silico SNP database (IBISS). Mamm Genome. 2004, 15 (10): 819-827. 10.1007/s00335-004-2382-4.
Bovine Genome Assembly. [ftp://ftp.hgsc.bcm.tmc.edu/pub/data/Btaurus/fasta/Btau20060815-freeze/]
The Wellcome Trust Centre for Human Genetics. [http://gscan.well.ox.ac.uk/]
Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21 (2): 263-265. 10.1093/bioinformatics/bth457.
Fung HC, Scholz S, Matarin M, Simon-Sanchez J, Hernandez D, Britton A, Gibbs JR, Langefeld C, Stiegert ML, Schymick J: Genome-wide genotyping in Parkinson's disease and neurologically normal controls: first stage analysis and public release of data. Lancet Neurol. 2006, 5 (11): 911-916. 10.1016/S1474-4422(06)70578-6.
NINDS Human Genetics DNA and Cell Line Repository. [http://ccr.coriell.org/ninds/index.html]
Hill WG, Robertson A: Linkage disequilibrium in finite populations. Theor Appl Genet. 1968, 38: 226-231. 10.1007/BF01245622.
Warnes G, Leisch F: genetics: Population Genetics. R package version 1.2.0. 2005, [http://cran.r-project.org/src/contrib/Descriptions/genetics.html]
R Development Core Team: R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0. 2005, [http://www.r-project.org]
Abecasis GR, Cookson WO: GOLD – graphical overview of linkage disequilibrium. Bioinformatics. 2000, 16 (2): 182-183. 10.1093/bioinformatics/16.2.182.
We thank Genetics Australia for semen samples and the Australian Dairy Herd Improvement Scheme for pedigree data. We thank Prof. Johann Sölkner for his editorial comments on an early version of the manuscript, and Mr Steve Christov and Prof. Johann Solkner for pointing out the web source for the human SNP data. This research is supported by the Cooperative Research Centre for Innovative Dairy Products, Victoria, Australia. JFT and RDS were supported by National Research Initiative grants 2005-35205-15448, 2005-35604-15615, 2006-35205-16701 and 2006-35616-16697 from the USDA Cooperative State Research, Education and Extension Service. We gratefully acknowledge that the genotype data on which this paper is based was able to be generated as a result of the considerable efforts of the Baylor College of Medicine Human Genome Sequencing Center and the Bovine Genome Sequencing Project Consortium.
MSK conceived the study, contributed in its design, data collection, analysis and was the primary author for assembling the manuscript. ARC participated in interpretation and manuscript preparation. FWN participated in data interpretation and manuscript preparation. KRZ contributed in generating SNP data and preliminary data analyses. JALC organized DNA sample preparation and participated in SNP discovery and SNP genotyping. WB contributed in SNP discovery and provided bioinformatic support in SNP positioning. RDS provided samples and edited the manuscript. JFT provided samples and edited the manuscript. HWR is project leader and contributed in project concept, design, interpretation and manuscript preparation.
Electronic supplementary material
Additional file 1: Figure S1. Frequency distribution of MAF for 9,195 SNPs used in the LD analysis. Figure S2. Comparison of Dvol and D' computed for non-syntenic SNP pairs. Figure S3. Relationship between mean minor allelic frequency (MAF) and D' for non-syntenic SNP pairs. Figure S4. Relationship between mean minor allelic frequency and r2 for non-syntenic SNP pairs. Figure S5. Relationship between mean minor allelic frequency and Dvol for non-syntenic SNP pairs. Figure S6. Mean D' at different physical distances pooled over all the autosomes and estimated using different sample sizes. The numbers with n_ in legend are the sample sizes randomly drawn from a total of 1000. The n_1546 are all the bulls genotyped in the present study. Figure S7. Mean r2 at different physical distances pooled over all autosomes and estimated using different sample sizes. The numbers with n_ in legend are the sample sizes randomly drawn from a total of 1000. The n_1546 are all the bulls genotyped in the present study. (DOC 118 KB)
Additional file 4: Figure S9. Distribution of D' between SNP pairs in relation to the physical distance between loci (Mb), for individual autosomes (1–29). The red line shows average D' in each 500 kb sliding window. The blue line shows the theoretical distribution from the Malécot model. (DOC 328 KB)
Additional file 6: Figure S11. Linkage disequilibrium (D') between SNPs on all the bovine autosomal chromosomes (BTA1-BTA29) presented in the form of heatmap of D' . All SNP with the MAF of less than 0.05 and showing deviation from HWE were excluded in the LD measurement. These figures were prepared by Haploview software. Here bright red color indicates D' = 1 with LOD ≥ 2, blue D' = 1 with LOD < 2, shades of pink/red D' < 1 with LOD ≥ 2 and white D' < 1 with LOD < 2. (DOC 2 MB)
Additional file 7: Figure S12. Linkage disequilibrium (r2) between SNPs on all the bovine autosomal chromosomes (BTA1-BTA29) presented in the form of heatmap of r2 . All SNP with the MAF of less than 0.05 and showing deviation from HWE were excluded in the LD measurement. These figures were prepared by Haploview software. Here black color indicates r2 = 1, shades of grey 0 <r2 < 1 and white r2 = 0. (DOC 3 MB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Khatkar, M.S., Nicholas, F.W., Collins, A.R. et al. Extent of genome-wide linkage disequilibrium in Australian Holstein-Friesian cattle based on a high-density SNP panel. BMC Genomics 9, 187 (2008) doi:10.1186/1471-2164-9-187
- Linkage Disequilibrium
- Minor Allelic Frequency
- Association Mapping
- Linkage Disequilibrium Analysis
- Quantitative Trait Nucleotide