A geographically explicit genetic model of worldwide human-settlement history

slendr: a framework for spatio-temporal population genomic simulations on geographic landscapes

One of the goals of population genetics is to understand how evolutionary forces shape patterns of genetic variation over time. However, because populations evolve across both time and space, most evolutionary processes also have an important spatial component, acting through phenomena such as isolation by distance, local mate choice, or uneven distribution of resources. This spatial dimension is often neglected, partly due to the lack of tools specifically designed for building and evaluating complex spatio-temporal population genetic models. To address this methodological gap, we present a new framework for simulating spatially-explicit genomic data, implemented in a new R package called slendr (www.slendr.net), which leverages a SLiM simulation back-end script bundled with the package. With this framework, the users can programmatically and visually encode spatial population ranges and their temporal dynamics (i.e., population displacements, expansions, and contractions) either on real Earth landscapes or on abstract custom maps, and schedule splits and gene-flow events between populations using a straightforward declarative language. Additionally, slendr can simulate data from traditional, non-spatial models, either with SLiM or using an alternative built-in coalescent msprime back end. Together with its R-idiomatic interface to the tskit library for tree-sequence processing and analysis, slendr opens up the possibility of performing efficient, reproducible simulations of spatio-temporal genomic data entirely within the R environment, leveraging its wealth of libraries for geospatial data analysis, statistics, and visualization. Here, we present the design of the slendr R package and demonstrate its features on several practical example workflows.

Historic Sampling of a Vanishing Beast: Population Structure and Diversity in the Black Rhinoceros

The black rhinoceros (Diceros bicornis L.) is a critically endangered species historically distributed across sub-Saharan Africa. Hunting and habitat disturbance have diminished both its numbers and distribution since the 19th century, but a poaching crisis in the late 20th century drove them to the brink of extinction. Genetic and genomic assessments can greatly increase our knowledge of the species and inform management strategies. However, when a species has been severely reduced, with the extirpation and artificial admixture of several populations, it is extremely challenging to obtain an accurate understanding of historic population structure and evolutionary history from extant samples. Therefore, we generated and analyzed whole genomes from 63 black rhinoceros museum specimens collected between 1775 and 1981. Results showed that the black rhinoceros could be genetically structured into six major historic populations (Central Africa, East Africa, Northwestern Africa, Northeastern Africa, Ruvuma, and Southern Africa) within which were nested four further subpopulations (Maasailand, southwestern, eastern rift, and northern rift), largely mirroring geography, with a punctuated north-south cline. However, we detected varying degrees of admixture among groups and found that several geographical barriers, most prominently the Zambezi River, drove population discontinuities. Genomic diversity was high in the middle of the range and decayed toward the periphery. This comprehensive historic portrait also allowed us to ascertain the ancestry of 20 resequenced genomes from extant populations. Lastly, using insights gained from this unique temporal data set, we suggest management strategies, some of which require urgent implementation, for the conservation of the remaining black rhinoceros diversity.

A Chinese indicine pangenome reveals a wealth of novel structural variants introgressed from other Bos species

Chinese indicine cattle harbor a much higher genetic diversity compared with other domestic cattle, but their genome architecture remains uninvestigated. Using PacBio HiFi sequencing data from 10 Chinese indicine cattle across southern China, we assembled 20 high-quality partially phased genomes and integrated them into a multiassembly graph containing 148.5 Mb (5.6%) of novel sequence. We identified 156,009 high-confidence nonredundant structural variants (SVs) and 206 SV hotspots spanning ∼195 Mb of gene-rich sequence. We detected 34,249 archaic introgressed fragments in Chinese indicine cattle covering 1.93 Gb (73.3%) of the genome. We inferred an average of 3.8%, 3.2%, 1.4%, and 0.5% of introgressed sequence originating, respectively, from banteng-like, kouprey-like, gayal-like, and gaur-like Bos species, as well as 0.6% of unknown origin. Introgression from multiple donors might have contributed to the genetic diversity of Chinese indicine cattle. Altogether, this study highlights the contribution of interspecies introgression to the genomic architecture of an important livestock population and shows how exotic genomic elements can contribute to the genetic variation available for selection.

Ultra-deep sequencing of Hadza hunter-gatherers recovers vanishing gut microbes

The gut microbiome modulates immune and metabolic health. Human microbiome data are biased toward industrialized populations, limiting our understanding of non-industrialized microbiomes. Here, we performed ultra-deep metagenomic sequencing on 351 fecal samples from the Hadza hunter-gatherers of Tanzania and comparative populations in Nepal and California. We recovered 91,662 genomes of bacteria, archaea, bacteriophages, and eukaryotes, 44% of which are absent from existing unified datasets. We identified 124 gut-resident species vanishing in industrialized populations and highlighted distinct aspects of the Hadza gut microbiome related to in situ replication rates, signatures of selection, and strain sharing. Industrialized gut microbes were found to be enriched in genes associated with oxidative stress, possibly a result of microbiome adaptation to inflammatory processes. This unparalleled view of the Hadza gut microbiome provides a valuable resource, expands our understanding of microbes capable of colonizing the human gut, and clarifies the extensive perturbation induced by the industrialized lifestyle.

Morphology, evolution, and the whole organism imperative: Why evolutionary questions need multi-trait evolutionary quantitative genetics

Since Washburn's New Physical Anthropology, researchers have sought to understand the complexities of morphological evolution among anatomical regions in human and non-human primates. Researchers continue, however, to preferentially use comparative and functional approaches to examine complex traits, but these methods cannot address questions about evolutionary process and often conflate function with fitness. Moreover, researchers also tend to examine anatomical elements in isolation, which implicitly assumes independent evolution among different body regions. In this paper, we argue that questions asked in primate evolution are best examined using multiple anatomical regions subjected to model-bound methods built from an understanding of evolutionary quantitative genetics. A nascent but expanding number of studies over the last two decades use this approach, examining morphological integration, evolvability, and selection modeling. To help readers learn how to use these methods, we review fundamentals of evolutionary processes within a quantitative genetic framework, explore the importance of neutral evolutionary theory, and explain the basics of evolutionary quantitative genetics, namely the calculation of evolutionary potential for multiple traits in response to selection. Leveraging these methods, we demonstrate their use to understand non-independence in possible evolutionary responses across the limbs, limb girdles, and basicranium of humans. Our results show that model-bound quantitative genetic methods can reveal unexpected genetic covariances among traits that create a novel but measurable understanding of evolutionary complexity among multiple traits. We advocate for evolutionary quantitative genetic methods to be a standard whenever appropriate to keep studies of primate morphological evolution relevant for the next seventy years and beyond.

