Molecular analysis reveals tighter social regulation of immigration in patrilocal populations than in matrilocal populations

Patrilineal segmentary systems provide a peaceful explanation for the post-Neolithic Y-chromosome bottleneck

Studies have found a pronounced decline in male effective population sizes worldwide around 3000-5000 years ago. This bottleneck was not observed for female effective population sizes, which continued to increase over time. Until now, this remarkable genetic pattern was interpreted as the result of an ancient structuring of human populations into patrilineal groups (gathering closely related males) violently competing with each other. In this scenario, violence is responsible for the repeated extinctions of patrilineal groups, leading to a significant reduction in male effective population size. Here, we propose an alternative hypothesis by modelling a segmentary patrilineal system based on anthropological literature. We show that variance in reproductive success between patrilineal groups, combined with lineal fission (i.e., the splitting of a group into two new groups of patrilineally related individuals), can lead to a substantial reduction in the male effective population size without resorting to the violence hypothesis. Thus, a peaceful explanation involving ancient changes in social structures, linked to global changes in subsistence systems, may be sufficient to explain the reported decline in Y-chromosome diversity.

Estimating parameters of a stochastic cell invasion model with fluorescent cell cycle labelling using approximate Bayesian computation

We develop a parameter estimation method based on approximate Bayesian computation (ABC) for a stochastic cell invasion model using fluorescent cell cycle labelling with proliferation, migration and crowding effects. Previously, inference has been performed on a deterministic version of the model fitted to cell density data, and not all parameters were identifiable. Considering the stochastic model allows us to harness more features of experimental data, including cell trajectories and cell count data, which we show overcomes the parameter identifiability problem. We demonstrate that, while difficult to collect, cell trajectory data can provide more information about the parameters of the cell invasion model. To handle the intractability of the likelihood function of the stochastic model, we use an efficient ABC algorithm based on sequential Monte Carlo. Rcpp and MATLAB implementations of the simulation model and ABC algorithm used in this study are available at https://github.com/michaelcarr-stats/FUCCI.

Coalescence times, life history traits and conservation concerns: An example from four coastal shark species from the Indo-Pacific

Dispersal abilities play a crucial role in shaping the extent of population genetic structure, with more mobile species being panmictic over large geographical ranges and less mobile ones organized in metapopulations exchanging migrants to different degrees. In turn, population structure directly influences the coalescence pattern of the sampled lineages, but the consequences on the estimated variation of the effective population size (N_e ) over time obtained by means of unstructured demographic models remain poorly understood. However, this knowledge is crucial for biologically interpreting the observed N_e trajectory and further devising conservation strategies in endangered species. Here we investigated the demographic history of four shark species (Carharhinus melanopterus, Carharhinus limbatus, Carharhinus amblyrhynchos, Galeocerdo cuvier) with different degrees of endangered status and life history traits related to dispersal distributed in the Indo-Pacific and sampled off New Caledonia. We compared several evolutionary scenarios representing both structured (metapopulation) and unstructured models and then inferred the N_e variation through time. By performing extensive coalescent simulations, we provided a general framework relating the underlying population structure and the observed N_e dynamics. On this basis, we concluded that the recent decline observed in three out of the four considered species when assuming unstructured demographic models can be explained by the presence of population structure. Furthermore, we also demonstrated the limits of the inferences based on the sole site frequency spectrum and warn that statistics based on linkage disequilibrium will be needed to exclude recent demographic events affecting meta-populations.

An Evolutionary Approach to the History of Barley (Hordeum vulgare) Cultivation in the Canary Islands

The Canary Islands are an archipelago that lies about 100 km west of North Africa. Barley (Hordeum vulgare) has been continuously cultivated since the colonization of the islands. To investigate the agricultural history of the islands, the DNA from multiple individuals of six extant landraces of barley was sequenced, and the resulting data were analyzed with ABC modeling. Estimates of separation times of barley populations on the different islands and the mainland were congruent with archaeological dating of the earliest settlements on the islands. The results of the genetic analyses were consistent with the continuous cultivation of barley on Lanzarote island since it was first colonized, but suggested cultivation was carried out at a smaller scale than on Gran Canaria and Tenerife. Contrary to archaeological evidence and early written historical sources, the genetic analyses suggest that barley was cultivated on a larger scale on Tenerife than on Gran Canaria. The genetic analysis of contemporary barley added support to the dating of the colonization of the islands and pointed to the need for more archaeological data concerning barley cultivation on Tenerife.

Sex-biased patterns shaped the genetic history of Roma

The Roma population is a European ethnic minority characterized by recent and multiple dispersals and founder effects. After their origin in South Asia around 1,500 years ago, they migrated West. In Europe, they diverged into ethnolinguistically distinct migrant groups that spread across the continent. Previous genetic studies based on genome-wide data and uniparental markers detected Roma founder events and West-Eurasian gene flow. However, to the best of our knowledge, it has not been assessed whether these demographic processes have equally affected both sexes in the population. The present study uses the largest and most comprehensive dataset of complete mitochondrial and Y chromosome Roma sequences to unravel the sex-biased patterns that have shaped their genetic history. The results show that the Roma maternal genetic pool carries a higher lineage diversity from South Asia, as opposed to a single paternal South Asian lineage. Nonetheless, the European gene flow events mainly occurred through the maternal lineages; however, a signal of this gene flow is also traceable in the paternal lineages. We also detect a higher female migration rate among European Roma groups. Altogether, these results suggest that sociocultural factors influenced the emergence of sex-biased genetic patterns at global and local scales in the Roma population through time.

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.

From matrimonial practices to genetic diversity in Southeast Asian populations: the signature of the matrilineal puzzle

