Berbers and Arabs: Tracing the genetic diversity and history of Southern Tunisia through genome wide analysis.

OBJECTIVES: Tunisia has been a crossroads for people from Africa, Europe, and the Middle East since prehistoric times. At present, it is inhabited by two main ethnic groups, Arabs and Berbers, and several minorities. This study aims to advance knowledge regarding their genetic structure using new population samplings and a genome-wide approach. MATERIALS AND METHODS: We investigated genomic variation, estimated ancestry components and dated admixture events in three Berber and two Arab populations from Southern Tunisia, mining a dataset including Middle Eastern, sub-Saharan, and European populations. RESULTS: Differences in the proportion of North African, Arabian, and European ancestries and the varying impact of admixture and isolation determined significant heterogeneity in the genetic structure of Southern Tunisian populations. Admixture time estimates show a multilayer pattern of admixture events, involving both ethno-linguistic groups, which started around the mid XI century and lasted for nearly five centuries. DISCUSSION: Our study provides evidence that the relationships between genetic and cultural diversity of old and new inhabitants of North Africa in southern Tunisia follow different patterns. The Berbers seem to have preserved a significant part of their common genomic heritage despite Islamization, Arab cultural influence, and linguistic diversity. Compared to Morocco and Algeria, southern Tunisian Arabs have retained a higher level of Arabian ancestry. This is more evident in the semi-nomad R'Baya, who have kept their original Bedouin lifestyle, than in the population from Douz, who have undergone multiple events of stratification and admixture.

Evaluating mtDNA patterns of genetic isolation using a re-sampling procedure: A case study on Italian populations.

BACKGROUND: A number of studies which have investigated isolation patterns in human populations rely on the analysis of intra- and inter-population genetic statistics of mtDNA polymorphisms. However, this approach makes it difficult to differentiate between the effects of long-term genetic isolation and the random fluctuations of statistics due to reduced sample size. AIM: To overcome the confounding effect of sample size when detecting signatures of genetic isolation. SUBJECTS AND METHODS: A re-sampling based procedure was employed to evaluate reduction in intra-population diversity, departure from surrounding genetic background and demographic stationarity in 34 Italian populations subject to isolation factors. RESULTS: Signatures of genetic isolation were detected for all three statistics in seven populations: Pusteria valley, Sappada, Sauris, Timau settled in the eastern Italian Alps and Cappadocia, Filettino and Vallepietra settled in the Appenines. However, this study was unable to find signals for any of the statistics analysed in 19 populations. Finally, eight populations showing signals of isolation were found for one or two statistics. CONCLUSION: The analysis revealed that the use of population genetic statistics combined with re-sampling procedure can help detect signatures of genetic isolation in human populations, even using a single, although highly informative, locus like mtDNA.

The Y-chromosome tree bursts into leaf: 13,000 high-confidence SNPs covering the majority of known clades.

Many studies of human populations have used the male-specific region of the Y chromosome (MSY) as a marker, but MSY sequence variants have traditionally been subject to ascertainment bias. Also, dating of haplogroups has relied on Y-specific short tandem repeats (STRs), involving problems of mutation rate choice, and possible long-term mutation saturation. Next-generation sequencing can ascertain single nucleotide polymorphisms (SNPs) in an unbiased way, leading to phylogenies in which branch-lengths are proportional to time, and allowing the times-to-most-recent-common-ancestor (TMRCAs) of nodes to be estimated directly. Here we describe the sequencing of 3.7 Mb of MSY in each of 448 human males at a mean coverage of 51×, yielding 13,261 high-confidence SNPs, 65.9% of which are previously unreported. The resulting phylogeny covers the majority of the known clades, provides date estimates of nodes, and constitutes a robust evolutionary framework for analyzing the history of other classes of mutation. Different clades within the tree show subtle but significant differences in branch lengths to the root. We also apply a set of 23 Y-STRs to the same samples, allowing SNP- and STR-based diversity and TMRCA estimates to be systematically compared. Ongoing purifying selection is suggested by our analysis of the phylogenetic distribution of nonsynonymous variants in 15 MSY single-copy genes.

A degradation-sensitive anionic trypsinogen (PRSS2) variant protects against chronic pancreatitis.

