Demic expansions and human evolution

Selective sweep probabilities in spatially expanding populations

Evolution during range expansions shapes biological systems from microbial communities and tumours up to invasive species. A fundamental question is whether, when a beneficial mutation arises during a range expansion, it will evade clonal interference and sweep through the population to fixation. However, most theoretical investigations of range expansions have been confined to regimes in which selective sweeps are effectively impossible, while studies of selective sweeps have either assumed constant population size or have ignored spatial structure. Here we use mathematical modelling and analysis to investigate selective sweep probabilities in the alternative yet biologically relevant scenario in which mutants can outcompete and displace a slowly spreading wildtype. Assuming constant radial expansion speed, we derive probability distributions for the arrival time and location of the first surviving mutant and hence find surprisingly simple approximate and exact expressions for selective sweep probabilities in one, two and three dimensions, which are independent of mutation rate. Namely, the selective sweep probability is approximately 1 - c w t / c m d , where c w t and c m are the wildtype and mutant radial expansion speeds, and d the spatial dimension. Using agent-based simulations, we show that our analytical results accurately predict selective sweep frequencies in the two-dimensional spatial Moran process. We further compare our results with those obtained for alternative growth laws. Parameterizing our model for human tumours, we find that selective sweeps are predicted to be rare except during very early solid tumour growth, thus providing a general, pan-cancer explanation for findings from recent sequencing studies.

Bayesian analyses of direct radiocarbon dates reveal geographic variations in the rate of rice farming dispersal in prehistoric Japan

The adoption of rice farming during the first millennium BC was a turning point in Japanese prehistory, defining the subsequent cultural, linguistic, and genetic variation in the archipelago. Here, we use a suite of novel Bayesian techniques to estimate the regional rates of dispersal and arrival time of rice farming using radiocarbon dates on charred rice remains. Our results indicate substantial variations in the rate of dispersal of rice within the Japanese islands, hinting at the presence of a mixture of demic and cultural diffusion, geographic variations in the suitability of its cultivation, and the possible role of existing social networks in facilitating or hindering the adoption of the new subsistence economy.

Spatially explicit paleogenomic simulations support cohabitation with limited admixture between Bronze Age Central European populations

The Bronze Age is a complex period of social, cultural and economic changes. Recent paleogenomic studies have documented a large and rapid genetic change in early Bronze Age populations from Central Europe. However, the detailed demographic and genetic processes involved in this change are still debated. Here we have used spatially explicit simulations of genomic components to better characterize the demographic and migratory conditions that may have led to this change. We investigated various scenarios representing the expansion of pastoralists from the Pontic steppe, potentially linked to the Yamnaya cultural complex, and their interactions with local populations in Central Europe, considering various eco-evolutionary factors, such as population admixture, competition and long-distance dispersal. Our results do not support direct competition but rather the cohabitation of pastoralists and farmers in Central Europe, with limited gene flow between populations. They also suggest occasional long-distance migrations accompanying the expansion of pastoralists and a demographic decline in both populations following their initial contact. These results link recent archaeological and paleogenomic observations and move further the debate of genomic changes during the early Bronze Age.

Population history of Brazilian south and southeast shellmound builders inferred through dental morphology

OBJECTIVE

The Southeast and South Coast of Brazil was inhabited during most of the Holocene by shellmound builders. Although there are cultural differences in the archaeological record between regions, it is still debatable how these differences may relate to different population histories. Here, we contribute to this discussion by exploring dental morphological affinities between several regional series.

MATERIALS AND METHODS

Dental morphology of 385 individuals from 14 archaeological sites was analyzed using the Arizona State University Dental Anthropology System. Fifteen traits were used to explore morphological affinities among series through Euclidean distance, Mean Measure of Divergence, and Principal Component Analysis. Mantel matrix correlation and partial correlation tests were used to examine the association between biological, geographic, and chronological distances.

RESULTS

Morphological affinities show that ceramic and nonceramic South Coast groups cluster and differ from most Southeast series. In contrast, Southeast coastal and riverine groups display high morphological variance, showing less biological coherence among them. These biological distances between regions are partially explained by geography, but not by chronology.

CONCLUSIONS

The results support that these coastal populations were low-mobility groups. Although interactions between individuals of different regions likely existed, gene flow occurred mostly among individuals from local or adjacent areas. The introduction of ceramic in the South Coast is not associated with changes in dental morphology patterns, suggesting its adoption is not exclusively associated with the arrival of different biological groups. Southeast coastal and riverine groups show high phenotypic diversity, suggesting a different history of human occupation and cultural development than observed in the South Coast.

Evaluating Causes of Current Genetic Gradients of Modern Humans of the Iberian Peninsula

The history of modern humans in the Iberian Peninsula includes a variety of population arrivals sometimes presenting admixture with resident populations. Genetic data from current Iberian populations revealed an overall east–west genetic gradient that some authors interpreted as a direct consequence of the Reconquista, where Catholic Kingdoms expanded their territories toward the south while displacing Muslims. However, this interpretation has not been formally evaluated. Here, we present a qualitative analysis of the causes of the current genetic gradient observed in the Iberian Peninsula using extensive spatially explicit computer simulations based on a variety of evolutionary scenarios. Our results indicate that the Neolithic range expansion clearly produces the orientation of the observed genetic gradient. Concerning the Reconquista (including political borders among Catholic Kingdoms and regions with different languages), if modeled upon a previous Neolithic expansion, it effectively favored the orientation of the observed genetic gradient and shows local isolation of certain regions (i.e., Basques and Galicia). Despite additional evolutionary scenarios could be evaluated to more accurately decipher the causes of the Iberian genetic gradient, here we show that this gradient has a more complex explanation than that previously hypothesized.

Principals about principal components in statistical genetics

Principal components (PCs) are widely used in statistics and refer to a relatively small number of uncorrelated variables derived from an initial pool of variables, while explaining as much of the total variance as possible. Also in statistical genetics, principal component analysis (PCA) is a popular technique. To achieve optimal results, a thorough understanding about the different implementations of PCA is required and their impact on study results, compared to alternative approaches. In this review, we focus on the possibilities, limitations and role of PCs in ancestry prediction, genome-wide association studies, rare variants analyses, imputation strategies, meta-analysis and epistasis detection. We also describe several variations of classic PCA that deserve increased attention in statistical genetics applications.

Whole Genome Sequencing of Four Representatives From the Admixed Population of the United Arab Emirates

