Mitochondrial and nuclear genes present conflicting portraits of human origins

Genetic variation in the Cytb gene of human cerebral Taenia solium cysticerci recovered from clinically and radiologically heterogeneous patients with neurocysticercosis

Neurocysticercosis (NC) is a clinically and radiologically heterogeneous parasitic disease caused by the establishment of larval Taenia solium in the human central nervous system. Host and/or parasite variations may be related to this observed heterogeneity. Genetic differences between pig and human-derived T. solium cysticerci have been reported previously. In this study, 28 cysticerci were surgically removed from 12 human NC patients, the mitochondrial gene that encodes cytochrome b was amplified from the cysticerci and genetic variations that may be related to NC heterogeneity were characterised. Nine different haplotypes (Ht), which were clustered in four haplogroups (Hg), were identified. Hg 3 and 4 exhibited a tendency to associate with age and gender, respectively. However, no significant associations were found between NC heterogeneity and the different T. solium cysticerci Ht or Hg. Parasite variants obtained from patients with similar NC clinical or radiological features were genetically closer than those found in groups of patients with a different NC profile when using the Mantel test. Overall, this study establishes the presence of genetic differences in the Cytb gene of T. solium isolated from human cysticerci and suggests that parasite variation could contribute to NC heterogeneity.

mtDNA sequence, phylogeny and evolution of laboratory mice

Laboratory mice are important tools for biomedical research. Aiming to investigate the phylogeny and evolution of laboratory mice, we investigated the mtDNA sequences of classic inbred strains, classic outbred stocks and wild-derived inbred strains. Our results showed that the most classic outbred stocks and classic inbred strains are descended from a single mtDNA ancestor. The phylogenic analysis supports the topology of M. m. castaneus/M. m. domesticus as sister subspecies, and the divergence time between the two sister subspecies and M. m. musculus was 493,000 (435,000-557,000) years ago. Furthermore, the mtDNA polymorphisms accumulated in the last 100years in the laboratory mice are under a relaxed purifying selection.

Autosomal resequence data reveal Late Stone Age signals of population expansion in sub-Saharan African foraging and farming populations

Background

A major unanswered question in the evolution of Homo sapiens is when anatomically modern human populations began to expand: was demographic growth associated with the invention of particular technologies or behavioral innovations by hunter-gatherers in the Late Pleistocene, or with the acquisition of farming in the Neolithic?

Methodology/Principal Findings

We investigate the timing of human population expansion by performing a multilocus analysis of≥20 unlinked autosomal noncoding regions, each consisting of ∼6 kilobases, resequenced in ∼184 individuals from 7 human populations. We test the hypothesis that the autosomal polymorphism data fit a simple two-phase growth model, and when the hypothesis is not rejected, we fit parameters of this model to our data using approximate Bayesian computation.

Conclusions/Significance

The data from the three surveyed non-African populations (French Basque, Chinese Han, and Melanesians) are inconsistent with the simple growth model, presumably because they reflect more complex demographic histories. In contrast, data from all four sub-Saharan African populations fit the two-phase growth model, and a range of onset times and growth rates is inferred for each population. Interestingly, both hunter-gatherers (San and Biaka) and food-producers (Mandenka and Yorubans) best fit models with population growth beginning in the Late Pleistocene. Moreover, our hunter-gatherer populations show a tendency towards slightly older and stronger growth (∼41 thousand years ago, ∼13-fold) than our food-producing populations (∼31 thousand years ago, ∼7-fold). These dates are concurrent with the appearance of the Late Stone Age in Africa, supporting the hypothesis that population growth played a significant role in the evolution of Late Pleistocene human cultures.

Culture, population structure, and low genetic diversity in Pleistocene hominins

Paleogenomic research has shown that modern humans, Neanderthals, and their most recent common ancestor have displayed less genetic diversity than living great apes. The traditional interpretation that low levels of genetic diversity in modern humans resulted from a relatively recent demographic bottleneck cannot account for similarly low levels of genetic diversity in Middle Pleistocene hominins. A more parsimonious hypothesis proposes that the effective population size of the human lineage has been low for more than 500,000 years, but the mechanism responsible for suppressing genetic diversity in Pleistocene hominin populations without similarly affecting that of their hominoid contemporaries remains unknown. Here we use agent-based simulation to study the effect of culturally mediated migration on neutral genetic diversity in structured populations. We show that, in populations structured by culturally mediated migration, selection can suppress neutral genetic diversity over thousands of generations, even in the absence of bottlenecks or expansions in census population size. In other words, selection could have suppressed the effective population size of Pleistocene hominins for as long as the degree of cultural similarity between regionally differentiated groups played an important role in mediating intraspecific gene flow.

Contrasting signatures of population growth for mitochondrial DNA and Y chromosomes among human populations in Africa

A history of Pleistocene population expansion has been inferred from the frequency spectrum of polymorphism in the mitochondrial DNA (mtDNA) of many human populations. Similar patterns are not typically observed for autosomal and X-linked loci. One explanation for this discrepancy is a recent population bottleneck, with different rates of recovery for haploid and autosomal loci as a result of their different effective population sizes. This hypothesis predicts that mitochondrial and Y chromosomal DNA will show a similar skew in the frequency spectrum in populations that have experienced a recent increase in effective population size. We test this hypothesis by resequencing 6.6 kb of noncoding Y chromosomal DNA and 780 basepairs of the mtDNA cytochrome c oxidase subunit III (COIII) gene in 172 males from 5 African populations. Four tests of population expansion are employed for each locus in each population: Fu's Fs statistic, the R(2) statistic, coalescent simulations, and the mismatch distribution. Consistent with previous results, patterns of mtDNA polymorphism better fit a model of constant population size for food-gathering populations and a model of population expansion for food-producing populations. In contrast, none of the tests reveal evidence of Y chromosome growth for either food-gatherers or food-producers. The distinct mtDNA and Y chromosome polymorphism patterns most likely reflect sex-biased demographic processes in the recent history of African populations. We hypothesize that males experienced smaller effective population sizes and/or lower rates of migration during the Bantu expansion, which occurred over the last 5,000 years.

