Neutral Demographic Processes Research Articles

Interpreting generative adversarial networks to infer natural selection from genetic data.

Understanding natural selection and other forms of non-neutrality is a major focus for the use of machine learning in population genetics. Existing methods rely on computationally intensive simulated training data. Unlike efficient neutral coalescent simulations for demographic inference, realistic simulations of selection typically require slow forward simulations. Because there are many possible modes of selection, a high dimensional parameter space must be explored, with no guarantee that the simulated models are close to the real processes. Finally, it is difficult to interpret trained neural networks, leading to a lack of understanding about what features contribute to classification. Here we develop a new approach to detect selection and other local evolutionary processes that requires relatively few selection simulations during training. We build upon a generative adversarial network trained to simulate realistic neutral data. This consists of a generator (fitted demographic model), and a discriminator (convolutional neural network) that predicts whether a genomic region is real or fake. As the generator can only generate data under neutral demographic processes, regions of real data that the discriminator recognizes as having a high probability of being "real" do not fit the neutral demographic model and are therefore candidates for targets of selection. To incentivize identification of a specific mode of selection, we fine-tune the discriminator with a small number of custom non-neutral simulations. We show that this approach has high power to detect various forms of selection in simulations, and that it finds regions under positive selection identified by state-of-the-art population genetic methods in three human populations. Finally, we show how to interpret the trained networks by clustering hidden units of the discriminator based on their correlation patterns with known summary statistics.

Island biogeography theory provides a plausible explanation for why larger vertebrates and taller humans have more diverse gut microbiomes.

Prior work has shown a positive scaling relationship between vertebrate body size, human height, and gut microbiome alpha diversity. This observation mirrors commonly observed species area relationships (SARs) in many other ecosystems. Here, we expand these observations to several large datasets, showing that this size-diversity scaling relationship is independent of relevant covariates, like diet, body mass index, age, sex, bowel movement frequency, antibiotic usage, and cardiometabolic health markers. Island biogeography theory (IBT), which predicts that larger islands tend to harbor greater species diversity through neutral demographic processes, provides a simple mechanism for positive SARs. Using a gut-adapted IBT model, we demonstrated that increasing the length of a flow-through ecosystem led to increased species diversity, closely matching our empirical observations. We delve into the possible clinical implications of these SARs in the American Gut cohort. Consistent with prior observations that lower alpha diversity is a risk factor for Clostridioides difficile infection (CDI), we found that individuals who reported a history of CDI were shorter than those who did not and that this relationship was mediated by alpha diversity. We observed that vegetable consumption had a much stronger association with CDI history, which was also partially mediated by alpha diversity. In summary, we find that the positive scaling observed between body size and gut alpha diversity can be plausibly explained by a gut-adapted IBT model, may be related to CDI risk, and vegetable intake appears to independently mitigate this risk, although additional work is needed to validate the potential disease risk implications.

Adaptive and Neutral Polymorphisms of the Onne-DAB Gene from the Major Histocompatibility Complex (MHC) in Sockeye Salmon Populations on the Asian Range

The variability of an MHC complex gene in sockeye salmon populations throughout the Asian range was studied to identify “footprints” indicative of pathogen-mediated selection and neutral demographic processes that have influenced these populations in both the recent and distant past. Genotype frequencies of a haplotype block consisting of two SNP loci (One_MHC2_109 and One_MHC2_190v2) in the Onne-DAB gene encoding the β-chain of the MHC class II molecule as well as allelic frequencies of 29 putative neutral SNPs have been traced in 27 sockeye salmon populations in the Asian Pacific coast. Differently directed clines of genetic diversity at the MHC2 loci were observed in sockeye salmon populations inhabiting the Sea of Okhotsk and Pacific coasts of the Kamchatka Peninsula. The formation of these clines can be attributed to a combination of historical processes associated with the colonization of the Asian range and the latitudinal gradient of abiotic and biotic factors influencing the variability of the Onne-DAB gene. In continental populations of sockeye salmon, balancing selection was not intense enough to conceal the impacts of demographic and historical processes associated with the fragmentation of the area in the late Pleistocene. In contrast, in island populations, balancing selection effectively maintained the diversity of the Onne-DAB gene despite a significant decrease in polymorphism observed in neutral regions of the genome.

