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.     

Patterns of linkage disequilibrium in mitochondrial DNA of 16 ruminant populations

Mitochondrial DNA (mtDNA) is a widely employed molecular tool in phylogeography, in the inference of human evolutionary history, in dating the domestication of livestock and in forensic science. In humans and other vertebrates the popularity of mtDNA can be partially attributed to an assumption of strict maternal inheritance, such that there is no recombination between mitochondrial lineages. The recent demonstration that linkage disequilibrium (LD) declines as a function of distance between polymorphic sites in hominid mitochondrial genomes has been interpreted as evidence of recombination between mtDNA haplotypes, and hence nonclonal inheritance. However, critics of mtDNA recombination have suggested that this association is an artefact of an inappropriate measure of LD or of sequencing error, and subsequent studies of other populations have failed to replicate the initial finding. Here we report the analysis of 16 ruminant populations and present evidence that LD significantly declines with distance in five of them. A meta-analysis of the data indicates a nonsignificant trend of LD declining with distance. Most of the earlier criticisms of patterns between LD and distance in hominid mtDNA are not applicable to this data set. Our results suggest that either ruminant mtDNA is not strictly clonal or that compensatory selection has influenced patterns of variation at closely linked sites within the mitochondrial control region. The potential impact of these processes should be considered when using mtDNA as a tool in vertebrate population genetic, phylogenetic and forensic studies.     

Physical linkage and mate preference generate linkage disequilibrium for behavioral isolation in two parapatric crickets

Behavioral isolation is a potent barrier to gene flow and a source of striking diversity in the animal kingdom. However, it remains unclear if the linkage disequilibrium (LD) between sex‐specific traits required for behavioral isolation results mostly from physical linkage between signal and preference loci or from directional mate preferences. Here, we test this in the field crickets Gryllus rubens and G. texensis. These closely related species diverged with gene flow and have strongly differentiated songs and preference functions for the mate calling song rhythm. We map quantitative trait loci for signal and preference traits (pQTL) as well as for gene expression associated with these traits (eQTL). We find strong, positive genetic covariance between song traits and between song and preference. Our results show that this is in part explained by incomplete physical linkage: although both linked pQTL and eQTL couple male and female traits, major effect loci for different traits were never on the same chromosome. We suggest that the finely tuned, highly divergent preference functions are likely an additional source of LD between male and female traits in this system. Furthermore, pleiotropy of gene expression presents an underappreciated mechanism to link sexually dimorphic phenotypes.     

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.     

An evaluation of the methods to estimate effective population size from measures of linkage disequilibrium

In 1971, John Sved derived an approximate relationship between linkage disequilibrium (LD) and effective population size for an ideal finite population. This seminal work was extended by Sved and Feldman (Theor Pop Biol 4, 129, 1973) and Weir and Hill (Genetics 95, 477, 1980) who derived additional equations with the same purpose. These equations yield useful estimates of effective population size, as they require a single sample in time. As these estimates of effective population size are now commonly used on a variety of genomic data, from arrays of single nucleotide polymorphisms to whole genome data, some authors have investigated their bias through simulation studies and proposed corrections for different mating systems. However, the cause of the bias remains elusive. Here, we show the problems of using LD as a statistical measure and, analogously, the problems in estimating effective population size from such measure. For that purpose, we compare three commonly used approaches with a transition probability‐based method that we develop here. It provides an exact computation of LD. We show here that the bias in the estimates of LD and effective population size are partly due to low‐frequency markers, tightly linked markers or to a small total number of crossovers per generation. These biases, however, do not decrease when increasing sample size or using unlinked markers. Our results show the issues of such measures of effective population based on LD and suggest which of the method here studied should be used in empirical studies as well as the optimal distance between markers for such estimates.     

Impact of polymorphic transposable elements on linkage disequilibrium along chromosomes