Effective population size and tests of neutrality at cytoplasmic genes inArabidopsis

Cytoplasmic genomes typically lack recombination, implying that genetic hitch-hiking could be a predominant force structuring nucleotide polymorphism in the chloroplast and mitochondria. We test this hypothesis by analysing nucleotide polymorphism data at 28 loci across the chloroplast and mitochondria of the outcrossing plantArabidopsis lyrata, and compare patterns with multiple nuclear loci, and the highly selfingArabidopsis thaliana. The maximum likelihood estimate of the ratio of effective population size at cytoplasmic relative to nuclear genes inA. lyratadoes not depart from the neutral expectation of 0·5. Similarly, the ratio of effective size inA. thalianais close to unity, the neutral expectation for a highly selfing species. The results are thus consistent with neutral organelle polymorphism in these species or with comparable effects of hitch-hiking in both cytoplasmic and nuclear genes, in contrast to the results of recent studies on gynodioecious taxa. The four-gamete test and composite likelihood estimation provide evidence for very low levels of recombination in the organelles ofA. lyrata, although permutation tests do not suggest that adjacent polymorphic sites are more closely linked than more distant sites across the two genomes, suggesting that mutation hotspots or very low rates of gene conversion could explain the data.

Rapid adaptive divergence in new world achillea, an autopolyploid complex of ecological races

Adaptive evolution is often associated with speciation. In plants, however, ecotypic differentiation is common within widespread species, suggesting that climatic and edaphic specialization can outpace cladogenesis and the evolution of postzygotic reproductive isolation. We used cpDNA sequence (5 noncoding regions, 3.5 kb) and amplified fragment length polymorphisms (AFLPs: 4 primer pairs, 1,013 loci) to evaluate the history of ecological differentiation in the North American Achillea millefolium, an autopolyploid complex of "ecological races" exhibiting morphological, physiological, and life-history adaptations to diverse environments. Phylogenetic analyses reveal North American A. millefolium to be a monophyletic group distinct from its European and Asian relatives. Based on patterns of sequence divergence, as well as fossil and paleoecological data, colonization of North America appears to have occurred via the Bering Land Bridge during the Pleistocene (1.8 MYA to 11,500 years ago). Population genetic analyses indicate negligible structure within North American A. millefolium associated with varietal identity, geographic distribution, or ploidy level. North American populations, moreover, exhibit the signature of demographic expansion. These results affirm the "ecotype" concept of the North American Achillea advocated by classical research and demonstrate the rapid rate of ecological differentiation that sometimes occurs in plants.

Population size changes reshape genomic patterns of diversity

Elucidating the forces responsible for genomic variation is critical for understanding evolution. Under standard conditions, X-linked diversity is expected to be three-quarters the level of autosomal diversity. Empirical data often deviate from this prediction, but the reasons for these departures are unclear. We demonstrate that population size changes can greatly alter relative levels of X-linked and autosomal variation: population size reductions lead to particularly low X-linked diversity, whereas growth elevates X-linked relative to autosomal diversity. Genetic variation from a diverse array of taxa supports an important role for this effect in accounting for population differences in the ratio of X-linked to autosomal diversity. Consideration of this effect may improve the inference of population history and other evolutionary processes.

A critical review of adaptive genetic variation in Atlantic salmon: implications for conservation

Here we critically review the scale and extent of adaptive genetic variation in Atlantic salmon (Salmo salar L.), an important model system in evolutionary and conservation biology that provides fundamental insights into population persistence, adaptive response and the effects of anthropogenic change. We consider the process of adaptation as the end product of natural selection, one that can best be viewed as the degree of matching between phenotype and environment. We recognise three potential sources of adaptive variation: heritable variation in phenotypic traits related to fitness, variation at the molecular level in genes influenced by selection, and variation in the way genes interact with the environment to produce phenotypes of varying plasticity. Of all phenotypic traits examined, variation in body size (or in correlated characters such as growth rates, age of seaward migration or age at sexual maturity) generally shows the highest heritability, as well as a strong effect on fitness. Thus, body size in Atlantic salmon tends to be positively correlated with freshwater and marine survival, as well as with fecundity, egg size, reproductive success, and offspring survival. By contrast, the fitness implications of variation in behavioural traits such as aggression, sheltering behaviour, or timing of migration are largely unknown. The adaptive significance of molecular variation in salmonids is also scant and largely circumstantial, despite extensive molecular screening on these species. Adaptive variation can result in local adaptations (LA) when, among other necessary conditions, populations live in patchy environments, exchange few or no migrants, and are subjected to differential selective pressures. Evidence for LA in Atlantic salmon is indirect and comes mostly from ecological correlates in fitness-related traits, the failure of many translocations, the poor performance of domesticated stocks, results of a few common-garden experiments (where different populations were raised in a common environment in an attempt to dissociate heritable from environmentally induced phenotypic variation), and the pattern of inherited resistance to some parasites and diseases. Genotype x environment interactions occurr for many fitness traits, suggesting that LA might be important. However, the scale and extent of adaptive variation remains poorly understood and probably varies, depending on habitat heterogeneity, environmental stability and the relative roles of selection and drift. As maladaptation often results from phenotype-environment mismatch, we argue that acting as if populations are not locally adapted carries a much greater risk of mismanagement than acting under the assumption for local adaptations when there are none. As such, an evolutionary approach to salmon conservation is required, aimed at maintaining the conditions necessary for natural selection to operate most efficiently and unhindered. This may require minimising alterations to native genotypes and habitats to which populations have likely become adapted, but also allowing for population size to reach or extend beyond carrying capacity to encourage competition and other sources of natural mortality.

Reconstructing human origins in the genomic era

Analyses of recently acquired genomic sequence data are leading to important insights into the early evolution of anatomically modern humans, as well as into the more recent demographic processes that accompanied the global radiation of Homo sapiens. Some of the new results contradict early, but still influential, conclusions that were based on analyses of gene trees from mitochondrial DNA and Y-chromosome sequences. In this review, we discuss the different genetic and statistical methods that are available for studying human population history, and identify the most plausible models of human evolution that can accommodate the contrasting patterns observed at different loci throughout the genome.

