The application of molecular genetic approaches to the study of human evolution.

The past decade of advances in molecular genetic technology has heralded a new era for all evolutionary studies, but especially the science of human evolution. Data on various kinds of DNA variation in human populations have rapidly accumulated. There is increasing recognition of the importance of this variation for medicine and developmental biology and for understanding the history of our species. Haploid markers from mitochondrial DNA and the Y chromosome have proven invaluable for generating a standard model for evolution of modern humans. Conclusions from earlier research on protein polymorphisms have been generally supported by more sophisticated DNA analysis. Co-evolution of genes with language and some slowly evolving cultural traits, together with the genetic evolution of commensals and parasites that have accompanied modern humans in their expansion from Africa to the other continents, supports and supplements the standard model of genetic evolution. The advances in our understanding of the evolutionary history of humans attests to the advantages of multidisciplinary research.

A serial founder effect model for human settlement out of Africa.

The increasing abundance of human genetic data has shown that the geographical patterns of worldwide genetic diversity are best explained by human expansion out of Africa. This expansion is modelled well by prolonged migration from a single origin in Africa with multiple subsequent serial founding events. We discuss a new simulation model for the serial founder effect out of Africa and compare it with results from previous studies. Unlike previous models, we distinguish colonization events from the continued exchange of people between occupied territories as a result of mating. We conduct a search through parameter space to estimate the range of parameter values that best explain key statistics from published data on worldwide variation in microsatellites. The range of parameters we use is chosen to be compatible with an out-of-Africa migration at 50-60Kyr ago and archaeo-ethno-demographic information. In addition to a colonization rate of 0.09-0.18, for an acceptable fit to the published microsatellite data, incorporation into existing models of exchange between neighbouring populations is essential, but at a very low rate. A linear decay of genetic diversity with geographical distance from the origin of expansion could apply to any species, especially if it moved recently into new geographical niches.

Paleolithic Y-haplogroup heritage predominates in a Cretan highland plateau.

The island of Crete, credited by some historical scholars as a central crucible of western civilization, has been under continuous archeological investigation since the second half of the nineteenth century. In the present work, the geographic stratification of the contemporary Cretan Y-chromosome gene pool was assessed by high-resolution haplotyping to investigate the potential imprints of past colonization episodes and the population substructure. In addition to analyzing the possible geographic origins of Y-chromosome lineages in relatively accessible areas of the island, this study includes samples from the isolated interior of the Lasithi Plateau--a mountain plain located in eastern Crete. The potential significance of the results from the latter region is underscored by the possibility that this region was used as a Minoan refugium. Comparisons of Y-haplogroup frequencies among three Cretan populations as well as with published data from additional Mediterranean locations revealed significant differences in the frequency distributions of Y-chromosome haplogroups within the island. The most outstanding differences were observed in haplogroups J2 and R1, with the predominance of haplogroup R lineages in the Lasithi Plateau and of haplogroup J lineages in the more accessible regions of the island. Y-STR-based analyses demonstrated the close affinity that R1a1 chromosomes from the Lasithi Plateau shared with those from the Balkans, but not with those from lowland eastern Crete. In contrast, Cretan R1b microsatellite-defined haplotypes displayed more resemblance to those from Northeast Italy than to those from Turkey and the Balkans.

Y-chromosome short tandem repeat DYS458.2 non-consensus alleles occur independently in both binary haplogroups J1-M267 and R1b3-M405.

AIM: To determine the human Y-chromosome haplogroup backgrounds of non-consensus DYS458.2 short tandem repeat alleles and evaluate their phylogenetic substructure and frequency in representative samples from the Middle East, Europe, and Pakistan. METHODS: Molecular characterization of lineages was achieved using a combination of Y-chromosome haplogroup defining binary polymorphisms and up to 37 short tandem repeat loci, including DYS388 to construct haplotypes. DNA sequencing of the DYS458 locus and median-joining network analyses were used to evaluate Y-chromosome lineages displaying the DYS458.2 motif. RESULTS: We showed that the DYS458.2 allelic innovation arose independently on at least two distinctive binary haplogroup backgrounds and possibly a third as well. The partial allele length pattern was fixed in all haplogroup J1 chromosomes examined, including its known rare sub-haplogroups. Within the alternative R1b3 associated M405 defined sub-haplogroup, both DYS458.0 and DYS458.2 allele classes occurred. A single chromosome also allocated to the R1b3-M269*(xM405) classification. The physical position of the partial insertion/deletion occurrence within the normal tetramer tract differed distinctly in each haplogroup context. CONCLUSIONS: While unusual DYS458.2 alleles are informative, additional information for other linked polymorphic loci is required when using such non-conforming alleles to infer haplogroup background and common ancestry.