Genomic Responses to Climate Change: Making the Most of the Drosophila Model.

It is pressing to understand how animal populations evolve in response to climate change. We argue that new sequencing technologies and the use of historical samples are opening unprecedented opportunities to investigate genome-wide responses to changing environments. However, there are important challenges in interpreting the emerging findings. First, it is essential to differentiate genetic adaptation from phenotypic plasticity. Second, it is extremely difficult to map genotype, phenotype, and fitness. Third, neutral demographic processes and natural selection affect genetic variation in similar ways. We argue that Drosophila melanogaster, a classical model organism with decades of climate adaptation research, is uniquely suited to overcome most of these challenges. In the near future, long-term time series genome-wide datasets of D. melanogaster natural populations will provide exciting opportunities to study adaptation to recent climate change and will lay the groundwork for related research in non-model systems.

Complex histories of gene flow and a mitochondrial capture event in a nonsister pair of birds.

Hybridization, introgression, and reciprocal gene flow during speciation, specifically the generation of mitonuclear discordance, are increasingly observed as parts of the speciation process. Genomic approaches provide insight into where, when, and how adaptation operates during and after speciation and can measure historical and modern introgression. Whether adaptive or neutral in origin, hybridization can cause mitonuclear discordance by placing the mitochondrial genome of one species (or population) in the nuclear background of another species. The latter, introgressed species may eventually have its own mtDNA replaced or “captured” by other species across its entire geographical range. Intermediate stages in the capture process should be observable. Two nonsister species of Australasian monarch‐flycatchers, Spectacled Monarch (Symposiachrus trivirgatus) mostly of Australia and Indonesia and Spot‐winged Monarch (S. guttula) of New Guinea, present an opportunity to observe this process. We analysed thousands of single nucleotide polymorphisms (SNPs) derived from ultraconserved elements of all subspecies of both species. Mitochondrial DNA sequences of Australian populations of S. trivirgatus form two paraphyletic clades, one being sister to and presumably introgressed by S. guttula despite little nuclear signal of introgression. Population genetic analyses (e.g., tests for modern and historical gene flow and selection) support at least one historical gene flow event between S. guttula and Australian S. trivirgatus. We also uncovered introgression from the Maluku Islands subspecies of S. trivirgatus into an island population of S. guttula, resulting in apparent nuclear paraphyly. We find that neutral demographic processes, not adaptive introgression, are the most likely cause of these complex population histories. We suggest that a Pleistocene extinction of S. guttula from mainland Australia resulted from range expansion by S. trivirgatus.

Patterns of introgression vary within an avian hybrid zone

BackgroundExploring hybrid zone dynamics at different spatial scales allows for better understanding of local factors that influence hybrid zone structure. In this study, we tested hypotheses about drivers of introgression at two spatial scales within the Saltmarsh Sparrow (Ammospiza caudacuta) and Nelson’s Sparrow (A. nelsoni) hybrid zone. Specifically, we evaluated the influence of neutral demographic processes (relative species abundance), natural selection (exogenous environmental factors and genetic incompatibilities), and sexual selection (assortative mating) in this mosaic hybrid zone. By intensively sampling adults (n = 218) and chicks (n = 326) at two geographically proximate locations in the center of the hybrid zone, we determined patterns of introgression on a fine scale across sites of differing habitat. We made broadscale comparisons of patterns from the center with those of prior studies in the southern edge of the hybrid zone.ResultsA panel of fixed SNPs (135) identified from ddRAD sequencing was used to calculate a hybrid index and determine genotypic composition/admixture level of the populations. Another panel of polymorphic SNPs (589) was used to assign paternity and reconstruct mating pairs to test for sexual selection. On a broad-scale, patterns of introgression were not explained by random mating within marshes. We found high rates of back-crossing and similarly low rates of recent-generation (F1/F2) hybrids in the center and south of the zone. Offspring genotypic proportions did not meet those expected from random mating within the parental genotypic distribution. Additionally, we observed half as many F1/F2 hybrid female adults than nestlings, while respective male groups showed no difference, in support of Haldane’s Rule. The observed proportion of interspecific mating was lower than expected when accounting for mate availability, indicating assortative mating was limiting widespread hybridization. On a fine spatial scale, we found variation in the relative influence of neutral and selective forces between inland and coastal habitats, with the smaller, inland marsh influenced primarily by neutral demographic processes, and the expansive, coastal marsh experiencing higher selective pressures in the form of natural (exogenous and endogenous) and sexual selection.ConclusionsMultiple drivers of introgression, including neutral and selective pressures (exogenous, endogenous, and sexual selection), are structuring this hybrid zone, and their relative influence is site and context-dependent.