Evolutionary relationships of Spirurina (Nematoda: Chromadorea: Rhabditida) with special emphasis on dracunculoid nematodes inferred from SSU rRNA gene sequences☆

The analysis of 26 new small subunit rRNA sequences obtained from helminths that primarily parasitize fishes sampled from five continents provided well-supported trees, allowing us to study the phylogenetic relationships among spirurid nematodes. The analyses have shown that Dracunculoidea is a paraphyletic taxon and Anguillicolidae and Gnathostomatidae constitute the basal branch of the suborder Spirurina. The genera Philometra and Philometroides appear to be paraphyletic, while on the higher taxonomic level, good correlation between the morphology-based system and molecular data was observed. Neither co-evolution of the studied helminths with their hosts, nor phylogeographic pattern, are apparent in our dataset.

HvrBase++: a phylogenetic database for primate species

HvrBase++ is the improved and extended version of HvrBase. Extensions are made by adding more population-based sequence samples from all primates including humans. The current collection comprises 13 873 hypervariable region I (HVRI) sequences and 4940 hypervariable region II (HVRII) sequences. In addition, we included 1376 complete mitochondrial genomes, 205 sequences from X-chromosomal loci and 202 sequences from autosomal chromosomes 1, 8, 11 and 16. In order to reduce the introduction of erroneous data into HvrBase++, we have developed a procedure that monitors GenBank for new versions of the current data in HvrBase++ and automatically updates the collection if necessary. For the stored sequences, supplementary information such as geographic origin, population affiliation and language of the sequence donor can be retrieved. HvrBase++ is Oracle based and easily accessible by a web interface (). As a new key feature, HvrBase++ provides an interactive graphical tool to easily access data from dynamically created geographical maps.

Interrogating multiple aspects of variation in a full resequencing data set to infer human population size changes

We present an expanded data set of 50 unlinked autosomal noncoding regions, resequenced in samples of Hausa from Cameroon, Italians, and Chinese. We use these data to make inferences about human demographic history by using a technique that combines multiple aspects of genetic data, including levels of polymorphism, the allele frequency spectrum, and linkage disequilibrium. We explore an extensive range of demographic parameters and demonstrate that our method of combining multiple aspects of the data results in a significant reduction of the compatible parameter space. In agreement with previous reports, we find that the Hausa data are compatible with demographic equilibrium as well as a set of recent population expansion models. In contrast to the Hausa, when multiple aspects of the data are considered jointly, the non-Africans depart from an equilibrium model of constant population size and are compatible with a range of simple bottleneck models, including a 50-90% reduction in effective population size occurring some time after the appearance of modern humans in Africa 160,000-120,000 years ago.

Heterogeneous patterns of variation among multiple human x-linked Loci: the possible role of diversity-reducing selection in non-africans

Studies of human DNA sequence polymorphism reveal a range of diversity patterns throughout the genome. This variation among loci may be due to natural selection, demographic influences, and/or different sampling strategies. Here we build on a continuing study of noncoding regions on the X chromosome in a panel of 41 globally sampled humans representing African and non-African populations by examining patterns of DNA sequence variation at four loci (APXL, AMELX, TNFSF5, and RRM2P4) and comparing these patterns with those previously reported at six loci in the same panel of 41 individuals. We also include comparisons with patterns of noncoding variation seen at five additional X-linked loci that were sequenced in similar global panels. We find that, while almost all loci show a reduction in non-African diversity, the magnitude of the reduction varies substantially across loci. The large observed variance in non-African levels of diversity results in the rejection of a neutral model of molecular evolution with a multi-locus HKA test under both a constant size and a bottleneck model. In non-Africans, some loci harbor an excess of rare mutations over neutral equilibrium predictions, while other loci show no such deviation in the distribution of mutation frequencies. We also observe a positive relationship between recombination rate and frequency spectra in our non-African, but not in our African, sample. These results indicate that a simple out-of-Africa bottleneck model is not sufficient to explain the observed patterns of sequence variation and that diversity-reducing selection acting at a subset of loci and/or a more complex neutral model must be invoked.

Nucleotide variation at Msn and Alas2, two genes flanking the centromere of the X chromosome in humans

The centromeric region of the X chromosome in humans experiences low rates of recombination over a considerable physical distance. In such a region, the effects of selection may extend to linked sites that are far away. To investigate the effects of this recombinational environment on patterns of nucleotide variability, we sequenced 4581 bp at Msn and 4697 bp at Alas2, two genes situated on either side of the X chromosome centromere, in a worldwide sample of 41 men, as well as in one common chimpanzee and one orangutan. To investigate patterns of linkage disequilibrium (LD) across the centromere, we also genotyped several informative sites from each gene in 120 men from sub-Saharan Africa. By studying X-linked loci in males, we were able to recover haplotypes and study long-range patterns of LD directly. Overall patterns of variability were remarkably similar at these two loci. Both loci exhibited (i) very low levels of nucleotide diversity (among the lowest seen in the human genome); (ii) a strong skew in the distribution of allele frequencies, with an excess of both very-low and very-high-frequency derived alleles in non-African populations; (iii) much less variation in the non-African than in the African samples; (iv) very high levels of population differentiation; and (v) complete LD among all sites within loci. We also observed significant LD between Msn and Alas2 in Africa, despite the fact that they are separated by approximately 10 Mb. These observations are difficult to reconcile with a simple demographic model but may be consistent with positive and/or purifying selection acting on loci within this large region of low recombination.

The allele frequency spectrum in genome-wide human variation data reveals signals of differential demographic history in three large world populations