Whole genome sequences (WGS) of four nationals of the United Arab Emirates (UAE) at an average coverage of 33X have been completed and described. The selection of suitable subpopulation representatives was informed by a preceding comprehensive population structure analysis. Representatives were chosen based on their central location within the subpopulation on a principal component analysis (PCA) and the degree to which they were admixed. Novel genomic variations among the different subgroups of the UAE population are reported here. Specifically, the WGS analysis identified 4,161,067-4,798,806 variants in the four individual samples, where approximately 80% were single nucleotide polymorphisms (SNPs) and 20% were insertions or deletions (indels). An average of 2.75% was found to be novel variants according to dbSNP (build 151). This is the first report of structural variants (SV) from WGS data from UAE nationals. There were 15,677-20,339 called SVs, of which around 13.5% were novel. The four samples shared 1,399,178 variants, each with distinct variants as follows: 1,085,524 (for the individual denoted as UAE S011), 1,228,559 (UAE S012), 791,072 (UAE S013), and 906,818 (UAE S014). These results show a previously unappreciated population diversity in the region. The synergy of WGS and genotype array data was demonstrated through variant annotation of the former using 2.3 million allele frequencies for the local population derived from the latter technology platform. This novel approach of combining breadth and depth of array and WGS technologies has guided the choice of population genetic representatives and provides complementary, regionalized allele frequency annotation to new genomes comprising millions of loci.

Influence of Paleolithic range contraction, admixture and long-distance dispersal on genetic gradients of modern humans in Asia

Cavalli-Sforza and coauthors originally explored the genetic variation of modern humans throughout the world and observed an overall east-west genetic gradient in Asia. However, the specific environmental and population genetics processes causing this gradient were not formally investigated and promoted discussion in recent studies. Here we studied the influence of diverse environmental and population genetics processes on Asian genetic gradients and identified which could have produced the observed gradient. To do so, we performed extensive spatially-explicit computer simulations of genetic data under the following scenarios: (a) variable levels of admixture between Paleolithic and Neolithic populations, (b) migration through long-distance dispersal (LDD), (c) Paleolithic range contraction induced by the last glacial maximum (LGM), and (d) Neolithic range expansions from one or two geographic origins (the Fertile Crescent and the Yangzi and Yellow River Basins). Next, we estimated genetic gradients from the simulated data and we found that they were sensible to the analysed processes, especially to the range contraction induced by LGM and to the number of Neolithic expansions. Some scenarios were compatible with the observed east-west genetic gradient, such as the Paleolithic expansion with a range contraction induced by the LGM or two Neolithic range expansions from both the east and the west. In general, LDD increased the variance of genetic gradients among simulations. We interpreted the obtained gradients as a consequence of both allele surfing caused by range expansions and isolation by distance along the vast east-west geographic axis of this continent.

Localizing Components of Shared Transethnic Genetic Architecture of Complex Traits from GWAS Summary Data

Despite strong transethnic genetic correlations reported in the literature for many complex traits, the non-transferability of polygenic risk scores across populations suggests the presence of population-specific components of genetic architecture. We propose an approach that models GWAS summary data for one trait in two populations to estimate genome-wide proportions of population-specific/shared causal SNPs. In simulations across various genetic architectures, we show that our approach yields approximately unbiased estimates with in-sample LD and slight upward-bias with out-of-sample LD. We analyze nine complex traits in individuals of East Asian and European ancestry, restricting to common SNPs (MAF > 5%), and find that most common causal SNPs are shared by both populations. Using the genome-wide estimates as priors in an empirical Bayes framework, we perform fine-mapping and observe that high-posterior SNPs (for both the population-specific and shared causal configurations) have highly correlated effects in East Asians and Europeans. In population-specific GWAS risk regions, we observe a 2.8× enrichment of shared high-posterior SNPs, suggesting that population-specific GWAS risk regions harbor shared causal SNPs that are undetected in the other GWASs due to differences in LD, allele frequencies, and/or sample size. Finally, we report enrichments of shared high-posterior SNPs in 53 tissue-specific functional categories and find evidence that SNP-heritability enrichments are driven largely by many low-effect common SNPs.

Analysis of forensic genetic parameters of 22 autosomal STR markers (PowerPlex® Fusion System) in a population sample from Bosnia and Herzegovina

Background: Bosnia and Herzegovina is a multinational and multireligious country, located in the western part of the Balkan Peninsula. Migrations through history were a key factor in the genetic identity of the Bosnian-Herzegovinian population.Aim: To analyse genetic polymorphisms of 22 autosomal short tandem repeat (STR) loci in the population of Bosnia and Herzegovina and to compare STR allele frequencies for STR loci with the reference data for European populations.Subjects and methods: The study was conducted among 600 unrelated individuals from all regions of Bosnia and Herzegovina. Genotyping was performed using the PowerPlex^® Fusion amplification kit. Allele frequencies and statistical parameters were calculated, as well as the genetic distance among analysed populations through the construction of a neighbor-joining dendrogram.Results: STR loci included in the PowerPlex^® Fusion amplification kit showed high discriminatory power indicating their reliability for human identification and paternity testing. The neighbor-joining dendrogram based on the results of genetic distance analysis showed that the Bosnian and Herzegovinian population has the greatest genetic distance from Turkish and Hungarian populations and greatest similarity with Croatian, Slovenian, and Serbian populations.Conclusion: The results of this study strongly support the application of 22 autosomal genetic markers for paternity testing and personal identity testing and are in agreement with most previous human studies in the investigated human populations.

The Early Peopling of the Philippines based on mtDNA

Despite the efforts made to reconstruct the history of modern humans, there are still poorly explored regions that are key for understanding the phylogeography of our species. One of them is the Philippines, which is crucial to unravel the colonization of Southeast Asia and Oceania but where little is known about when and how the first humans arrived. In order to shed light into this settlement, we collected samples from 157 individuals of the Philippines with the four grandparents belonging to the same region and mitochondrial variants older than 20,000 years. Next, we analyzed the hypervariable I mtDNA region by approximate Bayesian computation based on extensive spatially explicit computer simulations to select among several migration routes towards the Philippines and to estimate population genetic parameters of this colonization. We found that the colonization of the Philippines occurred more than 60,000 years ago, with long-distance dispersal and from both north and south migration routes. Our results also suggest an environmental scenario especially optimal for humans, with large carrying capacity and population growth, in comparison to other regions of Asia. In all, our study suggests a rapid expansion of modern humans towards the Philippines that could be associated with the establisment of maritime technologies and favorable environmental conditions.