Disentangling the determinants of transposable elements dynamics in vertebrate genomes using empirical evidences and simulations.

The interactions between transposable elements (TEs) and their hosts constitute one of the most profound co-evolutionary processes found in nature. The population dynamics of TEs depends on factors specific to each TE families, such as the rate of transposition and insertional preference, the demographic history of the host and the genomic landscape. How these factors interact has yet to be investigated holistically. Here we are addressing this question in the green anole (Anolis carolinensis) whose genome contains an extraordinary diversity of TEs (including non-LTR retrotransposons, SINEs, LTR-retrotransposons and DNA transposons). We observed a positive correlation between recombination rate and frequency of TEs and densities for LINEs, SINEs and DNA transposons. For these elements, there was a clear impact of demography on TE frequency and abundance, with a loss of polymorphic elements and skewed frequency spectra in recently expanded populations. On the other hand, some LTR-retrotransposons displayed patterns consistent with a very recent phase of intense amplification. To determine how demography, genomic features and intrinsic properties of TEs interact we ran simulations using SLiM3. We determined that i) short TE insertions are not strongly counter-selected, but long ones are, ii) neutral demographic processes, linked selection and preferential insertion may explain positive correlations between average TE frequency and recombination, iii) TE insertions are unlikely to have been massively recruited in recent adaptation. We demonstrate that deterministic and stochastic processes have different effects on categories of TEs and that a combination of empirical analyses and simulations can disentangle these mechanisms.

Genetic variation in the constitutive defensive metabolome and its inducibility are geographically structured and largely determined by demographic processes in maritime pine

Abstract Interspecific phenotypic variation in plant secondary metabolites (PSM) is often explained by biotic and abiotic factors. However, patterns of variation within species do not clearly fit the theoretical predictions. Exploring how genetics, environment and demographic processes shape such variation among and within populations is crucial for understanding evolution of PSM, particularly in long‐lived plants such as forest trees. Here, we quantified genetic variation in PSM among and within populations, and explored drivers of local adaptation by studying the role of climate as a source of population differentiation in PSM of maritime pine. Constitutive profile and concentrations of 63 PSM and their herbivory‐associated inducibility were determined in the bark of 130 clonally replicated genotypes with known familial structure from 10 populations covering the distribution range of the species. We compared neutral and quantitative population genetic differentiation of PSM (FST and QST). Also, we accounted for population genetic structure and kinship among individuals when exploring climate–trait relationships. We found large population differentiation and additive genetic variation in constitutive PSM. Many PSM were inducible, although very low genetic variation was observed with respect to their inducibility. QST–FST comparisons suggest that differentiation of most diterpenes, monoterpenes, and phenolics can be explained by neutral demographic processes. Spatially heterogeneous selection across populations leading to local adaptation was only found for total constitutive sesquiterpenes and a few individual PSM. After accounting for population genetic structure, only the constitutive concentration of two sesquiterpenes showing signs of diversifying selection was predicted by climate, with decreasing concentrations along a growth‐prone climatic gradient. Synthesis. Evolutionary patterns of plant secondary metabolites depended on their chemical nature, with neutral differentiation governing most plant secondary metabolites. Evidence of local adaptation was only found for total constitutive sesquiterpenes and a few individual plant secondary metabolites. The low genetic variation in the inducibility of plant secondary metabolites suggests a conserved model of defensive induction in this species. Since population differentiation linked to past demographic history could lead to false positives of adaptive differentiation signals, accounting for the genetic relatedness among populations is required to infer the environmental determinants of intraspecific genetic variation in putatively adaptive traits such as plant defences.

Positive selection and climatic effects on MHC class II gene diversity in hares (Lepus capensis) from a steep ecological gradient