In matrilineal populations, the descent group affiliation is transmitted by women whereas the socio-political power frequently remains in the hands of men. This situation, named the 'matrilineal puzzle', is expected to promote local endogamy as a coping mechanism allowing men to maintain their decision-making power over their natal descent group. In this paper, we revisit this 'matrilineal puzzle' from a population genetics' point of view. Indeed, such tendency for local endogamy in matrilineal populations is expected to increase their genetic inbreeding and generate isolation-by-distance patterns between villages. To test this hypothesis, we collected ethno-demographic data for 3261 couples and high-density genetic data for 675 individuals from 11 Southeast Asian populations with a wide range of social organizations: matrilineal and matrilocal populations (M), patrilineal and patrilocal populations (P) or cognatic populations with predominant matrilocal residence (C). We observed that M and C populations have higher levels of village endogamy than P populations, and that such higher village endogamy leads to higher genetic inbreeding. M populations also exhibit isolation-by-distance patterns between villages. We interpret such genetic patterns as the signature of the 'matrilineal puzzle'. Notably, our results suggest that any form of matrilocal marriage (whatever the descent rule is) increases village endogamy. These findings suggest that male dominance, when combined with matrilocality, constrains inter-village migrations, and constitutes an underexplored cultural process shaping genetic patterns in human populations. This article is part of the theme issue 'The evolution of female-biased kinship in humans and other mammals'.

A genetic perspective on Longobard-Era migrations

From the first century AD, Europe has been interested by population movements, commonly known as Barbarian migrations. Among these processes, the one involving the Longobard culture interested a vast region, but its dynamics and demographic impact remains largely unknown. Here we report 87 new complete mitochondrial sequences coming from nine early-medieval cemeteries located along the area interested by the Longobard migration (Czech Republic, Hungary and Italy). From the same areas, we sampled necropoleis characterized by cultural markers associated with the Longobard culture (LC) and coeval burials where no such markers were found, or with a chronology slightly preceding the presumed arrival of the Longobards in that region (NLC). Population genetics analysis and demographic modeling highlighted a similarity between LC individuals, as reflected by the sharing of quite rare haplogroups and by the degree of genetic resemblance between Hungarian and Italian LC necropoleis estimated via a Bayesian approach, ABC. The demographic model receiving the strongest statistical support also postulates a contact between LC and NLC communities, thus indicating a complex dynamics of admixture in medieval Europe.

Cultural Innovations Influence Patterns of Genetic Diversity in Northwestern Amazonia

Human populations often exhibit contrasting patterns of genetic diversity in the mtDNA and the nonrecombining portion of the Y-chromosome (NRY), which reflect sex-specific cultural behaviors and population histories. Here, we sequenced 2.3 Mb of the NRY from 284 individuals representing more than 30 Native American groups from Northwestern Amazonia (NWA) and compared these data to previously generated mtDNA genomes from the same groups, to investigate the impact of cultural practices on genetic diversity and gain new insights about NWA population history. Relevant cultural practices in NWA include postmarital residential rules and linguistic exogamy, a marital practice in which men are required to marry women speaking a different language. We identified 2,969 SNPs in the NRY sequences, only 925 of which were previously described. The NRY and mtDNA data showed different sex-specific demographic histories: female effective population size has been larger than that of males through time, which might reflect larger variance in male reproductive success. Both markers show an increase in lineage diversification beginning ∼5,000 years ago, which may reflect the intensification of agriculture, technological innovations, and the expansion of regional trade networks documented in the archaeological evidence. Furthermore, we find similar excesses of NRY versus mtDNA between-population divergence at both the local and continental scale, suggesting long-term stability of female versus male migration. We also find evidence of the impact of sociocultural practices on diversity patterns. Finally, our study highlights the importance of analyzing high-resolution mtDNA and NRY sequences to reconstruct demographic history, since this can differ considerably between sexes.

Comparing genetic diversity and demographic history in co-distributed wild South American camelids

Vicuñas and guanacos are two species of wild South American camelids that are key ruminants in the ecosystems where they occur. Although closely related, these species feature differing ecologies and life history characters, which are expected to influence both their genetic diversity and population differentiation at different spatial scales. Here, using mitochondrial and microsatellite genetic markers, we show that vicuña display lower genetic diversity within populations than guanaco but exhibit more structure across their Peruvian range, which may reflect a combination of natural genetic differentiation linked to geographic isolation and recent anthropogenic population declines. Coalescent-based demographic analyses indicate that both species have passed through a strong bottleneck, reducing their effective population sizes from over 20,000 to less than 1000 individuals. For vicuña, this bottleneck is inferred to have taken place ~3300 years ago, but to have occurred more recently for guanaco at ~2000 years ago. These inferred dates are considerably later than the onset of domestication (when the alpaca was domesticated from the vicuña while the llama was domesticated from the guanaco), coinciding instead with a major human population expansion following the mid-Holocene cold period. As importantly, they imply earlier declines than the well-documented Spanish conquest, where major mass mortality events were recorded for Andean human and camelid populations. We argue that underlying species' differences and recent demographic perturbations have influenced genetic diversity in modern vicuña and guanaco populations, and these processes should be carefully evaluated in the development and implementation of management strategies for these important genetic resources.

Residence rule flexibility and descent groups dynamics shape uniparental genetic diversities in South East Asia

OBJECTIVES

Social organization plays a major role in shaping human population genetic diversity. In particular, matrilocal populations tend to exhibit less mitochondrial diversity than patrilocal populations, and the other way around for Y chromosome diversity. However, several studies have not replicated such findings. The objective of this study is to understand the reasons for such inconsistencies and further evaluate the influence of social organization on genetic diversity.

MATERIALS AND METHODS

We explored uniparental diversity patterns using mitochondrial HV1 sequences and 17 Y-linked short tandem repeats (STRs) in 12 populations (n = 619) from mainland South-East Asia exhibiting a wide range of social organizations, along with quantitative ethno-demographic information sampled at the individual level.

RESULTS

MtDNA diversity was lower in matrilocal than in multilocal and patrilocal populations while Y chromosome diversity was similar among these social organizations. The reasons for such asymmetry at the genetic level were understood by quantifying sex-specific migration rates from our ethno-demographic data: while female migration rates varied between social organizations, male migration rates did not. This unexpected lack of difference in male migrations resulted from a higher flexibility in residence rule in patrilocal than in matrilocal populations. In addition, our data suggested an impact of clan fission process on uniparental genetic patterns.

CONCLUSIONS

The observed lack of signature of patrilocality on Y chromosome patterns might be attributed to the higher residence flexibility in the studied patrilocal populations, thus providing a potential explanation for the apparent discrepancies between social and genetic structures. Altogether, this study highlights the need to quantify the actual residence and descent patterns to fit social to genetic structures.