Chronic pancreatitis is a common inflammatory disease of the pancreas. Mutations in the genes encoding cationic trypsinogen (PRSS1) and the pancreatic secretory trypsin inhibitor (SPINK1) are associated with chronic pancreatitis. Because increased proteolytic activity owing to mutated PRSS1 enhances the risk for chronic pancreatitis, mutations in the gene encoding anionic trypsinogen (PRSS2) may also predispose to disease. Here we analyzed PRSS2 in individuals with chronic pancreatitis and controls and found, to our surprise, that a variant of codon 191 (G191R) is overrepresented in control subjects: G191R was present in 220/6,459 (3.4%) controls but in only 32/2,466 (1.3%) affected individuals (odds ratio 0.37; P = 1.1 x 10(-8)). Upon activation by enterokinase or trypsin, purified recombinant G191R protein showed a complete loss of trypsin activity owing to the introduction of a new tryptic cleavage site that renders the enzyme hypersensitive to autocatalytic proteolysis. In conclusion, the G191R variant of PRSS2 mitigates intrapancreatic trypsin activity and thereby protects against chronic pancreatitis.

A revised root for the human Y chromosomal phylogenetic tree: the origin of patrilineal diversity in Africa.

To shed light on the structure of the basal backbone of the human Y chromosome phylogeny, we sequenced about 200 kb of the male-specific region of the human Y chromosome (MSY) from each of seven Y chromosomes belonging to clades A1, A2, A3, and BT. We detected 146 biallelic variant sites through this analysis. We used these variants to construct a patrilineal tree, without taking into account any previously reported information regarding the phylogenetic relationships among the seven Y chromosomes here analyzed. There are several key changes at the basal nodes as compared with the most recent reference Y chromosome tree. A different position of the root was determined, with important implications for the origin of human Y chromosome diversity. An estimate of 142 KY was obtained for the coalescence time of the revised MSY tree, which is earlier than that obtained in previous studies and easier to reconcile with plausible scenarios of modern human origin. The number of deep branchings leading to African-specific clades has doubled, further strengthening the MSY-based evidence for a modern human origin in the African continent. An analysis of 2204 African DNA samples showed that the deepest clades of the revised MSY phylogeny are currently found in central and northwest Africa, opening new perspectives on early human presence in the continent.

The influence of habitats on female mobility in Central and Western Africa inferred from human mitochondrial variation.

BACKGROUND: When studying the genetic structure of human populations, the role of cultural factors may be difficult to ascertain due to a lack of formal models. Linguistic diversity is a typical example of such a situation. Patrilocality, on the other hand, can be integrated into a biological framework, allowing the formulation of explicit working hypotheses. The present study is based on the assumption that patrilocal traditions make the hypervariable region I of the mtDNA a valuable tool for the exploration of migratory dynamics, offering the opportunity to explore the relationships between genetic and linguistic diversity. We studied 85 Niger-Congo-speaking patrilocal populations that cover regions from Senegal to Central African Republic. A total of 4175 individuals were included in the study. RESULTS: By combining a multivariate analysis aimed at investigating the population genetic structure, with a Bayesian approach used to test models and extent of migration, we were able to detect a stepping-stone migration model as the best descriptor of gene flow across the region, with the main discontinuities corresponding to forested areas. CONCLUSIONS: Our analyses highlight an aspect of the influence of habitat variation on human genetic diversity that has yet to be understood. Rather than depending simply on geographic linear distances, patterns of female genetic variation vary substantially between savannah and rainforest environments. Our findings may be explained by the effects of recent gene flow constrained by environmental factors, which superimposes on a background shaped by pre-agricultural peopling.

An early divergence of KhoeSan ancestors from those of other modern humans is supported by an ABC-based analysis of autosomal resequencing data.

Sub-Saharan Africa has consistently been shown to be the most genetically diverse region in the world. Despite the fact that a substantial portion of this variation is partitioned between groups practicing a variety of subsistence strategies and speaking diverse languages, there is currently no consensus on the genetic relationships of sub-Saharan African populations. San (a subgroup of KhoeSan) and many Pygmy groups maintain hunter-gatherer lifestyles and cluster together in autosomal-based analysis, whereas non-Pygmy Niger-Kordofanian speakers (non-Pygmy NKs) predominantly practice agriculture and show substantial genetic homogeneity despite their wide geographic range throughout sub-Saharan Africa. However, KhoeSan, who speak a set of relatively unique click-based languages, have long been thought to be an early branch of anatomically modern humans based on phylogenetic analysis. To formally test models of divergence among the ancestors of modern African populations, we resequenced a sample of San, Eastern, and Western Pygmies and non-Pygmy NKs individuals at 40 nongenic (∼2 kb) regions and then analyzed these data within an Approximate Bayesian Computation (ABC) framework. We find substantial support for a model of an early divergence of KhoeSan ancestors from a proto-Pygmy-non-Pygmy NKs group ∼110 thousand years ago over a model incorporating a proto-KhoeSan-Pygmy hunter-gatherer divergence from the ancestors of non-Pygmy NKs. The results of our analyses are consistent with previously identified signals of a strong bottleneck in Mbuti Pygmies and a relatively recent expansion of non-Pygmy NKs. We also develop a number of methodologies that utilize "pseudo-observed" data sets to optimize our ABC-based inference. This approach is likely to prove to be an invaluable tool for demographic inference using genome-wide resequencing data.