In natural populations, allelic diversity of the major histocompatibility complex (MHC) is commonly interpreted as resulting from positive selection in varying spatiotemporal pathogenic landscapes. Composite pathogenic landscape data are, however, rarely available. We studied the spatial distribution of allelic diversity at two MHC class II loci (DQA, DQB) in hares, Lepus capensis, along a steep ecological gradient in North Africa and tested the role of climatic parameters for the spatial distribution of DQA and DQB proteins. Climatic parameters were considered to reflect to some extent pathogenic landscape variation. We investigated historical and contemporary forces that have shaped the variability at both genes, and tested for differential selective pressure across the ecological gradient by comparing allelic variation at MHC and neutral loci. We found positive selection on both MHC loci and significantly decreasing diversity from North to South Tunisia. Our multinomial linear models revealed significant effects of geographical positions that were correlated with mean annual temperature and precipitation on the occurrence of protein variants, but no effects of co-occurring DQA or DQB proteins, respectively. Diversifying selection, recombination, adaptation to local pathogenic landscapes (supposedly reflected by climate parameters) and neutral demographic processes have shaped the observed MHC diversity and differentiation patterns.

Sequence and Structural Diversity of Mouse Y Chromosomes.

Over the 180 My since their origin, the sex chromosomes of mammals have evolved a gene repertoire highly specialized for function in the male germline. The mouse Y chromosome is unique among mammalian Y chromosomes characterized to date in that it is large, gene-rich and euchromatic. Yet, little is known about its diversity in natural populations. Here, we take advantage of published whole-genome sequencing data to survey the diversity of sequence and copy number of sex-linked genes in three subspecies of house mice. Copy number of genes on the repetitive long arm of both sex chromosomes is highly variable, but sequence diversity in nonrepetitive regions is decreased relative to expectations based on autosomes. We use simulations and theory to show that this reduction in sex-linked diversity is incompatible with neutral demographic processes alone, but is consistent with recent positive selection on genes active during spermatogenesis. Our results support the hypothesis that the mouse sex chromosomes are engaged in ongoing intragenomic conflict.

Purifying Selection, Density Blocking and Unnoticed Mitochondrial DNA Diversity in the Red Deer, Cervus elaphus.

The trajectories of postglacial range expansions, the occurrence of lineage patches and the formation and maintenance of secondary contact between lineages may mostly reflect neutral demographic processes, including density blocking, that may leave long-lasting genetic signatures. However, a few studies have recently shown that climate may also play a role. We used red deer, a large, mobile herbivore that is assumed to be sensitive to climate change, to test hypotheses of possible selection on the mitochondrial DNA cytochrome b gene (mtDNA cytb) and competitive and/or density-blocking (using mtDNA control region). We searched for a possible link between the phylogeographic structure and abiotic climatic variables. Finally, we tested for isolation by distance and isolation by environment and assessed the impact of human-mediated translocations on the genetic structure of red deer. Our analysis of 30 red deer populations in Poland using the mtDNA control region (N = 357) and cytochrome b (N = 50) markers not only confirmed the presence of the Western and South-Eastern lineages of the species but also indicated the presence of a previously unnoticed, rare relic haplotype that grouped together C. e. italicus from Italy (the Mesola deer). No significant signs of positive selection were detected for the mtDNA cytb gene in the studied red deer. However, a significant signal for purifying selection was found in our study that may explain the narrowness of the contact zone because gene flow between the Western and South-Eastern lineages should drive relatively strong mito-nuclear incompatibilities. MtDNA control region differentiation among red deer populations in Poland correlated with different abiotic climatic variables. Strikingly, the southernmost ice sheet limits during the Elsterian was the most important factor, and it explained the largest amount of variation. However, neither isolation by distance (IBD) nor isolation by environment (IBE) were recorded, and a very limited impact of human translocations was evident. The above-mentioned results suggest that in contemporary red deer populations in Poland, the phylogeographic pattern is well preserved, and long-term processes (density and/or competitive blocking) still play a major role.

Current hypotheses to explain genetic chaos under the sea.

