Divergence between Human Populations Estimated from Linkage Disequilibrium

Observed linkage disequilibrium (LD) between genetic markers in different populations descended independently from a common ancestral population can be used to estimate their absolute time of divergence, because the correlation of LD between populations will be reduced each generation by an amount that, approximately, depends only on the recombination rate between markers. Although drift leads to divergence in allele frequencies, it has less effect on divergence in LD values. We derived the relationship between LD and time of divergence and verified it with coalescent simulations. We then used HapMap Phase II data to estimate time of divergence between human populations. Summed over large numbers of pairs of loci, we find a positive correlation of LD between African and non-African populations at levels of up to ~0.3 cM. We estimate that the observed correlation of LD is consistent with an effective separation time of approximately 1,000 generations or ~25,000 years before present. The most likely explanation for such relatively low separation times is the existence of substantial levels of migration between populations after the initial separation. Theory and results from coalescent simulations confirm that low levels of migration can lead to a downward bias in the estimate of separation time.     

Inferring Admixture Histories of Human Populations Using Linkage Disequilibrium

Long-range migrations and the resulting admixtures between populations have been important forces shaping human genetic diversity. Most existing methods for detecting and reconstructing historical admixture events are based on allele frequency divergences or patterns of ancestry segments in chromosomes of admixed individuals. An emerging new approach harnesses the exponential decay of admixture-induced linkage disequilibrium (LD) as a function of genetic distance. Here, we comprehensively develop LD-based inference into a versatile tool for investigating admixture. We present a new weighted LD statistic that can be used to infer mixture proportions as well as dates with fewer constraints on reference populations than previous methods. We define an LD-based three-population test for admixture and identify scenarios in which it can detect admixture events that previous formal tests cannot. We further show that we can uncover phylogenetic relationships among populations by comparing weighted LD curves obtained using a suite of references. Finally, we describe several improvements to the computation and fitting of weighted LD curves that greatly increase the robustness and speed of the calculations. We implement all of these advances in a software package, ALDER, which we validate in simulations and apply to test for admixture among all populations from the Human Genome Diversity Project (HGDP), highlighting insights into the admixture history of Central African Pygmies, Sardinians, and Japanese.     

Linkage Disequilibrium in Subdivided Populations

The linkage disequilibrium in a subdivided populaton is shown to be equal to the sum of the average linkage disequilibrium for all subpopulations and the covariance between gene frequencies of the loci concerned. Thus, in a subdivided population the linkage disequilibrium may not be 0 even if the linkage disequilibrium in each subpopulation is 0. If a population is divided into two subpopulations between which migration occurs, the asymptotic rate of approach to linkage equilibrium is equal to either r or 2(m(1) + m(2)) - (m(1) + m(2))(2), whichever is smaller, where r is the recombination value and m(1) and m(2) are the proportions of immigrants in subpopulations 1 and 2, respectively. Thus, if migration rate is high compared with recombination value, the change of linkage disequilibrium in subdivided populations is similar to that of a single random mating population. On the other hand, if migration rate is low, the approach to lnkage equilibrium may be retarded in subdivided populations. If isolated populations begin to exchange genes by migration, linkage disequilibrium may increase temporarily even for neutral loci. If overdominant selection operates and the equilibrium gene frequencies are different in the two subpopulations, a permanent linkage disequilibrium may be produced without epistasis in each subpopulation.     

Linkage Disequilibrium in Growing and Stable Populations

Nonrandom associations between alleles at different loci can be tested for using Fisher's exact test. Extensive simulations show that there is a substantial probability of obtaining significant nonrandom associations between closely or completely linked polymorphic neutral loci in a population of constant size at equilibrium under mutation and genetic drift. In a rapidly growing population, however, there will be little chance of finding significant nonrandom associations even between completely linked loci if the growth has been sufficiently rapid. This result is illustrated by the analysis of mitochondrial DNA sequence data from humans. In comparing all pairs of informative sites, fewer than 5% of the pairs show significant disequilibrium in Sardinians, which have apparently undergone rapid population growth, while 20% to 30% in !Kung and Pygmies, which apparently have not undergone rapid growth, show significance. The extent of linkage disequilibrium in a population is closely related to the gene genealogies of the loci examined, with ``starlike' genealogies making significant linkage disequilibrium unlikely.     

