Evolution of modern humans: evidence from nuclear DNA polymorphisms.

Previously we have described studies of the evolution of modern humans based upon data for classical genetic markers and for nuclear DNA polymorphisms. Such polymorphisms provide a different point of view regarding human evolution than do mitochondrial DNA sequences. Here we compare revised dates for major migrations of anatomically modern humans, estimated from archaeological data, with separations suggested by a genetic tree constructed from classical marker allele frequencies. Analyses of DNA polymorphisms have now been extended and compared with those of classical markers; genetic trees continue to support the hypothesis of an initial African and non-African divergence for modern humans. We have also begun testing non-human primates for a set of human DNA polymorphisms. For most polymorphisms tested so far, humans share a single allele with other primates; such shared alleles are likely to be ancestral. Populations living in humid tropical environments have significantly higher frequencies of ancestral alleles than do other populations, supporting the hypothesis that natural selection acts to maintain high frequencies of particular alleles in some environments.     

HLA DNA sequence variation among human populations: molecular signatures of demographic and selective events

Molecular differences between HLA alleles vary up to 57 nucleotides within the peptide binding coding region of human Major Histocompatibility Complex (MHC) genes, but it is still unclear whether this variation results from a stochastic process or from selective constraints related to functional differences among HLA molecules. Although HLA alleles are generally treated as equidistant molecular units in population genetic studies, DNA sequence diversity among populations is also crucial to interpret the observed HLA polymorphism. In this study, we used a large dataset of 2,062 DNA sequences defined for the different HLA alleles to analyze nucleotide diversity of seven HLA genes in 23,500 individuals of about 200 populations spread worldwide. We first analyzed the HLA molecular structure and diversity of these populations in relation to geographic variation and we further investigated possible departures from selective neutrality through Tajima's tests and mismatch distributions. All results were compared to those obtained by classical approaches applied to HLA allele frequencies.Our study shows that the global patterns of HLA nucleotide diversity among populations are significantly correlated to geography, although in some specific cases the molecular information reveals unexpected genetic relationships. At all loci except HLA-DPB1, populations have accumulated a high proportion of very divergent alleles, suggesting an advantage of heterozygotes expressing molecularly distant HLA molecules (asymmetric overdominant selection model). However, both different intensities of selection and unequal levels of gene conversion may explain the heterogeneous mismatch distributions observed among the loci. Also, distinctive patterns of sequence divergence observed at the HLA-DPB1 locus suggest current neutrality but old selective pressures on this gene. We conclude that HLA DNA sequences advantageously complement HLA allele frequencies as a source of data used to explore the genetic history of human populations, and that their analysis allows a more thorough investigation of human MHC molecular evolution.     

The timing of human adaptation from Neanderthal introgression

Admixture has the potential to facilitate adaptation by providing alleles that are immediately adaptive in a new environment or by simply increasing the long-term reservoir of genetic diversity for future adaptation. A growing number of cases of adaptive introgression are being identified in species across the tree of life, however the timing of selection, and therefore the importance of the different evolutionary roles of admixture, is typically unknown. Here, we investigate the spatio-temporal history of selection favoring Neanderthal-introgressed alleles in modern human populations. Using both ancient and present-day samples of modern humans, we integrate the known demographic history of populations, namely population divergence and migration, with tests for selection. We model how a sweep placed along different branches of an admixture graph acts to modify the variance and covariance in neutral allele frequencies among populations at linked loci. Using a method based on this model of allele frequencies, we study previously identified cases of adaptive Neanderthal introgression. From these, we identify cases in which Neanderthal-introgressed alleles were quickly beneficial and other cases in which they persisted at low frequency for some time. For some of the alleles that persisted at low frequency, we show that selection likely independently favored them later on in geographically separated populations. Our work highlights how admixture with ancient hominins has contributed to modern human adaptation and contextualizes observed levels of Neanderthal ancestry in present-day and ancient samples.     