Fruit extracts to control pathogenic Escherichia coli: A sweet solution

Escherichia coli is a major cause of diarrhea and is as well responsible for extraintestinal infections in humans and animals. Many pathotypes have been defined for this ubiquitous microorganism on the basis of the virulence attributes. For the last 70 years, antibiotics have been used to control infections caused by E. coli. However, with the resistance observed with many strains these drugs are less recommended. Plant extracts, in particular fruit, represent a source of bioactive compounds that could be beneficial in the control of infectious diseases caused by E. coli. These could have bacteriostatic or bactericidal potential or could be used as synergic agents to amplify the activity of antibiotics for which the germs present some level of resistance. Certain studies also revealed that fruit extracts could act directly on virulence characters to attenuate the pathogenic capacity of microorganisms. This review intent to expose the scant but rapidly growing information available that shows that fruit, used as crude extracts or purified molecules, should be considered to manage diverse types of infections caused by E. coli.

Post-LGM coastline evolution of the NW Sicilian Channel: Comparing high-resolution geophysical data with Glacial Isostatic Adjustment modeling

Since about 20,000 years ago, the geography of the Earth has been profoundly modified by the gradual sea-level rise caused by the melting of continental ice sheets. Flat areas and regions characterized by very low gradients experienced, more than others, rapid flooding, with the progressive disappearance of vast coastal territories. Here we present a reconstruction of the late Quaternary coastline evolution of the north-western sector of the Sicilian Channel, constrained by high-resolution seismic profiles where the marker of the post-Last Glacial Maximum (LGM) marine transgression has been clearly identified and mapped. The locations of the post-LGM seismic horizon have been compared with predictions of a Glacial Isostatic Adjustment (GIA) model, which accounts for the migration of the shorelines in response to sea-level rise and for Earth’s rotational and deformational effects associated with deglaciation. We have verified that most of the points mapped through seismic data interpretation fall along the palaeo-coastline that the GIA model predicts for the 21 kyrs B.P. time frame. However, the model shows a misfit in the marine sector between Mazara del Vallo and Sciacca, where the available data indicate a Quaternary tectonic uplift. The analysis of the seismic profiles provides useful constraints to current GIA models. These add on existing histories of relative sea level in the Mediterranean Sea, allowing to gain new insight into the evolution of the palaeo-geography of the region of study and of the whole Sicilian Channel since the LGM, even in areas where direct geophysical observations are not available yet. In this respect, one of the most attractive implications of the ancient coastline evolution is linked with the underwater archaeology. The sea-level rise heavily impacted the distribution of human settlements, possibly forcing site abandonment and migrations, and this is particularly relevant in the Mediterranean basin, the cradle of the western civilization. The underwater traces left by these ancient populations represent the fundamental proofs to reconstruct the early history of our precursors.

Population Structure and Implications on the Genetic Architecture of HIV-1 Phenotypes Within Southern Africa

The interesting history of Southern Africa has put the region in the spotlight for population medical genetics. Major events including the Bantu expansion and European colonialism have imprinted unique genetic signatures within autochthonous populations of Southern Africa, this resulting in differential allele frequencies across the region. This genetic structure has potential implications on susceptibility and resistance to infectious diseases such as human immunodeficiency virus (HIV) infection. Southern Africa is the region affected worst by HIV. Here, we discuss advances made in genome-wide association studies (GWAS) of HIV-1 in the past 12 years and dissect population diversity within Southern Africa. Our findings accentuate that a plethora of factors such as migration, language and culture, admixture, and natural selection have profiled the genetics of the people of Southern Africa. Genetic structure has been observed among the Khoe-San, among Bantu speakers, and between the Khoe-San, Coloureds, and Bantu speakers. Moreover, Southern African populations have complex admixture scenarios. Few GWAS of HIV-1 have been conducted in Southern Africa, with only one of these identifying two novel variants (HCG22rs2535307 and CCNG1kgp22385164) significantly associated with HIV-1 acquisition and progression. High genetic diversity, multi-wave genetic mixture and low linkage disequilibrium of Southern African populations constitute a challenge in identifying genetic variants with modest risk or protective effect against HIV-1. We therefore posit that it is compelling to assess genome-wide contribution of ancestry to HIV-1 infection. We further suggest robust methods that can pin-point population-specific variants that may contribute to the control of HIV-1 in Southern Africa.

Effect of population stratification on SNP-by-environment interaction

Proportions of false-positive rates in genome-wide association analysis are affected by population stratification, and if it is not correctly adjusted, the statistical analysis can produce the large false-negative finding. Therefore various approaches have been proposed to adjust such problems in genome-wide association studies. However, in spite of its importance, a few studies have been conducted in genome-wide single nucleotide polymorphism (SNP)-by-environment interaction studies. In this report, we illustrate in which scenarios can lead to the false-positive rates in association mapping and approach to maintaining the overall type-1 error rate.

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.

Selecting among Alternative Scenarios of Human Evolution by Simulated Genetic Gradients

Selecting among alternative scenarios of human evolution is nowadays a common methodology to investigate the history of our species. This strategy is usually based on computer simulations of genetic data under different evolutionary scenarios, followed by a fitting of the simulated data with the real data. A recent trend in the investigation of ancestral evolutionary processes of modern humans is the application of genetic gradients as a measure of fitting, since evolutionary processes such as range expansions, range contractions, and population admixture (among others) can lead to different genetic gradients. In addition, this strategy allows the analysis of the genetic causes of the observed genetic gradients. Here, we review recent findings on the selection among alternative scenarios of human evolution based on simulated genetic gradients, including pros and cons. First, we describe common methodologies to simulate genetic gradients and apply them to select among alternative scenarios of human evolution. Next, we review previous studies on the influence of range expansions, population admixture, last glacial period, and migration with long-distance dispersal on genetic gradients for some regions of the world. Finally, we discuss this analytical approach, including technical limitations, required improvements, and advice. Although here we focus on human evolution, this approach could be extended to study other species.

Consequences of diverse evolutionary processes on american genetic gradients of modern humans