Reevaluating a model of gender-biased gene flow among Sub-Saharan Hunter-gatherers and farmers.

In a previous study, we proposed a model for genetic admixture between African hunter-gatherers and food producers, in which we integrated demographic and genetic aspects together with ethnographic knowledge (Destro-Bisol et al. 2004b). In that study it was possible to test the model only using genetic information from widely dispersed and genetically heterogeneous populations. Here we reevaluate the congruence between the model and patterns of genetic variation using an anthropologically and geographically more homogeneous data set that includes Pygmies and farmers from Cameroon, Congo, and the Central African Republic. As implied by the model, the ratios of mtDNA to Y chromosome Nm estimates (effective population size, N, times the migration rate, m; 0.154 in Pygmies and 6.759 in farmers), support an asymmetric gene flow, with a higher Bantu-to-Pygmy gene flow for paternal than for maternal lineages, and vice versa for farmers. Analyses of intra- and interpopulation genetic variation further support the above observation, showing a prevailing effect of genetic drift on maternal lineages and gene flow on paternal lineages among Pygmies, and an opposite pattern among farmers. We also detected differences between patterns for classical and molecular measures of Y chromosome intrapopulation variation, which likely represent signatures of the introgression of Bantu lineages into the gene pool of Pygmy populations. On the whole, our results seem to reflect differences in the demographic history and the degree of patrilocality and polygyny between the two population groups, thus providing further support to our microevolutionary model in an anthropologically coherent framework.

The genomic echoes of the last Green Sahara on the Fulani and Sahelian people.

The population history of the Sahara/Sahelian belt is understudied, despite previous work highlighting complex dynamics.1,2,3,4,5,6,7 The Sahelian Fulani, i.e., the largest nomadic pastoral population in the world,8 represent an interesting case because they show a non-negligible proportion of an Eurasian genetic component, usually explained by recent admixture with northern Africans.1,2,5,6,7,9,10,11,12 Nevertheless, their origins are largely unknown, although several hypotheses have been proposed, including a possible link to ancient peoples settled in the Sahara during its last humid phase (Green Sahara, 12,000-5,000 years before present [BP]).13,14,15 To shed light about the Fulani ancient genetic roots, we produced 23 high-coverage (30×) whole genomes from Fulani individuals from 8 Sahelian countries, plus 17 samples from other African groups and 3 from Europeans as controls, for a total of 43 new whole genomes. These data have been compared with 814 published modern whole genomes2,16,17,18 and with relevant published ancient sequences (> 1,800 samples).19 These analyses showed some evidence that the non-sub-Saharan genetic ancestry component of the Fulani might have also been shaped by older events,1,5,6 possibly tracing the Fulani origins to unsampled ancient Green Saharan population(s). The joint analysis of modern and ancient samples allowed us to shed light on the genetic ancestry composition of such ancient Saharans, suggesting a similarity with Late Neolithic Moroccans and possibly pointing to a link with the spread of cattle herding. We also identified two different Fulani clusters whose admixture pattern may be informative about the historical Fulani movements and their later involvement in the western African empires.

Insights into the demographic history of African Pygmies from complete mitochondrial genomes.