Complete mitochondrial sequences from Mesolithic Sardinia

Little is known about the genetic prehistory of Sardinia because of the scarcity of pre-Neolithic human remains. From a genetic perspective, modern Sardinians are known as genetic outliers in Europe, showing unusually high levels of internal diversity and a close relationship to early European Neolithic farmers. However, how far this peculiar genetic structure extends and how it originated was to date impossible to test. Here we present the first and oldest complete mitochondrial sequences from Sardinia, dated back to 10,000 yBP. These two individuals, while confirming a Mesolithic occupation of the island, belong to rare mtDNA lineages, which have never been found before in Mesolithic samples and that are currently present at low frequencies not only in Sardinia, but in the whole Europe. Preliminary Approximate Bayesian Computations, restricted by biased reference samples for Mesolithic Sardinia (the two typed samples) and Neolithic Europe (limited to central and north European sequences), suggest that the first inhabitants of the island have had a small or negligible contribution to the present-day Sardinian population, which mainly derives its genetic diversity from continental migration into the island by Neolithic times.

Spatially Explicit Models to Investigate Geographic Patterns in the Distribution of Forensic STRs: Application to the North-Eastern Mediterranean

Human forensic STRs used for individual identification have been reported to have little power for inter-population analyses. Several methods have been developed which incorporate information on the spatial distribution of individuals to arrive at a description of the arrangement of diversity. We genotyped at 16 forensic STRs a large population sample obtained from many locations in Italy, Greece and Turkey, i.e. three countries crucial to the understanding of discontinuities at the European/Asian junction and the genetic legacy of ancient migrations, but seldom represented together in previous studies. Using spatial PCA on the full dataset, we detected patterns of population affinities in the area. Additionally, we devised objective criteria to reduce the overall complexity into reduced datasets. Independent spatially explicit methods applied to these latter datasets converged in showing that the extraction of information on long- to medium-range geographical trends and structuring from the overall diversity is possible. All analyses returned the picture of a background clinal variation, with regional discontinuities captured by each of the reduced datasets. Several aspects of our results are confirmed on external STR datasets and replicate those of genome-wide SNP typings. High levels of gene flow were inferred within the main continental areas by coalescent simulations. These results are promising from a microevolutionary perspective, in view of the fast pace at which forensic data are being accumulated for many locales. It is foreseeable that this will allow the exploitation of an invaluable genotypic resource, assembled for other (forensic) purposes, to clarify important aspects in the formation of local gene pools.

Estimating wildlife disease dynamics in complex systems using an Approximate Bayesian Computation framework

Emerging infectious diseases of wildlife are of increasing concern to managers and conservation policy makers, but are often difficult to study and predict due to the complexity of host-disease systems and a paucity of empirical data. We demonstrate the use of an Approximate Bayesian Computation statistical framework to reconstruct the disease dynamics of bovine tuberculosis in Kruger National Park's lion population, despite limited empirical data on the disease's effects in lions. The modeling results suggest that, while a large proportion of the lion population will become infected with bovine tuberculosis, lions are a spillover host and long disease latency is common. In the absence of future aggravating factors, bovine tuberculosis is projected to cause a lion population decline of ~3% over the next 50 years, with the population stabilizing at this new equilibrium. The Approximate Bayesian Computation framework is a new tool for wildlife managers. It allows emerging infectious diseases to be modeled in complex systems by incorporating disparate knowledge about host demographics, behavior, and heterogeneous disease transmission, while allowing inference of unknown system parameters.

Seed exchange networks, ethnicity, and sorghum diversity

Recent studies investigating the relationship between crop genetic diversity and human cultural diversity patterns showed that seed exchanges are embedded in farmers' social organization. However, our understanding of the social processes involved remains limited. We investigated how farmers' membership in three major social groups interacts in shaping sorghum seed exchange networks in a cultural contact zone on Mount Kenya. Farmers are members of residence groups at the local scale and of dialect groups clustered within larger ethnolinguistic units at a wider scale. The Chuka and Tharaka, who are allied in the same ethnolinguistic unit, coexist with the Mbeere dialect group in the study area. We assessed farmers' homophily, propensity to exchange seeds with members of the same group, using exponential random graph models. We showed that homophily is significant within both residence and ethnolinguistic groups. At these two levels, homophily is driven by the kinship system, particularly by the combination of patrilocal residence and ethnolinguistic endogamy, because most seeds are exchanged among relatives. Indeed, residential homophily in seed exchanges results from local interactions between women and their in-law family, whereas at a higher level, ethnolinguistic homophily is driven by marriage endogamy. Seed exchanges and marriage ties are interrelated, and both are limited between the Mbeere and the other groups, although frequent between the Chuka and Tharaka. The impact of these social homophily processes on crop diversity is discussed.

Genetic Evidence Supports the Multiethnic Character of Teopancazco, a Neighborhood Center of Teotihuacan, Mexico (AD 200-600)

Multiethnicity in Teopancazco, Teotihuacan, is supported by foreign individuals found in the neighborhood center as well as by the diversity observed in funerary rituals at the site. Studies of both stable and strontium isotopes as well as paleodietary analysis, suggest that the population of Teopancazco was composed by three population groups: people from Teotihuacan, people from nearby sites (Tlaxcala-Hidalgo-Puebla), and people from afar, including the coastal plains. In an attempt to understand the genetic dynamics in Teopancazco we conducted an ancient DNA (aDNA) analysis based on mtDNA. Our results show that the level of genetic diversity is consistent with the multiethnicity phenomenon at the neighborhood center. Levels of genetic diversity at different time periods of Teopancazco’s history show that multiethnicity was evident since the beginning and lasted until the collapse of the neighborhood center. However, a PCA and a Neighbor-Joining tree suggested the presence of a genetically differentiated group (buried at the Transitional phase) compared to the population from the initial phase (Tlamimilolpa) as well as the population from the final phase (Xolalpan) of the history of Teopancazco. Genetic studies showed no differences in genetic diversity between males and females in the adult population of Teopancazco, this data along with ample archaeological evidence, suggest a neolocal post-marital pattern of residence in Teopancazco. Nevertheless, genetic analyses on the infant population showed that the males are significantly more heterogeneous than the females suggesting a possible differential role in cultural practices by sex in the infant sector. Regarding interpopulation analysis, we found similar indices of genetic diversity between Teopancazco and heterogeneous native groups, which support the multiethnic character of Teopancazco. Finally, our data showed a close genetic relationship between Teopancazco and populations from the “Teotihuacan corridor” and from Oaxaca and the Maya region, in agreement with previous archaeological evidence.