Frequency Distribution of Esterase-5 Alleles in Two Populations of DROSOPHILA PSEUDOOBSCURA

Statistical tests comparing allele frequencies in natural populations with those predicted by various theories of genic variation depend critically on the accurate enumeration of alleles. This study used a series of five sequential electrophoretic conditions to characterize the allele frequency distributions of esterase-5 in two large population samples of Drosophila pseudoobscura from California. In Standard chromosome lines 12 electromorphs were discriminated using a single electrophoretic condition. When four additional criteria were used, the number of electromorphs increased to 41, 33 in one population and 22 in the other. Both populations had the same two alleles in high frequency, with other alleles present in frequencies of 6% or less. Although each population had a number of unique alleles, a χ² contingency test demonstrated no significant genetic divergence between them. A statistical comparison of allele frequencies in both populations with that predicted by neutral models suggests that the individual and combined distributions deviate from neutrality in the direction of purifying selection.—Sex-Ratio chromosomes differed markedly from Standard chromosomes in both allelic content and diversity. In 32 Sex-Ratio chromosomes from one population only three alleles were found, all of which were detected under the initial "standard" electrophoretic conditions. Moreover, none of these alleles was found in the Standard chromosome lines.     

The consequences of fluctuating selection for isozyme polymorphisms in Daphnia

Temporal sequences of allele frequencies in natural populations of Daphnia are analyzed to obtain the mean and variance of the selection coefficient for both asexual and sexual phases. In general, the alleles at enzyme loci appear to be quasi-neutral. Although significant variation exists for the estimated selection coefficients, the means are in all cases close to zero. Estimates of the variance of selection intensity are applied to existing models to demonstrate the implications of fluctuating selection for the spatial and temporal distribution of gene frequencies in Daphnia. The empirical and analytical results are shown to provide a possible solution to some previously puzzling aspects of Daphnia population genetic surveys. Neither genetic drift nor diversifying selection are necessary conditions for the local diversification of gene frequencies.     

Statistical Mechanics and the Evolution of Polygenic Quantitative Traits

The evolution of quantitative characters depends on the frequencies of the alleles involved, yet these frequencies cannot usually be measured. Previous groups have proposed an approximation to the dynamics of quantitative traits, based on an analogy with statistical mechanics. We present a modified version of that approach, which makes the analogy more precise and applies quite generally to describe the evolution of allele frequencies. We calculate explicitly how the macroscopic quantities (i.e., quantities that depend on the quantitative trait) depend on evolutionary forces, in a way that is independent of the microscopic details. We first show that the stationary distribution of allele frequencies under drift, selection, and mutation maximizes a certain measure of entropy, subject to constraints on the expectation of observable quantities. We then approximate the dynamical changes in these expectations, assuming that the distribution of allele frequencies always maximizes entropy, conditional on the expected values. When applied to directional selection on an additive trait, this gives a very good approximation to the evolution of the trait mean and the genetic variance, when the number of mutations per generation is sufficiently high (4Nμ > 1). We show how the method can be modified for small mutation rates (4Nμ → 0). We outline how this method describes epistatic interactions as, for example, with stabilizing selection.     

Frequency distributions for chloroplast genes in Chlamydomonas zygote clones: Evidence for random drift