Patterns of variation in canal and root number in human post-canine teeth

Descriptive morphology of tooth roots traditionally focuses on number of canals and roots. However, how or if canal and root number are related is poorly understood. While it is often assumed that canal number is concomitant with root number and morphology, in practice canal number and morphology do not always covary with external root features. To investigate the relationship between canal and root number, fully developed, adult post‐canine teeth were examined and quantified from computerized tomography scans from a global sample of 945 modern humans. We tested the hypotheses that root and canal number do not follow a 1:1 ratio, that canal to root ratios differ between teeth, and that canal to root ratios differ across major human geographical groups. Results indicate that not only is root number dependent on canal number, but that this relationship becomes more variable as canal number increases, varies between individual teeth and by major geographical group, and changes as these groups increase in geographical distance from Sub‐Saharan Africa. These results show that the ratio of canal number to root number is an important indicator of variation in dental phenotypes.

Global patterns of the cranial form of modern human populations described by analysis of a 3D surface homologous model

This study assessed the regional diversity of the human cranial form by using geometric homologous models based on scanned data from 148 ethnic groups worldwide. This method adopted a template-fitting technique for a nonrigid transformation via the iterative closest point algorithm to generate the homologous meshes. Through the application of principal component analysis to 342 sampled homologous models, the largest variation was detected in overall size, and small South Asian crania were clearly verified. The next greatest diversity was found in the length/breadth proportion of the neurocranium, which showed the contrast between the elongated crania of Africans and the globular crania of Northeast Asians. Notably, this component was slightly correlated with the facial profile. Well-known facial features, such as the forward projection of the cheek among Northeast Asians and compaction of the European maxilla, were reconfirmed. These facial variations were highly correlated with the calvarial outline, particularly the degree of frontal and occipital inclines. An allometric pattern was detected in facial proportions in relation to overall cranial size; in larger crania, the facial profiles tend to be longer and narrower, as demonstrated among many American natives and Northeast Asians. Although our study did not include data on environmental variables that are likely to affect cranial morphology, such as climate or dietary conditions, the large datasets of homologous cranial models will be usefully available for seeking various attributions to phenotypic skeletal characteristics.

Pan-Africanism vs. single-origin of Homo sapiens: Putting the debate in the light of evolutionary biology

The scenario of Homo sapiens origin/s within Africa has become increasingly complex, with a pan-African perspective currently challenging the long-established single-origin hypothesis. In this paper, we review the lines of evidence employed in support of each model, highlighting inferential limitations and possible terminological misunderstandings. We argue that the metapopulation scenario envisaged by pan-African proponents well describes a mosaic diversification among late Middle Pleistocene groups. However, this does not rule out a major contribution that emerged from a single population where crucial derived features-notably, a globular braincase-appeared as the result of a punctuated, cladogenetic event. Thus, we suggest that a synthesis is possible and propose a scenario that, in our view, better reconciles with consolidated expectations in evolutionary theory. These indicate cladogenesis in allopatry as an ordinary pattern for the origin of a new species, particularly during phases of marked climatic and environmental instability.

North Asian population relationships in a global context

Population genetic studies of North Asian ethnic groups have focused on genetic variation of sex chromosomes and mitochondria. Studies of the extensive variation available from autosomal variation have appeared infrequently. We focus on relationships among population samples using new North Asia microhaplotype data. We combined genotypes from our laboratory on 58 microhaplotypes, distributed across 18 autosomes, on 3945 individuals from 75 populations with corresponding data extracted for 26 populations from the Thousand Genomes consortium and for 22 populations from the GenomeAsia 100 K project. A total of 7107 individuals in 122 total populations are analyzed using STRUCTURE, Principal Component Analysis, and phylogenetic tree analyses. North Asia populations sampled in Mongolia include: Buryats, Mongolians, Altai Kazakhs, and Tsaatans. Available Siberians include samples of Yakut, Khanty, and Komi Zyriane. Analyses of all 122 populations confirm many known relationships and show that most populations from North Asia form a cluster distinct from all other groups. Refinement of analyses on smaller subsets of populations reinforces the distinctiveness of North Asia and shows that the North Asia cluster identifies a region that is ancestral to Native Americans.

Ultra-deep Sequencing of Hadza Hunter-Gatherers Recovers Vanishing Microbes.. [PMID: 36238714 PMCID: PMC9558438 DOI: 10.1101/2022.03.30.486478]

The gut microbiome is a key modulator of immune and metabolic health. Human microbiome data is biased towards industrialized populations, providing limited understanding of the distinct and diverse non-industrialized microbiomes. Here, we performed ultra-deep metagenomic sequencing and strain cultivation on 351 fecal samples from the Hadza, hunter-gatherers in Tanzania, and comparative populations in Nepal and California. We recover 94,971 total genomes of bacteria, archaea, bacteriophages, and eukaryotes, 43% of which are absent from existing unified datasets. Analysis of in situ growth rates, genetic pN/pS signatures, high-resolution strain tracking, and 124 gut-resident species vanishing in industrialized populations reveals differentiating dynamics of the Hadza gut microbiome. Industrialized gut microbes are enriched in genes associated with oxidative stress, possibly a result of microbiome adaptation to inflammatory processes. This unparalleled view of the Hadza gut microbiome provides a valuable resource that expands our understanding of microbes capable of colonizing the human gut and clarifies the extensive perturbation brought on by the industrialized lifestyle.

Female craniometrics support the 'two-layer model' of human dispersal in Eastern Eurasia

This study reports a cranio-morphometric analysis of female human remains from seven archaeological sites in China, Vietnam and Taiwan that date between 16,000 and 5300 BP. The aim of the analysis is to test the “two-layer” model of human dispersal in eastern Eurasia, using previously unanalysed female remains to balance the large sample of previously-analysed males. The resulting craniometric data indicate that the examined specimens all belong to the “first layer” of dispersal, and share a common ancestor with recent Australian and Papuan populations, and the ancient Jomon people of Japan. The analysed specimens pre-date the expansion of agricultural populations of East/Northeast Asian origin—that is, the “second layer” of human dispersal proposed by the model. As a result of this study, the two-layer model, which has hitherto rested on evidence only from male skeletons, is now strongly supported by female-derived data. Further comparisons reveal that the people of the first layer were closer in terms of their facial morphology to modern Africans and Sri Lankan Veddah than to modern Asians and Europeans, suggesting that the Late Pleistocene through Middle Holocene hunter-gatherers examined in this study were direct descendants of the anatomically modern humans who first migrated out of Africa through southern Eurasia.