European genetic gradients of modern humans were initially interpreted as a consequence of the demic diffusion of expanding Neolithic farmers. However, recent studies showed that these gradients may also be influenced by other evolutionary processes such as population admixture or range contractions. Genetic gradients were observed in the Americas, although their specific evolutionary causes were not investigated. Here we extended the approach used to study genetic gradients in Europe to analyze the influence of diverse evolutionary scenarios on American genetic gradients. Using extensive computer simulations, we evaluated the impact of (i) admixture between expansion waves of modern humans, (ii) the presence of ice-sheets during the last glacial maximum (LGM) and (iii) long-distance dispersal (LDD) events, on the genetic gradients (detected by principal component analysis) of the entire continent, North America and South America. The specific simulation of North and South America showed that genetic gradients are usually orthogonal to the direction of range expansions-either expansions from Bering or posterior re-expansions to recolonize northern regions after ice sheets melting-and we suggest that they result from allele surfing processes. Conversely, our results on the entire continent show a northwest-southeast gradient obtained with any scenario, which we interpreted as a consequence of isolation by distance along the long length of the continent. These findings suggest that distinct genetic gradients can be detected at different regions of the Americas and that subcontinent regions present gradients more sensible to evolutionary and environmental factors (such as LDD and the LGM) than the whole continent.

Robust inference of population structure from next-generation sequencing data with systematic differences in sequencing

Motivation

Inferring population structure is important for both population genetics and genetic epidemiology. Principal components analysis (PCA) has been effective in ascertaining population structure with array genotype data but can be difficult to use with sequencing data, especially when low depth leads to uncertainty in called genotypes. Because PCA is sensitive to differences in variability, PCA using sequencing data can result in components that correspond to differences in sequencing quality (read depth and error rate), rather than differences in population structure. We demonstrate that even existing methods for PCA specifically designed for sequencing data can still yield biased conclusions when used with data having sequencing properties that are systematically different across different groups of samples (i.e. sequencing groups). This situation can arise in population genetics when combining sequencing data from different studies, or in genetic epidemiology when using historical controls such as samples from the 1000 Genomes Project.

Results

To allow inference on population structure using PCA in these situations, we provide an approach that is based on using sequencing reads directly without calling genotypes. Our approach is to adjust the data from different sequencing groups to have the same read depth and error rate so that PCA does not generate spurious components representing sequencing quality. To accomplish this, we have developed a subsampling procedure to match the depth distributions in different sequencing groups, and a read-flipping procedure to match the error rates. We average over subsamples and read flips to minimize loss of information. We demonstrate the utility of our approach using two datasets from 1000 Genomes, and further evaluate it using simulation studies.

Availability and implementation

TASER-PC software is publicly available at http://web1.sph.emory.edu/users/yhu30/software.html.

Contact

yijuan.hu@emory.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

The origin and expansion of Pama-Nyungan languages across Australia

It remains a mystery how Pama-Nyungan, the world's largest hunter-gatherer language family, came to dominate the Australian continent. Some argue that social or technological advantages allowed rapid language replacement from the Gulf Plains region during the mid-Holocene. Others have proposed expansions from refugia linked to climatic changes after the last ice age or, more controversially, during the initial colonization of Australia. Here, we combine basic vocabulary data from 306 Pama-Nyungan languages with Bayesian phylogeographic methods to explicitly model the expansion of the family across Australia and test between these origin scenarios. We find strong and robust support for a Pama-Nyungan origin in the Gulf Plains region during the mid-Holocene, implying rapid replacement of non-Pama-Nyungan languages. Concomitant changes in the archaeological record, together with a lack of strong genetic evidence for Holocene population expansion, suggests that Pama-Nyungan languages were carried as part of an expanding package of cultural innovations that probably facilitated the absorption and assimilation of existing hunter-gatherer groups.

Reconciling evidence from ancient and contemporary genomes: a major source for the European Neolithic within Mediterranean Europe

Important gaps remain in our understanding of the spread of farming into Europe, due partly to apparent contradictions between studies of contemporary genetic variation and ancient DNA. It seems clear that farming was introduced into central, northern, and eastern Europe from the south by pioneer colonization. It is often argued that these dispersals originated in the Near East, where the potential source genetic pool resembles that of the early European farmers, but clear ancient DNA evidence from Mediterranean Europe is lacking, and there are suggestions that Mediterranean Europe may have resembled the Near East more than the rest of Europe in the Mesolithic. Here, we test this proposal by dating mitogenome founder lineages from the Near East in different regions of Europe. We find that whereas the lineages date mainly to the Neolithic in central Europe and Iberia, they largely date to the Late Glacial period in central/eastern Mediterranean Europe. This supports a scenario in which the genetic pool of Mediterranean Europe was partly a result of Late Glacial expansions from a Near Eastern refuge, and that this formed an important source pool for subsequent Neolithic expansions into the rest of Europe.

Multi-locus reassessment of a striking discord between mtDNA gene trees and taxonomy across two congeneric species complexes

Resolving relationships among members of the yellow and citrine wagtail species complexes is among the greatest challenges in avian systematics due to arguably the most dramatic disagreements between traditional taxonomy and mtDNA phylogeny. Each species complex is divided into three geographically cohesive mtDNA clades. Each clade from one species complex has a sister from the other complex. Furthermore, one cross-complex pair is more distantly related to the remaining two pairs than are several other wagtail species. To test mtDNA gene tree topology, we sequenced the mtDNA ND2 gene and 11 nuclear introns for seven wagtail species. Our mtDNA gene tree reconstruction supported the results of previous studies, thereby confirming the disagreement between mtDNA phylogeny and taxonomy. However, our multi-locus species tree which used mtDNA clades as "taxa" was consistent with traditional taxonomy regardless of whether mtDNA was included in the analysis or not. Our multi-locus data suggest that despite the presence of strongly supported, geographically structured mtDNA variation, the mtDNA gene tree misrepresents the evolutionary history of the yellow and citrine wagtail complexes. This mito-nuclear discord results from mtDNA representing the biogeographic, but not evolutionary history of these recently radiated Palearctic wagtails.

Successive range expansion promotes diversity and accelerates evolution in spatially structured microbial populations