Chaotic genetic patchiness (CGP) refers to surprising patterns of spatial and temporal genetic structure observed in some marine species at a scale where genetic variation should be efficiently homogenized by gene flow via larval dispersal. Here we review and discuss 4 mechanisms that could generate such unexpected patterns: selection, sweepstakes reproductive success, collective dispersal, and temporal shifts in local population dynamics. First, we review examples where genetic differentiation at specific loci was driven by diversifying selection, which was historically the first process invoked to explain CGP. Second, we turn to neutral demographic processes that may drive genome-wide effects, and whose effects on CGP may be enhanced when they act together. We discuss how sweepstakes reproductive success accelerates genetic drift and can thus generate genetic structure, provided that gene flow is not too strong. Collective dispersal is another mechanism whereby genetic structure can be maintained regardless of dispersal intensity, because it may prevent larval cohorts from becoming entirely mixed. Theoretical analyses of both the sweepstakes and the collective dispersal ideas are presented. Finally, we discuss an idea that has received less attention than the other ones just mentioned, namely temporal shifts in local population dynamics.

Association of transcriptome-wide sequence variation with climate gradients in valley oak (Quercus lobata)

A fundamental goal of evolutionary biology is to understand how environment shapes genetic variation through its effect on demographic processes and through natural selection. In non-model species, transcriptome sequencing generates large single nucleotide polymorphism (SNP) panels to disentangle these influences. Quercus lobata (valley oak) offers an excellent system for such analyses because it has stably occupied a climatically heterogeneous landscape throughout California. We used 220,427 diallelic SNPs from 22 individuals identified against a recently assembled reference transcriptome to (1) quantify transcriptome-wide associations of SNPs with climate indicative of demographic responses to climate, (2) identify SNPs especially associated with climate and thus potential targets of natural selection, and (3) test the hypothesis that genetic diversity is high in climate-adaptive candidate genes. Constrained ordinations (redundancy analysis) and variance partitioning showed that genetic structure in Q. lobata was explained by spatial location (49 %) and climate (24 %), especially minimum temperature and summer/spring precipitation balance, suggesting that climate influences neutral demographic processes and gene flow. After accounting for underlying structure, individual-based environmental association analyses identified 79 SNPs from 49 transcripts as candidates under natural selection by climate. These candidate genes had significantly higher SNP rates per base pair per locus (θ W), nucleotide diversity (π), and gene diversity (G) than non-candidate genes. These results provide preliminary support for the hypothesis that balancing selection maintains diversity in climate-adaptive genes. Climate has likely shaped both population demography and local adaptation in valley oak.

Craniometric analysis of European Upper Palaeolithic and Mesolithic samples supports discontinuity at the Last Glacial Maximum.

The Last Glacial Maximum (LGM) represents the most significant climatic event since the emergence of anatomically modern humans (AMH). In Europe, the LGM may have played a role in changing morphological features as a result of adaptive and stochastic processes. We use craniometric data to examine morphological diversity in pre- and post-LGM specimens. Craniometric variation is assessed across four periods—pre-LGM, late glacial, Early Holocene and Middle Holocene—using a large, well-dated, dataset. Our results show significant differences across the four periods, using a MANOVA on size-adjusted cranial measurements. A discriminant function analysis shows separation between pre-LGM and later groups. Analyses repeated on a subsample, controlled for time and location, yield similar results. The results are largely influenced by facial measurements and are most consistent with neutral demographic processes. These findings suggest that the LGM had a major impact on AMH populations in Europe prior to the Neolithic.

Sexual dimorphism in the size and shape of the os coxae and the effects of microevolutionary processes

Sexual dimorphism in the human pelvis has been studied widely for forensic purposes, but it is still unclear to what extent it varies among human populations. There is evidence that microevolutionary processes, both neutral (i.e., population history) and selective (e.g., thermoregulatory adaptation and size-related obstetrical constraints) contribute to explain pelvic variation among populations, but the extent to which these factors affect pelvic sexual dimorphism is unknown. In this study, I analyze sexual dimorphism of the os coxae in 20 globally distributed human populations, using 3D morphometric data to separate the size and shape components of sexual differences. After evaluating population differences in the degree and pattern of sexual dimorphism, I test for the effect of population history, climate, and body size in shaping global diversity. The results show that size and shape dimorphism follow different patterns. Coxal size dimorphism is generally quite consistent through populations, with males bigger than females, but it appears to be reduced in small-bodied populations, possibly in relation to obstetrically-related selective pressures for a spacious birth canal. Beyond a general species-wide pattern of shape dimorphism, commonly used for forensic sex determination, other aspects of sexual differences in coxal shape vary among human populations, reflecting the effects of neutral demographic processes and climatic adaptation.

