Go to Top Go to Bottom
Anim Biosci > Volume 31(12); 2018 > Article
Shin, Won, Kim, and Kim: Extent of linkage disequilibrium and effective population size of Korean Yorkshire swine

Abstract

Objective

We aimed to characterize linkage disequilibrium (LD) and effective population size (Ne) in a Korean Yorkshire population using genomic data from thousands of individuals.

Methods

We genotyped 2,470 Yorkshire individuals from four major Grand-Grand-Parent farms in Korea using the Illumina PorcineSNP60 version2 BeadChip, which covers >61,565 single nucleotide polymorphisms (SNPs) located across all chromosomes and mitochondria. We estimated the expected LD and inferred current Ne as well as ancestral Ne.

Results

We identified 61,565 SNP from autosomes, mitochondria, and sex chromosomes and characterized the LD of the Yorkshire population, which was relatively high between closely linked markers (>0.55 at 50 kb) and declined with increasing genetic distance. The current Ne of this Korean Yorkshire population was 122.87 (106.90; 138.84), while the historical Ne of Yorkshire pigs suggests that the ancestor Ne has decreased by 99.6% over the last 10,000 generations.

Conclusion

To maintain genetic diversity of a domesticated animal population, we must carefully consider appropriate breed management methods to avoid inbreeding. Although attenuated selection can affect short-term genetic gain, it is essential for maintaining the long-term genetic variability of the Korean Yorkshire population. Continuous and long-term monitoring would also be needed to maintain the pig population to avoid an unintended reduction of Ne. The best way to preserve a sustainable population is to maintain a sufficient Ne.

INTRODUCTION

Several important evolutionary processes in finite populations, including migration, artificial or natural selection, and genetic drift lead to non-random association of alleles between two different loci, or linkage disequilibrium (LD) [1]. Recent genomic methods in animal breeding, such as genome-wide association studies (GWAS) and genomic selection using single nucleotide polymorphism (SNP) data, depend on the extent of LD and its rate of decline with distance between loci of the population. Researchers have applied SNP chips for GWAS [24] and genomic selection studies [57] in pigs but have found that most traits of interest are complex and more suitable for genomic selection techniques than GWAS for studying significant associations in some genomic regions [8]. Methods used to study animal breeding genetics rely strongly on the quality of LD and sample size. Therefore, characterization of LD is very important for planning future genomic studies of complex traits relevant to animal breeding.
The LD among loci can provide insights into the evolutionary history of each population by using effective population size (Ne), where Ne is the number of individuals in an idealized population that would yield the same inbreeding degree as the real population [9]. Thus, we can monitor genetic diversity in each domesticated animal population based on Ne and explain the observed extent and pattern of genetic variation in population genetic terms. Using Ne, we can also predict loss of genetic variation from a prospective point of view and the accuracy of natural genomic selection before the emergence of artificial genomic selection in domesticated animal breeding. Additionally, we can infer ancestral Ne using the strength of LD at different genetic distances between markers. The pattern of historical Ne in each animal population can increase our understanding of the impact of recent strong artificial selection breeding methods on population-level genetic variation. If pedigrees are incomplete or unavailable, we can use inbreeding information of populations of interest with respect to Ne.
In Korea, several pig breeds are economically significant, including Yorkshire, Landrace, Duroc, and Berkshire, though Yorkshire is one of the most important breeds because has an excellent maternal line index and outnumbers all other maternal line pig breeds in Korean Grand-Grand-Parent (GGP) farms. Patterns of LD in Yorkshire populations in other countries have been characterized previously, and their Ne were predicted using SNP chip data. For example, the extent of LD of multiple pig breeds, including Yorkshire, from the United States, Denmark, and The Netherlands were investigated [1012]. In addition, Uimari determined the extent of LD and estimated LD-based actual Ne and ancestral Ne of 32 Finnish Yorkshire boars [13]. Although the Korean pig industry is active, Korea has no original Yorkshire breed and relies on significant pig imports. As a result, the Korean Yorkshire population has diverse genetic sources from several countries, so Ne could be an important measurement of both the Korean and global pig industries.
In our current study, LD and Ne in a Korean Yorkshire population were estimated using data generated with the Illumina PorcineSNP60 v2 BeadChip (Illumina Inc., San Diego, CA, USA). We also investigated the ancestral Ne of the population. Together with findings from other studies, our results can inform the establishment and implementation of the most effective animal breeding genomic methods for Korean Yorkshire swine.