Successive range expansions occur within all domains of life, where one population expands first (primary expansion) and one or more secondary populations then follow (secondary expansion). In general, genetic drift reduces diversity during range expansion. However, it is not clear whether the same effect applies during successive range expansion, mainly because the secondary population must expand into space occupied by the primary population. Here we used an experimental microbial model system to show that, in contrast to primary range expansion, successive range expansion promotes local population diversity. Because of mechanical constraints imposed by the presence of the primary population, the secondary population forms fractal-like dendritic structures. This divides the advancing secondary population into many small sub-populations and promotes intermixing between the primary and secondary populations. We further developed a mathematical model to simulate the formation of dendritic structures in the secondary population during succession. By introducing mutations in the primary or dendritic secondary populations, we found that mutations are more likely to accumulate in the dendritic secondary populations. Our results thus show that successive range expansion can promote intermixing over the short term and increase genetic diversity over the long term. Our results therefore have potentially important implications for predicting the ecological processes and evolutionary trajectories of microbial communities.

Genetic Heritage of the Balto-Slavic Speaking Populations: A Synthesis of Autosomal, Mitochondrial and Y-Chromosomal Data

The Slavic branch of the Balto-Slavic sub-family of Indo-European languages underwent rapid divergence as a result of the spatial expansion of its speakers from Central-East Europe, in early medieval times. This expansion–mainly to East Europe and the northern Balkans–resulted in the incorporation of genetic components from numerous autochthonous populations into the Slavic gene pools. Here, we characterize genetic variation in all extant ethnic groups speaking Balto-Slavic languages by analyzing mitochondrial DNA (n = 6,876), Y-chromosomes (n = 6,079) and genome-wide SNP profiles (n = 296), within the context of other European populations. We also reassess the phylogeny of Slavic languages within the Balto-Slavic branch of Indo-European. We find that genetic distances among Balto-Slavic populations, based on autosomal and Y-chromosomal loci, show a high correlation (0.9) both with each other and with geography, but a slightly lower correlation (0.7) with mitochondrial DNA and linguistic affiliation. The data suggest that genetic diversity of the present-day Slavs was predominantly shaped in situ, and we detect two different substrata: ‘central-east European’ for West and East Slavs, and ‘south-east European’ for South Slavs. A pattern of distribution of segments identical by descent between groups of East-West and South Slavs suggests shared ancestry or a modest gene flow between those two groups, which might derive from the historic spread of Slavic people.

Modeling the 3D geometry of the cortical surface with genetic ancestry

Knowing how the human brain is shaped by migration and admixture is a critical step in studying human evolution [1, 2], as well as in preventing the bias of hidden population structure in brain research [3, 4]. Yet, the neuroanatomical differences engendered by population history are still poorly understood. Most of the inference relies on craniometric measurements, because morphology of the brain is presumed to be the neurocranium's main shaping force before bones are fused and ossified [5]. Although studies have shown that the shape variations of cranial bones are consistent with population history [6-8], it is unknown how much human ancestry information is retained by the human cortical surface. In our group's previous study, we found that area measures of cortical surface and total brain volumes of individuals of European descent in the United States correlate significantly with their ancestral geographic locations in Europe [9]. Here, we demonstrate that the three-dimensional geometry of cortical surface is highly predictive of individuals' genetic ancestry in West Africa, Europe, East Asia, and America, even though their genetic background has been shaped by multiple waves of migratory and admixture events. The geometry of the cortical surface contains richer information about ancestry than the areal variability of the cortical surface, independent of total brain volumes. Besides explaining more ancestry variance than other brain imaging measurements, the 3D geometry of the cortical surface further characterizes distinct regional patterns in the folding and gyrification of the human brain associated with each ancestral lineage.

Vestiges of an Ancient Border in the Contemporary Genetic Diversity of North-Eastern Europe

It has previously been demonstrated that the advance of the Neolithic Revolution from the Near East through Europe was decelerated in the northernmost confines of the continent, possibly as a result of space and resource competition with lingering Mesolithic populations. Finland was among the last domains to adopt a farming lifestyle, and is characterized by substructuring in the form of a distinct genetic border dividing the northeastern and southwestern regions of the country. To explore the origins of this divergence, the geographical patterns of mitochondrial and Y-chromosomal haplogroups of Neolithic and Mesolithic ancestry were assessed in Finnish populations. The distribution of these uniparental markers revealed a northeastern bias for hunter-gatherer haplogroups, while haplogroups associated with the farming lifestyle clustered in the southwest. In addition, a correlation could be observed between more ancient mitochondrial haplogroup age and eastern concentration. These results coupled with prior archeological evidence suggest the genetic northeast/southwest division observed in contemporary Finland represents an ancient vestigial border between Mesolithic and Neolithic populations undetectable in most other regions of Europe.

Implications of Nubian-Like Core Reduction Systems in Southern Africa for the Identification of Early Modern Human Dispersals

Lithic technologies have been used to trace dispersals of early human populations within and beyond Africa. Convergence in lithic systems has the potential to confound such interpretations, implying connections between unrelated groups. Due to their reductive nature, stone artefacts are unusually prone to this chance appearance of similar forms in unrelated populations. Here we present data from the South African Middle Stone Age sites Uitpanskraal 7 and Mertenhof suggesting that Nubian core reduction systems associated with Late Pleistocene populations in North Africa and potentially with early human migrations out of Africa in MIS 5 also occur in southern Africa during early MIS 3 and with no clear connection to the North African occurrence. The timing and spatial distribution of their appearance in southern and northern Africa implies technological convergence, rather than diffusion or dispersal. While lithic technologies can be a critical guide to human population flux, their utility in tracing early human dispersals at large spatial and temporal scales with stone artefact types remains questionable.

Characterizing Race/Ethnicity and Genetic Ancestry for 100,000 Subjects in the Genetic Epidemiology Research on Adult Health and Aging (GERA) Cohort

Using genome-wide genotypes, we characterized the genetic structure of 103,006 participants in the Kaiser Permanente Northern California multi-ethnic Genetic Epidemiology Research on Adult Health and Aging Cohort and analyzed the relationship to self-reported race/ethnicity. Participants endorsed any of 23 race/ethnicity/nationality categories, which were collapsed into seven major race/ethnicity groups. By self-report the cohort is 80.8% white and 19.2% minority; 93.8% endorsed a single race/ethnicity group, while 6.2% endorsed two or more. Principal component (PC) and admixture analyses were generally consistent with prior studies. Approximately 17% of subjects had genetic ancestry from more than one continent, and 12% were genetically admixed, considering only nonadjacent geographical origins. Self-reported whites were spread on a continuum along the first two PCs, indicating extensive mixing among European nationalities. Self-identified East Asian nationalities correlated with genetic clustering, consistent with extensive endogamy. Individuals of mixed East Asian-European genetic ancestry were easily identified; we also observed a modest amount of European genetic ancestry in individuals self-identified as Filipinos. Self-reported African Americans and Latinos showed extensive European and African genetic ancestry, and Native American genetic ancestry for the latter. Among 3741 genetically identified parent-child pairs, 93% were concordant for self-reported race/ethnicity; among 2018 genetically identified full-sib pairs, 96% were concordant; the lower rate for parent-child pairs was largely due to intermarriage. The parent-child pairs revealed a trend toward increasing exogamy over time; the presence in the cohort of individuals endorsing multiple race/ethnicity categories creates interesting challenges and future opportunities for genetic epidemiologic studies.