Accelerated evolution of the pituitary adenylate cyclase-activating polypeptide precursor gene during human origin.

Pituitary adenylate cyclase-activating polypeptide (PACAP) is a neuropeptide abundantly expressed in the central nervous system and involved in regulating neurogenesis and neuronal signal transduction. The amino acid sequence of PACAP is extremely conserved across vertebrate species, indicating a strong functional constraint during the course of evolution. However, through comparative sequence analysis, we demonstrated that the PACAP precursor gene underwent an accelerated evolution in the human lineage since the divergence from chimpanzees, and the amino acid substitution rate in humans is at least seven times faster than that in other mammal species resulting from strong Darwinian positive selection. Eleven human-specific amino acid changes were identified in the PACAP precursors, which are conserved from murine to African apes. Protein structural analysis suggested that a putative novel neuropeptide might have originated during human evolution and functioned in the human brain. Our data suggested that the PACAP precursor gene underwent adaptive changes during human origin and may have contributed to the formation of human cognition.

Adaptive evolution of MRGX2, a human sensory neuron specific gene involved in nociception.

MRGX2, a G-protein-coupled receptor, is specifically expressed in the sensory neurons of the human peripheral nervous system and involved in nociception. Here, we studied DNA polymorphism patterns and evolution of the MRGX2 gene in world-wide human populations and the representative nonhuman primate species. Our results demonstrated that MRGX2 had undergone adaptive changes in the path of human evolution, which were likely caused by Darwinian positive selection. The patterns of DNA sequence polymorphisms in human populations showed an excess of derived substitutions, which against the expectation of neutral evolution, implying that the adaptive evolution of MRGX2 in humans was a relatively recent event. The reconstructed secondary structure of the human MRGX2 revealed that three of the four human-specific amino acid substitutions were located in the extra-cellular domains. Such critical substitutions may alter the interactions between MRGX2 protein and its ligand, thus, potentially led to adaptive changes of the pain-perception-related nervous system during human evolution.

Neutrality condition and response law for nonlinear reaction-diffusion equations, with application to population genetics.

We study a general class of nonlinear macroscopic evolution equations with "transport" and "reaction" terms which describe the dynamics of a species of moving individuals (atoms, molecules, quasiparticles, organisms, etc.). We consider that two types of individuals exist, "not marked" and "marked," respectively. We assume that the concentrations of both types of individuals are measurable and that they obey a neutrality condition, that is, the kinetic and transport properties of the "not marked" and "marked" individuals are identical. We suggest a response experiment, which consists in varying the fraction of "marked" individuals with the preservation of total fluxes, and show that the response of the system can be represented by a linear superposition law even though the underlying dynamics of the system is in general highly nonlinear. The linear response law is valid even for large perturbations and is not the result of a linearization procedure but rather a necessary consequence of the neutrality condition. First, we apply the response theorem to chemical kinetics, where the "marked species" is a molecule labeled with a radioactive isotope and there is no kinetic isotope effect. The susceptibility function of the response law can be related to the reaction mechanism of the process. Secondly we study the geographical distribution of the nonrecurrent, nonreversible neutral mutations of the nonrecombining portion of the Y chromosome from human populations and show that the fraction of mutants at a given point in space and time obeys a linear response law of the type introduced in this paper. The theory may be used for evaluating the geographic position and the moment in time where and when a mutation originated.

Y-chromosome distribution within the geo-linguistic landscape of northwestern Russia.