Linkage Disequilibrium in Natural Populations of DROSOPHILA MELANOGASTER

Two large, stable populations (Texas and Japan) of Drosophila melanogaster were surveyed at 21 allozyme loci on the second and third chromosomes and for chromosomal gene arrangements on those two chromosomes. Over 220 independent gametes were sampled from each population. The types and frequencies of the surveyed genetic variation are similar to those observed previously and suggest only slight differentiation among geographically distant populations. Linkage disequilibrium among linked allozymes loci is only slightly, if at all, detectable with these sample sizes. Linkage disequilibrium between linked inversions and allozymes loci is common especially when located in the same arm. These disequilibria appear to be in the same direction for most comparisons in the two population samples. This result is interpreted as evidence of similar selective environments (ecological and genetic) in the two populations. It is also noted that the direction of these linkage disequilibria appears to be oriented with respect to the gene frequencies at the component loci.     

Linkage Disequilibrium in Natural and Experimental Populations of Drosophila Melanogaster

We have studied linkage disequilibrium in Drosophila melanogaster in two samples from a wild population and in four large laboratory populations derived from the wild samples. We have assayed four polymorphic enzyme loci, fairly closely linked in the third chromosome: Sod Est-6, Pgm, and Odh. The assay method used allows us to identify the allele associations separately in each of the two homologous chromosomes from each male sampled. We have detected significant linkage disequilibrium between two loci in 16.7% of the cases in the wild samples and in 27.8% of the cases in the experimental populations, considerably more than would be expected by chance alone. We have also found three-locus disequilibria in more instances than would be expected by chance. Some disequilibria present in the wild samples disappear in the experimental populations derived from them, but new ones appear over the generations. The effective population sizes required to generate the observed disequilibria by randomness range from 40 to more than 60,000 individuals in the natural population, depending on which locus pair is considered, and from 100 to more than 60,000 in the experimental populations. These population sizes are unrealistic; the fact that different locus-pairs yield disparate estimates within the same population argues against the likelihood that the disequilibria may have arisen as a consequence of population bottlenecks. Migration, or population mixing, cannot be excluded as the process generating the disequilibria in the wild samples, but can in the experimental populations. We conclude that linkage disequilibrium in these populations is most likely due to natural selection acting on the allozymes, or on loci very tightly linked to them.     

A Study of Linkage Disequilibrium in Populations of DROSOPHILA MELANOGASTER

This paper presents the results of a study of linkage disequilibrium between five polymorphic enzyme genes located on chromosome 3 of D. melanogaster. Three sets of chromosomes were examined: two represented samples from successive years of the same natural population, and one came from a large laboratory population. Out of the thirty possible tests for linkage disequilibrium between pairs of loci, two were significant at the 5% level and two at the 1% level. This result cannot reasonably be ascribed to chance alone. The pairs of loci that had a significant correlation in one sample had higher than average correlations in the other samples (though not necessarily in the same direction); this effect was highly significant statistically. There was no tendency for the high correlations to be associated with tightness of linkage between the loci concerned. All five loci were involved in at least one significant effect. It was concluded that these results are difficult to explain on the neutral allele theory of protein polymorphism, but are consistent with the concept of selective control of allele frequencies.     

Linkage Disequilibrium between Allozymes in Natural Populations of Lodgepole Pine