Pygmy populations are among the few hunter-gatherers currently living in sub-Saharan Africa and are mainly represented by two groups, Eastern and Western, according to their current geographical distribution. They are scattered across the Central African belt and surrounded by Bantu-speaking farmers, with whom they have complex social and economic interactions. To investigate the demographic history of Pygmy groups, a population approach was applied to the analysis of 205 complete mitochondrial DNA (mtDNA) sequences from ten central African populations. No sharing of maternal lineages was observed between the two Pygmy groups, with haplogroup L1c being characteristic of the Western group but most of Eastern Pygmy lineages falling into subclades of L0a, L2a, and L5. Demographic inferences based on Bayesian coalescent simulations point to an early split among the maternal ancestors of Pygmies and those of Bantu-speaking farmers (∼ 70,000 years ago [ya]). Evidence for population growth in the ancestors of Bantu-speaking farmers has been observed, starting ∼ 65,000 ya, well before the diffusion of Bantu languages. Subsequently, the effective population size of the ancestors of Pygmies remained constant over time and ∼ 27,000 ya, coincident with the Last Glacial Maximum, Eastern and Western Pygmies diverged, with evidence of subsequent migration only among the Western group and the Bantu-speaking farmers. Western Pygmies show signs of a recent bottleneck 4,000-650 ya, coincident with the diffusion of Bantu languages, whereas Eastern Pygmies seem to have experienced a more ancient decrease in population size (20,000-4,000 ya). In conclusion, the results of this first attempt at analyzing complete mtDNA sequences at the population level in sub-Saharan Africa not only support previous findings but also offer new insights into the demographic history of Pygmy populations, shedding new light on the ancient peopling of the African continent.

Signatures of the preagricultural peopling processes in sub-Saharan Africa as revealed by the phylogeography of early Y chromosome lineages.

The study of Y chromosome variation has helped reconstruct demographic events associated with the spread of languages, agriculture, and pastoralism in sub-Saharan Africa, but little attention has been given to the early history of the continent. In order to overcome this lack of knowledge, we carried out a phylogeographic analysis of haplogroups A and B in a broad data set of sub-Saharan populations. These two lineages are particularly suitable for this objective because they are the two most deeply rooted branches of the Y chromosome genealogy. Their distribution is almost exclusively restricted to sub-Saharan Africa where their frequency peaks at 65% in groups of foragers. The combined high-resolution single nucleotide polymorphism analysis with short tandem repeats variation of their subclades reveals strong geographic and population structure for both haplogroups. This has allowed us to identify specific lineages related to regional preagricultural dynamics in different areas of sub-Saharan Africa. In addition, we observed signatures of relatively recent contact, both among Pygmies and between them and Khoisan speaker groups from southern Africa, thus contributing to the understanding of the complex evolutionary relationships among African hunter-gatherers. Finally, by revising the phylogeography of the very early human Y chromosome lineages, we have obtained support for the role of southern Africa as a sink, rather than a source, of the first migrations of modern humans from eastern and central parts of the continent. These results open new perspectives on the early history of Homo sapiens in Africa, with particular attention to areas of the continent where human fossil remains and archaeological data are scant.

Evidence of high genetic variation among linguistically diverse populations on a micro-geographic scale: a case study of the Italian Alps.

Although essential for the fine-scale reconstruction of genetic structure, only a few micro-geographic studies have been carried out in European populations. This study analyzes mitochondrial variation (651 bp of the hypervariable region plus 17 single-nucleotide polymorphisms) in 393 samples from nine populations from Trentino (Eastern Italian Alps), a small area characterized by a complex geography and high linguistic diversity. A high level of genetic variation, comparable to geographically dispersed European groups, was observed. We found a difference in the intensity of peopling processes between two longitudinal areas, as populations from the west-central part of the region show stronger signatures of expansion, whereas those from the eastern area are closer to the expectations of a stationary demographic state. This may be explained by geomorphological factors and is also supported by archeological data. Finally, our results reveal a striking difference in the way in which the two linguistically isolated populations are genetically related to the neighboring groups. The Ladin speakers were found to be genetically close to the Italian-speaking populations and differentiated from the other Dolomitic Ladins, whereas the German-speaking Cimbri behave as an outlier, showing signatures of founder effects and low growth rate.

When gender matters: new insights into the relationships between social systems and the genetic structure of human populations.