Mitochondria and chloroplasts of eucaryotic cells contain populations of DNA molecules. In certain cases, e.g., the chloroplasts of Chlamydomonas reinhardtii and the mitochondria of Saccharomyces cerevisiae, organelles contributed by the two parents are known to fuse in the zygote, creating a single population of DNA molecules. In a cross, this population will include molecules of both parental genotypes. There is reason to suspect that organelle DNA molecules in this population are selected randomly for replication and recombination. This would result in random changes in the frequency of a particular allele or genotype within the organelle gene pool of a single zygote and also within its clone of progeny cells. A given gene frequency would increase in some zygote clones and decrease in others, analogous to random drift of gene frequencies in small Mendelian populations. To test this, we have examined the distribution of chloroplast gene frequencies among the zygote clones produced in each of a number of crosses of Chlamydomonas. These distributions are typically U or L shaped as predicted by the random drift hypothesis. They include uniparental zygote clones, in which a chloroplast allele from one parent has been fixed (frequency 100%) and the alternative allele from the other parent has been lost (frequency 0%). Among the remaining (biparental) zygote clones, there is a linear distribution of allele frequencies, showing a great increase in variance over the input frequencies. In these experiments both biparental and uniparental zygotes show a bias favoring chloroplast alleles from the mt⁺ (maternal) parent, and there is no statistically significant mode at the allele frequency of 0.5 corresponding to the equal input of alleles from the maternal and paternal (mt^?) parents. The observed distributions support the hypothesis that both uniparental inheritance and the high variance of allele frequencies among zygote clones are due to random drift of allele frequencies, coupled with a directional force which favors fixation of the maternal allele. In addition, statistical analysis of the data shows a strong but incomplete tendency for linked chloroplast markers to be fixed or lost together in uniparental zygotes. Possible cellular and molecular mechanisms for these observations are discussed.     

The Maintenance of Single-Locus Polymorphism. IV. Models with Mutation from Existing Alleles

The ability of viability selection to maintain allelic polymorphism is investigated using a constructionist approach. In extensions to the models we have previously proposed, a population is bombarded with a series of mutations whose fitnesses in conjunction with other alleles are functions of the corresponding fitnesses with a particular allele, the parent allele, already in the population. Allele frequencies are iterated simultaneously, thus allowing alleles to be driven to extinction by selection. Such models allow very high levels of polymorphism to evolve: up to 38 alleles in one case. Alleles that are lethal as homozygotes can evolve to surprisingly high frequencies. The joint evolution of allele frequencies and viabilities highlights the necessity to consider more than the current morphology of a population. Comparisons are made with the neutral theory of evolution and it is suggested that failure to reject neutrality using the Ewens-Watterson test cannot be regarded as evidence for the neutral theory.     

A principal target of human immunity to malaria identified by molecular population genetic and immunological analyses

New strategies are required to identify the most important targets of protective immunity in complex eukaryotic pathogens. Natural selection maintains allelic variation in some antigens of the malaria parasite Plasmodium falciparum. Analysis of allele frequency distributions could identify the loci under most intense selection. The merozoite surface protein 1 (Msp1) is the most-abundant surface component on the erythrocyte-invading stage of P. falciparum. Immunization with whole Msp1 has protected monkeys completely against homologous and partially against non-homologous parasite strains. The single-copy msp1 gene, of about 5 kilobases, has highly divergent alleles with stable frequencies in endemic populations. To identify the region of msp1 under strongest selection to maintain alleles within populations, we studied multiple intragenic sequence loci in populations in different regions of Africa and Southeast Asia. On both continents, the locus with the lowest inter-population variance in allele frequencies was block 2, indicating selection in this part of the gene. To test the hypothesis of immune selection, we undertook a large prospective longitudinal cohort study. This demonstrated that serum IgG antibodies against each of the two most frequent allelic types of block 2 of the protein were strongly associated with protection from P. falciparum malaria.     

New Statistical Tests of Neutrality for DNA Samples from a Population

The purpose of this paper is to develop statistical tests of the neutral model of evolution against a class of alternative models with the common characteristic of having an excess of mutations that occurred a long time ago or a reduction of recent mutations compared to the neutral model. This class of population genetics models include models for structured populations, models with decreasing effective population size and models of selection and mutation balance. Four statistical tests were proposed in this paper for DNA samples from a population. Two of these tests, one new and another a modification of an existing test, are based on EWENS'' sampling formula, and the other two new tests make use of the frequencies of mutations of various classes. Using simulated samples and regression analyses, the critical values of these tests can be computed from regression equations. This approach for computing the critical values of a test was found to be appropriate and quite effective. We examined the powers of these four tests using simulated samples from structured populations, populations with linearly decreasing sizes and models of selection and mutation balance and found that they are more powerful than existing statistical tests of the neutral model of evolution.     