Understanding population structure in an evolutionary context: population-specific FST and pairwise FST

Populations are shaped by their history. It is crucial to interpret population structure in an evolutionary context. Pairwise F_ST measures population structure, whereas population-specific F_ST measures deviation from the ancestral population. To understand the current population structure and a population’s history of range expansion, we propose a representation method that overlays population-specific F_ST estimates on a sampling location map, and on an unrooted neighbor-joining tree and a multi-dimensional scaling plot inferred from a pairwise F_ST distance matrix. We examined the usefulness of our procedure using simulations that mimicked population colonization from an ancestral population and by analyzing published human, Atlantic cod, and wild poplar data. Our results demonstrated that population-specific F_ST values identify the source population and trace the evolutionary history of its derived populations. Conversely, pairwise F_ST values represent the current population structure. By integrating the results of both estimators, we obtained a new picture of the population structure that incorporates evolutionary history. The generalized least squares estimate of genome-wide population-specific F_ST indicated that the wild poplar population expanded its distribution to the north, where daylight hours are long in summer, to coastal areas with abundant rainfall, and to the south where summers are dry. Genomic data highlight the power of the bias-corrected moment estimators of F_ST, whether global, pairwise, or population-specific, that provide unbiased estimates of F_ST. All F_ST moment estimators described in this paper have reasonable processing times and are useful in population genomics studies.

Ancient DNA and multimethod dating confirm the late arrival of anatomically modern humans in southern China

The expansion of anatomically modern humans (AMHs) from Africa around 65,000 to 45,000 y ago (ca. 65 to 45 ka) led to the establishment of present-day non-African populations. Some paleoanthropologists have argued that fossil discoveries from Huanglong, Zhiren, Luna, and Fuyan caves in southern China indicate one or more prior dispersals, perhaps as early as ca. 120 ka. We investigated the age of the human remains from three of these localities and two additional early AMH sites (Yangjiapo and Sanyou caves, Hubei) by combining ancient DNA (aDNA) analysis with a multimethod geological dating strategy. Although U-Th dating of capping flowstones suggested they lie within the range ca. 168 to 70 ka, analyses of aDNA and direct AMS ¹⁴C dating on human teeth from Fuyan and Yangjiapo caves showed they derive from the Holocene. OSL dating of sediments and AMS ¹⁴C analysis of mammal teeth and charcoal also demonstrated major discrepancies from the flowstone ages; the difference between them being an order of magnitude or more at most of these localities. Our work highlights the surprisingly complex depositional history recorded at these subtropical caves which involved one or more episodes of erosion and redeposition or intrusion as recently as the late Holocene. In light of our findings, the first appearance datum for AMHs in southern China should probably lie within the timeframe set by molecular data of ca. 50 to 45 ka.

On the Unfounded Enthusiasm for Soft Selective Sweeps III: The Supervised Machine Learning Algorithm That Isn't

In the last 15 years or so, soft selective sweep mechanisms have been catapulted from a curiosity of little evolutionary importance to a ubiquitous mechanism claimed to explain most adaptive evolution and, in some cases, most evolution. This transformation was aided by a series of articles by Daniel Schrider and Andrew Kern. Within this series, a paper entitled “Soft sweeps are the dominant mode of adaptation in the human genome” (Schrider and Kern, Mol. Biol. Evolut. 2017, 34(8), 1863–1877) attracted a great deal of attention, in particular in conjunction with another paper (Kern and Hahn, Mol. Biol. Evolut. 2018, 35(6), 1366–1371), for purporting to discredit the Neutral Theory of Molecular Evolution (Kimura 1968). Here, we address an alleged novelty in Schrider and Kern’s paper, i.e., the claim that their study involved an artificial intelligence technique called supervised machine learning (SML). SML is predicated upon the existence of a training dataset in which the correspondence between the input and output is known empirically to be true. Curiously, Schrider and Kern did not possess a training dataset of genomic segments known a priori to have evolved either neutrally or through soft or hard selective sweeps. Thus, their claim of using SML is thoroughly and utterly misleading. In the absence of legitimate training datasets, Schrider and Kern used: (1) simulations that employ many manipulatable variables and (2) a system of data cherry-picking rivaling the worst excesses in the literature. These two factors, in addition to the lack of negative controls and the irreproducibility of their results due to incomplete methodological detail, lead us to conclude that all evolutionary inferences derived from so-called SML algorithms (e.g., S/HIC) should be taken with a huge shovel of salt.

A Revised Model of Anatomically Modern Human Expansions Out of Africa through a Machine Learning Approximate Bayesian Computation Approach

There is a wide consensus in considering Africa as the birthplace of anatomically modern humans (AMH), but the dispersal pattern and the main routes followed by our ancestors to colonize the world are still matters of debate. It is still an open question whether AMH left Africa through a single process, dispersing almost simultaneously over Asia and Europe, or in two main waves, first through the Arab Peninsula into southern Asia and Australo-Melanesia, and later through a northern route crossing the Levant. The development of new methodologies for inferring population history and the availability of worldwide high-coverage whole-genome sequences did not resolve this debate. In this work, we test the two main out-of-Africa hypotheses through an Approximate Bayesian Computation approach, based on the Random-Forest algorithm. We evaluated the ability of the method to discriminate between the alternative models of AMH out-of-Africa, using simulated data. Once assessed that the models are distinguishable, we compared simulated data with real genomic variation, from modern and archaic populations. This analysis showed that a model of multiple dispersals is four-fold as likely as the alternative single-dispersal model. According to our estimates, the two dispersal processes may be placed, respectively, around 74,000 and around 46,000 years ago.

Sometimes hidden but always there: the assumptions underlying genetic inference of demographic histories

Demographic processes directly affect patterns of genetic variation within contemporary populations as well as future generations, allowing for demographic inference from patterns of both present-day and past genetic variation. Advances in laboratory procedures, sequencing and genotyping technologies in the past decades have resulted in massive increases in high-quality genome-wide genetic data from present-day populations and allowed retrieval of genetic data from archaeological material, also known as ancient DNA. This has resulted in an explosion of work exploring past changes in population size, structure, continuity and movement. However, as genetic processes are highly stochastic, patterns of genetic variation only indirectly reflect demographic histories. As a result, past demographic processes need to be reconstructed using an inferential approach. This usually involves comparing observed patterns of variation with model expectations from theoretical population genetics. A large number of approaches have been developed based on different population genetic models that each come with assumptions about the data and underlying demography. In this article I review some of the key models and assumptions underlying the most commonly used approaches for past demographic inference and their consequences for our ability to link the inferred demographic processes to the archaeological and climate records. This article is part of the theme issue 'Cross-disciplinary approaches to prehistoric demography'.