Skull and limb morphology differentially track population history and environmental factors in the transition to agriculture in Europe

The Neolithic transition in Europe was a complex mosaic spatio-temporal process, involving both demic diffusion from the Near East and the cultural adoption of farming practices by indigenous hunter-gatherers. Previous analyses of Mesolithic hunter-gatherers and Early Neolithic farmers suggest that cranial shape variation preserves the population history signature of the Neolithic transition. However, the extent to which these same demographic processes are discernible in the postcranium is poorly understood. Here, for the first time, crania and postcranial elements from the same 11 prehistoric populations are analysed together in an internally consistent theoretical and methodological framework. Results show that while cranial shape reflects the population history differences between Mesolithic and Neolithic lineages, relative limb dimensions exhibit significant congruence with environmental variables such as latitude and temperature, even after controlling for geography and time. Also, overall limb size is found to be consistently larger in hunter-gatherers than farmers, suggesting a reduction in size related to factors other than thermoregulatory adaptation. Therefore, our results suggest that relative limb dimensions are not tracking the same demographic population history as the cranium, and point to the strong influence of climatic, dietary and behavioural factors in determining limb morphology, irrespective of underlying neutral demographic processes.

Human Pelvis and Long Bones Reveal Differential Preservation of Ancient Population History and Migration Out of Africa

One of the main events in the history of our species has been our expansion out of Africa. A clear signature of this expansion has been found on global patterns of neutral genetic variation, whereby a serial founder effect accompanied the colonization of new regions, in turn creating a wilhin-pupulation decrease in neutral genetic diversity with increasing distance from Africa. This same distinctive pattern has also been described for cranial and dental morphological variation in human populations distributed across the globe. Here, we used a data set of postcranial linear measurements for 30 globally distributed human populations, and a climatic data set of minimum annual temperature, maximum annual temperature, and precipitation in order to separate for the first time the relative effect of neutral demographic processes and climatic selection on four long (limb) bones (femur, tibia, radius, and humerus) versus the pelvic bones of the human appendicular skeleton. We implemented a stepwise regression procedure in which phenotypic variance is assumed to be affected by the iterative founder events that accompanied human expansion from Africa, as well as by climate. This model included, as independent factors, geographic distance from central Africa, the three climatic variables, and all possible interactions between the three climatic variables. We excluded all nonsignificant factors by backward stepwise elimination with the aim of identifying the minimal model significantly explaining variation in the phenotypic data. Our results indicate a sharp difference in the way the pelvis and the limb bones reflect the neutral signature of the out-of-Africa expansion. Consistent with previous analyses of the cranium and dentition, pelvic shape variation shows a significant within-population decrease with increasing distance from Africa. However, no such pattern could be found in the long bones. Rather, in the case of both the tibia and the femur, a significant relationship between population-level variance and minimum temperature was demonstrated. Hence, in the case of these limb bones, it is probable that the effects of climatic selection have obliterated the demographic signature of human dispersal from Africa. Our finding mat pelvic variation exhibits the neutral effects of demographic history suggests that consideration of this skeletal element might be used to shed light on factors of human population history, just as the cranium has done.

Human Pelvis and Long Bones Reveal Differential Preservation of Ancient Population History and Migration Out of Africa