Due to its important effects on the ecological dynamics and the genetic structure of species, biologists have long been interested in gender-biased dispersal, a condition where one gender is more prone to move from the natal site. More recently, this topic has attracted a great attention from human evolutionary geneticists. Considering the close relations between residential rules and social structure, gender-biased dispersal is, in fact, regarded as an important case study concerning the effects of socio-cultural factors on human genetic variation. It all started with the seminal paper by Mark Seielstad, Erich Minch and Luigi Luca Cavalli Sforza from Stanford University (Seielstad et al. 1998). They observed a larger differentiation for Y-chromosome than mitochondrial DNA between extant human populations, purportedly a consequence of the prevalence of long-term patrilocality in human societies. Subsequent studies, however, have highlighted the need to consider geographically close and culturally homogeneous groups, disentangle signals due to different peopling events and obtain unbiased estimates of genetic diversity. In this issue of Molecular Ecology, not only do Marks et al. (2012) adopt an experimental design which addresses these concerns, but they also take a further and important step forward by integrating the genetic analysis of two distant populations, the Basotho and Spanish, with data regarding migration rates and matrimonial distances. Using both empirical evidence and simulations, the authors show that female-biased migration due to patrilocality might shape the genetic structure of human populations only at short ranges and under substantial differences in migration rates between genders. Providing a quantitative framework for future studies of the effects of residential rules on the human genome, this study paves the way for further developments in the field. On a wider perspective, Marks et al.'s work demonstrates the power of approaches which integrate biological, cultural and demographic lines of evidence in the study of relations between social and genetic structures of human populations.

How can we get more open access to medical studies? Simple, let's take the green road.

Open Access (OA) to research publications is a fundamental resource for the advancement of scientific research. To facilitate the transition to OA, publishers and institutions have begun negotiating the so-called transformative agreements, contracts combining access to subscription journals with the ability to publish OA. While the debate on transformative agreements is very much alive, little attention is being paid to the "green road," the practice of openly self-archiving manuscripts that have been accepted by journals but not yet typeset. Here we focus on medical literature, showing how the green road could outperform transformative agreements as a means of increasing the full and free availability of peer-reviewed scientific papers.

From the Alps to the Mediterranean and beyond: genetics, environment, culture and the "impossible beauty" of Italy.

Since prehistoric times, Italy has represented a bridge between peoples, genes and cultures. Its peculiar geographical position explains why: it is located in the center of the Mediterranean Sea, flanked by the Balkans and the Hellenic Peninsula to the east, Iberia to the west and surrounded by North Africa to the south and central Europe to the north. This makes Italy of extraordinary interest for the study of some different aspects of human diversity. Here we overview current knowledge regarding the relationships between the structure of the genetic variation of Italian populations and the geographical, ecological and cultural factors that have characterized their evolutionary history. Human presence in Italian territory is deeply rooted in the past. Lithic artifacts produced by the genus Homo and remains of Homo sapiens are among the earliest to have been found on the continent, as shown by the lithic industry of Pirro Nord (between 1.3 and 1.6 Mya) and the dental remains of the "Grotta del Cavallo" (between 45 and 43 Kya). Genetic and genomic studies relating to existing and extinct human groups have shed light on the migrations from Europe, Africa and Asia that created the ancient layers of the genetic structure of today's Italian populations, especially before the Iron Age. The important role of isolation (genetic and cultural) in shaping genetic structure is clearly visible in the patterns of intra- and inter-population diversity observed among Italian ethno-linguistic minorities that settled on the peninsula and on the major islands until the 19th century. Finally, selective pressures have likely driven the distribution of originally adaptive variants and haplotypes that now confer protection or susceptibility to major diseases such as diabetes and cardiovascular disease (in northern Italy) and tuberculosis and leprosy (in the south). What emerges is a picture where the combined effects of migration, isolation and natural selection generated by the interplay of geography, environment and culture have shaped a complex pattern of human diversity that is unique in Europe and which goes hand in hand with today's rich animal and plant biodiversity. In a nutshell, scientific evidence and cultural heritage paint Italy as a place with extremely diverse environments where distant peoples have met since the deep past, bringing and sharing genes and ideas.

HLA class-I diversity in Cameroon: evidence for a north-south structure of genetic variation and relationships with African populations.

HLA class I diversity (loci A, B and C) was analysed in four populations, two from North Cameroon (Podokwo and Uldeme) and two from South Cameroon (Ewondo and Bamileke). Northern and southern Cameroon populations show a substantial genetic diversity in terms of haplotype sharing and genetic distances, even despite the low percentage of variance due to differences among populations evidenced by analysis of molecular variance. The signals of differentiation among populations are consistent with their linguistic affiliation, and support previous evidence, based on autosomal microsatellites and protein loci, which has shown that the complex pattern of genetic variation of Cameroon can in part be described by contrasting the northern and southern part of the country. Looking at our results in the more general framework of HLA diversity in sub-Saharan Africa, it turns out that the Podokwo and Uldeme show some genetic links to populations of the southern western branch of the Sahel corridor, while their high frequency of A*02 and C*04 alleles is congruent with previously hypothesised introgression of non-sub-Saharan alleles. On the other hand, signals of shared ancestry between the Bamileke and Ewondo and the Bantu speakers from central and southern Africa were detected.