We have studied a genome-wide set of single-nucleotide polymorphism (SNP) allele frequency measures for African-American, East Asian, and European-American samples. For this analysis we derived a simple, closed mathematical formulation for the spectrum of expected allele frequencies when the sampled populations have experienced nonstationary demographic histories. The direct calculation generates the spectrum orders of magnitude faster than coalescent simulations do and allows us to generate spectra for a large number of alternative histories on a multidimensional parameter grid. Model-fitting experiments using this grid reveal significant population-specific differences among the demographic histories that best describe the observed allele frequency spectra. European and Asian spectra show a bottleneck-shaped history: a reduction of effective population size in the past followed by a recent phase of size recovery. In contrast, the African-American spectrum shows a history of moderate but uninterrupted population expansion. These differences are expected to have profound consequences for the design of medical association studies. The analytical methods developed for this study, i.e., a closed mathematical formulation for the allele frequency spectrum, correcting the ascertainment bias introduced by shallow SNP sampling, and dealing with variable sample sizes provide a general framework for the analysis of public variation data.

Natural selection and molecular evolution in PTC, a bitter-taste receptor gene

The ability to taste phenylthiocarbamide (PTC) is a classic phenotype that has long been known to vary in human populations. This phenotype is of genetic, epidemiologic, and evolutionary interest because the ability to taste PTC is correlated with the ability to taste other bitter substances, many of which are toxic. Thus, variation in PTC perception may reflect variation in dietary preferences throughout human history and could correlate with susceptibility to diet-related diseases in modern populations. To test R. A. Fisher's long-standing hypothesis that variability in PTC perception has been maintained by balancing natural selection, we examined patterns of DNA sequence variation in the recently identified PTC gene, which accounts for up to 85% of phenotypic variance in the trait. We analyzed the entire coding region of PTC (1,002 bp) in a sample of 330 chromosomes collected from African (n=62), Asian (n=138), European (n=110), and North American (n=20) populations by use of new statistical tests for natural selection that take into account the potentially confounding effects of human population growth. Two intermediate-frequency haplotypes corresponding to "taster" and "nontaster" phenotypes were found. These haplotypes had similar frequencies across Africa, Asia, and Europe. Genetic differentiation between the continental population samples was low (FST=0.056) in comparison with estimates based on other genes. In addition, Tajima's D and Fu and Li's D and F statistics demonstrated a significant deviation from neutrality because of an excess of intermediate-frequency variants when human population growth was taken into account (P<.01). These results combine to suggest that balancing natural selection has acted to maintain "taster" and "nontaster" alleles at the PTC locus in humans.

Reduced polymorphism in the chimpanzee semen coagulating protein, semenogelin I

The semen of many primate species coagulates into a mating plug believed to prevent the sperm of subsequent mating events from accessing the ova. The texture of the coagulum varies among species: from a semisoft mass in humans to a firm plug in chimpanzees. In humans, a component of the coagulum, semenogelin I, also inhibits sperm motility. We tested the hypothesis that polymorphism and divergence at semenogelin I differ among hominoid species with different mating systems. Sequence data for the semenogelin I locus were obtained from 12 humans, 10 chimpanzees, 7 gorillas, and 1 bonobo. Mitochondrial D-loop data were collected from a subset of individuals to assess levels of variation at an unlinked locus. HKA tests using D-loop sequence data revealed a significant reduction of polymorphism at semenogelin I in chimpanzees, consistent with predictions of a selective sweep at this locus. This result was supported by independent HKA tests using polymorphism data from a putatively neutral locus from the literature. Humans show a similar trend toward reduced polymorphism, although HKA tests were only marginally significant. Gorilla sequence data show evidence of functional loss at the semenogelin I locus, indicated by stop codons within the putative open reading frame as well as high levels of polymorphism. Elevated Ka/Ks ratios within the Pan-Homo clade suggest a history of positive selection at semenogelin I. Our results suggest that there is a positive relationship between the intensity of sperm competition in a species and the strength of positive Darwinian selection on the seminal protein semenogelin I.

Patterns of human genetic diversity: implications for human evolutionary history and disease

Since the completion of the human genome sequencing project, the discovery and characterization of human genetic variation is a principal focus for future research. Comparative studies across ethnically diverse human populations and across human and nonhuman primate species is important for reconstructing human evolutionary history and for understanding the genetic basis of human disease. In this review, we summarize data on patterns of human genetic diversity and the evolutionary forces (mutation, genetic drift, migration, and selection) that have shaped these patterns of variation across both human populations and the genome. African population samples typically have higher levels of genetic diversity, a complex population substructure, and low levels of linkage disequilibrium (LD) relative to non-African populations. We discuss these differences and their implications for mapping disease genes and for understanding how population and genomic diversity have been important in the evolution, differentiation, and adaptation of humans.

Low genetic divergence obscures phylogeny among populations of Sphenodon, remnant of an ancient reptile lineage

Tuatara (two species of Sphenodon) are the last representatives of a branch of an ancient reptilian lineage, Sphenodontia, that have been isolated on the New Zealand landmass for 82 million years. We present analyses of geographic variation in allozymes, mitochondrial DNA, nuclear DNA sequences, and one-way albumin immunological comparisons. These all confirm a surprisingly low level of genetic diversity within Sphenodon for such an ancient lineage. We hypothesise a recent extended population bottleneck, probably during the Pliocene/Pleistocene glaciation cycles, to explain the current paucity of variation. All data sets reveal clear genetic differentiation between the northern populations and those in Cook Strait, but offer conflicting views of the history and taxonomic relationships of the Cook Strait population on North Brother Island, currently recognised as Sphenodon guntheri. Allozymes show this population to be the most divergent of all tuatara populations, but preliminary mitochondrial DNA data indicate few differences between S. guntheri and Cook Strait Sphenodon punctatus. Interpretation of the trees is confounded by the lack of a suitable outgroup. As in other cases of conflicting nuclear and mitochondrial data sets, the different data sets likely reveal different aspects of the animals' evolutionary history, and introgression is not uncommon between species pairs.

Human population structure and its effects on sampling Y chromosome sequence variation