Large-scale recent expansion of European patrilineages shown by population resequencing

The proportion of Europeans descending from Neolithic farmers ∼ 10 thousand years ago (KYA) or Palaeolithic hunter-gatherers has been much debated. The male-specific region of the Y chromosome (MSY) has been widely applied to this question, but unbiased estimates of diversity and time depth have been lacking. Here we show that European patrilineages underwent a recent continent-wide expansion. Resequencing of 3.7 Mb of MSY DNA in 334 males, comprising 17 European and Middle Eastern populations, defines a phylogeny containing 5,996 single-nucleotide polymorphisms. Dating indicates that three major lineages (I1, R1a and R1b), accounting for 64% of our sample, have very recent coalescent times, ranging between 3.5 and 7.3 KYA. A continuous swathe of 13/17 populations share similar histories featuring a demographic expansion starting ∼ 2.1-4.2 KYA. Our results are compatible with ancient MSY DNA data, and contrast with data on mitochondrial DNA, indicating a widespread male-specific phenomenon that focuses interest on the social structure of Bronze Age Europe.

Genealogical relationships between early medieval and modern inhabitants of Piedmont

In the period between 400 to 800 AD, also known as the period of the Barbarian invasions, intense migration is documented in the historical record of Europe. However, little is known about the demographic impact of these historical movements, potentially ranging from negligible to substantial. As a pilot study in a broader project on Medieval Europe, we sampled 102 specimens from 5 burial sites in Northwestern Italy, archaeologically classified as belonging to Lombards or Longobards, a Germanic people ruling over a vast section of the Italian peninsula from 568 to 774. We successfully amplified and typed the mitochondrial hypervariable region I (HVR-I) of 28 individuals. Comparisons of genetic diversity with other ancient populations and haplotype networks did not suggest that these samples are heterogeneous, and hence allowed us to jointly compare them with three isolated contemporary populations, and with a modern sample of a large city, representing a control for the effects of recent immigration. We then generated by serial coalescent simulations 16 millions of genealogies, contrasting a model of genealogical continuity with one in which the contemporary samples are genealogically independent from the medieval sample. Analyses by Approximate Bayesian Computation showed that the latter model fits the data in most cases, with one exception, Trino Vercellese, in which the evidence was compatible with persistence up to the present time of genetic features observed among this early medieval population. We conclude that it is possible, in general, to detect evidence of genealogical ties between medieval and specific modern populations. However, only seldom did mitochondrial DNA data allow us to reject with confidence either model tested, which indicates that broader analyses, based on larger assemblages of samples and genetic markers, are needed to understand in detail the effects of medieval migration.

Genetic variation reveals large-scale population expansion and migration during the expansion of Bantu-speaking peoples

The majority of sub-Saharan Africans today speak a number of closely related languages collectively referred to as 'Bantu' languages. The current distribution of Bantu-speaking populations has been found to largely be a consequence of the movement of people rather than a diffusion of language alone. Linguistic and single marker genetic studies have generated various hypotheses regarding the timing and the routes of the Bantu expansion, but these hypotheses have not been thoroughly investigated. In this study, we re-analysed microsatellite markers typed for large number of African populations that-owing to their fast mutation rates-capture signatures of recent population history. We confirm the spread of west African people across most of sub-Saharan Africa and estimated the expansion of Bantu-speaking groups, using a Bayesian approach, to around 5600 years ago. We tested four different divergence models for Bantu-speaking populations with a distribution comprising three geographical regions in Africa. We found that the most likely model for the movement of the eastern branch of Bantu-speakers involves migration of Bantu-speaking groups to the east followed by migration to the south. This model, however, is only marginally more likely than other models, which might indicate direct movement from the west and/or significant gene flow with the western Branch of Bantu-speakers. Our study use multi-loci genetic data to explicitly investigate the timing and mode of the Bantu expansion and it demonstrates that west African groups rapidly expanded both in numbers and over a large geographical area, affirming the fact that the Bantu expansion was one of the most dramatic demographic events in human history.

Spatial and temporal simulation of human evolution. Methods, frameworks and applications

Analyses of human evolution are fundamental to understand the current gradients of human diversity. In this concern, genetic samples collected from current populations together with archaeological data are the most important resources to study human evolution. However, they are often insufficient to properly evaluate a variety of evolutionary scenarios, leading to continuous debates and discussions. A commonly applied strategy consists of the use of computer simulations based on, as realistic as possible, evolutionary models, to evaluate alternative evolutionary scenarios through statistical correlations with the real data. Computer simulations can also be applied to estimate evolutionary parameters or to study the role of each parameter on the evolutionary process. Here we review the mainly used methods and evolutionary frameworks to perform realistic spatially explicit computer simulations of human evolution. Although we focus on human evolution, most of the methods and software we describe can also be used to study other species. We also describe the importance of considering spatially explicit models to better mimic human evolutionary scenarios based on a variety of phenomena such as range expansions, range shifts, range contractions, sex-biased dispersal, long-distance dispersal or admixtures of populations. We finally discuss future implementations to improve current spatially explicit simulations and their derived applications in human evolution.

Gene-flow in a mosaic hybrid zone: is local introgression adaptive?