Evolvability in human postcranial traits across ecogeographic regions

OBJECTIVES

Though recent quantitative genetic analyses have indicated that directional selection appears to be acting on limb lengths and measures of body size in modern humans, these studies assume equal evolvability across modern human groups. However, differences in trait covariance structure due to ancient migration patterns and/or selection may limit the evolvability of populations further from Africa. This study therefore explores patterns of human evolvability across ecogeographic regions.

MATERIALS AND METHODS

Mean evolvability, respondability, conditional evolvability, and autonomy were calculated from variance-covariance matrices of limb length and body size measures representing 14 human groups spanning four ecogeographic regions. Measures of evolvability were compared across groups and regions, and the minimum sample size, inaccuracy, and bias were calculated for each.

RESULTS

When compared between regions, humans demonstrate significant differences between indices of evolvability across regions. Despite the relatively recent evolution of modern humans, several measures of evolvability show a strong negative correlation with latitude across regions, demonstrating a reduction in genetic variance that is potentially reflective of human migration and/or response to selection.

CONCLUSIONS

These results demonstrate the importance of establishing patterns of evolvability prior to additional quantitative genetic analyses, and emphasize the influence of sample size on the accuracy of estimated evolvability measures. These findings also suggest that while modern human groups share similar covariance structures, there is evidence for emergent differentiation in evolvability and respondability between human groups across ecogeographic regions, further complicating our ability to apply results derived from modern human groups to ancient hominin lineages.

Counting Oceanians of Non-European, Non-Asian Descent (ONENA) in the South Pacific to Make Them Count in Global Health

Several diseases and vulnerabilities associated with genetic or microbial factors are more frequent among populations of Oceanian, Non-European, Non-Asian descent (ONENA). ONENA are specific and have long been isolated geographically. To our knowledge, there are no published official, quantitative, aggregated data on the populations impacted by these excess vulnerabilities in Oceania. We searched official census reports for updated estimates of the total population for each of the Pacific Island Countries and Territories (including Australia) and the US State of Hawaii, privileging local official statistical or censual sources. We multiplied the most recent total population estimate by the cumulative percentage of the ONENA population as determined in official reports. Including Australia and the US State of Hawaii, Oceania counts 27 countries and territories, populated in 2016 by approximately 41 M inhabitants (17 M not counting Australia) among which approximately 12.5 M (11.6 M not counting Australia) consider themselves of entire or partial ONENA ancestry. Specific genetic and microbiome traits of ONENA may be unique and need further investigation to adjust risk estimates, risk prevention, diagnostic and therapeutic strategies, to the benefit of populations in the Pacific and beyond.

Minimum founding populations for the first peopling of Sahul

The timing, context and nature of the first people to enter Sahul is still poorly understood owing to a fragmented archaeological record. However, quantifying the plausible demographic context of this founding population is essential to determine how and why the initial peopling of Sahul occurred. We developed a stochastic, age-structured model using demographic rates from hunter-gatherer societies, and relative carrying capacity hindcasted with LOVECLIM's net primary productivity for northern Sahul. We projected these populations to determine the resilience and minimum sizes required to avoid extinction. A census founding population of between 1,300 and 1,550 individuals was necessary to maintain a quasi-extinction threshold of ≲0.1. This minimum founding population could have arrived at a single point in time, or through multiple voyages of ≥130 people over ~700-900 years. This result shows that substantial population amalgamation in Sunda and Wallacea in Marine Isotope Stages 3-4 provided the conditions for the successful, large-scale and probably planned peopling of Sahul.

Climate change impact on migration, travel, travel destinations and the tourism industry

Background: Climate change is not only increasing ambient temperature but also accelerating the frequency, duration and intensity of extreme weather and climate events, such as heavy precipitation and droughts, and causing sea level rise, which can lead to population displacement. Climate change-related reductions in land productivity and habitability and in food and water security can also interact with demographic, economic and social factors to increase migration. In addition to migration, climate change has also implications for travel and the risk of disease. This article discusses the impact of climate change on migration and travel with implications for public health practice. Methods: Literature review. Results: Migrants may be at increased risk of communicable and non-communicable diseases, due to factors in their country of origin and their country of destination or conditions that they experience during migration. Although migration has not been a significant driver of communicable disease outbreaks to date, public health authorities need to ensure that effective screening and vaccination programmes for priority communicable diseases are in place.Population growth coupled with socio-economic development is increasing travel and tourism, and advances in technology have increased global connectivity and reduced the time required to cover long distances. At the same time, as a result of climate change, many temperate regions, including high-income countries, are now suitable for vector-borne disease transmission. This is providing opportunities for importation of vectors and pathogens from endemic areas that can lead to cases or outbreaks of communicable diseases with which health professionals may be unfamiliar. Conclusion: Health systems need to be prepared for the potential population health consequences of migration, travel and tourism and the impact of climate change on these. Integrated surveillance, early detection of cases and other public health interventions are critical to protect population health and prevent and control communicabledisease outbreaks.

Out of Africa by spontaneous migration waves

Hominin evolution is characterized by progressive regional differentiation, as well as migration waves, leading to anatomically modern humans that are assumed to have emerged in Africa and spread over the whole world. Why or whether Africa was the source region of modern humans and what caused their spread remains subject of ongoing debate. We present a spatially explicit, stochastic numerical model that includes ongoing mutations, demic diffusion, assortative mating and migration waves. Diffusion and assortative mating alone result in a structured population with relatively homogeneous regions bound by sharp clines. The addition of migration waves results in a power-law distribution of wave areas: for every large wave, many more small waves are expected to occur. This suggests that one or more out-of-Africa migrations would probably have been accompanied by numerous smaller migration waves across the world. The migration waves are considered "spontaneous", as the current model excludes environmental or other extrinsic factors. Large waves preferentially emanate from the central areas of large, compact inhabited areas. During the Pleistocene, Africa was the largest such area most of the time, making Africa the statistically most likely origin of anatomically modern humans, without a need to invoke additional environmental or ecological drivers.

Genetic diversity of Mycobacterium tuberculosis in south coastal Karnataka, India, using spoligotyping