Perturbation-Reperturbation Test of Selection Vs. Hitchhiking of the Two Major Alleles of Esterase-5 in Drosophila Pseudoobscura

A perturbation-reperturbation tests selective neutrality of 100/100/100/100/100 and 106/100/100/100/100, the two most common alleles at the highly polymorphic X-linked locus Esterase-5 in Drosophila pseudoobscura. A total of 22 replicate populations are set up in cages, 11 start at a high frequency of 76% (U) and 11 at a low frequency of 21% (N) of the 106 allele. Allele frequencies change directionally and decrease in both U and N populations as groups and reach equilibria of 60 and 14%, respectively, after 200-300 days. These changes suggest natural selection. A hypothesis of balancing selection accounts for the pattern and predicts a dynamic equilibrium. A rival neutral hypothesis accounts for the pattern equally well by postulating hitchhiking and breakup of linkage leaving the Est-5 variants to drift at neutral equilibria. A reperturbation of allele frequencies in each population, creating 22 additional reperturbed populations EN and EU, with the original populations as controls, directly addresses the question of balancing selection or hitchhiking and breakup of linkage effects. Allele frequencies do not change directionally among the reperturbed populations as a group. The hypothesis of balancing selection is rejected in favor of the hypothesis of initial hitchhiking and dissipated linkage effects. The power of the experimental design to detect selection is studied by simulation. Within the limits of power set by the design, it is concluded that the 100 and 106 are iso-fitness alleles of Est-5 under the environmental conditions of the laboratory populations. The requirements of a method of perturbation and reperturbation are discussed.     

Parallel evolution of adaptive mutations in Plasmodium falciparum mitochondrial DNA during atovaquone-proguanil treatment

Here we provide direct evidence that two adaptive nucleotide changes in the same codon (268) of the cytochrome b gene (pfcytb) each occurred repeatedly in independent Plasmodium falciparum lineages exposed to the antimalarial drug atovaquone-proguanil (AP). We analyzed the history of 7 AP resistance alleles from clinical isolates by sequencing the mitochondrial (mt) genome that encodes the pfcytb gene and found that a distinct mt haplotype was associated with each AP resistance allele. By comparing mt sequences and microsatellite genotypes of the isolates both before treatment initiation and at the day of failure for each uncured patient, we observed that the AP resistance alleles occurred and spread within the patients. These data demonstrate that identical AP resistance alleles have multiple independent origins and provide an example of parallel evolution driven by drug treatment selection in P. falciparum.     

Usefulness of molecular markers for detecting population bottlenecks via monitoring genetic change

It is important to detect population bottlenecks in threatened and managed species because bottlenecks can increase the risk of population extinction. Early detection is critical and can be facilitated by statistically powerful monitoring programs for detecting bottleneck-induced genetic change. We used Monte Carlo computer simulations to evaluate the power of the following tests for detecting genetic changes caused by a severe reduction in a population's effective size ( N _e): a test for loss of heterozygosity, two tests for loss of alleles, two tests for change in the distribution of allele frequencies, and a test for small N _e based on variance in allele frequencies (the 'variance test'). The variance test was most powerful; it provided an 85% probability of detecting a bottleneck of size N _e = 10 when monitoring five microsatellite loci and sampling 30 individuals both before and one generation after the bottleneck. The variance test was almost 10-times more powerful than a commonly used test for loss of heterozygosity, and it allowed for detection of bottlenecks before 5% of a population's heterozygosity had been lost. The second most powerful tests were generally the tests for loss of alleles. However, these tests had reduced power for detecting genetic bottlenecks caused by skewed sex ratios. We provide guidelines for the number of loci and individuals needed to achieve high-power tests when monitoring via the variance test. We also illustrate how the variance test performs when monitoring loci that have widely different allele frequency distributions as observed in five wild populations of mountain sheep ( Ovis canadensis ).     

