Measuring the extent of linkage disequilibrium in commercial pig populations

To evaluate the extent of linkage disequilibrium in domestic pigs, we genotyped 33 and 44 unrelated individuals from two commercial populations for 29 and five microsatellite markers located on chromosomes 15 and 2 respectively. A high proportion of marker pairs up to 40 cM apart exhibited significant linkage disequilibrium in both populations. Pair-wise r(2) values averaged between 0.15 and 0.50 (depending on chromosome and population) for markers <1 cM apart and declined to values of 0.05 for more distant syntenic markers. Our results suggest that both populations underwent a bottleneck approximately 20 generations ago, which reduced the effective population size from thousands to <200 animals.     

Estimating effective population size from linkage disequilibrium: severe bias in small samples

Effective population size (N _e) is a central concept in evolutionary biology and conservation genetics. It predicts rates of loss of neutral genetic variation, fixation of deleterious and favourable alleles, and the increase of inbreeding experienced by a population. A method exists for the estimation of N _e from the observed linkage disequilibrium between unlinked loci in a population sample. While an increasing number of studies have applied this method in natural and managed populations, its reliability has not yet been evaluated. We developed a computer program to calculate this estimator of N _e using the most widely used linkage disequilibrium algorithm and used simulations to show that this estimator is strongly biased when the sample size is small (<‰100) and below the true N _e. This is probably due to the linkage disequilibrium generated by the sampling process itself and the inadequate correction for this phenomenon in the method. Results suggest that N _e estimates derived using this method should be regarded with caution in many cases. To improve the method’s reliability and usefulness we propose a way to determine whether a given sample size exceeds the population N _e and can therefore be used for the computation of an unbiased estimate.     

The extent of linkage disequilibrium in a large cattle population of western Africa and its consequences for association studies

Several previous studies concluded that linkage disequilibrium (LD) in livestock populations from developed countries originated from the impact of strong selection. Here, we assessed the extent of LD in a cattle population from western Africa that was bred in an extensive farming system. The analyses were performed on 363 individuals in a Bos indicus x Bos taurus population using 42 microsatellite markers on BTA04, BTA07 and BTA13. A high level of expected heterozygosity (0.71), a high mean number of alleles per locus (9.7) and a mild shift in Hardy-Weinberg equilibrium were found. Linkage disequilibrium extended over shorter distances than what has been observed in cattle from developed countries. Effective population size was assessed using two methods; both methods produced large values: 1388 when considering heterozygosity (assuming a mutation rate of 10(-3)) and 2344 when considering LD on whole linkage groups (assuming a constant population size over generations). However, analysing the decay of LD as a function of marker spacing indicated a decreasing trend in effective population size over generations. This decrease could be explained by increasing selective pressure and/or by an admixture process. Finally, LD extended over small distances, which suggested that whole-genome scans will require a large number of markers. However, association studies using such populations will be effective.     

利用全基因组连锁不平衡估计中国荷斯坦牛有效群体大小

有效群体大小是群体遗传学研究的一个重要内容,有助于我们更清楚地了解群体的遗传变异、进化和复杂性状的遗传机制等。随着高密度SNP标记的出现,越来越多的研究利用SNP标记间连锁不平衡估计有效群体大小。文章采集北京地区中国荷斯坦牛2 093个样本,并利用牛SNP芯片(Illumina BovineSNP50,含5 4001 SNPs)进行基因型测定,估计不同世代中国荷斯坦牛的有效群体大小。质量控制标准设定为SNP检出率0.95,最小等位基因频率>0.05,样本检出率0.95,哈代温伯格平衡检验显著性水平P<0.0001。经过质量控制,共1 968个样本和38 796个SNPs用于连锁不平衡分析。文章选取SNP间距0.1、0.2、0.5、1、2、5、10、15(Mb),估计中国荷斯坦牛在4世代之前有效群体大小。结果表明,中国荷斯坦牛的有效群体呈逐代下降趋势,至4世代前,中国荷斯坦牛平均有效群体为45头左右。     