Populations of northeastern Europe and the Uralic mountain range are found in close geographic proximity, but they have been subject to different demographic histories. The current study attempts to better understand the genetic paternal relationships of ethnic groups residing in these regions. We have performed high-resolution haplotyping of 236 Y-chromosomes from populations in northwestern Russia and the Uralic mountains, and compared them to relevant previously published data. Haplotype variation and age estimation analyses using 15 Y-STR loci were conducted for samples within the N1b, N1c1 and R1a1 single-nucleotide polymorphism backgrounds. Our results suggest that although most genetic relationships throughout Eurasia are dependent on geographic proximity, members of the Uralic and Slavic linguistic families and subfamilies, yield significant correlations at both levels of comparison making it difficult to denote either linguistics or geographic proximity as the basis for their genetic substrata. Expansion times for haplogroup R1a1 date approximately to 18,000 YBP, and age estimates along with Network topology of populations found at opposite poles of its range (Eastern Europe and South Asia) indicate that two separate haplotypic foci exist within this haplogroup. Data based on haplogroup N1b challenge earlier findings and suggest that the mutation may have occurred in the Uralic range rather than in Siberia and much earlier than has been proposed (12.9+/-4.1 instead of 5.2+/-2.7 kya). In addition, age and variance estimates for haplogroup N1c1 suggest that populations from the western Urals may have been genetically influenced by a dispersal from northeastern Europe (eg, eastern Slavs) rather than the converse.

Mitochondrial DNA variation of modern Tuscans supports the near eastern origin of Etruscans.

The origin of the Etruscan people has been a source of major controversy for the past 2,500 years, and several hypotheses have been proposed to explain their language and sophisticated culture, including an Aegean/Anatolian origin. To address this issue, we analyzed the mitochondrial DNA (mtDNA) of 322 subjects from three well-defined areas of Tuscany and compared their sequence variation with that of 55 western Eurasian populations. Interpopulation comparisons reveal that the modern population of Murlo, a small town of Etruscan origin, is characterized by an unusually high frequency (17.5%) of Near Eastern mtDNA haplogroups. Each of these haplogroups is represented by different haplotypes, thus dismissing the possibility that the genetic allocation of the Murlo people is due to drift. Other Tuscan populations do not show the same striking feature; however, overall, ~5% of mtDNA haplotypes in Tuscany are shared exclusively between Tuscans and Near Easterners and occupy terminal positions in the phylogeny. These findings support a direct and rather recent genetic input from the Near East--a scenario in agreement with the Lydian origin of Etruscans. Such a genetic contribution has been extensively diluted by admixture, but it appears that there are still locations in Tuscany, such as Murlo, where traces of its arrival are easily detectable.

Polarity and temporality of high-resolution y-chromosome distributions in India identify both indigenous and exogenous expansions and reveal minor genetic influence of Central Asian pastoralists.

Although considerable cultural impact on social hierarchy and language in South Asia is attributable to the arrival of nomadic Central Asian pastoralists, genetic data (mitochondrial and Y chromosomal) have yielded dramatically conflicting inferences on the genetic origins of tribes and castes of South Asia. We sought to resolve this conflict, using high-resolution data on 69 informative Y-chromosome binary markers and 10 microsatellite markers from a large set of geographically, socially, and linguistically representative ethnic groups of South Asia. We found that the influence of Central Asia on the pre-existing gene pool was minor. The ages of accumulated microsatellite variation in the majority of Indian haplogroups exceed 10,000-15,000 years, which attests to the antiquity of regional differentiation. Therefore, our data do not support models that invoke a pronounced recent genetic input from Central Asia to explain the observed genetic variation in South Asia. R1a1 and R2 haplogroups indicate demographic complexity that is inconsistent with a recent single history. Associated microsatellite analyses of the high-frequency R1a1 haplogroup chromosomes indicate independent recent histories of the Indus Valley and the peninsular Indian region. Our data are also more consistent with a peninsular origin of Dravidian speakers than a source with proximity to the Indus and with significant genetic input resulting from demic diffusion associated with agriculture. Our results underscore the importance of marker ascertainment for distinguishing phylogenetic terminal branches from basal nodes when attributing ancestral composition and temporality to either indigenous or exogenous sources. Our reappraisal indicates that pre-Holocene and Holocene-era--not Indo-European--expansions have shaped the distinctive South Asian Y-chromosome landscape.

The Human Genome Diversity Project: past, present and future.

The Human Genome Project, in accomplishing its goal of sequencing one human genome, heralded a new era of research, a component of which is the systematic study of human genetic variation. Despite delays, the Human Genome Diversity Project has started to make progress in understanding the patterns of this variation and its causes, and also promises to provide important information for biomedical studies.

Support from the relationship of genetic and geographic distance in human populations for a serial founder effect originating in Africa.