Effects of asymmetric nuclear introgression, introgressive mitochondrial sweep, and purifying selection on phylogenetic reconstruction and divergence estimates in the Pacific clade of Locustella warblers

When isolated but reproductively compatible populations expand geographically and meet, simulations predict asymmetric introgression of neutral loci from a local to invading taxon. Genetic introgression may affect phylogenetic reconstruction by obscuring topology and divergence estimates. We combined phylogenetic analysis of sequences from one mtDNA and 12 nuDNA loci with analysis of gene flow among 5 species of Pacific Locustella warblers to test for presence of genetic introgression and its effects on tree topology and divergence estimates. Our data showed that nuDNA introgression was substantial and asymmetrical among all members of superspecies groups whereas mtDNA showed no introgression except a single species pair where the invader's mtDNA was swept by mtDNA of the local species. This introgressive sweep of mtDNA had the opposite direction of the nuDNA introgression and resulted in the paraphyly of the local species' mtDNA haplotypes with respect to those of the invader. Тhe multilocus nuDNA species tree resolved all inter- and intraspecific relationships despite substantial introgression. However, the node ages on the species tree may be underestimated as suggested by the differences in node age estimates based on non-introgressing mtDNA and introgressing nuDNA. In turn, the introgressive sweep and strong purifying selection appear to elongate internal branches in the mtDNA gene tree.

Aspects of ancient mitochondrial DNA analysis in different populations for understanding human evolution

The evolution of modern humans is a long and difficult process which started from their first appearance and continues to the present day. The study of the genetic origin of populations can help to determine population kinship and to better understand the gradual changes of the gene pool in space and time. Mitochondrial DNA (mtDNA) is a proper tool for the determination of the origin of populations due to its high evolutionary importance. Ancient mitochondrial DNA retrieved from museum specimens, archaeological finds and fossil remains can provide direct evidence for population origins and migration processes. Despite the problems with contaminations and authenticity of ancient mitochondrial DNA, there is a developed set of criteria and platforms for obtaining authentic ancient DNA. During the last two decades, the application of different methods and techniques for analysis of ancient mitochondrial DNA gave promising results. Still, the literature is relatively poor with information for the origin of human populations. Using comprehensive phylogeographic and population analyses we can observe the development and formation of the contemporary populations. The aim of this study was to shed light on human migratory processes and the formation of populations based on available ancient mtDNA data.

Natural selection in a postglacial range expansion: the case of the colour cline in the European barn owl

Gradients of variation--or clines--have always intrigued biologists. Classically, they have been interpreted as the outcomes of antagonistic interactions between selection and gene flow. Alternatively, clines may also establish neutrally with isolation by distance (IBD) or secondary contact between previously isolated populations. The relative importance of natural selection and these two neutral processes in the establishment of clinal variation can be tested by comparing genetic differentiation at neutral genetic markers and at the studied trait. A third neutral process, surfing of a newly arisen mutation during the colonization of a new habitat, is more difficult to test. Here, we designed a spatially explicit approximate Bayesian computation (ABC) simulation framework to evaluate whether the strong cline in the genetically based reddish coloration observed in the European barn owl (Tyto alba) arose as a by-product of a range expansion or whether selection has to be invoked to explain this colour cline, for which we have previously ruled out the actions of IBD or secondary contact. Using ABC simulations and genetic data on 390 individuals from 20 locations genotyped at 22 microsatellites loci, we first determined how barn owls colonized Europe after the last glaciation. Using these results in new simulations on the evolution of the colour phenotype, and assuming various genetic architectures for the colour trait, we demonstrate that the observed colour cline cannot be due to the surfing of a neutral mutation. Taking advantage of spatially explicit ABC, which proved to be a powerful method to disentangle the respective roles of selection and drift in range expansions, we conclude that the formation of the colour cline observed in the barn owl must be due to natural selection.

Patterns and dynamics of genetic diversity in Plasmodium falciparum: what past human migrations tell us about malaria

Plasmodium falciparum is the main agent of malaria, one of the major human infectious diseases affecting millions of people worldwide. The genetic diversity of P. falciparum populations is an essential factor in the parasite's ability to adapt to changes in its environment, enabling the development of drug resistance and the evasion from the host immune system through antigenic variation. Therefore, characterizing these patterns and understanding the main drivers of the pathogen's genetic diversity can provide useful inputs for informing control strategies. In this paper, we review the pioneering work led by Professor Kazuyuki Tanabe on the genetic diversity of P. falciparum populations. In a first part, we recall basic results from population genetics for quantifying within-population genetic diversity, and discuss the main mechanisms driving this diversity. Then, we show how these approaches have been used for reconstructing the historical spread of malaria worldwide, and how current patterns of genetic diversity suggest that the pathogen followed our ancestors in their journey out of Africa. Because these results are robust to different types of genetic markers, they provide a baseline for predicting the pathogen's diversity in unsampled populations, and some useful elements for predicting vaccine efficacy and informing malaria control strategies.

Out of Africa, but how and when? The case of hamadryas baboons (Papio hamadryas)