Genome-wide signatures of population bottlenecks and diversifying selection in European wolves

Genomic resources developed for domesticated species provide powerful tools for studying the evolutionary history of their wild relatives. Here we use 61K single-nucleotide polymorphisms (SNPs) evenly spaced throughout the canine nuclear genome to analyse evolutionary relationships among the three largest European populations of grey wolves in comparison with other populations worldwide, and investigate genome-wide effects of demographic bottlenecks and signatures of selection. European wolves have a discontinuous range, with large and connected populations in Eastern Europe and relatively smaller, isolated populations in Italy and the Iberian Peninsula. Our results suggest a continuous decline in wolf numbers in Europe since the Late Pleistocene, and long-term isolation and bottlenecks in the Italian and Iberian populations following their divergence from the Eastern European population. The Italian and Iberian populations have low genetic variability and high linkage disequilibrium, but relatively few autozygous segments across the genome. This last characteristic clearly distinguishes them from populations that underwent recent drastic demographic declines or founder events, and implies long-term bottlenecks in these two populations. Although genetic drift due to spatial isolation and bottlenecks seems to be a major evolutionary force diversifying the European populations, we detected 35 loci that are putatively under diversifying selection. Two of these loci flank the canine platelet-derived growth factor gene, which affects bone growth and may influence differences in body size between wolf populations. This study demonstrates the power of population genomics for identifying genetic signals of demographic bottlenecks and detecting signatures of directional selection in bottlenecked populations, despite their low background variability.     

Analysis of microsatellite variation in Pinus radiata reveals effects of genetic drift but no recent bottlenecks

Most conifer species occur in large continuous populations, but radiata pine, Pinus radiata, occurs only in five disjunctive natural populations in California and Mexico. The Mexican island populations were presumably colonized from the mainland millions of years ago. According to Axelrod (1981), the mainland populations are relicts of an earlier much wider distribution, reduced some 8,000 years ago, whereas according to Millar (1997, 2000), the patchy metapopulation-like structure is typical of the long-term population demography of the species. We used 19 highly polymorphic microsatellite loci to describe population structure and to search for signs of the dynamics of population demography over space and time. Frequencies of null alleles at microsatellite loci were estimated using an approach based on the probability of identity by descent. Microsatellite genetic diversities were high in all populations [expected heterozygosity (H(e)) = 0.68-0.77], but the island populations had significantly lower estimates. Variation between loci in genetic differentiation (F(ST)) was high, but no locus deviated statistically significantly from the rest at an experiment wide level of 0.05. Thus, all loci were included in subsequent analysis. The average differentiation was measured as F(ST) = 0.14 (SD 0.012), comparable with earlier allozyme results. The island populations were more diverged from the other populations and from an inferred common ancestral gene pool than the mainland ones. All populations showed a deficiency of expected heterozygosity given the number of alleles, the mainland populations more so than the island ones. The results thus do not support a recent important contraction in the mainland range of radiata pine.     

Genome-Wide Estimation of Linkage Disequilibrium from Population-Level High-Throughput Sequencing Data

Rapidly improving sequencing technologies provide unprecedented opportunities for analyzing genome-wide patterns of polymorphisms. In particular, they have great potential for linkage-disequilibrium analyses on both global and local genetic scales, which will substantially improve our ability to derive evolutionary inferences. However, there are some difficulties with analyzing high-throughput sequencing data, including high error rates associated with base reads and complications from the random sampling of sequenced chromosomes in diploid organisms. To overcome these difficulties, we developed a maximum-likelihood estimator of linkage disequilibrium for use with error-prone sampling data. Computer simulations indicate that the estimator is nearly unbiased with a sampling variance at high coverage asymptotically approaching the value expected when all relevant information is accurately estimated. The estimator does not require phasing of haplotypes and enables the estimation of linkage disequilibrium even when all individual reads cover just single polymorphic sites.     