The excess of rare variants in global sequencing studies of the nonrecombining portion of the Y chromosome (NRY) has been interpreted as evidence for the effects of human demographic expansion. However, many NRY polymorphisms are geographically localized and the effect of different geographical sampling on patterns of NRY variation is unknown. We use two sampling designs to detect population structure and its effects on patterns of human NRY polymorphism. First, we sequence 26.5 kb of noncoding Y chromosome DNA from 92 globally distributed males representing 35 populations. We find that the number of polymorphisms with singleton variants is positively correlated with the number of populations sampled and that there is a significant negative correlation of Tajima's D (TD) and Fu and Li's D (FD) statistics with the number of pooled populations. We then sequence the same region in a total of 73 males sampled from 3 distinct populations and find that TD and FD values for the 3 pooled and individual population samples were much less negative than those in the aforementioned global sample. Coalescent simulations show that a simple splitting model of population structure, with no changes in population size, is sufficient to produce the negative values of TD seen in our pooled samples. These empirical and simulation results suggest that observed levels of NRY population structure may lead to an upward bias in the number of singleton variants in global surveys and call into question inferences of population expansion based on global sampling strategies.

The study of structured populations--new hope for a difficult and divided science

Natural populations, including those of humans, have complex geographies and histories. Studying how they evolve is difficult, but it is possible with population-based DNA sequence data. However, the study of structured populations is divided by two distinct schools of thought and analysis. The phylogeographic approach is fundamentally graphical and begins with a gene-tree estimate. By contrast, the more traditional approach of using summary statistics is fundamentally mathematical. Both approaches have limitations, but there is promise in newer probabilistic methods that offer the flexibility and data exploitation of the phylogeographic approach in an explicitly model-based mathematical framework.

Linkage disequilibrium in human populations

Whereas the human linkage map appears on limited evidence to be constant over populations, maps of linkage disequilibrium (LD) vary among populations that differ in gene history. The greatest difference is between populations of sub-Saharan origin and populations remotely derived from Africa after a major bottleneck that reduced their heterozygosity and altered their Malecot parameters, increasing the intercept M that reflects association in founders and decreasing the exponential decline epsilon. Variation among populations within this ethnic dichotomy is much smaller. These observations validate use of a cosmopolitan LD map based on a sizeable sample representing a large population reliably typed for markers at high density. Then an LD map for a region or isolate within an ethnic group may be created by fitting the sample LD to the cosmopolitan map, estimating Malecot parameters simultaneously. The cosmopolitan map scaled by epsilon recovers 95% of the information that a local map at the same density gives and therefore more than the information in a low-resolution local map. Relative to a Eurasian cosmopolitan map the scaling factors are estimated to be 0.82 for isolates of European descent, 1.53 for Yorubans, and 1.74 for African Americans. These observations are consistent with a common bottleneck (perhaps but not necessarily speciation) approximately 173,500 years ago, if the bottleneck associated with migration out of Africa was 100,000 years ago. Eurasian populations (especially isolates with numerous cases) are efficient for genome scans, and populations of recent African origin (such as African Americans) are efficient for identification of causal polymorphisms within a candidate sequence.

Evidence for balancing selection from nucleotide sequence analyses of human G6PD

Glucose-6-phosphate dehydrogenase (G6PD) mutations that result in reduced enzyme activity have been implicated in malarial resistance and constitute one of the best examples of selection in the human genome. In the present study, we characterize the nucleotide diversity across a 5.2-kb region of G6PD in a sample of 160 Africans and 56 non-Africans, to determine how selection has shaped patterns of DNA variation at this gene. Our global sample of enzymatically normal B alleles and A, A-, and Med alleles with reduced enzyme activities reveals many previously uncharacterized silent-site polymorphisms. In comparison with the absence of amino acid divergence between human and chimpanzee G6PD sequences, we find that the number of G6PD amino acid polymorphisms in human populations is significantly high. Unlike many other G6PD-activity alleles with reduced activity, we find that the age of the A variant, which is common in Africa, may not be consistent with the recent emergence of severe malaria and therefore may have originally had a historically different adaptive function. Overall, our observations strongly support previous genotype-phenotype association studies that proposed that balancing selection maintains G6PD deficiencies within human populations. The present study demonstrates that nucleotide sequence analyses can reveal signatures of both historical and recent selection in the genome and may elucidate the impact that infectious disease has had during human evolution.

Genetic evidence for long-term population decline in a savannah-dwelling primate: inferences from a hierarchical bayesian model

The purpose of this study was to test for evidence that savannah baboons (Papio cynocephalus) underwent a population expansion in concert with a hypothesized expansion of African human and chimpanzee populations during the late Pleistocene. The rationale is that any type of environmental event sufficient to cause simultaneous population expansions in African humans and chimpanzees would also be expected to affect other codistributed mammals. To test for genetic evidence of population expansion or contraction, we performed a coalescent analysis of multilocus microsatellite data using a hierarchical Bayesian model. Markov chain Monte Carlo (MCMC) simulations were used to estimate the posterior probability density of demographic and genealogical parameters. The model was designed to allow interlocus variation in mutational and demographic parameters, which made it possible to detect aberrant patterns of variation at individual loci that could result from heterogeneity in mutational dynamics or from the effects of selection at linked sites. Results of the MCMC simulations were consistent with zero variance in demographic parameters among loci, but there was evidence for a 10- to 20-fold difference in mutation rate between the most slowly and most rapidly evolving loci. Results of the model provided strong evidence that savannah baboons have undergone a long-term historical decline in population size. The mode of the highest posterior density for the joint distribution of current and ancestral population size indicated a roughly eightfold contraction over the past 1,000 to 250,000 years. These results indicate that savannah baboons apparently did not share a common demographic history with other codistributed primate species.

Assessing DNA sequence variations in human ESTs in a phylogenetic context using high-density oligonucleotide arrays