Pairwise linkage disequilibrium values (D) were estimated for 14 allozyme loci in two natural populations of lodgepole pine (Pinus contorta ssp. latifolia). Maternal multilocus genotypes were inferred from samples of (haploid) megagametophytic seed-endosperms. Coupling/repulsion double heterozygotes were distinguished for closely linked pairs of loci. Assays of seven of the loci in seed embryos allowed estimates of D for these loci in the outcross pollen pool (estimates of outcrossing rates indicate no significant departures from random mating in either population). No disequilibrium was observed between unlinked loci in either maternal genotypes or outcross pollen. However, significant disequilibrium was observed within and between gametes for some allelic combinations of four tightly linked loci; the assumption of random association of gamete types within individuals is thus invalid for some loci in lodgepole pine. Possible causes of the observed D were examined using the noncentrality parameter of the general noncentral chi square distribution. We concluded, from estimates of population size, linkage and measurements of population substructure, that neither drift nor population subdivision was responsible for the significant values of D which were observed and that epistatic selection was the most likely cause of the disequilibrium observed.     

A Study of Linkage Disequilibrium in British Populations of DROSOPHILA SUBOBSCURA

Data have been obtained concerning the genetic content of samples of O chromosomes from three British populations, and J chromosomes from one population, of Drosophila subobscura. Some improvements to the genetic map of the O and J chromosomes have been made. Allele frequencies at the loci studied do not show much geographical variation, except where associations with geographically varying gene arrangements distort the picture. Striking nonrandom associations between alleles at three enzyme loci and closely linked O chromosome gene arrangements are present. Some historical explanation for these associations cannot at present be ruled out, but it is clear that a very high degree of genetic differentiation must exist between different gene arrangements in this species. There is no convincing direct evidence for linkage disequilibrium between pairs of enzyme loci, although there is a significant association between close linkage and a high value of the linkage disequilibrium measure. This suggests that there may be disequilibria between closely linked enzyme loci that are too small to be individually detectable. These results are in broad agreement with those reported by workers on other Drosophila species. At present there appears to be no evidence to support the concept that selection is sufficiently strong at individual enzyme loci to produce a high degree of nonrandom associations. (Franklin and Lewontin 1970).     

Linkage Disequilibrium in Natural Populations of DROSOPHILA MELANOGASTER. Seasonal Variation

Linkage disequilibrium among ten polymorphic allozyme loci and polymorphic inversions on chromosomes 2 and 3 in a natural population of Drosophila melanogaster was examined early and late in the annual season. Similar to previous studies, little linkage disequilibrium was observed among allozymes. The two significant cases that were observed in the first sample behaved in a contradictory way. One declined much more rapidly than expected due simply to recombination; the other declined slowly as expected. There was little change in allozyme or inversion frequencies during the season.     

Extensive Long-Range and Nonsyntenic Linkage Disequilibrium in Livestock Populations: Deconstruction of a Conundrum

Great interest was aroused by reports, based on microsatellite markers, of high levels of statistically significant long-range and nonsyntenic linkage disequilibrium (LD) in livestock. Simulation studies showed that this could result from population family structure. In contrast, recent SNP-based studies of livestock populations report much lower levels of LD. In this study we show, on the basis of microsatellite data from four cattle populations, that high levels of long-range LD are indeed obtained when using the multi-allelic D′ measure of LD. Long-range and nonsyntenic LD are exceedingly low, however, when evaluated by the standardized chi-square measure of LD, which stands in relation to the predictive ability of LD. Furthermore, specially constructed study populations provided no evidence for appreciable LD resulting from family structure at the grandparent level. We propose that the high statistical significance and family structure effects observed in the earlier studies are due to the use of large sample sizes, which accord high statistical significance to even slight deviations from asymptotic expectations under the null hypothesis. Nevertheless, even after taking sample size into account, our results indicate that microsatellites testify to the presence of usable LD at considerably wider separation distances than SNPs, suggesting that use of SNP haplotypes may considerably increase the usefulness of a given fixed SNP array.     

Allozymic Variation and Linkage Disequilibrium in Some Laboratory Populations of DROSOPHILA MELANOGASTER