Background & objectives:

Despite high occurrence of tuberculosis in India very little information is available about the genetic diversity of Mycobacterium tuberculosis isolates prevailing in coastal Karnataka, India. Thus, the present study was undertaken to explore the genetic biodiversity of M. tuberculosis isolates prevailing in south coastal region of Karnataka (Udupi District), India.

Methods:

A total of 111 Mycobacterial isolates were cultured in Lowenstein Jensen (LJ) medium and after obtaining growth, DNA was extracted and spoligotyping was performed. SITVIT WEB database was used to locate families of spoligotypes.

Results:

On analyzing the hybridization results of all 111 isolates on SITVIT WEB database 57 (51.35%) isolates were clustered into 11 Spoligotype International Types (SIT). The largest cluster of 14 (12.61%) isolates was SIT-48 (EAI1-SOM), followed by SIT-1942 (CAS1-Delhi) with 11 isolates (9.9%) and SIT-11 with seven (6.30%). Moreover, 23 isolates (20.72%) had unique spoligotypes and 31 (27.92%) were orphans. Spotclust analysis revealed that majority (67%) of orphan isolates were variants of CAS (37%) and EAI-5 (34%).

Interpretation & conclusions:

The present study revealed high biodiversity among the circulating isolates of M. tuberculosis in this region with the presence of mixed genotypes earlier reported from north and south India along with certain new genotypes with unique SITs. The study highlights the need for further longitudinal studies to explore the genetic diversity and to understand the transmission dynamics of prevailing isolates.

Hominoid intraspecific cranial variation mirrors neutral genetic diversity

In humans, patterns of cranial variation mirror genetic diversity globally, implicating population history as a key driver of cranial disparity. Here, we demonstrate that the magnitude of genetic diversity within 12 extant ape taxa explains a large proportion of cranial shape variation. Taxa that are more genetically diverse tend to be more cranially diverse also. Our results suggest that neutral evolutionary processes such as mutation, genetic drift, and gene flow are reflected in both genetic and cranial diversity in apes. This work provides a perspective on intraspecific cranial variation in apes which has important implications for interpreting selective and developmental pressures on the cranium and for understanding shape variation in fossil hominin crania.

Natural selection, developmental constraint, and plasticity have all been invoked as explanations for intraspecific cranial variation in humans and apes. However, global patterns of human cranial variation are congruent with patterns of genetic variation, demonstrating that population history has influenced cranial variation in humans. Here we show that this finding is not unique to Homo sapiens but is also broadly evident across extant ape species. Specifically, taxa that exhibit greater intraspecific cranial shape variation also exhibit greater genetic diversity at neutral autosomal loci. Thus, cranial shape variation within hominoid taxa reflects the population history of each species. Our results suggest that neutral evolutionary processes such as mutation, gene flow, and genetic drift have played an important role in generating cranial variation within species. These findings are consistent with previous work on human cranial morphology and improve our understanding of the evolutionary processes that generate intraspecific cranial shape diversity within hominoids. This work has implications for the analysis of selective and developmental pressures on the cranium and for interpreting shape variation in fossil hominin crania.

Ancient Genomics of Modern Humans: The First Decade

The first decade of ancient genomics has revolutionized the study of human prehistory and evolution. We review new insights based on prehistoric modern human genomes, including greatly increased resolution of the timing and structure of the out-of-Africa expansion, the diversification of present-day non-African populations, and the earliest expansions of those populations into Eurasia and America. Prehistoric genomes now document population transformations on every inhabited continent—in particular the effect of agricultural expansions in Africa, Europe, and Oceania—and record a history of natural selection that shapes present-day phenotypic diversity. Despite these advances, much remains unknown, in particular about the genomic histories of Asia (the most populous continent) and Africa (the continent that contains the most genetic diversity). Ancient genomes from these and other regions, integrated with a growing understanding of the genomic basis of human phenotypic diversity, will be in focus during the next decade of research in the field.

Who were the Nataruk people? Mandibular morphology among late Pleistocene and early Holocene fisher-forager populations of West Turkana (Kenya)