Recombination and selection drive the extent of linkage disequilibrium (LD) among loci and therefore affect the reshuffling of adaptive genetic variation. However, it is poorly known to what extent the enrichment of transposable elements (TEs) in recombinationally‐inert regions reflects their inefficient removal by purifying selection and whether the presence of polymorphic TEs can modify the local recombination rate. In this study, we investigate how TEs and recombination interact at fine scale along chromosomes and possibly support linked selection in natural populations. Whole‐genome sequencing data of 304 individuals from nearby alpine populations of Arabis alpina were used to show that the density of polymorphic TEs is specifically correlated with local LD along chromosomes. Consistent with TEs modifying recombination, the characterization of 28 such LD blocks of up to 5.5 Mb in length revealed strong evidence of selective sweeps at a few loci through either site frequency spectrum or haplotype structure. A majority of these blocks were enriched in genes related to ecologically relevant functions such as responses to cold, salt stress or photoperiodism. In particular, the S‐locus (i.e., supergene responsible for strict outcrossing) was identified in a LD block with high levels of polymorphic TEs and evidence of selection. Another such LD block was enriched in cold‐responding genes and presented evidence of adaptive loci related to photoperiodism and flowering being increasingly linked by polymorphic TEs. These results are consistent with the hypothesis that TEs modify recombination landscapes and thus interact with selection in driving blocks of linked adaptive loci in natural populations.     

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.     

An evaluation of a novel estimator of linkage disequilibrium

The analysis of systems involving many loci is important in population and quantitativegenetics. An important problem is the study of linkage disequilibrium (LD), a conceptrelevant in genome-enabled prediction of quantitative traits and in exploration ofmarker–phenotype associations. This article introduces a new estimator of a LDparameter (ρ²) that is much easier to compute than a maximumlikelihood (or Bayesian) estimate of a tetra-choric correlation. We examined theconjecture that the sampling distribution of the estimator of ρ²could be less frequency dependent than that of the estimator ofr², a widely used metric for assessing LD. This was donevia an empirical evaluation of LD in 806 Holstein–Friesian cattle using 771single-nucleotide polymorphism (SNP) markers and of HapMap III data on 21 991 SNPs(chromosome 3) observed in 88 unrelated individuals from Tuscany. Also, 1600 haplotypesover a region of 1 Mb simulated under the coalescent were used to estimate LD usingthe two measures. Subsequently, a simulation study compared the new estimator with that ofr² using several scenarios of LD and allelic frequencies.From these studies, it is concluded that ρ² provides a usefulmetric for the study of LD as the distribution of its estimator is less frequencydependent than that of the standard estimator of r².     

Associative overdominance: Evaluating the effects of inbreeding and linkage disequilibrium

Expressions are obtained for the expected phenotypic values of homozygous and heterozygous genotypes for a neutral marker locus linked to a locus segregating for a recessive deleterious gene. The phenotypic values are functions of the allele frequencies at the marker locus, the inbreeding coefficient and the degree of association of the deleterious gene with the marker alleles. The analysis is extended to more than two alleles at the marker locus. Either linkage disequilibrium or inbreeding alone can produce an apparent superiority of heterozygotes for the marker locus (unless specified otherwise, the terms ‘homozygote’ and ‘heterozygote’ will refer to the marker locus). The effect of linkage disequilibrium on the difference between the heterozygote and homozygote values can be positive (associative overdominance) or negative (associative underdominance), depending on the frequencies of the marker alleles and the degree of their association with the deleterious gene. Inbreeding has always a positive effect. In general, the expected value of a homozygote is a positive function of its allele frequency. When the various homozygous genotypes are combined into one class and the various heterozygous genotypes into another, the phenotypic difference of the two classes is a function of the evenness of the allelic frequency distribution. Inbreeding is a more likely explanation of associative overdominance if the frequency of the deleterious gene is low, but its effect on the character high. Conversely, linkage disequilibrium is more likely if the frequency is high and the effect low. The degrees of association between marker alleles and the deleterious gene can, in principle, be estimated from the observed phenotypic scores and used to calculate expected multi-locus genotype scores. This could provide the basis for statistical tests of the associative overdominance hypothesis as an explanation of observed correlations between multi-locus heterozygosity and phenotypic traits.     

speed‐ne: Software to simulate and estimate genetic effective population size (Ne) from linkage disequilibrium observed in single samples