The Bantu expansion revisited: a new analysis of Y chromosome variation in Central Western Africa.

The current distribution of Bantu languages is commonly considered to be a consequence of a relatively recent population expansion (3-5kya) in Central Western Africa. While there is a substantial consensus regarding the centre of origin of Bantu languages (the Benue River Valley, between South East Nigeria and Western Cameroon), the identification of the area from where the population expansion actually started, the relation between the processes leading to the spread of languages and peoples and the relevance of local migratory events remain controversial. In order to shed new light on these aspects, we studied Y chromosome variation in a broad dataset of populations encompassing Nigeria, Cameroon, Gabon and Congo. Our results evidence an evolutionary scenario which is more complex than had been previously thought, pointing to a marked differentiation of Cameroonian populations from the rest of the dataset. In fact, in contrast with the current view of Bantu speakers as a homogeneous group of populations, we observed an unexpectedly high level of interpopulation genetic heterogeneity and highlighted previously undetected diversity for lineages associated with the diffusion of Bantu languages (E1b1a (M2) sub-branches). We also detected substantial differences in local demographic histories, which concord with the hypotheses regarding an early diffusion of Bantu languages into the forest area and a subsequent demographic expansion and migration towards eastern and western Africa.

The emerging complexity of Open Science: assessing Intelligent Data Openness in Genomic Anthropology and Human Genomics.

In recent decades, the scientific community has become aware of the importance of science being effectively open in order to speed up scientific and technological progress. In this context, the achievement of a robust, effective and responsible form of data sharing is now widely acknowledged as a fundamental part of the research process. The production and resolution of human genomic data has steadily increased in recent years, mainly due to technological advances and decreasing costs of DNA genotyping and sequencing. There is, however, a downside to this process due to the huge increase in the complexity of the data and related metadata. This means it is advisable to go beyond traditional forms of sharing analysis, which have focused on data availability only. Here we present a pilot study that aims to complement a survey on the availability of data related to peer-reviewed publications with an analysis of their findability, accessibility, useability and assessability (according to the "intelligent data openness" scheme). Sharing rates in genomic anthropology (73.0%) were found to be higher than human genomics (32.4%), but lower than closely related research fields (from 96.8% to 79.2% for paleogenetics and evolutionary genetics, respectively). We discuss the privacy and methodological issues that could be linked to this finding. Comparisons of sharing rates across a wide range of disciplines has suggested that the idea of human genomics as a forerunner for the open data movement should be questioned. Finally, both in genomic anthropology and human genomics, findability and useability were found to be compliant with the expectations of an intelligent data openness, whereas only a minor part of studies met the need to make the data completely assessable.

Evolution of detoxifying systems: the role of environment and population history in shaping genetic diversity at human CYP2D6 locus.

OBJECTIVE: The transition from food collection to food production (FP) modified the nature of selective pressures, and several studies illustrate that genetic adaptation to new lifestyle has occurred in humans since the agricultural revolution. Here we test the hypothesis that high levels of genetic variation at CYP2D6, a locus coding for a detoxifying enzyme of the cytochrome P450 complex, reflect this change. METHODS: We compared DNA sequences and predicted the levels of enzyme activity across 10 African, Asian and European populations, six of which currently rely on hunting and gathering (HG) while four on food production (FP). RESULTS AND CONCLUSION: HG and FP showed similar levels of CYP2D6 diversity, but displayed different substitution patterns at coding DNA sites possibly related to selective differences. Comparison with variation at presumably neutral independent loci confirmed this finding, despite the confounding effects of population history, resulting in higher overall variation in Africans than in Eurasians. The differences between HG and FP populations suggest that new lifestyle and dietary habits acquired in the transition to agriculture affected the variation pattern at CYP2D6, leading to an increase in FP populations of the frequency of alleles that are associated with a slower rate of metabolism. These alleles reached a balanced co-existence with other important and previously selected variants. We suggest that the pronounced substrate-dependent activity of most of these enzymes expanded the spectrum of the metabolic response.