MATERIALS AND METHODS

Study sampling and selection of genotypic data

This study included genotypic data from 2,470 Yorkshire pigs from four representative major GGP farms (1,473 from Farm 1, 542 from Farm 2, 312 from Farm 3, and 143 from Farm 4) in Korean, constituting a representative sample of the Korean Yorkshire population. All pigs were pigs from each major GGP farms in 2015 through 2016 (Supplementary Figure S1). The Illumina PorcineSNP60 v2 BeadChip (Illumina Inc., USA), which contains 61,565 SNPs located across all autosomes, mitochondria, and sex chromosomes, was used in this study. These informative SNPs were selected from the recent porcine reference genome (Build 9) from the Ensembl database. The generated genotypic data were analyzed using PLINK (version 1.90) for quality control. First, we extracted 52,257 autosomal SNPs from the total 61,565 SNPs on all autosomes and then removed certain SNPs from analysis due to poor genotyping qualify; 9,303 SNPs were excluded based on Hardy-Weinberg equilibrium testing (p≤0.000001). The total genotyping rate of remaining sequences was 0.996, but 480 SNPs failed the missingness test (missing genotype>0.1), and 10,160 SNPs with a minor allele frequency (MAF) of <0.05 were excluded. The criteria for genotyping quality used were similar to those described previously [13] and yielded 33,418 SNPs for downstream analysis. We performed imputation to replace missing data with substituted genotypes [14] before LD and Ne were characterized.

Characterization of linkage disequilibrium in Korean Yorkshire pigs

We estimated LD using the LDcorSV R package (version 1.3.1). Pairwise haplotype frequencies were estimated from the genotype data and used to calculate the squared correlation coefficient between two loci (r2). The r2 are respectively equivalent to the covariance and correlation between alleles at two different loci computed as:
(1)
r2=D2PAPaPBPb
Where PA, Pa, PB, and Pb are the respective frequencies of alleles A, a, B, and b, and D is PAB - PAPB.
For each SNP, pairwise LD was calculated for its adjacent 100 SNPs. The SNP quality and distance requirements resulted in approximately 3.2 million SNP pairs distributed over the total genome. This r2 calculation was performed on a chromosome-by-chromosome basis and illustrated a relationship of physical distance and r2 between the two target loci per chromosome (Supplementary Figure S2). Details of the SNPs’ physical position can be found in the Illumina product literature. To determine LD in relation to the physical distance between SNPs, SNP pairs were divided into distance bins. We established two classes (0 to 0.5 Mb and 0 to 5 Mb), and applicable SNP pairs in each class were placed into 1 of 50 distance bins with class-dependent bin ranges (Supplementary Table S1). The two types of mean r2 for each distance bin were plotted against the median of the distance bin range, which are presented in Figure 1.

Construction model of linkage disequilibrium with distance