Equilibrium models of isolation by distance predict an increase in genetic differentiation with geographic distance. Here we find a linear relationship between genetic and geographic distance in a worldwide sample of human populations, with major deviations from the fitted line explicable by admixture or extreme isolation. A close relationship is shown to exist between the correlation of geographic distance and genetic differentiation (as measured by F(ST)) and the geographic pattern of heterozygosity across populations. Considering a worldwide set of geographic locations as possible sources of the human expansion, we find that heterozygosities in the globally distributed populations of the data set are best explained by an expansion originating in Africa and that no geographic origin outside of Africa accounts as well for the observed patterns of genetic diversity. Although the relationship between F(ST) and geographic distance has been interpreted in the past as the result of an equilibrium model of drift and dispersal, simulation shows that the geographic pattern of heterozygosities in this data set is consistent with a model of a serial founder effect starting at a single origin. Given this serial-founder scenario, the relationship between genetic and geographic distance allows us to derive bounds for the effects of drift and natural selection on human genetic variation.

Fisher's theorems for multivariable, time- and space-dependent systems, with applications in population genetics and chemical kinetics.

We study different physical, chemical, or biological processes involving replication, transformation, and disappearance processes, as well as transport processes, and assume that the time and space dependence of the species densities are known. We derive two types of Fisher equations. The first type relates the average value of the time derivative of the relative time-specific rates of growth of the different species to the variance of the relative, time-specific rates of growth. A second type relates the average value of the gradient or the divergence of the relative, space-specific rates of growth to the space correlation matrix of the relative, space-specific rates of growth. These Fisher equations are exact results, which are independent of the detailed kinetics of the process: they are valid whether the evolution equations are linear or nonlinear, local or nonlocal in space and/or time and can be applied for the study of a large class of physical, chemical, and biological systems described in terms of time- and/or space-dependent density fields. We examine the implications of our generalized Fisher relations in population genetics, biochemistry, and chemical kinetics (reaction-diffusion systems). We show that there is a connection between the enhanced (hydrodynamic) transport of mutations induced by population growth and space-specific rate vectors: the velocity of enhanced transport is proportional to the product of the diffusion coefficient of the species and the space rate vector; this relation is similar to a fluctuation-dissipation relation in statistical mechanics.

Enhanced (hydrodynamic) transport induced by population growth in reaction-diffusion systems with application to population genetics.

We consider a system made up of different physical, chemical, or biological species undergoing replication, transformation, and disappearance processes, as well as slow diffusive motion. We show that for systems with net growth the balance between kinetics and the diffusion process may lead to fast, enhanced hydrodynamic transport. Solitary waves in the system, if they exist, stabilize the enhanced transport, leading to constant transport speeds. We apply our theory to the problem of determining the original mutation position from the current geographic distribution of a given mutation. We show that our theory is in good agreement with a simulation study of the mutation problem presented in the literature. It is possible to evaluate migratory trajectories from measured data related to the current distribution of mutations in human populations.

Mutations arising in the wave front of an expanding population.

The ability to infer the time and place of origin of a mutation can be very useful when reconstructing the evolutionary histories of populations and species. We use forward computer simulations of population growth, migration, and mutation in an analysis of an expanding population with a wave front that advances at a constant slow rate. A pronounced founder effect can be observed among mutations arising in this wave front where extreme population bottlenecks arise and are followed by major population growth. A fraction of mutations travel with the wave front and generate mutant populations that are on average much larger than those that remain stationary. Analysis of the diffusion of these mutants makes it possible to reconstruct migratory trajectories during population expansions, thus helping us better understand observed patterns in the evolution of species such as modern humans. Examination of some historical data supports our model.

Ethiopians and Khoisan share the deepest clades of the human Y-chromosome phylogeny.

The genetic structure of 126 Ethiopian and 139 Senegalese Y chromosomes was investigated by a hierarchical analysis of 30 diagnostic biallelic markers selected from the worldwide Y-chromosome genealogy. The present study reveals that (1) only the Ethiopians share with the Khoisan the deepest human Y-chromosome clades (the African-specific Groups I and II) but with a repertoire of very different haplotypes; (2) most of the Ethiopians and virtually all the Senegalese belong to Group III, whose precursor is believed to be involved in the first migration out of Africa; and (3) the Ethiopian Y chromosomes that fall into Groups VI, VIII, and IX may be explained by back migrations from Asia. The first observation confirms the ancestral affinity between the Ethiopians and the Khoisan, which has previously been suggested by both archaeological and genetic findings.