Genome-wide scans of genetic differentiation between hybridizing taxa can identify genome regions with unusual rates of introgression. Regions of high differentiation might represent barriers to gene flow, while regions of low differentiation might indicate adaptive introgression-the spread of selectively beneficial alleles between reproductively isolated genetic backgrounds. Here we conduct a scan for unusual patterns of differentiation in a mosaic hybrid zone between two mussel species, Mytilus edulis and M. galloprovincialis. One outlying locus, mac-1, showed a characteristic footprint of local introgression, with abnormally high frequency of edulis-derived alleles in a patch of M. galloprovincialis enclosed within the mosaic zone, but low frequencies outside of the zone. Further analysis of DNA sequences showed that almost all of the edulis allelic diversity had introgressed into the M. galloprovincialis background in this patch. We then used a variety of approaches to test the hypothesis that there had been adaptive introgression at mac-1. Simulations and model fitting with maximum-likelihood and approximate Bayesian computation approaches suggested that adaptive introgression could generate a "soft sweep," which was qualitatively consistent with our data. Although the migration rate required was high, it was compatible with the functioning of an effective barrier to gene flow as revealed by demographic inferences. As such, adaptive introgression could explain both the reduced intraspecific differentiation around mac-1 and the high diversity of introgressed alleles, although a localized change in barrier strength may also be invoked. Together, our results emphasize the need to account for the complex history of secondary contacts in interpreting outlier loci.

Demographic processes shaping genetic variation of the solitarious phase of the desert locust

Between plagues, the solitarious desert locust (Schistocerca gregaria) is generally thought to exist as small populations, which are particularly prone to extinction events in arid regions of Africa and Asia. Given the high genetic structuring observed in one geographical area (the Eritrean coast) by former authors, a metapopulation dynamics model involving repeated extinction and colonization events was favoured. In this study, we assessed the validity of a demographic scenario involving temporary populations of the solitarious phase of the desert locust by analysing large-scale population genetic data. We scored 24 microsatellites in 23 solitarious population samples collected over most of the species range during remission. We found very little genetic structuring and little evidence of declining genetic diversity. A Bayesian clustering method distinguished four genetically differentiated units. Three groups were largely consistent with three population samples which had undergone recent bottleneck events. Nevertheless, the last genetically homogeneous unit included all individuals from the remaining 18 population samples and did not show evidence of demographic disequilibrium. An approximate Bayesian computation treatment indicated a large population size for this main genetic group, moderately reduced between plague and remission but still containing tens of thousands of individuals. Our results diverge from the hypothesis of a classical metapopulation dynamics model. They instead support the scenario in which large populations persist in the solitarious phase of the desert locust.

Association mapping, patterns of linkage disequilibrium and selection in the vicinity of the PHYTOCHROME C gene in pearl millet

Linkage analysis confirmed the association in the region of PHYC in pearl millet. The comparison of genes found in this region suggests that PHYC is the best candidate. Major efforts are currently underway to dissect the phenotype-genotype relationship in plants and animals using existing populations. This method exploits historical recombinations accumulated in these populations. However, linkage disequilibrium sometimes extends over a relatively long distance, particularly in genomic regions containing polymorphisms that have been targets for selection. In this case, many genes in the region could be statistically associated with the trait shaped by the selected polymorphism. Statistical analyses could help in identifying the best candidate genes into such a region where an association is found. In a previous study, we proposed that a fragment of the PHYTOCHROME C gene (PHYC) is associated with flowering time and morphological variations in pearl millet. In the present study, we first performed linkage analyses using three pearl millet F2 families to confirm the presence of a QTL in the vicinity of PHYC. We then analyzed a wider genomic region of ~100 kb around PHYC to pinpoint the gene that best explains the association with the trait in this region. A panel of 90 pearl millet inbred lines was used to assess the association. We used a Markov chain Monte Carlo approach to compare 75 markers distributed along this 100-kb region. We found the best candidate markers on the PHYC gene. Signatures of selection in this region were assessed in an independent data set and pointed to the same gene. These results foster confidence in the likely role of PHYC in phenotypic variation and encourage the development of functional studies.

Genetic evidence does not support an Etruscan origin in Anatolia

The debate on the origins of Etruscans, documented in central Italy between the eighth century BC and the first century AD, dates back to antiquity. Herodotus described them as a group of immigrants from Lydia, in Western Anatolia, whereas for Dionysius of Halicarnassus they were an indigenous population. Dionysius' view is shared by most modern archeologists, but the observation of similarities between the (modern) mitochondrial DNAs (mtDNAs) of Turks and Tuscans was interpreted as supporting an Anatolian origin of the Etruscans. However, ancient DNA evidence shows that only some isolates, and not the bulk of the modern Tuscan population, are genetically related to the Etruscans. In this study, we tested alternative models of Etruscan origins by Approximate Bayesian Computation methods, comparing levels of genetic diversity in the mtDNAs of modern and ancient populations with those obtained by millions of computer simulations. The results show that the observed genetic similarities between modern Tuscans and Anatolians cannot be attributed to an immigration wave from the East leading to the onset of the Etruscan culture in Italy. Genetic links between Tuscany and Anatolia do exist, but date back to a remote stage of prehistory, possibly but not necessarily to the spread of farmers during the Neolithic period.

Indian Ocean crossroads: human genetic origin and population structure in the Maldives

The Maldives are an 850 km-long string of atolls located centrally in the northern Indian Ocean basin. Because of this geographic situation, the present-day Maldivian population has potential for uncovering genetic signatures of historic migration events in the region. We therefore studied autosomal DNA-, mitochondrial DNA-, and Y-chromosomal DNA markers in a representative sample of 141 unrelated Maldivians, with 119 from six major settlements. We found a total of 63 different mtDNA haplotypes that could be allocated to 29 mtDNA haplogroups, mostly within the M, R, and U clades. We found 66 different Y-STR haplotypes in 10 Y-chromosome haplogroups, predominantly H1, J2, L, R1a1a, and R2. Parental admixture analysis for mtDNA- and Y-haplogroup data indicates a strong genetic link between the Maldive Islands and mainland South Asia, and excludes significant gene flow from Southeast Asia. Paternal admixture from West Asia is detected, but cannot be distinguished from admixture from South Asia. Maternal admixture from West Asia is excluded. Within the Maldives, we find a subtle genetic substructure in all marker systems that is not directly related to geographic distance or linguistic dialect. We found reduced Y-STR diversity and reduced male-mediated gene flow between atolls, suggesting independent male founder effects for each atoll. Detected reduced female-mediated gene flow between atolls confirms a Maldives-specific history of matrilocality. In conclusion, our new genetic data agree with the commonly reported Maldivian ancestry in South Asia, but furthermore suggest multiple, independent immigration events and asymmetrical migration of females and males across the archipelago.