Under the assumption of an isolated population with random mating, Sved derived an approximate expression for the expectation of r2 [15]:
(2)
E(r2)=11+4Nc
Where N is Ne, and c is the recombination frequency. We replaced c with linkage distance in morgans [1618], and our calculations were further supported by approximating the more precise equation for E(r2) given by Sved [15]. Based on this formula, a non-linear least-squares approach to statistically model the observed r2 was implemented within R using this model:
(3)
yi=1a+4bdi+ei
Where yi is the value of r2 for the SNP pair i at linkage distance di in morgans. Parameters a and b were estimated iteratively using the least-squares method. Chromosome-specific Mb-to-centimorgan (cM) conversion rates were calculated based on total physical chromosome length stated on the UCSC Web site (https://genome.ucsc.edu/) and each chromosome’s length from a porcine linkage map from maps of four pedigrees (ILL, UIUC, USDA, ROS) [19]. Because UIUC consisted of Meishan and Yorkshire, we selected maps of this pedigree for use in our study. This model was applied to data for each chromosome and estimated parameter. Similar to Corbin and Shin, we combined the estimated parameters into a meta-analysis using an inverse variance method for pooling and a random-effects approach based on the DerSimonian-Laird method [18,20].

Ancestral Ne estimation

Equation (2) can predict the Ne at a given point in time expressed as the past generation [16]:
(4)
NT(t)=14c*(1rc2-1)
Where NT(t) is the Ne at t generations ago, c is the distance between markers in morgans, rc2 is the mean value of r2 for SNP pairs c morgans apart, and c = 1/2t when assuming linear growth [16]. To estimate NT(t), the number of prior generations was selected, and a suitable range for c was calculated. The binning process was designed to ensure sufficient SNP pairs within each bin to obtain a representative r2 mean. This process was performed for SNPs pooled across autosomes. Bin information used for estimating ancestral Ne is presented in Supplementary Table S3.

RESULTS

Characteristics of genotypic data

To investigate genetic variations in 2,470 Yorkshire pigs from four representative major Korean major GPP farms, chip analysis was performed on 33,418 SNPs after quality control processing. The number of SNPs per autosome remaining after filtering and imputation ranged from 831 to 3,834 and was closely related to chromosome length and total number of SNPs (Figure 2). The MAFs of remaining SNPs followed a uniform distribution and averaged (±standard deviation [SD]) 0.287±0.127. The average distance between SNP pairs (±SD) was 3,610±2,831 kb, with the distance between SNPs ranging from 0.009 kb to 24,564 kb (Supplementary Figure S3).

Estimation of linkage disequilibrium

The average LD (r2) of 3.2 million SNP pairs (±SD) analyzed in this study was 0.103 (±0.179). Distance of 28,969 pairs from all studied pairs was less than 50 kb. Overall, 51.22% of these 28,969 SNP pairs had r2>0.3, and 60.75% had r2>0.2. The average LD of distinct autosomes for SNPs at least 50 kb apart varied from r2 = 0.361 to r2 = 0.508, and the average LD for SNPs at least 5 Mb apart ranged from r2 = 0.071 to r2 = 0.185 (Supplementary Table S4). Although the aim of this study was not to compare LD in different chromosomes, we did observe some variation in the extent of LD by chromosome. We also determined that chromosomes 1, 13, and 14 had the highest average LD values, while chromosomes 10 and 12 showed the lowest average LD values. These results corroborate findings from a previous Yorkshire LD study.
Average LD values across all autosomes showed that the most rapid decrease was seen over the first ten bins (distance between SNPs ranged from 0 to 0.1Mb), with the mean r2 decreasing by approximately 40% (Figure 1a). The most rapid decrease was seen over the first five bins with the mean r2 decreasing by approximately 50% (Figure 1b). Mean r2 decreased much more slowly as distanced increased and was constant after 3 Mb of distance. According to our r2 calculations, 3,965 of the 3.2 million SNP pairs were in complete LD.

Determination of the relationship between linkage disequilibrium and single nucleotide polymorphism distance

Application of the non-linear regression model of declining LD in accordance with distance determined both parameters a and b in Equation (2) are significantly different from zero. For parameter estimation using r2, the mean estimate and 95% confidence interval by meta-analysis across autosomes for parameters a and b (±SD) were 2.71 (2.61; 2.82) and 122.87 (106.90; 138.84), respectively. Values for estimated parameters a and b in Equation (2) per chromosome are shown in Figure 3. Parameter b showed greater variability between chromosomes than parameter a. No such relationship was observed between each estimated parameter and chromosome length in centimorgans. We estimated predicted r2 by distance of SNP pairs using our estimated parameters a and b in Equation (2) and compared predicted r2 with observed r2 as performed in other studies [17,18]. We determined that the predicted r2 from the non-linear regression equation was similar to the mean observed r2 (Figure 4), indicating our estimated parameters in Equation (2) accurately represent Korean Yorkshire population history.

Ancestral Ne estimation

Figure 5 and Supplementary Figure S4 display estimated Ne at t generations ago. Based on the genomic data, the current Ne of Korean Yorkshire is approximately 122.87 (106.90; 138.84) individuals. Supplementary Figure S4 also shows that the reduction of Ne in Korean Yorkshire populations was continuous and gradual over the last 100 generations, ranging from 208.4 individuals to 122.87 individuals as determined in our study. Additionally, we found that Korean Yorkshire Ne has decreased by 99.6% over the last 10,000 generations (30,380.28 initial individuals) to the present estimate (Figure 5).

DISCUSSION

Four GGP farms in this study were operated by each different farm owners, but they were connected by sharing some semen. To show the genetic background of these population, we made principal component analysis plot including PC1 and PC2 (Figure 6). We could found that four groups were clearly distinguished, but the difference is not large, and the explanation variance of PC1 and PC2 is also very small. So we regarded four representative GGP farms as one Yorkshire population in Korea. We removed 9,303 SNPs by the quality control of Hardy-Weinberg equilibrium test for 2,470 Yorkshire pigs. We thought that the reason why so many SNPs removed by the quality control was related to number of heterozygous SNP alleles. A strong degree of selection in GGP would have reduced the effective population and had an impact on the number of heterozygous alleles, which could increase the degree of loss of heterozygosity. For this reason, we thought that many SNPs have been removed. Additionally, becasue core objective was estimation of Ne of Yorkshire population in Korea, we need one Yorkshire population in Korea and high quality SNPs of this population. And we thought that number of SNPs (33,418 SNPs) after quality control was sufficient for estimation of Ne when compared to finnish Yorkshire (Uimari [13]). So we thought that regarding 4 GGP farms as one Korean population certainly was appropriate.
We investigated the extent of LD and changes of the Ne of the Korean Yorkshire population based on whole-genome SNP data. The observed LD, measured as r2, extended for a long distance based on the adjacent 100 SNPs of each SNP studied in the genome. A previous study used large pedigree datasets and small genomic datasets [13], but instead of pedigree data, we used large-scale genomic data from swine in GGP farms to characterize LD and estimate Ne. Because domesticated pigs have long been strongly and artificially selected, observed LD in the Korean Yorkshire population was higher with shorter genomic distances and more extensive compared to human populations [21]. Declining LD in the Korean Yorkshire population is consistent with previous studied pigs [13,22] as well as other domesticated animals [17,18].
We estimated Ne of Korean Yorkshire swine based on a non-linear regression model that describes the relationship between linkage distance and LD. Estimating Ne using this equation raises difficulties in handling values in the limits of the parameter space, because if r2 = 0.0, the estimated Ne is infinite, and if r2 = 1.0, Ne is zero. In this study, we calculated r2 for the adjacent 100 SNPs of each SNP to decrease bias in handling these values. The results from this simplified approach yielded quite similar estimates of Ne as other previous studies. Another concern related to the relationship between estimated LD and distance between SNPs lies in the accuracy of the porcine reference genome used in this study. The order and distance between SNPs in the commercial Illumina PorcineSNP60 v2 BeadChip will likely be refined because the reference genome version will also be updated. However, the bias from incorrectly ordered SNPs or wrong SNP distances between SNPs may be minimized by the large number of SNP pairs used in this study, so some overestimated or underestimated LD can be overlooked. Because the relationship between genetic and physical distances varies across chromosomes and chromosomal regions, we inferred the Mb/cM ratio per chromosome using physical map position information from the porcine reference genome and from a previous Yorkshire pig study [19]. We then used genetic distances based on physical distances to estimate Ne and could obtain more reliable Ne estimates with such detailed estimates of genetic distances between SNPs. Finally, one study reported that a limited sample size can bias the estimates of r2 and recommends correcting the estimates of r2 for sample size n (r2–1/2n) and using the equation of Sved [15]. However, our sample size was sufficient enough to correct estimates of r2. To estimate Ne of Korean Yorkshire pigs, we used alternative version of this equation further derived by Tenesa [23], which adds a new parameter a to account for mutations. Based on the new formula, the initial value of parameter a = 2 in estimations of parameters using a non-linear regression model in R. Regarding variance heterogeneity of the observed r2, the variance declined with increased distances between SNPs, which may have impacted our results when estimating parameter b in Equation (2) (Supplementary Figure S2). A significant, negative relationship between chromosome length and estimates of parameter b obtained from a non-linear model have been observed [17], while others have noted a positive relationship in domestic livestock species [24,25] or did not investigate directionality of the relationship [18]. Because the evolutionary history of each species and breed is different, the relationship between chromosome length and parameter b is also different per population. In this study, all marker pairs were calculated in each bin so that r2 was not affected by chromosome length. These results are in agreement with the Yorkshire LD characterization findings of Uimari [13]. We did not observe a significant relationship between chromosome length and estimates of b in our study population.
Our estimate of b represents an estimated Ne assuming a constant present population size, because we used genetic data from the Korean Yorkshire population consisting of pigs in major GGP farms. In the calculation of Ne, b in Equation (2) represents a conceptual average of Ne over the period inferred from the SNP pair distance ranges per chromosome [26]. We regarded parameter b combined by meta-analysis as reflecting the current Ne of the Korean Yorkshire population.
As Table 1 showed, we produced SNP chip data of 2,470 individuals in between 2011 and 2015. And approximately 87% (2,149 individuals) of total SNP data was produced in between 2014 and 2015. Because sampling period was short, we regarded sampled individuals in this study as “current generation population” in Korean Yorksrhie population. So we thought that we need not to sort 2,470 individuals data according to accurate generation. Instead, we infered long time generation-related change of Ne using another method in this study. Hayes [16] reported that the degree of linkage equilibrium according to genetic distance had reflected genetic diversity of past generation. After r2 estimation using SNP data, we divided r2 by distance and inferred Ne of past generation. The LD over greater genetic distances reflects a population’s recent history, whereas LD over shorter distances depends on the Ne many generations ago [16,27] (Supplementary Table 3). Historical Ne estimation suggests a linear population as reported in a previous study (Figure 5; Supplementary Figure S4) [16]. The observed pattern displayed a consistent decrease in Ne from 100 generations ago to the present, decreasing by 99.6% from 10,000 to 100 generation ago. Several explanations exist for this pattern, including bottlenecks associated with domestication, selection, and breed administration, business strategy, and endangerment of the breed. Therefore, our results should be considered in context of the demographic history of the Yorkshire population in Korea. The reliability of predicting changes in Ne over time depends both on technical implementation and proper iteration based on previous studies using this approach [17,18].
We aimed to characterize LD and Ne in a Korean Yorkshire population using genomic data from thousands of individuals. Our observed LD patterns are similar to the average value presented by Du for six commercial pig lines and Uimari for Finnish Yorkshire pig breeds [13,22]. The overall LD in Finnish Yorkshire breed appears to be stronger than in Korean Yorkshire pigs. Because the Korean Yorkshire population consists of seed pigs from several original Yorkshire breeds, the genetic diversity of Korean Yorkshire pigs is greater than that of the single Finnish Yorkshire breed, and the LD of Finnish Yorkshire pigs is higher than that of Korean Yorkshires.
The minimum number of breeding animals recommended by the UN Food and Agriculture Organization is 50, although Meuwissen suggested this number is the lower limit for a critical population size, proposing that the actual size should be between 50 and 100 [28]. The current Ne of the Korean Yorkshire population is 122.87, which is sufficient to maintain the population’s viability. The population’s genetic variation enables an acceptable inbreeding rate, including compromising genetic gain in commercially important traits. This genetic variation is necessary to apply methods that maximize selection efficacy with a fixed rate of inbreeding or optimize the use of genetic resources from the parental generation [29].
When we apply a new genetic method, such as genomic selection, for estimating breeding values, the Ne may be very small or continually decrease [8]. Therefore, one must carefully consider appropriate breed management methods to avoid inbreeding. Although attenuated selection can affect short-term genetic gain, it is essential for maintaining the long-term genetic variability of the Korean Yorkshire population. Long-term continuous monitoring would also be needed to maintain the pig population to avoid an unintended reduction of Ne. The best way to preserve a sustainable population is to ensure sure its production populations maintain a sufficient Ne.

Notes

CONFLICT OF INTEREST

We certify that there is no conflict of interest with any financial organization regarding the material discussed in the manuscript. Kim SH is an employee of PigGene Korea Inc.

ACKNOWLEDGMENTS

This study was supported by a grant from the KRIBB initiative program and the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (NRF-2017R1C1B3007144) and a grant from the Next-Generation BioGreen 21 Program (PJ01337702), Rural Development Administration, Republic of Korea.

Figure 1
Average linkage disequilibrium (LD) plotted against the median of the distance bin range (Mb). (a) Distances ranged from 0 to 0.5 Mb, and r2 values were averaged using bins of 0.01 Mb and pooled over autosomes. (b) Distances ranged from 0 to 5 Mb, and r2 values were averaged using bins of 0.1 Mb and pooled over autosomes.
ajas-31-12-1843f1.jpg
Figure 2
Number of single nucleotide polymorphisms (SNPs) per chromosome after quality control processing and imputation for estimation of effective population size. The x- and y-axes indicate SNP numbers and chromosomes, respectively.
ajas-31-12-1843f2.jpg
Figure 3
Parameter estimates from Equation (3) plotted against chromosome length in centimorgans (cM) according to the referenced porcine linkage map r2 used in our study. (a) Estimates of parameter a in Equation (3) plotted against chromosome length (cM). (b) Estimates of parameter b in Equation (3) plotted against chromosome length (cM).
ajas-31-12-1843f3.jpg
Figure 4
Predicted r2 (Black x) vs observed r2 (sky blue circle) against mean distance between markers (cM expressed on a log scale). Predicted r2 calculated using Equation (3) with a = 2.71 and b = 122.87.
ajas-31-12-1843f4.jpg
Figure 5
Average estimated Ne (effective population size) plotted against number of past generations truncated at 100,000 generations using linkage disequilibrium (r2) plotted on a log scale.
ajas-31-12-1843f5.jpg
Figure 6
PCA (Principal component analysis) plot of Yorkshire population (2,470 Individuals) in this study. The 2,470 individuals Yorkshire population belonged to 4 GGPs and each color indicated each farms.
ajas-31-12-1843f6.jpg
Table 1
Sample information of Yorkshire population (2,470 Individuals) in this study
Year Number of individual Farm1 Farm2 Farm3 Farm4 Female Male
2011 5 0 0 1 4 1 4
2012 137 4 0 59 74 134 3
2013 179 23 0 124 32 167 12
2014 687 122 430 109 26 674 13
2015 1,462 1,324 112 19 7 1,166 296
Total 2,470 1,473 542 312 143 2,142 328

REFERENCES

1. Wang J. Estimation of effective population sizes from data on genetic markers. Philos Trans R Soc Lond B Biol Sci 2005; 360:1395–409.
crossref pmid pmc
2. Qian D, Huiying L, Xinyun L, et al. A genome-wide association study of five meat quality traits in Yorkshire pigs. Front Agric Sci Eng 2014; 1:137–43.
crossref
3. Meng Q, Wang K, Liu X, et al. Identification of growth trait related genes in a Yorkshire purebred pig population by genome-wide association studies. Asian-Australas J Anim Sci 2017; 30:462–9.
crossref pmid pdf
4. Lee T, Shin DH, Cho S, et al. Genome-wide Association Study of Integrated Meat Quality-related Traits of the Duroc Pig Breed. Asian-Australas J Anim Sci 2014; 27:303–9.
crossref pmid pmc pdf
5. Tribout T, Larzul C, Phocas F. Efficiency of genomic selection in a purebred pig male line. J Anim Sci 2012; 90:4164–76.
crossref pmid pdf
6. Lillehammer M, Meuwissen TH, Sonesson AK. Genomic selection for two traits in a maternal pig breeding scheme. J Anim Sci 2013; 91:3079–87.
crossref pmid pdf
7. Christensen OF, Madsen P, Nielsen B, Ostersen T, Su G. Single-step methods for genomic evaluation in pigs. Animal 2012; 6:1565–71.
crossref pmid
8. Hayes B, Goddard M. Prediction of total genetic value using genome-wide dense marker maps. Genetics 2001; 157:1819–29.
crossref pmid pmc pdf
9. Falconer DS. Introduction to quantitative genetics. Edinburgh and London, UK: Oliver and Boyd; 1960.

10. Badke YM, Bates RO, Ernst CW, Schwab C, Steibel JP. Estimation of linkage disequilibrium in four US pig breeds. BMC Genomics 2012; 13:24
crossref pmid pmc
11. Veroneze R, Lopes PS, Guimaraes SE, et al. Linkage disequilibrium and haplotype block structure in six commercial pig lines. J Anim Sci 2013; 91:3493–501.
crossref pmid pdf
12. Wang L, Sorensen P, Janss L, Ostersen T, Edwards D. Genome-wide and local pattern of linkage disequilibrium and persistence of phase for 3 Danish pig breeds. BMC Genet 2013; 14:115
crossref pmid pmc
13. Uimari P, Tapio M. Extent of linkage disequilibrium and effective population size in Finnish Landrace and Finnish Yorkshire pig breeds. J Anim Sci 2011; 89:609–14.
crossref pmid pdf
14. Howie BN, Donnelly P, Marchini J. A flexible and accurate genotype imputation method for the next generation of genome-wide association studies. PLoS Genet 2009; 5:e1000529
crossref pmid pmc
15. Sved JA. Linkage disequilibrium and homozygosity of chromosome segments in finite populations. Theor Popul Biol 1971; 2:125–41.
crossref pmid
16. Hayes BJ, Visscher PM, McPartlan HC, Goddard ME. Novel multilocus measure of linkage disequilibrium to estimate past effective population size. Genome Res 2003; 13:635–43.
crossref pmid pmc
17. Corbin LJ, Blott SC, Swinburne JE, Vaudin M, Bishop SC, Woolliams JA. Linkage disequilibrium and historical effective population size in the Thoroughbred horse. Anim Genet 2010; 41:Suppl 28–15.
crossref pmid
18. Shin DH, Cho KH, Park KD, Lee HJ, Kim H. Accurate estimation of effective population size in the korean dairy cattle based on linkage disequilibrium corrected by genomic relationship matrix. Asian-Australas J Anim Sci 2013; 26:1672–9.
crossref pmid pmc pdf
19. Tortereau F, Servin B, Frantz L, et al. A high density recombination map of the pig reveals a correlation between sex-specific recombination and GC content. BMC Genomics 2012; 13:586
crossref pmid pmc
20. DerSimonian R, Laird N. Meta-analysis in clinical trials. Controlled clinical trials 1986; 7:177–88.
crossref pmid
21. Shifman S, Kuypers J, Kokoris M, Yakir B, Darvasi A. Linkage disequilibrium patterns of the human genome across populations. Hum Mol Genet 2003; 12:771–6.
crossref pmid pdf
22. Du FX, Clutter AC, Lohuis MM. Characterizing linkage disequilibrium in pig populations. Int J Biol Sci 2007; 3:166–78.
crossref pmid pmc
23. Tenesa A, Navarro P, Hayes BJ, et al. Recent human effective population size estimated from linkage disequilibrium. Genome Res 2007; 17:520–6.
crossref pmid pmc
24. Khatkar MS, Nicholas FW, Collins AR, et al. Extent of genome-wide linkage disequilibrium in Australian Holstein-Friesian cattle based on a high-density SNP panel. BMC Genomics 2008; 9:187
crossref pmid pmc
25. Muir W, Wong G, Zhang Y, et al. Review of the initial validation and characterization of a 3K chicken SNP array. World’s Poult Sci J 2008; 64:219–26.
crossref
26. Toosi A, Fernando RL, Dekkers JC. Genomic selection in admixed and crossbred populations. J Anim Sci 2010; 88:32–46.
crossref pmid pdf
27. Hill WG. Estimation of effective population size from data on linkage disequilibrium. Genet Res 1981; 38:209–16.
crossref
28. Meuwissen T. Genetic management of small populations: a review. Acta Agric Scand A Anim Sci 2009; 59:71–9.
crossref
29. Colleau JJ, Tribout T. Optimized management of genetic variability in selected pig populations. J Anim Breed Genet 2008; 125:291–300.
crossref pmid


Editorial Office
Asian-Australasian Association of Animal Production Societies(AAAP)
Room 708 Sammo Sporex, 23, Sillim-ro 59-gil, Gwanak-gu, Seoul 08776, Korea   
TEL : +82-2-888-6558    FAX : +82-2-888-6559   
E-mail : editor@animbiosci.org               

Copyright © 2024 by Asian-Australasian Association of Animal Production Societies.

Developed in M2PI

Close layer
prev next