Africa is the birthplace of the species Homo sapiens, and Africans today are genetically more diverse than other populations of the world. However, the processes that underpinned the evolution of African populations remain largely obscure. Only a handful of late Pleistocene African fossils (∼50-12 Ka) are known, while the more numerous sites with human fossils of early Holocene age are patchily distributed. In particular, late Pleistocene and early Holocene human diversity in Eastern Africa remains little studied, precluding any analysis of the potential factors that shaped human diversity in the region, and more broadly throughout the continent. These periods include the Last Glacial Maximum (LGM), a moment of extreme aridity in Africa that caused the fragmentation of population ranges and localised extinctions, as well as the 'African Humid Period', a moment of abrupt climate change and enhanced connectivity throughout Africa. East Africa, with its range of environments, may have acted as a refugium during the LGM, and may have played a critical biogeographic role during the heterogene`ous environmental recovery that followed. This environmental context raises a number of questions about the relationships among early Holocene African populations, and about the role played by East Africa in shaping late hunter-gatherer biological diversity. Here, we describe eight mandibles from Nataruk, an early Holocene site (∼10 Ka) in West Turkana, offering the opportunity of exploring population diversity in Africa at the height of the 'African Humid Period'. We use 3D geometric morphometric techniques to analyze the phenotypic variation of a large mandibular sample. Our results show that (i) the Nataruk mandibles are most similar to other African hunter-fisher-gatherer populations, especially to the fossils from Lothagam, another West Turkana locality, and to other early Holocene fossils from the Central Rift Valley (Kenya); and (ii) a phylogenetic connection may have existed between these Eastern African populations and some Nile Valley and Maghrebian groups, who lived at a time when a Green Sahara may have allowed substantial contact, and potential gene flow, across a vast expanse of Northern and Eastern Africa.

A common copy-number variant within SIRPB1 correlates with human Out-of-Africa migration after genetic drift correction

Previous reports have proposed that personality may have played a role on human Out-Of-Africa migration, pinpointing some genetic variants that were positively selected in the migrating populations. In this work, we discuss the role of a common copy-number variant within the SIRPB1 gene, recently associated with impulsive behavior, in the human Out-Of-Africa migration. With the analysis of the variant distribution across forty-two different populations, we found that the SIRPB1 haplotype containing duplicated allele significantly correlated with human migratory distance, being one of the few examples of positively selected loci found across the human world colonization. Circular Chromosome Conformation Capture (4C-seq) experiments from the SIRPB1 promoter revealed important 3D modifications in the locus depending on the presence or absence of the duplication variant. In addition, a 3’ enhancer showed neural activity in transgenic models, suggesting that the presence of the CNV may compromise the expression of SIRPB1 in the central nervous system, paving the way to construct a molecular explanation of the SIRPB1 variants role in human migration.

Genetic and phenotypic variability of iris color in Buenos Aires population

The aim of this work was to describe the phenotypic and genotypic variability related to iris color for the population of Buenos Aires province (Argentina), and to assess the usefulness of current methods of analysis for this country. We studied five Single Nucleotide Polymorphisms (SNPs) included in the IrisPlex kit, in 118 individuals, and we quantified eye color with Digital Iris Analysis Tool. The markers fit Hardy-Weinberg equilibrium for the whole sample, but not for rs12913832 within the group of brown eyes (LR=8.429; p=0.004). We found a remarkable association of HERC2 rs12913832 GG with blue color (p < 0.01) but the other markers did not show any association with iris color. The results for the Buenos Aires population differ from those of other populations of the world for these polymorphisms (p < 0,01). The differences we found might respond to the admixed ethnic composition of Argentina; therefore, methods of analysis used in European populations should be carefully applied when studying the population of Argentina. These findings reaffirm the importance of this investigation in the Argentinian population for people identification based on iris color.

Novel polymorphisms in TICAM2 and NOD1 associated with tuberculosis progression phenotypes in Ethiopian populations

BACKGROUND

Infection by Mycobacterium tuberculosis (Mtb) is a necessary but not sufficient cause for tuberculosis (TB). Although numerous studies suggest human genetic variation may influence TB pathogenesis, there is a conspicuous lack of replication, likely due to imprecise phenotype definition. We aimed to replicate novel findings from a Ugandan cohort in Ethiopian populations.

METHOD

We ascertained TB cases and household controls (n = 292) from three different ethnic groups. Latent Mtb infection was determined using Quantiferon to develop reliable TB progression phenotypes. We sequenced exonic regions of TICAM2 and NOD1.

RESULT

Significant novel associations were observed between two variants in NOD1 and TB: rs751770147 [unadjusted p = 7.28 × 10-5] and chr7:30477156(T), a novel variant, [unadjusted p = 1.04 × 10-4]. Two SNPs in TICAM2 were nominally associated with TB, including rs2288384 [unadjusted p = 0.003]. Haplotype-based association tests supported the SNP-based results.

CONCLUSION

We replicated the association of TICAM2 and NOD1 with TB and identified novel genetic associations with TB in Ethiopian populations.

Global skin colour prediction from DNA

Human skin colour is highly heritable and externally visible with relevance in medical, forensic, and anthropological genetics. Although eye and hair colour can already be predicted with high accuracies from small sets of carefully selected DNA markers, knowledge about the genetic predictability of skin colour is limited. Here, we investigate the skin colour predictive value of 77 single-nucleotide polymorphisms (SNPs) from 37 genetic loci previously associated with human pigmentation using 2025 individuals from 31 global populations. We identified a minimal set of 36 highly informative skin colour predictive SNPs and developed a statistical prediction model capable of skin colour prediction on a global scale. Average cross-validated prediction accuracies expressed as area under the receiver-operating characteristic curve (AUC) ± standard deviation were 0.97 ± 0.02 for Light, 0.83 ± 0.11 for Dark, and 0.96 ± 0.03 for Dark-Black. When using a 5-category, this resulted in 0.74 ± 0.05 for Very Pale, 0.72 ± 0.03 for Pale, 0.73 ± 0.03 for Intermediate, 0.87±0.1 for Dark, and 0.97 ± 0.03 for Dark-Black. A comparative analysis in 194 independent samples from 17 populations demonstrated that our model outperformed a previously proposed 10-SNP-classifier approach with AUCs rising from 0.79 to 0.82 for White, comparable at the intermediate level of 0.63 and 0.62, respectively, and a large increase from 0.64 to 0.92 for Black. Overall, this study demonstrates that the chosen DNA markers and prediction model, particularly the 5-category level; allow skin colour predictions within and between continental regions for the first time, which will serve as a valuable resource for future applications in forensic and anthropologic genetics.

β-Thalassemia Distribution in the Old World: an Ancient Disease Seen from a Historical Standpoint

Background

Haemoglobinopathies constitute the commonest recessive monogenic disorders worldwide, and the treatment of affected individuals presents a substantial global disease burden. β-thalassaemia is characterised by the reduced synthesis (β⁺) or absence (β^o) of the β-globin chains in the HbA molecule, resulting in accumulation of excess unbound α-globin chains that precipitate in erythroid precursors in the bone marrow and in the mature erythrocytes, leading to ineffective erythropoiesis and peripheral haemolysis. Approximately 1.5% of the global population are heterozygotes (carriers) of the β-thalassemias; there is a high incidence in populations from the Mediterranean basin, throughout the Middle East, the Indian subcontinent, Southeast Asia, and Melanesia to the Pacific Islands.

Aim

The principal aim of this paper is to review, from a historical standpoint, our knowledge about an ancient disease, the β-thalassemias, and in particular, when, how and in what way β-thalassemia spread worldwide to reach such high incidences in certain populations.

Results

Mutations involving the β-globin gene are the most common cause of genetic disorders in humans. To date, more than 350 β-thalassaemia mutations have been reported. Considering the current distribution of β- thalassemia, the wide diversity of mutations and the small number of specific mutations in individual populations, it seems unlikely that β-thalassemia originated in a single place and time.

Conclusions

Various processes are known to determine the frequency of genetic disease in human populations. However, it is almost impossible to decide to what extent each process is responsible for the presence of a particular genetic disease. The wide spectrum of β-thalassemia mutations could well be explained by looking at their geographical distribution, the history of malaria, wars, invasions, mass migrations, consanguinity, and settlements. An analysis of the distribution of the molecular spectrum of haemoglobinopathies allows for the development and improvement of diagnostic tests and management of these disorders.

Genetic polymorphism and evolutionary differentiation of Eastern Chinese Han: a comprehensive and comparative analysis on KIRs

Killer cell immunoglobulin-like receptor genes, namely KIRs, cluster together within the 160 kb genomic DNA region. In this study, we used PCR-SSP approach and successfully identified the genotype of 17 KIR genes in 123 independent healthy donors residing in the Jiangsu province, China. All individuals were positive at the 7 genes. The observed carrier gene frequencies (OFs) of remaining 10 KIRs ranged from 14.63% (KIR2DS3) to 95.93% (KIR3DL1). We found 27 distinct genotypes excluding KIR1D. The most frequent occurred in 63 individuals (51.22%). The linkage disequilibrium analysis signified 29 positive and 6 negative relations in 45 pairwise comparisons. To study population differentiation, we drew a Heatmap based on the data of KIRs from 59 populations and conducted Hierarchical Clustering by Euclidean distances. We next validated our results by estimating pairwise D_A distances and illustrating a Neighbor-Joining tree, as well as a MDS plot covering 3 additional Chinese Han groups. The phylogenetic reconstruction and cluster analysis strongly indicated a genetically close relationship between Eastern and Jilin Hans. In conclusion, the present study provided a meritorious resource of KIR genotyping for population genetics, and could be helpful to uncover the genetic mechanism of KIRs in immune disease in the future.

Integrating the signatures of demic expansion and archaic introgression in studies of human population genomics

Human population genomic studies have repeatedly observed a decrease in heterozygosity and an increase in linkage disequilibrium with geographic distance from Africa. While multiple demographic models can generate these patterns, many studies invoke the serial founder effect model, in which populations expand from a single origin and each new population's founders represent a subset of genetic variation in the previous population. The model assumes no admixture with archaic hominins, however, recent studies have identified loci in Homo sapiens bearing signatures of archaic introgression. These results appear to contradict the validity of analyses invoking the serial founder effect model, but we show these two perspectives are compatible. We also propose using the serial founder effect model as a null model for determining the signature of archaic admixture in modern human genomes at different geographic and genomic scales.

Reconstructing Demography and Social Behavior During the Neolithic Expansion from Genomic Diversity Across Island Southeast Asia

Archaeology, linguistics, and increasingly genetics are clarifying how populations moved from mainland Asia, through Island Southeast Asia, and out into the Pacific during the farming revolution. Yet key features of this process remain poorly understood, particularly how social behaviors intersected with demographic drivers to create the patterns of genomic diversity observed across Island Southeast Asia today. Such questions are ripe for computer modeling. Here, we construct an agent-based model to simulate human mobility across Island Southeast Asia from the Neolithic period to the present, with a special focus on interactions between individuals with Asian, Papuan, and mixed Asian-Papuan ancestry. Incorporating key features of the region, including its complex geography (islands and sea), demographic drivers (fecundity and migration), and social behaviors (marriage preferences), the model simultaneously tracks a full suite of genomic markers (autosomes, X chromosome, mitochondrial DNA, and Y chromosome). Using Bayesian inference, model parameters were determined that produce simulations that closely resemble the admixture profiles of 2299 individuals from 84 populations across Island Southeast Asia. The results highlight that greater propensity to migrate and elevated birth rates are related drivers behind the expansion of individuals with Asian ancestry relative to individuals with Papuan ancestry, that offspring preferentially resulted from marriages between Asian women and Papuan men, and that in contrast to current thinking, individuals with Asian ancestry were likely distributed across large parts of western Island Southeast Asia before the Neolithic expansion.

Genetic Diversity and Societally Important Disparities

The magnitude of genetic diversity within human populations varies in a way that reflects the sequence of migrations by which people spread throughout the world. Beyond its use in human evolutionary genetics, worldwide variation in genetic diversity sometimes can interact with social processes to produce differences among populations in their relationship to modern societal problems. We review the consequences of genetic diversity differences in the settings of familial identification in forensic genetic testing, match probabilities in bone marrow transplantation, and representation in genome-wide association studies of disease. In each of these three cases, the contribution of genetic diversity to social differences follows from population-genetic principles. For a fourth setting that is not similarly grounded, we reanalyze with expanded genetic data a report that genetic diversity differences influence global patterns of human economic development, finding no support for the claim. The four examples describe a limit to the importance of genetic diversity for explaining societal differences while illustrating a distinction that certain biologically based scenarios do require consideration of genetic diversity for solving problems to which populations have been differentially predisposed by the unique history of human migrations.

Human Dispersal Out of Africa: A Lasting Debate

Unraveling the first migrations of anatomically modern humans out of Africa has invoked great interest among researchers from a wide range of disciplines. Available fossil, archeological, and climatic data offer many hypotheses, and as such genetics, with the advent of genome-wide genotyping and sequencing techniques and an increase in the availability of ancient samples, offers another important tool for testing theories relating to our own history. In this review, we report the ongoing debates regarding how and when our ancestors left Africa, how many waves of dispersal there were and what geographical routes were taken. We explore the validity of each, using current genetic literature coupled with some of the key archeological findings.

A three-dimensional geometric morphometrics view of the cranial shape variation and population history in the New World

OBJECTIVES

Craniofacial variation in past and present Amerindians has been attributed to the effect of multiple founder events, or to one major migration followed by in situ differentiation and possibly recurrent contacts among Circum-Arctic groups. Our study aims to: (i) detect morphological differences that may indicate several migrations; (ii) test for the presence of genetic isolation; and (iii) test the correlation between shape data and competing settlement hypotheses by taking into account geography, chronology, climate effects, the presence of genetic isolation and recurrent gene flow.

METHODS

We analyzed a large sample of three-dimensional (3D) cranial surface scans (803 specimens) including past and modern groups from America and Australasia. Shape variation was investigated using geometric morphometrics. Differential external gene flow was evaluated by applying genetic concepts to morphometric data (Relethford-Blangero approach). Settlement hypotheses were tested using a matrix correlation approach (Mantel tests).

RESULTS

Our results highlight the strong dichotomy between Circum-Arctic and continental Amerindians as well as the impact of climate adaptation, and possibly recurrent gene flow in the Circum-Arctic area. There is also evidence for the impact of genetic isolation on phenetic variation in Baja California. Several settlement hypotheses are correlated with our data.

CONCLUSIONS

The three approaches used in this study highlight the importance of local processes especially in Baja California, and caution against the use of overly simplistic models when searching for the number of migration events. The results stress the complexity of the settlement of the Americas as well as the mosaic nature of the processes involved in this process. Am. J. Hum. Biol. 28:646-661, 2016. © 2016 Wiley Periodicals, Inc.

Genetics and psychotic disorders: A fresh look at consanguinity

Consanguineous unions refer to marriages between related individuals who share a common ancestor. These unions are still commonplace in certain regions of the world such as the southern coast of the Mediterranean, throughout the Middle East and South-East Asia. According to available data, couples of second cousins or closer and their offspring currently represent 10.4% of the world's population, thus resulting in increased frequencies of autosomal recessive disorders. Furthermore, consanguinity may be implicated in the increased frequency of multifactorial pathologies such as mental disorders. The few existing epidemiological studies in consanguineous and/or geographically isolated populations confirm that there is a significant association between consanguinity and mental disorders and a higher risk of schizophrenia or bipolar disorders among offspring from consanguineous couples. There exists a strong and complex genetic component in the predisposition to psychotic disorders that has been confirmed in numerous studies. However, the genetic basis of these disorders remains poorly understood. GWAS studies (Genome Wide Association Studies) over the past 10 years have identified a few weak associations, thus refuting the "common diseases-common variants" hypothesis. A model implicating numerous rare variants has been supported by the recent discovery of CNVs (Copy Number Variants) and their statistically significant association with psychiatric disorders such as schizophrenia, bipolar disorders and autism. The study of consanguineous families may contribute to identifying rare variants in homogenous populations who have conserved certain alleles. Major developments in molecular biology techniques would facilitate these studies as well as contributing to identifying major genes. These results emphasize the need for genetic counseling in high-risk communities and the importance of implementing preventive actions and raising awareness concerning the risk of consanguineous unions.

Long-Distance Dispersal Shaped Patterns of Human Genetic Diversity in Eurasia

Most previous attempts at reconstructing the past history of human populations did not explicitly take geography into account or considered very simple scenarios of migration and ignored environmental information. However, it is likely that the last glacial maximum (LGM) affected the demography and the range of many species, including our own. Moreover, long-distance dispersal (LDD) may have been an important component of human migrations, allowing fast colonization of new territories and preserving high levels of genetic diversity. Here, we use a high-quality microsatellite data set genotyped in 22 populations to estimate the posterior probabilities of several scenarios for the settlement of the Old World by modern humans. We considered models ranging from a simple spatial expansion to others including LDD and a LGM-induced range contraction, as well as Neolithic demographic expansions. We find that scenarios with LDD are much better supported by data than models without LDD. Nevertheless, we show evidence that LDD events to empty habitats were strongly prevented during the settlement of Eurasia. This unexpected absence of LDD ahead of the colonization wave front could have been caused by an Allee effect, either due to intrinsic causes such as an inbreeding depression built during the expansion or due to extrinsic causes such as direct competition with archaic humans. Overall, our results suggest only a relatively limited effect of the LGM contraction on current patterns of human diversity. This is in clear contrast with the major role of LDD migrations, which have potentially contributed to the intermingled genetic structure of Eurasian populations.

Helicobacter pylori Eradication in Patients with Immune Thrombocytopenic Purpura: A Review and the Role of Biogeography

Idiopathic thrombocytopenic purpura (ITP) is typically a diagnosis of exclusion, assigned by clinicians after ruling out other identifiable etiologies. Since a report by Gasbarrini et al. in 1998, an accumulating body of evidence has proposed a pathophysiological link between ITP and chronic Helicobacter pylori (H. pylori) infection. Clinical reports have described a spontaneous resolution of ITP symptoms in about 50% of chronic ITP patients following empirical treatment of H. pylori infection, but response appears to be geography dependent. Studies have also documented that ITP patients in East Asian countries are more likely to express positive antibody titers against H. pylori-specific cytotoxic-associated gene A (CagA), a virulence factor that is associated with an increased risk for gastric diseases including carcinoma. While a definitive mechanism by which H. pylori may induce thrombocytopenia remains elusive, proposed pathways include molecular mimicry of CagA by host autoantibodies against platelet surface glycoproteins, as well as perturbations in the phagocytic activity of monocytes. Traditional treatments of ITP have been largely empirical, involving the use of immunosuppressive agents and immunoglobulin therapy. However, based on the findings of clinical reports emerging over the past 20 years, health organizations around the world increasingly suggest the detection and eradication of H. pylori as a treatment for ITP. Elucidating the exact molecular mechanisms of platelet activation in H. pylori-positive ITP patients, while considering biogeographical differences in response rates, could offer insight into how best to use clinical H. pylori eradication to treat ITP, but will require well-designed studies to confirm the suggested causative relationship between bacterial infection and an autoimmune disease state.

Genetic diversity and drug susceptibility profile of Mycobacterium tuberculosis isolated from different regions of India

OBJECTIVES

Molecular genotyping profiles of Mycobacterium tuberculosis (MTB) provide a valuable insight into the evolution and transmission of the bacilli. Due to the lack of comprehensive national level data from India on this subject, we performed this study to determine the recent trends and distribution of various MTB lineages circulating in India.

METHODS

A total of 628 MTB isolates were obtained from North, West, South, Central and Eastern India. Spoligotyping and drug susceptibility testing was performed by using manufacturer's instructions.

RESULTS

Spoligotyping detected 102 distinct spoligo-patterns. A total of 536 (85.3%) isolates were distributed into 85 SITs which matched the pre-existing database, whereas 17 SITs were newly created for 34 (5.4%) isolates. Overall, CAS family genotype was predominant, comprising 222 (35.4%) isolates, followed by EAI in 152 (24.2%), Beijing in 108 (17.2%), Manu in 41 (6.5%), T in 30 (4.8%), H in 6 (0.9%), X in 3 (0.5%) and one (0.2%) each in Ural and AFRI. Drug susceptibility testing identified 134 (21.3%) isolates as multi drug resistant (MDR).

CONCLUSIONS

The CAS lineage had a pan India presence but EAI lineage was confined to southern parts of India. Beijing genotype of MTB was significantly associated (p-value <0.0001) with MDR.

POPULATION GENETICS. Genomic evidence for the Pleistocene and recent population history of Native Americans

How and when the Americas were populated remains contentious. Using ancient and modern genome-wide data, we found that the ancestors of all present-day Native Americans, including Athabascans and Amerindians, entered the Americas as a single migration wave from Siberia no earlier than 23 thousand years ago (ka) and after no more than an 8000-year isolation period in Beringia. After their arrival to the Americas, ancestral Native Americans diversified into two basal genetic branches around 13 ka, one that is now dispersed across North and South America and the other restricted to North America. Subsequent gene flow resulted in some Native Americans sharing ancestry with present-day East Asians (including Siberians) and, more distantly, Australo-Melanesians. Putative "Paleoamerican" relict populations, including the historical Mexican Pericúes and South American Fuego-Patagonians, are not directly related to modern Australo-Melanesians as suggested by the Paleoamerican Model.