Many species of Arabian mammals are considered to be of Afrotropical origin and for most of them the Red Sea has constituted an obstacle for dispersal since the Miocene-Pliocene transition. There are two possible routes, the 'northern' and the 'southern', for terrestrial mammals (including humans) to move between Africa and Arabia. The 'northern route', crossing the Sinai Peninsula, is confirmed for several taxa by an extensive fossil record, especially from northern Egypt and the Levant, whereas the 'southern route', across the Bab-el-Mandab Strait, which links the Red Sea with the Gulf of Aden, is more controversial, although post-Pliocene terrestrial crossings of the Red Sea might have been possible during glacial maxima when sea levels were low. Hamadryas baboons (Papio hamadryas) are the only baboon taxon to disperse out of Africa and still inhabit Arabia. In this study, we investigate the origin of Arabian hamadryas baboons using mitochondrial sequence data from 294 samples collected in Arabia and Northeast Africa. Through the analysis of the geographic distribution of genetic diversity, the timing of population expansions, and divergence time estimates combined with palaeoecological data, we test: (i) if Arabian and African hamadryas baboons are genetically distinct; (ii) if Arabian baboons exhibit population substructure; and (iii) when, and via which route, baboons colonized Arabia. Our results suggest that hamadryas baboons colonized Arabia during the Late Pleistocene (130-12 kya [thousands of years ago]) and also moved back to Africa. We reject the hypothesis that hamadryas baboons were introduced to Arabia by humans, because the initial colonization considerably predates the earliest records of human seafaring in this region. Our results strongly suggest that the 'southern route' from Africa to Arabia could have been used by hamadryas baboons during the same time period as proposed for modern humans.

Impact of range expansions on current human genomic diversity

The patterns of population genetic diversity depend to a large extent on past demographic history. Most human populations are known to have gone recently through a series of range expansions within and out of Africa, but these spatial expansions are rarely taken into account when interpreting observed genomic diversity, possibly because they are difficult to model. Here we review available evidence in favour of range expansions out of Africa, and we discuss several of their consequences on neutral and selected diversity, including some recent observations on an excess of rare neutral and selected variants in large samples. We further show that in spatially subdivided populations, the sampling strategy can severely impact the resulting genetic diversity and be confounded by past demography. We conclude that ignoring the spatial structure of human population can lead to some misinterpretations of extant genetic diversity.

Relationship between angiotensin-converting enzyme gene polymorphism and respiratory distress syndrome in premature neonates

OBJECTIVE

The aim of this study was to investigate the possible relationship between angiotensin-converting enzyme (ACE) gene polymorphism (D/D and I/D genotypes) and respiratory distress syndrome (RDS) in preterm neonates.

STUDY DESIGN

Our study included 120 preterm neonates (<37 weeks of gestation) with RDS (the patient group) and 120 preterm neonates without RDS (the control group). Blood samples were obtained from patients and control groups, and ACE gene polymorphism was analysed using the polymerase chain reaction method.

RESULTS

D/D genotype was highly significant in the patient group compared with the control group (48.3% of RDS group vs 20% of the control group, P < 0.001). Meanwhile, I/D and I/I genotypes were significantly higher in the control group (75% and 5% of the control group vs 50% and 1.7% of the patient group, P < 0.001). D/D genotype was highly significant in neonates with bronchopulmonary dysplasia (BPD) compared with I/D genotype (P = 0.001).

CONCLUSION

Our results may suggest that D/D genotype is associated with increased risk of RDS and BPD development in preterm neonates.

Stronger spatial genetic structure in recolonized areas than in refugia in the European beech

Extant rear-edge populations located in former glacial refugia remain understudied despite their high conservation value. These populations should have experienced strong genetic drift due to their small size and long isolation. Moreover, the prolonged action of isolation by distance in refugial areas should result in stronger regional spatial genetic structure (SGS) than in recolonized areas, but empirical tests of this prediction are scarce. To fill this gap, we first used a set of 16 microsatellite markers to investigate the genetic structure of European beech in France in 65 populations from three refugial areas and one control recolonized (nonrefugial) area. Then, using the same approach, we reanalysed published isozyme data from 375 populations distributed across the entire species range. We found stronger genetic differentiation among populations in refugia than in recolonized areas. However, contrary to expectations, regional SGS was lower within refugia than within recolonized areas. Published studies presenting similar analyses suggest that our results could have generality across different biogeographical settings and types of organisms. Strong and prolonged genetic drift in refugial areas could have erased the signature of range expansions that is still visible in recolonized areas. Our results therefore suggest that Pleistocene population isolation has played a key role in increasing the genetic complexity of extant rear-edge populations.

Genetic association studies of endothelial nitric oxide synthase gene polymorphisms in women with unexplained recurrent pregnancy loss: a systematic and meta-analysis

Recurrent pregnancy loss (RPL) is a multifactorial disorder, both genetic and environmental factors contribute to the development of RPL. Recently, the association between endothelial nitric oxide synthase (eNOS) gene polymorphisms and unexplained recurrent pregnancy loss (URPL) has been reported, and the results were inconsistent. Hence, we performed the meta-analysis to drive a more precise estimation of association between eNOS polymorphisms and URPL. Odds ratio (OR) and its 95% confidence interval were calculated under co-dominant (AA vs. BB, TT vs. GG) and additive (A vs. B, T vs. G) genetic models. Studies of eNOS intron 4 VNTR and Glu298Asp were separated by ethnicities. 13 studies included 1,769 URPL cases and 1,376 healthy controls on eNOS intron 4 VNTR polymorphism, and 11 studies were involved in Glu298Asp polymorphism with 1,498 URPL cases and 1,123 healthy controls. The integrated results showed that eNOS Glu298Asp polymorphism was associated with URPL [ORs were 1.91 (1.42-2.56), P < 0.001; 1.67 (1.36-2.04), P < 0.001, respectively]. When analyses were separated by ethnic subgroups, the association between eNOS Glu298Asp polymorphism and URPL was only observed in East Asians [OR = 1.88 (1.52-2.33), P < 0.001 under additive model], and there was no association between eNOS intron 4 VNTR polymorphism and URPL in Caucasians and East Asians. The results indicated a significant association between eNOS Glu298Asp polymorphism and URPL in East Asians. No association was observed between eNOS intron 4 VNTR polymorphism and URPL in Caucasians and East Asians.

Assessing the maximum contribution from ancient populations