Signatures of environmental genetic adaptation pinpoint pathogens as the main selective pressure through human evolution

Previous genome-wide scans of positive natural selection in humans have identified a number of non-neutrally evolving genes that play important roles in skin pigmentation, metabolism, or immune function. Recent studies have also shown that a genome-wide pattern of local adaptation can be detected by identifying correlations between patterns of allele frequencies and environmental variables. Despite these observations, the degree to which natural selection is primarily driven by adaptation to local environments, and the role of pathogens or other ecological factors as selective agents, is still under debate. To address this issue, we correlated the spatial allele frequency distribution of a large sample of SNPs from 55 distinct human populations to a set of environmental factors that describe local geographical features such as climate, diet regimes, and pathogen loads. In concordance with previous studies, we detected a significant enrichment of genic SNPs, and particularly non-synonymous SNPs associated with local adaptation. Furthermore, we show that the diversity of the local pathogenic environment is the predominant driver of local adaptation, and that climate, at least as measured here, only plays a relatively minor role. While background demography by far makes the strongest contribution in explaining the genetic variance among populations, we detected about 100 genes which show an unexpectedly strong correlation between allele frequencies and pathogenic environment, after correcting for demography. Conversely, for diet regimes and climatic conditions, no genes show a similar correlation between the environmental factor and allele frequencies. This result is validated using low-coverage sequencing data for multiple populations. Among the loci targeted by pathogen-driven selection, we found an enrichment of genes associated to autoimmune diseases, such as celiac disease, type 1 diabetes, and multiples sclerosis, which lends credence to the hypothesis that some susceptibility alleles for autoimmune diseases may be maintained in human population due to past selective processes.     

CCR2-64Ⅰ在中国南方14个少数民族群体中的分布

为调查HIV-1感染相关等位基因CCR2-64Ⅰ在我国南方14个少数民族群体的频率和多态性分布, 从上述人群外周血中抽提基因组DNA, 采用PCR和PCR-RFLP等方法进行基因分型。在791例调查对象中, 636例是野生纯合子基因型, 104例为杂合子基因型, 51例为突变纯合子基因型。上述各群体等位基因型的分布符合Hardy-Weinberg平衡。14个民族群体的平均突变基因频率为13.6％, 等位基因频率范围分布在1.6％~30.3％之间, 14个民族群体之间突变基因频率具有显著差异(P<0.05)。广西壮族群体CCR2-64Ⅰ突变基因频率最低, 为1.6％, 云南的六库傈僳族频率最高, 为30.3％。12个群体的突变基因频率均低于中国汉族健康群体, 南方3个少数民族群体基因突变频率显著低于西南11个少数民族群体, 该突变基因在艾滋病发病过程中的影响值得进一步深入研究。     

Nearly identical allelic distributions of xanthine dehydrogenase in two populations of Drosophila pseudoobscura

In a previous study, Keith (1983) showed by sequential gel electrophoresis of the esterase-5 protein in Drosophila pseudoobscura that a highly polymorphic locus with many alleles can have very similar frequency distributions in populations separated by 500 km. The present work studies another highly polymorphic locus, xanthine dehydrogenase, in the same California population samples, using the same technique to distinguish allelic classes. Twelve electromorphs were found in one population and 15 in the other. Both populations shared a single very frequent (approximately 60%) allele, as well as five other alleles in low but similar frequencies. In addition, each population had an array of unique alleles present only once in one population sample but absent in the other. A statistical test against the stationary distribution for neutral alleles shows that, if the populations are at equilibrium, then purifying selection is operating on xanthine dehydrogenase. The extremely close similarity in frequency distributions of the alleles between populations for both the xanthine dehydrogenase and esterase-5 loci, despite differences in allele frequency distribution between loci, strongly emphasizes the importance of migration in influencing genic diversity in these populations.     

A method for estimating the intensity of overdominant selection from the distribution of allele frequencies