Nine laboratory populations of D. melanogaster were surveyed by starch gel electrophoresis for variation at 17 enzyme loci. A single-fly extract could be assayed for all 17 enzymes, so that the data consist of 17-locus genotypes.--Pairwise linkage disequilibria were estimated from the multilocus genotypic frequencies, using both Burrows' and Hill's methods. Large amounts of linkage disequilibrium were found, in contrast to the results reported for natural populations.-Knowledge of the approximate sizes of these populations was used to compare the observed heterozygosities and linkage disequilibria with predictions of the neutral allele hypothesis. The relatively large amount of linkage disequilibrium is consistent with the small sizes of the populations. However, the levels of heterozygosity in at least some populations suggest that some mechanism has been operating to retard the rate of decay by random drift. Several examples of significant deviation from Hardy-Weinberg frequencies and the large amount of linkage disequilibrum present in these populations indicate that a likely mechanism is selective effects associated with neutral alleles because of linkage disequilibrium with selected loci (e.g., "associative overdominance"). The results are therefore consistent with both neutralist, and selectionist hypotheses, but suggest the importance of considering linkage disequilibrium between neutral and selected loci when attempting to explain the dynamics of enzyme polymorphisms.     

Long Range Linkage Disequilibrium across the Human Genome

Long-range linkage disequilibria (LRLD) between sites that are widely separated on chromosomes may suggest that population admixture, epistatic selection, or other evolutionary forces are at work. We quantified patterns of LRLD on a chromosome-wide level in the YRI population of the HapMap dataset of single nucleotide polymorphisms (SNPs). We calculated the disequilibrium between all pairs of SNPs on each chromosome (a total of >2×10¹¹ values) and evaluated significance of overall disequilibrium using randomization. The results show an excess of associations between pairs of distant sites (separated by >0.25 cM) on all of the 22 autosomes. We discuss possible explanations for this observation.     

A Model of Functional Epistasis and Linkage Disequilibrium in Populations with Overlapping Generations

A model of functional epistasis is proposed in which it is assumed that coupling and repulsion genotypes differ in metabolic efficiency and thus in development time and net fecundity. The implications of this model are investigated for iteroparous populations with fluctuating rates of increase. It is found that the fluctuations in rate of increase can lead to large fluctuations in gamete frequency and D, the coefficient of linkage disequilibrium, but that D will almost always have a value of zero at some point during the populations' demographic cycle. Some of the model populations would be expected to be in a state of linkage disequilibrium only fleetingly: others would exhibit D-cycles interpretable as random fluctuation. Implications of the model for interpretations of existing data on linkage disequilibrium among enzyme loci in Drosophila are discussed.     

The Genetics of DROSOPHILA SUBOBSCURA Populations. IX. Studies on Linkage Disequilibrium in Four Natural Populations

Gametic frequencies were obtained in four natural populations of D. sub-obscura by extracting wild chromosomes and subsequently analyzing them for inversions and allozymes. The genes Lap and Pept-1, both located within the same inversions of chromosome O, were found in striking nonrandom associations with them of the same kind and degree in all populations studied. On the contrary, the gene Acph, also located within the previously mentioned inversions, was found in linkage disequilibrium with them only in two populations and of opposite directions. This is also the case for the genes Est-9 and Hk, both located within chromosome E inversions. While the gene Est-9 was in strong linkage disequilibrium with the inversions, of the same kind and degree in all populations studied, Hk was found to be in linkage equilibrium. Allele frequencies for the 29 genes studied do not show geographical variation except for the genes Lap, Pept-1 and Est-9, the ones found in linkage disequilibria with the geographically varying gene arrangements. Although mechanical or historical explanations for these equilibria cannot be ruled out, these data cannot be explained satisfactorily by the "middle gene explanation," which states that loci displaying such linkage disequilibria are the ones located near the break points of inversions, while the ones displaying linkage equilibria with them are located in the middle of them. There is no evidence for consistent linkage disequilibria between pairs of loci, except for the closely linked genes of the complex locus, Est-9. This would imply, if it is not a peculiarity of the Est-9 complex, that the linkage disequilibria are found only between very closely linked loci or that, for less closely linked genes, the associations are too weak to be detected by the usual samples sizes.     