Ancestral relationships between populations separated by time represent an often neglected dimension in population genetics, a field which historically has focused on analysis of spatially distributed samples from the same point in time. Models are usually straightforward when two time-separated populations are assumed to be completely isolated from all other populations. However, this is usually an unrealistically stringent assumption when there is gene flow with other populations. Here, we investigate continuity in the presence of gene flow from unknown populations. This setup allows a more nuanced treatment of questions regarding population continuity in terms of "level of contribution" from a particular ancient population to a more recent population. We propose a statistical framework which makes use of a biallelic marker sampled at two different points in time to assess population contribution, and present two different interpretations of the concept. We apply the approach to published data from a prehistoric human population in Scandinavia (Malmström H, Gilbert MTP, Thomas MG, Brandström M, Storå J, Molnar P, Andersen PK, Bendixen C, Holmlund G, Götherström A, et al. 2009. Ancient DNA reveals lack of continuity between Neolithic hunter-gatherers and contemporary Scandinavians. Curr Biol. 19:1758-1762) and Pleistocene woolly mammoth (Barnes I, Shapiro B, Lister A, Kuznetsova T, Sher A, Guthrie D, Thomas MG. 2007. Genetic structure and extinction of the woolly mammoth, Mammuthus primigenius. Curr Biol. 17:1072-1075; Debruyne R, Chu G, King CE, Bos K, Kuch M, Schwarz C, Szpak P, Gröcke DR, Matheus P, Zazula G, et al. 2008. Out of America: ancient DNA evidence for a new world origin of late quaternary woolly mammoths. Curr Biol. 18:1320-1326).

Separation of the largest eigenvalues in eigenanalysis of genotype data from discrete subpopulations

We present a mathematical model, and the corresponding mathematical analysis, that justifies and quantifies the use of principal component analysis of biallelic genetic marker data for a set of individuals to detect the number of subpopulations represented in the data. We indicate that the power of the technique relies more on the number of individuals genotyped than on the number of markers.

Reconciling migration models to the Americas with the variation of North American native mitogenomes

In this study we evaluated migration models to the Americas by using the information contained in native mitochondrial genomes (mitogenomes) from North America. Molecular and phylogeographic analyses of B2a mitogenomes, which are absent in Eskimo-Aleut and northern Na-Dene speakers, revealed that this haplogroup arose in North America ∼11-13 ka from one of the founder Paleo-Indian B2 mitogenomes. In contrast, haplogroup A2a, which is typical of Eskimo-Aleuts and Na-Dene, but also present in the easternmost Siberian groups, originated only 4-7 ka in Alaska, led to the first Paleo-Eskimo settlement of northern Canada and Greenland, and contributed to the formation of the Na-Dene gene pool. However, mitogenomes also show that Amerindians from northern North America, without any distinction between Na-Dene and non-Na-Dene, were heavily affected by an additional and distinctive Beringian genetic input. In conclusion, most mtDNA variation (along the double-continent) stems from the first wave from Beringia, which followed the Pacific coastal route. This was accompanied or followed by a second inland migratory event, marked by haplogroups X2a and C4c, which affected all Amerindian groups of Northern North America. Much later, the ancestral A2a carriers spread from Alaska, undertaking both a westward migration to Asia and an eastward expansion into the circumpolar regions of Canada. Thus, the first American founders left the greatest genetic mark but the original maternal makeup of North American Natives was subsequently reshaped by additional streams of gene flow and local population dynamics, making a three-wave view too simplistic.

Next generation sequencing and rare genetic variants: from human population studies to medical genetics

The allelic frequency spectrum emerging from several Next Generation Sequencing (NGS) projects is revealing important details about evolutionary and demographic forces that shaped the human genome. Herein, we discuss some of the achievements of the use of low-frequency and rare variants from NGS studies. The majority of variants that affect protein-coding regions are recent and rare. Often, the novel rare variants are enriched for deleterious alleles and are population-specific, making them suitable for the study of disease susceptibility. To investigate this kind of variation and its effects in association studies, very large sample sizes will be necessary to achieve sufficient statistical power. Moreover, as these variants are typically population-specific, the replication of disease associations across populations could be very difficult due to population stratification. Therefore, the design of experiments focusing on the identification of rare variants and their effects should be carefully planned. Although several successes have already been achieved through NGS for genetic epidemiology, pharmacogenetic and clinical purposes, with improvements of the sequencing technology and decreased costs, further advances are expected in the near future.

Localized versus generalist phenotypes in a broadly distributed tropical mammal: how is intraspecific variation distributed across disparate environments?

BACKGROUND

The extent of phenotypic differentiation in response to local environmental conditions is a key component of species adaptation and persistence. Understanding the structuring of phenotypic diversity in response to local environmental pressures can provide important insights into species evolutionary dynamics and responses to environmental change. This work examines the influence of steep environmental gradients on intraspecific phenotypic variation and tests two hypotheses about how the tropical soft grass mouse, Akodon mollis (Cricetidae, Rodentia), contends with the disparate environmental conditions encompassed by its broad distribution. Specifically, we test if the species expresses a geographically unstructured, or generalist, phenotype throughout its range or if it shows geographically localized morphological differentiation across disparate environments.

RESULTS

Using geometric morphometric and ecomorphological analyses of skull shape variation we found that despite distinct environmental conditions, geographically structured morphological variation is limited, with the notable exception of a distinct morphological disjunction at the high-elevation forest-grassland transition in the southern portion of A. mollis distribution. Based on genetic analyses, geographic isolation alone does not explain this localized phenotype, given that similar levels of genetic differentiation were also observed among individuals inhabiting other ecosystems that are nonetheless not distinct morphologically.

CONCLUSIONS

Instead of phenotypic specialization across environments in these tropical mountains, there was limited differentiation of skull shape and size across the broad range of A. mollis, with the exception of individuals from the puna, the highest-elevation ecosystem. The high morphological variance among individuals, together with a weak association with local environmental conditions, not only highlights the flexibility of A. mollis' skull, but also highlights the need for further study to understand what maintains the observed morphological patterns. The work also indicates that mechanisms other than processes linked to local ecological specialization as a driver of diversification may contribute to the high diversity of this tropical region.

A novel multi-scale modeling approach to infer whole genome divergence

We propose a novel and simple approach to elucidate genomic patterns of divergence using principal component analysis (PCA). We applied this methodology to the metric space generated by M. musculus genome-wide SNPs. Distance profiles were computed between M. musculus and its closely related species, M. spretus, which was used as external reference. While the speciation dynamics were apparent in the first principal component, the within M. musculus differentiation dimensions gave rise to three minor components. We were unable to obtain a clear divergence signature discriminating laboratory strains, suggesting a stronger effect of genetic drift. These results were at odds with wild strains which exhibit defined deterministic signals of divergence. Finally, we were able to rank novel and previously known genes according to their likelihood to be under selective pressure. In conclusion, we posit PCA as a robust methodology to unravel diverging DNA regions without any a priori forcing.