We have analyzed human genomic diversity in 32 individuals representing four continental populations of Homo sapiens in the context of four ape species. We used DNA resequencing chips covering 898 expressed sequence tags (ESTs), corresponding to 109 kb of sequence. Based on the intra-species data, the neutral hypothesis could not be rejected. However, the mutation rate was two times lower than typically observed in functionally unconstrained genomic segments, suggesting a certain level of selection. The worldwide diversity (297 segregating sites and nucleotide diversity of 0.054%) was partitioned among continents, with the greatest amount of variation observed in the African sample. The long-term effective population size of the human population was estimated at 13,000; a similar figure was obtained for the African sample and a 20% lower estimate was obtained for the other continents. Africans also differed in having a higher number of continental-specific polymorphisms contributing to the higher average nucleotide diversity. These results are consistent with the existence of two distinct lineages of modern humans: amalgamation of these lineages in Africa led to the higher present-day diversity on that continent, whereas colonization of other continents by one of them gave the effect of a population bottleneck.

Distinguishing between selection and population expansion in an experimental lineage of bacteriophage T7

Experimental evolution of short-lived organisms offers the opportunity to study the dynamics of polymorphism over time in a controlled environment. Here, we characterize DNA polymorphism data over time for four genes in bacteriophage T7. Our experiment ran for 2500 generations and populations were sampled after 500, 2000, and 2500 generations. We detect positive selection, purifying ("negative") selection, and population expansion in our experiment. We also present a statistical test that is able to distinguish demographic from selective events, processes that are hard to identify individually because both often produce an excess of rare mutations. Our "heterogeneity test" modifies common statistics measuring the frequency spectrum of polymorphism (e.g., Fu and Li's D) by looking for processes producing different patterns on nonsynonymous and synonymous mutations. Test results agree with the known conditions of the experiment, and we are therefore confident that this test offers a tool to evaluate natural populations. Our results suggest that instances of segregating deleterious mutations may be common, but as yet undetected, in nature.

Microsatellite data resolve phylogeographic patterns in European grayling, Thymallus thymallus, Salmonidae

The phylogeography of an endangered salmonid, European grayling (Thymallus thymallus), was studied based on analysis of 17 nuclear microsatellite DNA loci. In agreement with earlier mitochondrial DNA (mtDNA) studies, phylogenetic relationships of the populations suggested that northern Europe was colonized from two distinct Pleistocene refugia. Furthermore, microsatellites revealed highly supported grouping of mainland Swedish, Norwegian, Danish, German and Slovenian populations, suggesting that grayling from the northwestern and central Europe have descended from their southern conspecifics. The level of divergence between populations was substantial, even across short geographical distances. Although this was in part due to postglacial colonization patterns and contemporary barriers for gene flow, the high divergence estimates between hydrologically connected sampling locations implied efficient interpopulation reproductive isolation. Microsatellites revealed that the populations exhibited, on average, only 3.5 (+/-2.2) alleles per locus, indicating that T. thymallus has strikingly low levels of intrapopulation genetic diversity as compared with other freshwater fish species. Accordingly, as indicated by analysis of molecular variance (AMOVA), only 49.1-58.0% of the total grayling microsatellite diversity resided within populations. A latitudinal genetic diversity gradient, potentially resulting from glaciation-mediated founder events, was not evident. Alternatively, it is possible that grayling display limited dispersal behaviour/capability, leading to low long-term effective population sizes and, consequently, depauperate intrapopulation polymorphism. These findings have implications for conservation of T. thymallus. Importantly, they exemplify that microsatellites can be highly informative for intraspecific phylogeography studies dealing with substantial divergence scales.

Divergence of mitochondrial dna is not corroborated by nuclear dna, morphology, or behavior in Drosophila simulans

We ask whether the observed mitochondrial DNA (mtDNA) population subdivision of Drosophila simulans is indicative of organismal structure or of specific processes acting on the mitochondrial genome. Factors either intrinsic or extrinsic to the host genome may influence the evolutionary dynamics of mtDNA. Potential intrinsic factors include adaptation of the mitochondrial genome and of nucleomitochondrial gene complexes specific to the local environment. An extrinsic force that has been shown to influence mtDNA evolution in invertebrates is the bacterial endosymbiont Wolbachia. Evidence presented in this study suggests that mtDNA is not a good indicator of organismal subdivision in D. simulans. Furthermore, there is no evidence to suggest that Wolbachia causes any reduction in nuclear gene flow in this species. The observed differentiation in mtDNA is not corroborated by data from NADH: ubiquinone reductase 75kD subunit precursor or the Alcohol dehydrogenase-related loci, from the shape or size of the male genital arch, or from assortative premating behavior. We discuss these results in relation to a mitochondrial genetic species concept and the potential for Wolbachia-induced incompatibility to be a mechanism of speciation in insects. We conclude with an iterated appeal to include phylogenetic and statistical tests of neutrality as a supplement to phylogenetic and population genetic analyses when using mtDNA as an evolutionary marker.

Patterns of genetic variation in Mendelian and complex traits

This review discusses the prospects for understanding the genetic basis of complex traits in humans. We take the view that work done on Drosophila melanogaster can serve as a model for understanding complex traits in humans, and the literature on this model system, as well as on humans, is reviewed. The prospects for success in understanding the genetic basis of complex traits depend, in part, on the nature of the forces acting on genetic variation. We suggest that different experimental approaches should be undertaken for traits caused by common genetic variants versus those arising from rare genetic variants.

Methods to detect selection in populations with applications to the human

The development of statistical tests of natural selection at the DNA level in population samples has been ongoing for the past 13 years. The current state of the field is reviewed, and the available tests of selection are described. All tests use predictions from the theory of neutrally evolving sites as a null hypothesis. Departures from equilibrium-neutral expectations can indicate the presence of natural selection acting either at one or more of the sites under investigation or at a sufficiently tightly linked site. Complications can arise in the interpretation of departures from neutrality if populations are not at equilibrium for mutation and genetic drift or if populations are subdivided, both of which are likely scenarios for humans. Attempts to understand the nonequilibrium configuration of silent polymorphism in human mitochondrial DNA illustrate the difficulty of distinguishing between selection and alternative demographic hypotheses. The range of plausible alternatives to selection will become better defined, however, as additional population genetic data sets become available, allowing better null models to be constructed.

Genetic perspectives on human origins and differentiation

This is a review of genetic evidence about the ancient demography of the ancestors of our species and about the genesis of worldwide human diversity. The issue of whether or not a population size bottleneck occurred among our ancestors is under debate among geneticists as well as among anthropologists. The bottleneck, if it occurred, would confirm the Garden of Eden (GOE) model of the origin of modern humans. The competing model, multiregional evolution (MRE), posits that the number of human ancestors has been large, occupying much of the temperate Old World for the last two million years. While several classes of genetic marker seem to contain a strong signal of demographic recovery from a small number of ancestors, other nuclear loci show no such signal. The pattern at these loci is compatible with the existence of widespread balancing selection in humans. The study of human diversity at (putatively) neutral genetic marker loci has been hampered since the beginning by ascertainment bias since they were discovered in Europeans. The high levels of polymorphism at microsatellite loci means that they are free of this bias. Microsatellites exhibit a clear almost linear diversity gradient away from Africa, so that New World populations are approximately 15% less diverse than African populations. This pattern is not compatible with a model of a single large population expansion and colonization of most of the Earth by our ancestors but suggests, instead, gradual loss of diversity in successive colonization bottlenecks as our species grew and spread.

The discovery of single-nucleotide polymorphisms--and inferences about human demographic history

A method of historical inference that accounts for ascertainment bias is developed and applied to single-nucleotide polymorphism (SNP) data in humans. The data consist of 84 short fragments of the genome that were selected, from three recent SNP surveys, to contain at least two polymorphisms in their respective ascertainment samples and that were then fully resequenced in 47 globally distributed individuals. Ascertainment bias is the deviation, from what would be observed in a random sample, caused either by discovery of polymorphisms in small samples or by locus selection based on levels or patterns of polymorphism. The three SNP surveys from which the present data were derived differ both in their protocols for ascertainment and in the size of the samples used for discovery. We implemented a Monte Carlo maximum-likelihood method to fit a subdivided-population model that includes a possible change in effective size at some time in the past. Incorrectly assuming that ascertainment bias does not exist causes errors in inference, affecting both estimates of migration rates and historical changes in size. Migration rates are overestimated when ascertainment bias is ignored. However, the direction of error in inferences about changes in effective population size (whether the population is inferred to be shrinking or growing) depends on whether either the numbers of SNPs per fragment or the SNP-allele frequencies are analyzed. We use the abbreviation "SDL," for "SNP-discovered locus," in recognition of the genomic-discovery context of SNPs. When ascertainment bias is modeled fully, both the number of SNPs per SDL and their allele frequencies support a scenario of growth in effective size in the context of a subdivided population. If subdivision is ignored, however, the hypothesis of constant effective population size cannot be rejected. An important conclusion of this work is that, in demographic or other studies, SNP data are useful only to the extent that their ascertainment can be modeled.

Determinants of effective population size for loci with different modes of inheritance

Here we report an assessment of the determinants of effective population size (N(e)) in species with overlapping generations. Specifically, we used a stochastic demographic model to investigate the influence of different life-history variables on N(e)/N (where N = population census number) and the influence of sex differences in life-history variables on N(e) for loci with different modes of inheritance. We applied an individual-based modeling approach to two datasets: one from a natural population of savannah baboons (Papio cynocephalus) in the Amboseli basin of southern Kenya and one from a human tribal population (the Gainj of Papua New Guinea). Simulation-based estimates of N(e)/N averaged 0.329 for the Amboseli baboon population (SD = 0.116, 95% CI = 0.172 - 0.537) and 0.786 for the Gainj (SD = 0.184, 95% CI = 0.498 - 1.115). Although variance in male fitness had a substantial impact on N(e)/N in each of the two primate populations, ratios of N(e) values for autosomal and sex-linked loci exhibited no significant departures from Poisson-expected values. In each case, similarities in sex-specific N(e) values were attributable to the unexpectedly high variance in female fitness. Variance in male fitness resulted primarily from age-dependent variance in reproductive success, whereas variance in female fitness resulted primarily from stochastic variance in survival during the reproductive phase.

Hierarchical patterns of global human Y-chromosome diversity

We examined 43 biallelic polymorphisms on the nonrecombining portion of the Y chromosome (NRY) in 50 human populations encompassing a total of 2,858 males to study the geographic structure of Y-chromosome variation. Patterns of NRY diversity varied according to geographic region and method/level of comparison. For example, populations from Central Asia had the highest levels of heterozygosity, while African populations exhibited a higher level of mean pairwise differences among haplotypes. At the global level, 36% of the total variance of NRY haplotypes was attributable to differences among populations (i.e., Phi(ST) = 0.36). When a series of AMOVA analyses was performed on different groupings of the 50 populations, high levels of among-groups variance (Phi(CT)) were found between Africans, Native Americans, and a single group containing all 36 remaining populations. The same three population groupings formed distinct clusters in multidimensional scaling plots. A nested cladistic analysis (NCA) demonstrated that both population structure processes (recurrent gene flow restricted by isolation by distance and long-distance dispersals) and population history events (contiguous range expansions and long-distance colonizations) were instrumental in explaining this tripartite division of global NRY diversity. As in our previous analyses of smaller NRY data sets, the NCA detected a global contiguous range expansion out of Africa at the level of the total cladogram. Our new results support a general scenario in which, after an early out-of-Africa range expansion, global-scale patterns of NRY variation were mainly influenced by migrations out of Asia. Two other notable findings of the NCA were (1) Europe as a "receiver" of intercontinental signals primarily from Asia, and (2) the large number of intracontinental signals within Africa. Our AMOVA analyses also supported the hypothesis that patrilocality effects are evident at local and regional scales, rather than at intercontinental and global levels. Finally, our results underscore the importance of subdivision of the human paternal gene pool and imply that caution should be exercised when using models and experimental strategies based on the assumption of panmixia.

Y-chromosome mismatch distributions in Europe

Ancient demographic events can be inferred from the distribution of pairwise sequence differences (or mismatches) among individuals. We analyzed a database of 3,677 Y chromosomes typed for 11 biallelic markers in 48 human populations from Europe and the Mediterranean area. Contrary to what is observed in the analysis of mitochondrial polymorphisms, Tajima's test was insignificant for most Y-chromosome samples, and in 47 populations the mismatch distributions had multiple peaks. Taken at face value, these results would suggest either (1) that the size of the male population stayed essentially constant over time, while the female population size increased, or (2) that different selective regimes have shaped mitochondrial and Y-chromosome diversity, leading to an excess of rare alleles only in the mitochondrial genome. An alternative explanation would be that the 11 variable sites of the Y chromosome do not provide sufficient statistical power, so a comparison with mitochondrial data (where more than 200 variable sites are studied in Europe) is impossible at present. To discriminate between these possibilities, we repeatedly analyzed a European mitochondrial database, each time considering only 11 variable sites, and we estimated mismatch distributions in stable and growing populations, generated by simulating coalescent processes. Along with theoretical considerations, these tests suggest that the difference between the mismatch distributions inferred from mitochondrial and Y-chromosome data are not a statistical artifact. Therefore, the observed mismatch distributions appear to reflect different underlying demographic histories and/or selective pressures for maternally and paternally transmitted loci.

Contrasting patterns of polymorphisms at the ABO-secretor gene (FUT2) and plasma alpha(1,3)fucosyltransferase gene (FUT6) in human populations

The coding sequences ( approximately 1 kb) of FUT2 [ABO-Secretor type alpha(1,2)fucosyltransferase] and of FUT6 [plasma alpha(1,3)fucosyltransferase] were analyzed for allelic polymorphism by direct sequencing in five populations. The nucleotide diversities of FUT2 estimated from pairwise sequence differences were 0.0045, 0.0042, 0.0042, 0.0009, and 0.0008 in Africans, European-Africans, Iranians, Chinese, and Japanese, respectively. The nucleotide diversities of FUT6 were 0.0024, 0.0016, 0.0015, 0.0017, and 0.0020 in Africans, European-Africans, Iranians, Chinese, and Japanese, respectively. At FUT2, excesses in pairwise sequence differences compared to the number of polymorphic sites as indicated by a significantly positive Tajima's D were observed in European-Africans and in Iranians. The data do not fit expectations of the equilibrium neutral model with an infinite number of sites. On the other hand, Tajima's D's at FUT6 in each of the five populations and at FUT2 in Africans, Chinese, and Japanese were not significantly different from zero. F(ST) between the Asians and the others measured at FUT2 was higher than at FUT6. These results suggest that natural selection was responsible for the generation of the FUT2 polymorphism in European-Africans and in Iranians.

Sequence variability of a human pseudogene

We have obtained haplotypes from the autosomal glucocerebrosidase pseudogene (psGBA) for 100 human chromosomes from worldwide populations, as well as for four chimpanzee and four gorilla chromosomes. In humans, in a 5420-nucleotide stretch analyzed, variation comprises 17 substitutions, a 3-bp deletion, and a length polymorphism at a polyadenine tract. The substitution rate on the pseudogene (1.23 +/- 0.22 x 10(-9) per nucleotide and year) is within the range of previous estimates considering phylogenetic estimations. Recombination within the pseudogene was recognized, although the low variability of this locus prevented an accurate measure of recombination rates. At least 13% of the psGBA sequence could be attributed to gene conversion from the contiguous GBA gene, whereas the reciprocal event has been shown to lead to Gaucher disease. Human psGBA sequences showed a recent coalescence time (approximately 200,000 yr ago), and the most ancestral haplotype was found only in Africans; both observations are compatible with the replacement hypothesis of human origins. In a deeper timeframe, phylogenetic analysis showed that the duplication event that created psGBA could be dated at approximately 27 million years ago, in agreement with previous estimates.

Human populations show reduced DNA sequence variation at the Factor IX locus

Levels and patterns of human DNA sequence variation vary widely among loci. However, some of this variation may be due to the different populations used in different studies. So far, few studies of diverse human populations have compared different genetic loci for the same samples of populations and individuals. Here, we present new polymorphism data from intron 4 of the Factor IX gene (FIX) sequenced in diverse Old World populations. An explicit comparison is made with another X-linked gene, PDHA1, for which the sampling of individuals was very similar. Despite having a similar amount of divergence from chimpanzees, as do other nuclear genes, FIX has comparatively much less DNA sequence variation among humans. Nucleotide diversity at FIX is the lowest among the existing non-Y chromosome nuclear gene datasets and is less than 10% of the diversity found at PDHA1. Estimates of effective population size based on FIX are 8,558, about half of the value obtained for PDHA1, and the time to the most recent common ancestry among human FIX gene copies (282,000 years) is one of the most recent estimates reported for human genes. Analyses presented here suggest a history for the FIX region that includes recent positive directional selection, or background, selection. The general conclusion emerging is that very large variations can exist between the histories of similar genomic regions, even when sampling differences are minimized.

Brief communication: paleoanthropology and the population genetics of ancient genes

The Mezmaiskaya cave mtDNA is similar in many ways to the Feldhofer cave Neandertal sequence and the more recently obtained Vindija cave sequence. If we accept the contention that the Mezmaiskaya cave specimen is a Neandertal infant, its mtDNA provides no new information about the fate of the European Neandertals. However, there is reason to believe that the Mezmaiskaya cave infant is not a Neandertal, and this places its importance in another light, because it delimits the possible hypotheses of Neandertal and recent human genetic relationships. One possibility is a that the pattern found in ancient mtDNA results from the replacement of an isolated gene pool (Neandertals) by one of its contemporaries (modern humans). A second possibility is natural selection expressed as the substitution of an advantageous mtDNA variant within a single large species, including both Neandertals and modern humans. The geologic, archaeological, and dating evidence shows the Mezmaiskaya cave infant to be a burial from a level even more recent than the Upper Paleolithic preserved at the site, and its anatomy does not contradict the assessment that the Mezmaiskaya cave infant is not a Neandertal. Therefore, the second pattern can be favored over the first.