Genetic monitoring reveals temporal stability over 30 years in a small, lake-resident brown trout population

Knowledge of the degree of temporal stability of population genetic structure and composition is important for understanding microevolutionary processes and addressing issues of human impact of natural populations. We know little about how representative single samples in time are to reflect population genetic constitution, and we explore the temporal genetic variability patterns over a 30-year period of annual sampling of a lake-resident brown trout (Salmo trutta) population, covering 37 consecutive cohorts and five generations. Levels of variation remain largely stable over this period, with no indication of substructuring within the lake. We detect genetic drift, however, and the genetically effective population size (N(e)) was assessed from allele-frequency shifts between consecutive cohorts using an unbiased estimator that accounts for the effect of overlapping generation. The overall mean N(e) is estimated as 74. We find indications that N(e) varies over time, but there is no obvious temporal trend. We also estimated N(e) using a one-sample approach based on linkage disequilibrium (LD) that does not account for the effect of overlapping generations. Combining one-sample estimates for all years gives an N(e) estimate of 76. This similarity between estimates may be coincidental or reflecting a general robustness of the LD approach to violations of the discrete generations assumption. In contrast to the observed genetic stability, body size and catch per effort have increased over the study period. Estimates of annual effective number of breeders (N(b)) correlated with catch per effort, suggesting that genetic monitoring can be used for detecting fluctuations in abundance.     

Protein electrophoresis as a management tool: Detection of hybridization between banteng (Bos javanicus d'Alton) and domestic cattle

The ability of protein electrophoresis to detect hybridization was evaluated by a detailed examination of 30 banteng-domestic cattle hybrids. These hybrids were run as unknowns in a sample of 53 individuals that included both banteng and domestic cattle. Twenty-six of the 30 hybrid animals could be unambiguously identified as such, while the remaining hybrids were indistinguishable from domestic cattle. Theoretical analysis showed that the sensitivity of the protein assays at the individual level increased slowly with the number of available markers. The probability of detecting limited amounts of hybrid influence (5–10%) in an individual was less than 70% for any biologically realistic number of protein markers. Hybridization was much easier to detect at the population level, where the sensitivity of an electrophoretic assay increased rapidly with increased sample size. The combined results of the protein electrophoresis and the theoretical analysis suggest that there is an 94.8% probability that there are essentially no genes from domestic cattle in the North American banteng herd at the present time. Zoos interested in obtaining new animals should examine the source population as well as the individuals in question whenever it is feasible. An electrophoretic survey would also provide useful information when a species survival plan is organized.     

Effective population size of natural populations of Drosophila buzzatii, with a comparative evaluation of nine methods of estimation

Allozyme and microsatellite data from numerous populations of Drosophila buzzatii have been used (i) to determine to what degree N(e) varies among generations within populations, and among populations, and (ii) to evaluate the congruence of four temporal and five single-sample estimators of N(e) . Effective size of different populations varied over two orders of magnitude, most populations are not temporally stable in genetic composition, and N(e) showed large variation over generations in some populations. Short-term N(e) estimates from the temporal methods were highly correlated, but the smallest estimates were the most precise for all four methods, and the most consistent across methods. Except for one population, N(e) estimates were lower when assuming gene flow than when assuming populations that were closed. However, attempts to jointly estimate N(e) and immigration rate were of little value because the source of migrants was unknown. Correlations among the estimates from the single-sample methods generally were not significant although, as for the temporal methods, estimates were most consistent when they were small. These single-sample estimates of current N(e) are generally smaller than the short-term temporal estimates. Nevertheless, population genetic variation is not being depleted, presumably because of past or ongoing migration. A clearer picture of current and short-term effective population sizes will only follow with better knowledge of migration rates between populations. Different methods are not necessarily estimating the same N(e) , they are subject to different bias, and the biology, demography and history of the population(s) may affect different estimators differently.