Origins and evolution of the Etruscans' mtDNA

The Etruscan culture is documented in Etruria, Central Italy, from the 8^th to the 1^st century BC. For more than 2,000 years there has been disagreement on the Etruscans’ biological origins, whether local or in Anatolia. Genetic affinities with both Tuscan and Anatolian populations have been reported, but so far all attempts have failed to fit the Etruscans’ and modern populations in the same genealogy. We extracted and typed the hypervariable region of mitochondrial DNA of 14 individuals buried in two Etruscan necropoleis, analyzing them along with other Etruscan and Medieval samples, and 4,910 contemporary individuals from the Mediterranean basin. Comparing ancient (30 Etruscans, 27 Medieval individuals) and modern DNA sequences (370 Tuscans), with the results of millions of computer simulations, we show that the Etruscans can be considered ancestral, with a high degree of confidence, to the current inhabitants of Casentino and Volterra, but not to the general contemporary population of the former Etruscan homeland. By further considering two Anatolian samples (35 and 123 individuals) we could estimate that the genetic links between Tuscany and Anatolia date back to at least 5,000 years ago, strongly suggesting that the Etruscan culture developed locally, and not as an immediate consequence of immigration from the Eastern Mediterranean shores.

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

BACKGROUND

When studying the genetic structure of human populations, the role of cultural factors may be difficult to ascertain due to a lack of formal models. Linguistic diversity is a typical example of such a situation. Patrilocality, on the other hand, can be integrated into a biological framework, allowing the formulation of explicit working hypotheses. The present study is based on the assumption that patrilocal traditions make the hypervariable region I of the mtDNA a valuable tool for the exploration of migratory dynamics, offering the opportunity to explore the relationships between genetic and linguistic diversity. We studied 85 Niger-Congo-speaking patrilocal populations that cover regions from Senegal to Central African Republic. A total of 4175 individuals were included in the study.

RESULTS

By combining a multivariate analysis aimed at investigating the population genetic structure, with a Bayesian approach used to test models and extent of migration, we were able to detect a stepping-stone migration model as the best descriptor of gene flow across the region, with the main discontinuities corresponding to forested areas.

CONCLUSIONS

Our analyses highlight an aspect of the influence of habitat variation on human genetic diversity that has yet to be understood. Rather than depending simply on geographic linear distances, patterns of female genetic variation vary substantially between savannah and rainforest environments. Our findings may be explained by the effects of recent gene flow constrained by environmental factors, which superimposes on a background shaped by pre-agricultural peopling.

Sociocultural behavior, sex-biased admixture, and effective population sizes in Central African Pygmies and non-Pygmies

Sociocultural phenomena, such as exogamy or phylopatry, can largely determine human sex-specific demography. In Central Africa, diverging patterns of sex-specific genetic variation have been observed between mobile hunter-gatherer Pygmies and sedentary agricultural non-Pygmies. However, their sex-specific demography remains largely unknown. Using population genetics and approximate Bayesian computation approaches, we inferred male and female effective population sizes, sex-specific migration, and admixture rates in 23 Central African Pygmy and non-Pygmy populations, genotyped for autosomal, X-linked, Y-linked, and mitochondrial markers. We found much larger effective population sizes and migration rates among non-Pygmy populations than among Pygmies, in agreement with the recent expansions and migrations of non-Pygmies and, conversely, the isolation and stationary demography of Pygmy groups. We found larger effective sizes and migration rates for males than for females for Pygmies, and vice versa for non-Pygmies. Thus, although most Pygmy populations have patrilocal customs, their sex-specific genetic patterns resemble those of matrilocal populations. In fact, our results are consistent with a lower prevalence of polygyny and patrilocality in Pygmies compared with non-Pygmies and a potential female transmission of reproductive success in Pygmies. Finally, Pygmy populations showed variable admixture levels with the non-Pygmies, with often much larger introgression from male than from female lineages. Social discrimination against Pygmies triggering complex movements of spouses in intermarriages can explain these male-biased admixture patterns in a patrilocal context. We show how gender-related sociocultural phenomena can determine highly variable sex-specific demography among populations, and how population genetic approaches contrasting chromosomal types allow inferring detailed human sex-specific demographic history.

Testing hypotheses for genealogical discordance in a rainforest lizard

Genealogical discordance, or when different genes tell distinct stories although they evolved under a shared history, often emerges from either coalescent stochasticity or introgression. In this study, we present a strong case of mito-nuclear genealogical discordance in the Australian rainforest lizard species complex of Saproscincus basiliscus and S. lewisi. One of the lineages that comprises this complex, the Southern S. basiliscus lineage, is deeply divergent at the mitochondrial genome but shows markedly less divergence at the nuclear genome. By placing our results in a comparative context and reconstructing the lineages' demography via multilocus and coalescent-based approximate Bayesian computation methods, we test hypotheses for how coalescent variance and introgression contribute to this pattern. These analyses suggest that the observed genealogical discordance likely results from introgression. Further, to generate such strong discordance, introgression probably acted in concert with other factors promoting asymmetric gene flow between the mitochondrial and nuclear genomes, such as selection or sex-biased dispersal. This study offers a framework for testing sources of genealogical discordance and suggests that historical introgression can be an important force shaping the genetic diversity of species and their populations.

Estimating demographic parameters from large-scale population genomic data using Approximate Bayesian Computation

BACKGROUND

The Approximate Bayesian Computation (ABC) approach has been used to infer demographic parameters for numerous species, including humans. However, most applications of ABC still use limited amounts of data, from a small number of loci, compared to the large amount of genome-wide population-genetic data which have become available in the last few years.

RESULTS

We evaluated the performance of the ABC approach for three 'population divergence' models - similar to the 'isolation with migration' model - when the data consists of several hundred thousand SNPs typed for multiple individuals by simulating data from known demographic models. The ABC approach was used to infer demographic parameters of interest and we compared the inferred values to the true parameter values that was used to generate hypothetical "observed" data. For all three case models, the ABC approach inferred most demographic parameters quite well with narrow credible intervals, for example, population divergence times and past population sizes, but some parameters were more difficult to infer, such as population sizes at present and migration rates. We compared the ability of different summary statistics to infer demographic parameters, including haplotype and LD based statistics, and found that the accuracy of the parameter estimates can be improved by combining summary statistics that capture different parts of information in the data. Furthermore, our results suggest that poor choices of prior distributions can in some circumstances be detected using ABC. Finally, increasing the amount of data beyond some hundred loci will substantially improve the accuracy of many parameter estimates using ABC.