Linkage Disequilibrium in Human Ribosomal Genes: Implications for Multigene Family Evolution

Members of the rDNA multigene family within a species do not evolve independently, rather, they evolve together in a concerted fashion. Between species, however, each multigene family does evolve independently indicating that mechanisms exist which will amplify and fix new mutations both within populations and within species. In order to evaluate the possible mechanisms by which mutation, amplification and fixation occur we have determined the level of linkage disequilibrium between two polymorphic sites in human ribosomal genes in five racial groups and among individuals within two of these groups. The marked linkage disequilibrium we observe within individuals suggests that sister chromatid exchanges are much more important than homologous or nonhomologous recombination events in the concerted evolution of the rDNA family and further that recent models of molecular drive may not apply to the evolution of the rDNA multigene family.     

Evidence for Linkage Disequilibrium Maintained by Selection in Two Natural Populations of DROSOPHILA SUBOBSCURA

One island and one mainland population of Drosophila subobscura were found polymorphic at the XDH (xanthine dehydrogenase) and the AO (aldehyde oxidase) loci. It was observed that one allele at the XDH locus, which has a low frequency in both populations, is nonrandomly associated with the alleles at the AO locus. Two lines of evidence support the thesis that this linkage disequilibrium is due to epistasis rather than random drift: (1) D or r, measures of the disequilibrium, have the same sign and magnitude in both populations. (2) The linkage disequilibrium is not due to inversions. Inversions segregating on the chromosome carrying XDH and AO have been separated into two classes, between which exchange of alleles at the two loci is suppressed. Linkage disequilibrium for XDH and AO was observed within each class. In the absence of any exchange of alleles, these disequilibria must have arisen and been maintained independently. The suggestion is made that the epistatic disequilibrium results from the close structural and physiological relationship which exists between the two enzymes.     

Linkage Disequilibrium Decay and Past Population History in the Human Genome

The fluctuation of population size has not been well studied in the previous studies of theoretical linkage disequilibrium (LD) expectation. In this study, an improved theoretical prediction of LD decay was derived to account for the effects of changes in effective population sizes. The equation was used to estimate effective population size (N_e) assuming a constant N_e and LD at equilibrium, and these N_e estimates implied the past changes of N_e for a certain number of generations until equilibrium, which differed based on recombination rate. As the influence of recent population history on the N_e estimates is larger than old population history, recent changes in population size can be inferred more accurately than old changes. The theoretical predictions based on this improved expression showed accurate agreement with the simulated values. When applied to human genome data, the detailed recent history of human populations was obtained. The inferred past population history of each population showed good correspondence with historical studies. Specifically, four populations (three African ancestries and one Mexican ancestry) showed population growth that was significantly less than that of other populations, and two populations originated from China showed prominent exponential growth. During the examination of overall LD decay in the human genome, a selection pressure on chromosome 14, the gephyrin gene, was observed in all populations.     

人类混血群体连锁不平衡的遗传学模型与基因定位

人类混血群体可以说是混合群体的一种特例．在无选择、无突变、无限随机交配群体的假定前提下,研究了亲本群体的基因频率对混血群体及其衍生后代群体连锁不平衡结构的影响,导出了各群体连锁不平衡值的表达式,建立了一个估计基因间重组率的简便方法;同时, 采用估算分子标记与QTL之间连锁不平衡系数的统计分析方法,分析了人类混血群体及其衍生后代群体QTL检测与估计的关系,建立了该关系的系列理论公式．研究结果表明,本方法不仅适用于人类疾病(包括复杂遗传疾病)基因定位,而且适合于人类正常基因的定位,同时也适用于人类普通多基因性状的QTL分析．