Origin, diffusion, and differentiation of Y-chromosome haplogroups E and J: inferences on the neolithization of Europe and later migratory events in the Mediterranean area.

The phylogeography of Y-chromosome haplogroups E (Hg E) and J (Hg J) was investigated in >2400 subjects from 29 populations, mainly from Europe and the Mediterranean area but also from Africa and Asia. The observed 501 Hg E and 445 Hg J samples were subtyped using 36 binary markers and eight microsatellite loci. Spatial patterns reveal that (1). the two sister clades, J-M267 and J-M172, are distributed differentially within the Near East, North Africa, and Europe; (2). J-M267 was spread by two temporally distinct migratory episodes, the most recent one probably associated with the diffusion of Arab people; (3). E-M81 is typical of Berbers, and its presence in Iberia and Sicily is due to recent gene flow from North Africa; (4). J-M172(xM12) distribution is consistent with a Levantine/Anatolian dispersal route to southeastern Europe and may reflect the spread of Anatolian farmers; and (5). E-M78 (for which microsatellite data suggest an eastern African origin) and, to a lesser extent, J-M12(M102) lineages would trace the subsequent diffusion of people from the southern Balkans to the west. A 7%-22% contribution of Y chromosomes from Greece to southern Italy was estimated by admixture analysis.

Excavating Y-chromosome haplotype strata in Anatolia.

Analysis of 89 biallelic polymorphisms in 523 Turkish Y chromosomes revealed 52 distinct haplotypes with considerable haplogroup substructure, as exemplified by their respective levels of accumulated diversity at ten short tandem repeat (STR) loci. The major components (haplogroups E3b, G, J, I, L, N, K2, and R1; 94.1%) are shared with European and neighboring Near Eastern populations and contrast with only a minor share of haplogroups related to Central Asian (C, Q and O; 3.4%), Indian (H, R2; 1.5%) and African (A, E3*, E3a; 1%) affinity. The expansion times for 20 haplogroup assemblages was estimated from associated STR diversity. This comprehensive characterization of Y-chromosome heritage addresses many multifaceted aspects of Anatolian prehistory, including: (1) the most frequent haplogroup, J, splits into two sub-clades, one of which (J2) shows decreasing variances with increasing latitude, compatible with a northward expansion; (2) haplogroups G1 and L show affinities with south Caucasus populations in their geographic distribution as well as STR motifs; (3) frequency of haplogroup I, which originated in Europe, declines with increasing longitude, indicating gene flow arriving from Europe; (4) conversely, haplogroup G2 radiates towards Europe; (5) haplogroup E3b3 displays a latitudinal correlation with decreasing frequency northward; (6) haplogroup R1b3 emanates from Turkey towards Southeast Europe and Caucasia and; (7) high resolution SNP analysis provides evidence of a detectable yet weak signal (<9%) of recent paternal gene flow from Central Asia. The variety of Turkish haplotypes is witness to Turkey being both an important source and recipient of gene flow.

A back migration from Asia to sub-Saharan Africa is supported by high-resolution analysis of human Y-chromosome haplotypes.

The variation of 77 biallelic sites located in the nonrecombining portion of the Y chromosome was examined in 608 male subjects from 22 African populations. This survey revealed a total of 37 binary haplotypes, which were combined with microsatellite polymorphism data to evaluate internal diversities and to estimate coalescence ages of the binary haplotypes. The majority of binary haplotypes showed a nonuniform distribution across the continent. Analysis of molecular variance detected a high level of interpopulation diversity (PhiST=0.342), which appears to be partially related to the geography (PhiCT=0.230). In sub-Saharan Africa, the recent spread of a set of haplotypes partially erased pre-existing diversity, but a high level of population (PhiST=0.332) and geographic (PhiCT=0.179) structuring persists. Correspondence analysis shows that three main clusters of populations can be identified: northern, eastern, and sub-Saharan Africans. Among the latter, the Khoisan, the Pygmies, and the northern Cameroonians are clearly distinct from a tight cluster formed by the Niger-Congo-speaking populations from western, central western, and southern Africa. Phylogeographic analyses suggest that a large component of the present Khoisan gene pool is eastern African in origin and that Asia was the source of a back migration to sub-Saharan Africa. Haplogroup IX Y chromosomes appear to have been involved in such a migration, the traces of which can now be observed mostly in northern Cameroon.