The genetic effective population size, N_e, can be estimated from the average gametic disequilibrium () between pairs of loci, but such estimates require evaluation of assumptions and currently have few methods to estimate confidence intervals. speed‐ne is a suite of matlab computer code functions to estimate from with a graphical user interface and a rich set of outputs that aid in understanding data patterns and comparing multiple estimators. speed‐ne includes functions to either generate or input simulated genotype data to facilitate comparative studies of estimators under various population genetic scenarios. speed‐ne was validated with data simulated under both time‐forward and time‐backward coalescent models of genetic drift. Three classes of estimators were compared with simulated data to examine several general questions: what are the impacts of microsatellite null alleles on , how should missing data be treated, and does disequilibrium contributed by reduced recombination among some loci in a sample impact . Estimators differed greatly in precision in the scenarios examined, and a widely employed estimator exhibited the largest variances among replicate data sets. speed‐ne implements several jackknife approaches to estimate confidence intervals, and simulated data showed that jackknifing over loci and jackknifing over individuals provided ~95% confidence interval coverage for some estimators and should be useful for empirical studies. speed‐ne provides an open‐source extensible tool for estimation of from empirical genotype data and to conduct simulations of both microsatellite and single nucleotide polymorphism (SNP) data types to develop expectations and to compare estimators.     

Properties of 2-locus genealogies and linkage disequilibrium in temporally structured samples

Archeogenetics has been revolutionary, revealing insights into demographic history and recent positive selection. However, most studies to date have ignored the nonrandom association of genetic variants at different loci (i.e. linkage disequilibrium). This may be in part because basic properties of linkage disequilibrium in samples from different times are still not well understood. Here, we derive several results for summary statistics of haplotypic variation under a model with time-stratified sampling: (1) The correlation between the number of pairwise differences observed between time-staggered samples (πΔt) in models with and without strict population continuity; (2) The product of the linkage disequilibrium coefficient, D, between ancient and modern samples, which is a measure of haplotypic similarity between modern and ancient samples; and (3) The expected switch rate in the Li and Stephens haplotype copying model. The latter has implications for genotype imputation and phasing in ancient samples with modern reference panels. Overall, these results provide a characterization of how haplotype patterns are affected by sample age, recombination rates, and population sizes. We expect these results will help guide the interpretation and analysis of haplotype data from ancient and modern samples.     

Spatial genetic patterns generated by two admixing genetic lineages: a simulation study

Two formerly geographically separated lineages of the zebra mussel Dreissena polymorpha had been given the opportunity to mix extensively across the newly built German Main–Danube canal. We had monitored this admixture of mussel lineages and had described spatial patterns in different genetic measures [Müller et al. 2001 (Heredity, 86: 103); 2002 (Proc. R. Soc. Lond., 269: 1139)]. Here, we present an individual-based model to assess the potential of spatial genetic patterns of detecting and quantifying admixture of mussel lineages. Genetic measures studied are (1) allele frequencies, (2) deviations from Hardy–Weinberg expectations of loci (deficit of heterozygotes, HWD) and (3) linkage disequilibria between unlinked loci (LD). For allele frequencies, we observed a cline over the zone of admixture in all simulations of mixing mussel lineages suggesting that these are appropriate for verification of their mixture. The point of the first contact between lineages was always detectable from their intermediate allele frequencies. LD and HWD were only spatially informative for diagnostic loci or loci with very strong differences in allele frequencies of lineages. For such loci, the probability of disequilibria was highest where lineages had met and decreased towards both sources of lineages Main and Danube. The overall probability of detecting any disequilibrium was higher for LD than for HWD and increased with an increasing rate of genetic interchange. Our simulation results are corroborated by our zebra mussel data and studies from literature. They are applicable to any case of two known linearly mixing genetic lineages.     

The heterogeneous levels of linkage disequilibrium in white spruce genes and comparative analysis with other conifers