CONCLUSIONS

We conclude that the ABC approach can accommodate realistic genome-wide population genetic data, which may be difficult to analyze with full likelihood approaches, and that the ABC can provide accurate and precise inference of demographic parameters from these data, suggesting that the ABC approach will be a useful tool for analyzing large genome-wide datasets.

Investigating sex-biased migration during the Neolithic transition in Europe, using an explicit spatial simulation framework

Cultural practices can deeply influence genetic diversity patterns. The Neolithic transitions that took place at different times and locations around the world led to major cultural and demographic changes that influenced and therefore left their marks on human genetic diversity patterns. Several studies on the European Neolithic transition suggest that mitochondrial DNA (mtDNA) and Y-chromosome data can exhibit different patterns, which could be owing to different demographic histories for females and males. Archaeological and anthropological data suggest that the transition from hunter-gatherers (HGs) to farmers' societies is probably associated with changes in social organization, particularly in post-marital residence (PMR) rules (i.e. patrilocality, matrilocality or bilocality). The movements of humans and genes associated with these rules can be seen as sex-biased short-range migrations. We developed a new individual-based simulation approach to explore the genetic consequences of 45 different scenarios, where we varied the patterns of PMR and admixture between HGs and farmers. We recorded mtDNA and Y-chromosome data and analysed their diversity patterns within and between populations, through time and space. We also collected published mtDNA and Y-chromosome data from European and Near-Eastern populations in order to identify the scenarios that would best explain them. We show that: (i) different PMR systems can lead to different patterns of genetic diversity and differentiation, (ii) asymmetries between mtDNA and Y-chromosome can be owing to different behaviours between males and females, but also to different mutations rates, and (iii) patrilocality in farmers explains the present patterns of genetic diversity better than matrilocality or bilocality. Moreover, we found that (iv) the genetic diversity of farmers change depending on the HGs PMR rules even though they are assumed to disappear more than 5000 years ago in our simulations.

Migration in Afro-Brazilian rural communities: crossing demographic and genetic data

Many studies have used genetic markers to understand global migration patterns of our species. However, there are only few studies of human migration on a local scale. We, therefore, researched migration dynamics in three Afro-Brazilian rural communities, using demographic data and ten Ancestry Informative Markers. In addition to the description of migration and marriage structures, we carried out genetic comparisons between the three populations, as well as between locals and migrants from each community. Genetic admixture analyses were conducted according to the gene-identity method, with Sub-Saharan Africans, Amerindians, and Europeans as parental populations. The three analyzed Afro-Brazilian rural communities consisted of 16% to 30% of migrants, most of them women. The age pyramid revealed a gap in the segment of men aged between 20 to 30 yrs. While endogamous marriages predominated, exogamous marriages were mainly patrilocal. Migration dynamics are apparently associated with matrimonial customs and other social practices of such communities. The impact of migration upon the populations' genetic composition was low but showed an increase in European alleles with a concomitant decrease in the Amerindian contribution. Admixture analysis evidenced a higher African contribution to the gene pool of the studied populations, followed by the contribution of Europeans and Amerindians, respectively.

Evolutionary history of pearl millet (Pennisetum glaucum [L.] R. Br.) and selection on flowering genes since its domestication

The plant domestication process is associated with considerable modifications of plant phenotype. The identification of the genetic basis of this adaptation is of great interest for evolutionary biology. One of the methods used to identify such genes is the detection of signatures of selection. However, domestication is generally associated with major demographic effects. It is therefore crucial to disentangle the effects of demography and selection on diversity. In this study, we investigated selection in a flowering time pathway during domestication of pearl millet. We first used a random set of 20 genes to model pearl millet domestication using approximate Bayesian computation. This analysis showed that a model with exponential growth and wild-cultivated gene flow was well supported by our data set. Under this model, the domestication date of pearl millet is estimated at around 4,800 years ago. We assessed selection in 15 pearl millet DNA sequences homologous to flowering time genes and showed that these genes underwent selection more frequently than expected. We highlighted significant signatures of selection in six pearl millet flowering time genes associated with domestication or improvement of pearl millet. Moreover, higher deviations from neutrality were found for circadian clock-associated genes. Our study provides new insights into the domestication process of pearl millet and shows that a category of genes of the flowering pathway were preferentially selected during pearl millet domestication.

Does speciation between Arabidopsis halleri and Arabidopsis lyrata coincide with major changes in a molecular target of adaptation?

Ever since Darwin proposed natural selection as the driving force for the origin of species, the role of adaptive processes in speciation has remained controversial. In particular, a largely unsolved issue is whether key divergent ecological adaptations are associated with speciation events or evolve secondarily within sister species after the split. The plant Arabidopsis halleri is one of the few species able to colonize soils highly enriched in zinc and cadmium. Recent advances in the molecular genetics of adaptation show that the physiology of this derived ecological trait involves copy number expansions of the AhHMA4 gene, for which orthologs are found in single copy in the closely related A. lyrata and the outgroup A. thaliana. To gain insight into the speciation process, we ask whether adaptive molecular changes at this candidate gene were contemporary with important stages of the speciation process. We first inferred the scenario and timescale of speciation by comparing patterns of variation across the genomic backgrounds of A. halleri and A. lyrata. Then, we estimated the timing of the first duplication of AhHMA4 in A. halleri. Our analysis suggests that the historical split between the two species closely coincides with major changes in this molecular target of adaptation in the A. halleri lineage. These results clearly indicate that these changes evolved in A. halleri well before industrial activities fostered the spread of Zn- and Cd-polluted areas, and suggest that adaptive processes related to heavy-metal homeostasis played a major role in the speciation process.

Marriage exchanges, seed exchanges, and the dynamics of manioc diversity