One of the main events in the history of our species has been our expansion out of Africa. A clear signature of this expansion has been found on global patterns of neutral genetic variation, whereby a serial founder effect accompanied the colonization of new regions, in turn creating a within-population decrease in neutral genetic diversity with increasing distance from Africa. This same distinctive pattern has also been described for cranial and dental morphological variation in human populations distributed across the globe. Here, we used a data set of postcranial linear measurements for 30 globally distributed human populations, and a climatic data set of minimum annual temperature, maximum annual temperature, and precipitation in order to separate for the first time the relative effect of neutral demographic processes and climatic selection on four long (limb) bones (femur, tibia, radius, and humerus) versus the pelvic bones of the human appendicular skeleton. We implemented a stepwise regression procedure in which phenotypic variance is assumed to be affected by the iterative founder events that accompanied human expansion from Africa, as well as by climate. This model included, as independent factors, geographic distance from central Africa, the three climatic variables, and all possible interactions between the three climatic variables. We excluded all nonsignificant factors by backward stepwise elimination with the aim of identifying the minimal model significantly explaining variation in the phenotypic data. Our results indicate a sharp difference in the way the pelvis and the limb bones reflect the neutral signature of the out-of-Africa expansion. Consistent with previous analyses of the cranium and dentition, pelvic shape variation shows a significant within-population decrease with increasing distance from Africa. However, no such pattern could be found in the long bones. Rather, in the case of both the tibia and the femur, a significant relationship between population-level variance and minimum temperature was demonstrated. Hence, in the case of these limb bones, it is probable that the effects of climatic selection have obliterated the demographic signature of human dispersal from Africa. Our finding that pelvic variation exhibits the neutral effects of demographic history suggests that consideration of this skeletal element might be used to shed light on factors of human population history, just as the cranium has done.

Contrasting epidemic histories reveal pathogen-mediated balancing selection on class II MHC diversity in a wild songbird.

The extent to which pathogens maintain the extraordinary polymorphism at vertebrate Major Histocompatibility Complex (MHC) genes via balancing selection has intrigued evolutionary biologists for over half a century, but direct tests remain challenging. Here we examine whether a well-characterized epidemic of Mycoplasmal conjunctivitis resulted in balancing selection on class II MHC in a wild songbird host, the house finch (Carpodacus mexicanus). First, we confirmed the potential for pathogen-mediated balancing selection by experimentally demonstrating that house finches with intermediate to high multi-locus MHC diversity are more resistant to challenge with Mycoplasma gallisepticum. Second, we documented sequence and diversity-based signatures of pathogen-mediated balancing selection at class II MHC in exposed host populations that were absent in unexposed, control populations across an equivalent time period. Multi-locus MHC diversity significantly increased in exposed host populations following the epidemic despite initial compromised diversity levels from a recent introduction bottleneck in the exposed host range. We did not observe equivalent changes in allelic diversity or heterozygosity across eight neutral microsatellite loci, suggesting that the observations reflect selection rather than neutral demographic processes. Our results indicate that a virulent pathogen can exert sufficient balancing selection on class II MHC to rescue compromised levels of genetic variation for host resistance in a recently bottlenecked population. These results provide evidence for Haldane's long-standing hypothesis that pathogens directly contribute to the maintenance of the tremendous levels of genetic variation detected in natural populations of vertebrates.

Evolutionary forces shaping genomic islands of population differentiation in humans.

BackgroundLevels of differentiation among populations depend both on demographic and selective factors: genetic drift and local adaptation increase population differentiation, which is eroded by gene flow and balancing selection. We describe here the genomic distribution and the properties of genomic regions with unusually high and low levels of population differentiation in humans to assess the influence of selective and neutral processes on human genetic structure.MethodsIndividual SNPs of the Human Genome Diversity Panel (HGDP) showing significantly high or low levels of population differentiation were detected under a hierarchical-island model (HIM). A Hidden Markov Model allowed us to detect genomic regions or islands of high or low population differentiation.ResultsUnder the HIM, only 1.5% of all SNPs are significant at the 1% level, but their genomic spatial distribution is significantly non-random. We find evidence that local adaptation shaped high-differentiation islands, as they are enriched for non-synonymous SNPs and overlap with previously identified candidate regions for positive selection. Moreover there is a negative relationship between the size of islands and recombination rate, which is stronger for islands overlapping with genes. Gene ontology analysis supports the role of diet as a major selective pressure in those highly differentiated islands. Low-differentiation islands are also enriched for non-synonymous SNPs, and contain an overly high proportion of genes belonging to the 'Oncogenesis' biological process.ConclusionsEven though selection seems to be acting in shaping islands of high population differentiation, neutral demographic processes might have promoted the appearance of some genomic islands since i) as much as 20% of islands are in non-genic regions ii) these non-genic islands are on average two times shorter than genic islands, suggesting a more rapid erosion by recombination, and iii) most loci are strongly differentiated between Africans and non-Africans, a result consistent with known human demographic history.