In plants, knowledge about linkage disequilibrium (LD) is relevant for the design of efficient single-nucleotide polymorphism arrays in relation to their use in population and association genomics studies. Previous studies of conifer genes have shown LD to decay rapidly within gene limits, but exceptions have been reported. To evaluate the extent of heterogeneity of LD among conifer genes and its potential causes, we examined LD in 105 genes of white spruce (Picea glauca) by sequencing a panel of 48 haploid megagametophytes from natural populations and further compared it with LD in other conifer species. The average pairwise r(2) value was 0.19 (s.d.=0.19), and LD dropped quickly with a half-decay being reached at a distance of 65 nucleotides between sites. However, LD was significantly heterogeneous among genes. A first group of 29 genes had stronger LD (mean r(2)=0.28), and a second group of 38 genes had weaker LD (mean r(2)=0.12). While a strong relationship was found with the recombination rate, there was no obvious relationship between LD and functional classification. The level of nucleotide diversity, which was highly heterogeneous across genes, was also not significantly correlated with LD. A search for selection signatures highlighted significant deviations from the standard neutral model, which could be mostly attributed to recent demographic changes. Little evidence was seen for hitchhiking and clear relationships with LD. When compared among conifer species, on average, levels of LD were similar in genes from white spruce, Norway spruce and Scots pine, whereas loblolly pine and Douglas fir genes exhibited a significantly higher LD.     

Validation of whole genome linkage‐linkage disequilibrium and association results,and identification of markers to predict genetic merit for twinning

A previous genome‐wide search with a moderate density 10K marker set identified many marker associations with twinning rate, either through single‐marker analysis or combined linkage‐linkage disequilibrium (LLD; haplotype) analysis. The objective of the current study was to validate putative marker associations using an independent set of phenotypic data. Holstein bulls (n = 921) from 100 paternal half‐sib families were genotyped. Twinning rate predicted transmitting abilities were calculated using calving records from 1994 to 1998 (Data I) and 1999 to 2006 (Data II), and the underlying liability scores from threshold model analysis were used as the trait in marker association analyses. The previous analysis used 201 bulls with daughter records in Data I. In the current analysis, this was increased to 434, providing a revised estimate of effect and significance. Bulls with daughter records in Data II totaled 851, and analysis of this data provided the validation of results from analysis of Data I. Single nucleotide polymorphisms (SNPs) were selected to validate previously significant single‐marker associations and LLD results. Bulls were genotyped for a total of 306 markers. Nine of 13 LLD regions located on chromosomes 1, 2, 3, 6, 9, 22, 23(2) and 26 were validated, showing significant results for both Data I and II. Association analysis revealed 55 of 174 markers validated, equating to a single‐marker validation rate of 31%. Stepwise backward elimination and cross‐validation analyses identified 18 SNPs for use in a final reduced marker panel explaining 34% of the genetic variation, and to allow prediction of genetic merit for twinning rate.     

Evaluation of Bovine High-Density SNP Genotyping Array in Indigenous Dairy Cattle Breeds