The conservation of crop genetic resources requires understanding the different variables-cultural, social, and economic-that impinge on crop diversity. In small-scale farming systems, seed exchanges represent a key mechanism in the dynamics of crop genetic diversity, and analyzing the rules that structure social networks of seed exchange between farmer communities can help decipher patterns of crop genetic diversity. Using a combination of ethnobotanical and molecular genetic approaches, we investigated the relationships between regional patterns of manioc genetic diversity in Gabon and local networks of seed exchange. Spatially explicit Bayesian clustering methods showed that geographical discontinuities of manioc genetic diversity mirror major ethnolinguistic boundaries, with a southern matrilineal domain characterized by high levels of varietal diversity and a northern patrilineal domain characterized by low varietal diversity. Borrowing concepts from anthropology--kinship, bridewealth, and filiation--we analyzed the relationships between marriage exchanges and seed exchange networks in patrilineal and matrilineal societies. We demonstrate that, by defining marriage prohibitions, kinship systems structure social networks of exchange between farmer communities and influence the movement of seeds in metapopulations, shaping crop diversity at local and regional levels.

Phylogenetic inconsistency caused by ancient sex-biased gene migration

Inferences of ancient sex-biased migration patterns using sex-linked genetic markers are usually difficult because of a stochastic process of allele fixation. Nevertheless, incongruent phylogenetic trees between different sex-linked markers and between sex-linked and autosomal markers are frequently interpreted as a signature of sex-biased migration without further statistical evaluation. I investigated the types of incongruent phylogenetic trees from which past sex-biased migration events can be statistically supported under the coalescent model. In the case of mammals, detecting a sex-biased migration pattern is not guaranteed by comparing the phylogenetic pattern of mitochondrial and Y-chromosomal loci. Likewise, evidence of introgression at a mitochondrial locus, but not at autosomal loci, does not support the hypothesis of an ancient female-biased migration pattern with statistical significance. In contrast, evidence of introgression at ≥ 5 unlinked autosomal loci, but not at a Y-chromosomal locus, would reject the null hypothesis of a sexually equal migration rate with statistical significance. A similar argument can be made to infer a male-biased migration pattern. Furthermore, the investigation of many recombining sex-biased markers such as X-chromosomal loci in mammals has the potential to efficiently detect ancient sex-biased demographic patterns.

Biological ancestries, kinship connections, and projected identities in four central Anatolian settlements: insights from culturally contextualized genetic anthropology

Previous population genetics studies in Turkey failed to delineate recent historical and social factors that shaped Anatolian cultural and genetic diversity at the local level. To address this shortcoming, we conducted focused ethnohistorical fieldwork and screened biological samples collected from the Yuksekyer region for mitochondrial, Y chromosome, and autosomal markers and then analyzed the data within an ethnohistorical context. Our results revealed that, at the village level, paternal genetic diversity is structured among settlements, whereas maternal genetic diversity is distributed more homogenously, reflecting the strong patrilineal cultural traditions that transcend larger ethnic and religious structures. Local ancestries and origin myths, rather than ethnic or religious affiliations, delineate the social boundaries and projected identities among the villages. Therefore, we conclude that broad, ethnicity-based sampling is inadequate to capture the genetic signatures of recent social and historical dynamics, which have had a profound influence on contemporary genetic and cultural regional diversity.

Estimators of the human effective sex ratio detect sex biases on different timescales

Determining historical sex ratios throughout human evolution can provide insight into patterns of genomic variation, the structure and composition of ancient populations, and the cultural factors that influence the sex ratio (e.g., sex-specific migration rates). Although numerous studies have suggested that unequal sex ratios have existed in human evolutionary history, a coherent picture of sex-biased processes has yet to emerge. For example, two recent studies compared human X chromosome to autosomal variation to make inferences about historical sex ratios but reached seemingly contradictory conclusions, with one study finding evidence for a male bias and the other study identifying a female bias. Here, we show that a large part of this discrepancy can be explained by methodological differences. Specifically, through reanalysis of empirical data, derivation of explicit analytical formulae, and extensive simulations we demonstrate that two estimators of the effective sex ratio based on population structure and nucleotide diversity preferentially detect biases that have occurred on different timescales. Our results clarify apparently contradictory evidence on the role of sex-biased processes in human evolutionary history and show that extant patterns of human genomic variation are consistent with both a recent male bias and an earlier, persistent female bias.

Testing founder effect speciation: divergence population genetics of the spoonbills Platalea regia and Pl. minor (Threskiornithidae, Aves)

Although founder effect speciation has been a popular theoretical model for the speciation of geographically isolated taxa, its empirical importance has remained difficult to evaluate due to the intractability of past demography, which in a founder effect speciation scenario would involve a speciational bottleneck in the emergent species and the complete cessation of gene flow following divergence. Using regression-weighted approximate Bayesian computation, we tested the validity of these two fundamental conditions of founder effect speciation in a pair of sister species with disjunct distributions: the royal spoonbill Platalea regia in Australasia and the black-faced spoonbill Pl. minor in eastern Asia. When compared with genetic polymorphism observed at 20 nuclear loci in the two species, simulations showed that the founder effect speciation model had an extremely low posterior probability (1.55 × 10(-8)) of producing the extant genetic pattern. In contrast, speciation models that allowed for postdivergence gene flow were much more probable (posterior probabilities were 0.37 and 0.50 for the bottleneck with gene flow and the gene flow models, respectively) and postdivergence gene flow persisted for a considerable period of time (more than 80% of the divergence history in both models) following initial divergence (median = 197,000 generations, 95% credible interval [CI]: 50,000-478,000, for the bottleneck with gene flow model; and 186,000 generations, 95% CI: 45,000-477,000, for the gene flow model). Furthermore, the estimated population size reduction in Pl. regia to 7,000 individuals (median, 95% CI: 487-12,000, according to the bottleneck with gene flow model) was unlikely to have been severe enough to be considered a bottleneck. Therefore, these results do not support founder effect speciation in Pl. regia but indicate instead that the divergence between Pl. regia and Pl. minor was probably driven by selection despite continuous gene flow. In this light, we discuss the potential importance of evolutionarily labile traits with significant fitness consequences, such as migratory behavior and habitat preference, in facilitating divergence of the spoonbills.