A method is proposed for estimating the intensity of overdominant selection scaled by the effective population size, S = 2Ns, from allele frequencies. The method is based on the assumption that, with strong overdominant selection, allele frequencies are nearly at their deterministic equilibrium values and that, to a first approximation, deviations depend only on S. Simulations verify that reasonably accurate estimates of S can be obtained for realistic sample sizes. The method is applied to data from several loci in the major histocompatibility complex (Mhc) in numerous human populations. For alleles distinguished by both serological typing and the sequence of the peptide-binding region, our estimates of S are comparable to those obtained by analysis of DNA sequences in showing that selection is strongest on HLA-B and weaker on HLA-A, HLA-DRB1, and HLA-DQA1. The intensity of selection on HLA-B varied considerably among populations. Two populations, Native American and Inuit, showed an excess rather than a deficiency in homozygosity. Comparable estimates of S were obtained for alleles at Mhc class II loci distinguished by serological reactions (serotyping) and by differences in the amino acid sequences of the peptide-binding region (molecular typing). A comparison of two types of data for DQA1 and DRB1 showed that serotyping led to generally lower estimates of S.     

What paths do advantageous alleles take during short‐term evolutionary change?

Combining experimental evolution with whole‐genome resequencing is a promising new strategy for investigating the dynamics of evolutionary change. Published studies that have resequenced laboratory‐selected populations of sexual organisms have typically focused on populations sampled at the end of an evolution experiment. These studies have attempted to associate particular alleles with phenotypic change and attempted to distinguish between different theoretical models of adaptation. However, neither the population used to initiate the experiment nor multiple time points sampled during the evolutionary trajectory are generally available for examination. In this issue of Molecular Ecology, Orozco‐terWengel et al. (2012) take a significant step forward by estimating genome‐wide allele frequencies at the start, 15 generations into and at the end of a 37‐generation Drosophila experimental evolution study. The authors identify regions of the genome that have responded to laboratory selection and describe the temporal dynamics of allele frequency change. They identify two common trajectories for putatively adaptive alleles: alleles either gradually increase in frequency throughout the entire 37 generations or alleles plateau at a new frequency by generation 15. The identification of complex trajectories of alleles under selection contributes to a growing body of literature suggesting that simple models of adaptation, whereby beneficial alleles arise and increase in frequency unimpeded until they become fixed, may not adequately describe short‐term response to selection.     

Experimental evolution and the dynamics of genomic mutation rate modifiers

Because genes that affect mutation rates are themselves subject to mutation, mutation rates can be influenced by natural selection and other evolutionary forces. The population genetics of mutation rate modifier alleles has been a subject of theoretical interest for many decades. Here, we review experimental contributions to our understanding of mutation rate modifier dynamics. Numerous evolution experiments have shown that mutator alleles (modifiers that elevate the genomic mutation rate) can readily rise to high frequencies via genetic hitchhiking in non-recombining microbial populations. Whereas these results certainly provide an explanatory framework for observations of sporadically high mutation rates in pathogenic microbes and in cancer lineages, it is nonetheless true that most natural populations have very low mutation rates. This raises the interesting question of how mutator hitchhiking is suppressed or its phenotypic effect reversed in natural populations. Very little experimental work has addressed this question; with this in mind, we identify some promising areas for future experimental investigation.     

The relative importance of meiotic gene conversion,selection and mutation pressure,in population genetics and evolution

Disparity in the direction of meiotic gene conversion can change allele frequencies, favouring one allele of a pair in heterozygotes. Equilibrium allele frequencies for large diploid populations are examined by means of equations relating them to meiotic gene conversion, selection and mutation for deleterious recessives, deleterious dominants, and deleterious alleles with no dominance. Using observed conversion parameters from various fungi,Zea mays andDrosophila, it is shown that conversion is generally much more important than mutation pressure and may be of greater or lesser importance than selection, depending on dominance and the strength of selection and conversion forces for the alleles involved.