In total 52 samples of Sahiwal (19 Excoffier L, Laval G, Schneider S. Arlequin 3.01: An integrated software package for population genetics data analysis. Evol Bioinform Online 2005; 1:47–50.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), Tharparkar (17 Hayes BJ, Bowman PJ, Chamberlain, AJ, Goddard ME. Invited review: genomic selection in dairy cattle: progress and challenges. J Dairy Sci 2009; 92:433–443.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), and Gir (16 Melka HD, Jeon EK, Kim SW, Han JB, Yoon D, Kim KS. Identification of genomic differences between Hanwoo and Holstein breeds using the Illumina Bovine SNP50 BeadChip. Genomics Inform 2011; 9:69–73.[Crossref] , [Google Scholar]) were genotyped by using BovineHD SNP chip to analyze minor allele frequency (MAF), genetic diversity, and linkage disequilibrium among these cattle. The common SNPs of BovineHD and 54K SNP Chips were also extracted and evaluated for their performance. Only 40%?50% SNPs of these arrays was found informative for genetic analysis in these cattle breeds. The overall mean of MAF for SNPs of BovineHD SNPChip was 0.248?±?0.006, 0.241?±?0.007, and 0.242?±?0.009 in Sahiwal, Tharparkar and Gir, respectively, while that for 54K SNPs was on lower side. The average Reynold’s genetic distance between breeds ranged from 0.042 to 0.055 based on BovineHD Beadchip, and from 0.052 to 0.084 based on 54K SNP Chip. The estimates of genetic diversity based on HD and 54K chips were almost same and, hence, low density chip seems to be good enough to decipher genetic diversity of these cattle breeds. The linkage disequilibrium started decaying (r²?相似文献   

SHEsis, a powerful software platform for analyses of linkage disequilibrium, haplotype construction, and genetic association at polymorphism loci

In multiloci-based genetic association studies of complex diseases, a powerful and high efficient tool for analyses of linkage disequilibrium (LD) between markers, haplotype distributions and many chi-square/p values with a large number of samples has been sought for long. In order to achieve the goal of obtaining meaningful results directly from raw data,we developed a robust and user-friendly software platform with a series of tools for analysis in association study with high efficiency. The platform has been well evaluated by several sets of real data.     

The signature of positive selection on standing genetic variation

Considerable interest is focused on the use of polymorphism data to identify regions of the genome that underlie recent adaptations. These searches are guided by a simple model of positive selection, in which a mutation is favored as soon as it arises. This assumption may not be realistic, as environmental changes and range expansions may lead previously neutral or deleterious alleles to become beneficial. We examine what effect this mode of selection has on patterns of variation at linked neutral sites by implementing a new coalescent model of positive directional selection on standing variation. In this model, a neutral allele arises and drifts in the population, then at frequency f becomes beneficial, and eventually reaches fixation. Depending on the value of f, this scenario can lead to a large variance in allele frequency spectra and in levels of linkage disequilibrium at linked, neutral sites. In particular, for intermediate f, the beneficial substitution often leads to a loss of rare alleles--a pattern that differs markedly from the signature of directional selection currently relied on by researchers. These findings highlight the importance of an accurate characterization of the effects of positive selection, if we are to reliably identify recent adaptations from polymorphism data.     

Dissecting linkage disequilibrium in African-American genomes: roles of markers and individuals

Substantial increases of linkage disequilibrium (LD) both in magnitude and in range have been observed in recently admixed populations such as African-American (AfA). On the other hand, it has also been shown that LD in AfAs was very similar to that of African. In this study, we attempted to resolve these contradicting observations by conducting a systematic examination of the LD structure in AfAs by genotyping a sample of AfA individuals at 24,341 single nucleotide polymorphisms (SNPs) spanning almost the entire chromosome 21, with an average density of 1.5 kb/SNP. The overall LD in AfAs is similar to that in African populations and much less than that in European populations. Even when the ancestry-informative markers (AIMs) were used, extended LD in AfA was found to be limited to certain magnitude range (0.2 < or = r(2) < or = 0.8) and certain distance range, that is, between-marker distance more than 200 kb. Furthermore, the inclusion of AfA individuals with predominant African ancestry was found to reduce the overall magnitude of LD. Elevation of LD in the AfA population, compared with its parental populations, can only be observed at the markers with large allele frequency differences between 2 parental populations at limited scenario. AfA individuals of wholly African ancestry contribute little to the extended LD in the AfA population, and further genotyping or association analysis conducted using only admixed individuals may lead to higher statistical power and possibly reduced cost.     

Patterns of cyto-nuclear linkage disequilibrium in Silene latifolia: genomic heterogeneity and temporal stability

Non-random association of alleles in the nucleus and cytoplasmic organelles, or cyto-nuclear linkage disequilibrium (LD), is both an important component of a number of evolutionary processes and a statistical indicator of others. The evolutionary significance of cyto-nuclear LD will depend on both its magnitude and how stable those associations are through time. Here, we use a longitudinal population genetic data set to explore the magnitude and temporal dynamics of cyto-nuclear disequilibria through time. We genotyped 135 and 170 individuals from 16 and 17 patches of the plant species Silene latifolia in Southwestern VA, sampled in 1993 and 2008, respectively. Individuals were genotyped at 14 highly polymorphic microsatellite markers and a single-nucleotide polymorphism (SNP) in the mitochondrial gene, atp1. Normalized LD (D′) between nuclear and cytoplasmic loci varied considerably depending on which nuclear locus was considered (ranging from 0.005–0.632). Four of the 14 cyto-nuclear associations showed a statistically significant shift over approximately seven generations. However, the overall magnitude of this disequilibrium was largely stable over time. The observed origin and stability of cyto-nuclear LD is most likely caused by the slow admixture between anciently diverged lineages within the species'' newly invaded range, and the local spatial structure and metapopulation dynamics that are known to structure genetic variation in this system.