Vicariance divergence and gene flow among islet populations of an endemic lizard

Small Island Effects on the Thermal Biology of the Endemic Mediterranean Lizard Podarcis gaigeae

Ectotherms are vastly affected by climatic conditions as they rely on external sources of heat to regulate their body temperature, and changes in their habitat thermal quality could seriously affect their overall biology. To overcome the problems of a thermally unfavorable habitat, lizards need to either adjust their thermoregulatory behavior or respond to directional selection and shift their preferred body temperatures. To assess the impact of habitat thermal quality on the thermoregulatory profile, we studied multiple islet and 'mainland' populations of the Skyros wall lizard Podarcis gaigeae, an endemic lacertid to Skyros Archipelago, Greece. We evaluated the effectiveness of thermoregulation (E) using the three main thermal parameters: body (Tb), operative (Te), and preferred (Tpref) temperatures. We first hypothesized that the spatial homogeneity, the scarcity of thermal shelters, and the exposure to higher winds on islets would result in more demanding climate conditions. Second, we anticipated that islet lizards would achieve higher E in response to the lower thermal quality therein. As hypothesized, thermal parameters differed between populations but not in the expected manner. Skyros 'mainland' habitats reached higher temperatures, had more intense fluctuations, and were of lower thermal quality. As a result, lizards showed higher accuracy, precision, and effectiveness of thermoregulation. Noteworthy, we found that lizards from different populations have shifted their thermal profile and preferred body temperatures to cope with the particular conditions prevailing in their habitats. The latter supports the labile view on the evolution of thermoregulation.

Pleistocene divergence in the absence of gene flow among populations of a viviparous reptile with intraspecific variation in sex determination

Polymorphisms can lead to genetic isolation if there is differential mating success among conspecifics divergent for a trait. Polymorphism for sex-determining system may fall into this category, given strong selection for the production of viable males and females and the low success of heterogametic hybrids when sex chromosomes differ (Haldane's rule). Here we investigated whether populations exhibiting polymorphism for sex determination are genetically isolated, using the viviparous snow skink Carinascincus ocellatus. While a comparatively high elevation population has genotypic sex determination, in a lower elevation population there is an additional temperature component to sex determination. Based on 11,107 SNP markers, these populations appear genetically isolated. "Isolation with Migration" analysis also suggests these populations diverged in the absence of gene flow, across a period encompassing multiple Pleistocene glaciations and likely greater geographic proximity of populations. However, further experiments are required to establish whether genetic isolation may be a cause or consequence of differences in sex determination. Given the influence of temperature on sex in one lineage, we also discuss the implications for the persistence of this polymorphism under climate change.

A Genomewide Scan for Genetic Structure and Demographic History of Two Closely Related Species, Rhododendron dauricum and R. mucronulatum (Rhododendron, Ericaceae)

Understanding the processes of divergence and speciation is an important task for evolutionary research, and climate oscillations play a pivotal role. We estimated the genetic structure and demographic history of two closely related species of Rhododendron, R. dauricum, and R. mucronulatum, distributed in northeastern China using 664,406 single nucleotide polymorphic loci of specific-locus amplified fragment sequencing (SLAF-seq) and 4 chloroplast DNA (cpDNA) fragments, sampling 376 individuals from 39 populations of these two species across their geographic distributions. The geographical distribution of cpDNA haplotypes revealed that R. dauricum and R. mucronulatum have different spatial genetic structures and haplotype diversity. Analysis of molecular variance (AMOVA) results showed that these two species have significant genetic differentiation and that the phylogeny demonstrates that these two species clustered a monophyletic group based on SLAF data, respectively, but not in cpDNA data. The evidence of significant gene flow was also detected from R. mucronulatum to R. dauricum. A deep divergence between the two species was observed and occurred during the early Oligocene. The niche models showed that the two species have different demographic histories. Thus, our results imply that geography and climate changes played important roles in the evolutionary process of R. dauricum and R. mucronulatum, and although there was an interspecific gene flow, the divergence was maintained by natural selection.

Intraspecific diversification of the crop wild relative Brassica cretica Lam. using demographic model selection

BACKGROUND

Crop wild relatives (CWRs) contain genetic diversity, representing an invaluable resource for crop improvement. Many of their traits have the potential to help crops to adapt to changing conditions that they experience due to climate change. An impressive global effort for the conservation of various CWR will facilitate their use in crop breeding for food security. The genus Brassica is listed in Annex I of the International Treaty on Plant Genetic Resources for Food and Agriculture. Brassica oleracea (or wild cabbage), a species native to southern and western Europe, has become established as an important human food crop plant because of its large reserves stored over the winter in its leaves. Brassica cretica Lam. (Bc) is a CWR in the brassica group and B. cretica subsp. nivea (Bcn) has been suggested as a separate subspecies. The species Bc has been proposed as a potential gene donor to brassica crops, including broccoli, cabbage, cauliflower, oilseed rape, etc. RESULTS: We sequenced genomes of four Bc individuals, including two Bcn and two Bc. Demographic analysis based on our whole-genome sequence data suggests that populations of Bc are not isolated. Classification of the Bc into distinct subspecies is not supported by the data. Using only the non-coding part of the data (thus, the parts of the genome that has evolved nearly neutrally), we find the gene flow between different Bc population is recent and its genomic diversity is high.

CONCLUSIONS

Despite predictions on the disruptive effect of gene flow in adaptation, when selection is not strong enough to prevent the loss of locally adapted alleles, studies show that gene flow can promote adaptation, that local adaptations can be maintained despite high gene flow, and that genetic architecture plays a fundamental role in the origin and maintenance of local adaptation with gene flow. Thus, in the genomic era it is important to link the selected demographic models with the underlying processes of genomic variation because, if this variation is largely selectively neutral, we cannot assume that a diverse population of crop wild relatives will necessarily exhibit the wide-ranging adaptive diversity required for further crop improvement.

Elucidating the drivers of genetic differentiation in Malaysian torrent frogs (Anura: Ranidae: Amolops): a landscape genomics approach

The interplay between environmental attributes and evolutionary processes can provide valuable insights into how biodiversity is generated, partitioned and distributed. This study investigates the role of spatial, environmental and historical factors that could potentially drive diversification and shape genetic variation in Malaysian torrent frogs. Torrent frogs are ecologically conserved, and we hypothesize that this could impose tight constraints on dispersal routes, gene flow and consequently genetic structure. Moreover, levels of gene flow were shown to vary among populations from separate mountain ranges, indicating that genetic differentiation could be influenced by landscape features. Using genome-wide single nucleotide polymorphisms, in conjunction with landscape variables derived from Geographic Information Systems, we performed distance-based redundancy analyses and variance partitioning to disentangle the effects of isolation-by-distance (IBD), isolation-by-resistance (IBR) and isolation-by-colonization (IBC). Our results demonstrated that IBR contributed minimally to genetic variation. Intraspecific population structure can be largely attributed to IBD, whereas interspecific diversification was primarily driven by IBC. We also detected two distinct population bottlenecks, indicating that speciation events were likely driven by vicariance or founder events.

Spotlight on islands: on the origin and diversification of an ancient lineage of the Italian wall lizard Podarcis siculus in the western Pontine Islands

Groups of proximate continental islands may conceal more tangled phylogeographic patterns than oceanic archipelagos as a consequence of repeated sea level changes, which allow populations to experience gene flow during periods of low sea level stands and isolation by vicariant mechanisms during periods of high sea level stands. Here, we describe for the first time an ancient and diverging lineage of the Italian wall lizard Podarcis siculus from the western Pontine Islands. We used nuclear and mitochondrial DNA sequences of 156 individuals with the aim of unraveling their phylogenetic position, while microsatellite loci were used to test several a priori insular biogeographic models of migration with empirical data. Our results suggest that the western Pontine populations colonized the islands early during their Pliocene volcanic formation, while populations from the eastern Pontine Islands seem to have been introduced recently. The inter-island genetic makeup indicates an important role of historical migration, probably due to glacial land bridges connecting islands followed by a recent vicariant mechanism of isolation. Moreover, the most supported migration model predicted higher gene flow among islands which are geographically arranged in parallel. Considering the threatened status of small insular endemic populations, we suggest this new evolutionarily independent unit be given priority in conservation efforts.

Evidence for late Pleistocene origin of Astyanax mexicanus cavefish

BACKGROUND

Cavefish populations belonging to the Mexican tetra species Astyanax mexicanus are outstanding models to study the tempo and mode of adaptation to a radical environmental change. They are currently assigned to two main groups, the so-called "old" and "new" lineages, which would have populated several caves independently and at different times. However, we do not have yet accurate estimations of the time frames of evolution of these populations.

RESULTS

We reanalyzed the geographic distribution of mitochondrial and nuclear DNA polymorphisms and we found that these data do not support the existence of two cavefish lineages. Using IMa2, a program that allows dating population divergence in addition to demographic parameters, we found that microsatellite polymorphism strongly supports a very recent origin of cave populations (< 20,000 years). We identified a large number of single-nucleotide polymorphisms (SNPs) in transcript sequences of pools of embryos (Pool-seq) belonging to Pachón cave population and a surface population from Texas. Based on summary statistics that can be computed with this SNP data set together with simulations of evolution of SNP polymorphisms in two recently isolated populations, we looked for sets of demographic parameters that allow the computation of summary statistics with simulated populations that are similar to the ones with the sampled populations. In most simulations for which we could find a good fit between the summary statistics of observed and simulated data, the best fit occurred when the divergence between simulated populations was less than 30,000 years.

CONCLUSIONS

Although it is often assumed that some cave populations have a very ancient origin, a recent origin of these populations is strongly supported by our analyses of independent sets of nuclear DNA polymorphism. Moreover, the observation of two divergent haplogroups of mitochondrial and nuclear genes with different geographic distributions support a recent admixture of two divergent surface populations, before the isolation of cave populations. If cave populations are indeed only several thousand years old, many phenotypic changes observed in cavefish would thus have mainly involved the fixation of genetic variants present in surface fish populations and within a very short period of time.

Allopatric speciation despite historical gene flow: Divergence and hybridization in Carex furva and C. lucennoiberica (Cyperaceae) inferred from plastid and nuclear RAD-seq data

Gene flow among incipient species can act as a creative or destructive force in the speciation process, generating variation on which natural selection can act while, potentially, undermining population divergence. The flowering plant genus Carex exhibits a rapid and relatively recent radiation with many species limits still unclear. This is the case with the Iberian Peninsula (Spain and Portugal)-endemic C. lucennoiberica, which lay unrecognized within Carex furva until its recent description as a new species. In this study, we test how these species were impacted by interspecific gene flow during speciation. We sampled the full range of distribution of C. furva (15 individuals sampled) and C. lucennoiberica (88 individuals), sequenced two cpDNA regions (atpI-atpH, psbA-trnH) and performed genomic sequencing of 45,100 SNPs using restriction site-associated DNA sequencing (RAD-seq). We utilized a set of partitioned D-statistic tests and demographic analyses to study the degree and direction of introgression. Additionally, we modelled species distributions to reconstruct changes in range distribution during glacial and interglacial periods. Plastid, nuclear and morphological data strongly support divergence between species with subsequent gene flow. Combined with species distribution modelling, these data support a scenario of allopatry leading to species divergence, followed by secondary contact and gene flow due to long-distance dispersal and/or range expansions and contractions in response to Quaternary glacial cycles. We conclude that this is a case of allopatric speciation despite historical secondary contacts, which could have temporally influenced the speciation process, contributing to the knowledge of forces that are driving or counteracting speciation.

Closing the gap: Avian lineage splits at a young, narrow seaway imply a protracted history of mixed population response

The evolutionary significance of spatial habitat gaps has been well recognized since Alfred Russel Wallace compared the faunas of Bali and Lombok. Gaps between islands influence population structuring of some species, and flightless birds are expected to show strong partitioning even where habitat gaps are narrow. We examined the population structure of the most numerous living flightless land bird in New Zealand, Weka (Gallirallus australis). We surveyed Weka and their feather lice in native and introduced populations using genetic data gathered from DNA sequences of mitochondrial genes and nuclear β-fibrinogen and five microsatellite loci. We found low genetic diversity among extant Weka population samples. Two genetic clusters were evident in the mtDNA from Weka and their lice, but partitioning at nuclear loci was less abrupt. Many formerly recognized subspecies/species were not supported; instead, we infer one subspecies for each of the two main New Zealand islands. Although currently range restricted, North Island Weka have higher mtDNA diversity than the more wide-ranging southern Weka. Mismatch and neutrality statistics indicate North Island Weka experienced rapid and recent population reduction, while South Island Weka display the signature of recent expansion. Similar haplotype data from a widespread flying relative of Weka and other New Zealand birds revealed instances of North Island-South Island partitioning associated with a narrow habitat gap (Cook Strait). However, contrasting patterns indicate priority effects and other ecological factors have a strong influence on spatial exchange at this scale.

Can genomic data alone tell us whether speciation happened with gene flow?

The allopatric model, which requires a period of geographical isolation for speciation to complete, has been the standard model in the modern era. Recently, "speciation with gene flow" has been widely discussed in relation to the model of "strict allopatry" and the level of DNA divergence across genomic regions. We wish to caution that genomic data by themselves may only permit the rejection of the simplest form of allopatry. Even a slightly more complex and realistic model that starts with subdivided populations would be impossible to reject by the genomic data alone. To resolve this central issue of speciation, other forms of observations such as the sequencing of reproductive isolation genes or the identification of geographical barrier(s) will be necessary.

The Aegean archipelago: a natural laboratory of evolution, ecology and civilisations

The Aegean archipelago, comprising numerous islands and islets with great heterogeneity in topographic, geological, historical and environmental properties, offers an ideal natural laboratory for ecological and evolutionary research, and has been the stage for a very long interaction between human civilizations and local ecosystems. This work presents insights that have been gained from past and current relevant research in the area, highlighting also the importance of the Aegean archipelago as a useful model to address many major questions in biogeography, ecology and evolutionary processes. Among the most interesting findings from such studies concern the role of habitat heterogeneity as the most important determinant of species richness, the development of a new model (Choros) for the species–area–habitats relationship, the mechanistic aspects of the Small Island Effect, the very high rates of species turnover, the lack of a role for interspecific competition in shaping species co-occurrence patterns in most cases, the importance of non adaptive radiation in diversification of several taxa, the insights into the relative roles of vicariance and dispersal in speciation, the understanding of the interplay between human presence and the establishment of exotic species and extinction of indigenous biotas. Concluding, the Aegean archipelago is an ideal stage for research in evolution, ecology and biogeography, and has the potential to become a model study area at a global level, especially for land-bridge, continental islands.

Phylogeographic model selection leads to insight into the evolutionary history of four-eyed frogs

Phylogeographic research investigates biodiversity at the interface between populations and species, in a temporal and geographic context. Phylogeography has benefited from analytical approaches that allow empiricists to estimate parameters of interest from the genetic data (e.g., θ = 4Neμ, population divergence, gene flow), and the widespread availability of genomic data allow such parameters to be estimated with greater precision. However, the actual inferences made by phylogeographers remain dependent on qualitative interpretations derived from these parameters' values and as such may be subject to overinterpretation and confirmation bias. Here we argue in favor of using an objective approach to phylogeographic inference that proceeds by calculating the probability of multiple demographic models given the data and the subsequent ranking of these models using information theory. We illustrate this approach by investigating the diversification of two sister species of four-eyed frogs of northeastern Brazil using single nucleotide polymorphisms obtained via restriction-associated digest sequencing. We estimate the composite likelihood of the observed data given nine demographic models and then rank these models using Akaike information criterion. We demonstrate that estimating parameters under a model that is a poor fit to the data is likely to produce values that lead to spurious phylogeographic inferences. Our results strongly imply that identifying which parameters to estimate from a given system is a key step in the process of phylogeographic inference and is at least as important as being able to generate precise estimates of these parameters. They also illustrate that the incorporation of model uncertainty should be a component of phylogeographic hypothesis tests.

Secondary contact seeds phenotypic novelty in cichlid fishes

Theory proposes that genomic admixture between formerly reproductively isolated populations can generate phenotypic novelty for selection to act upon. Secondary contact may therefore be a significant promoter of phenotypic novelty that allows species to overcome environmental challenges and adapt to novel environments, including during adaptive radiation. To date, this has largely been considered from the perspective of interspecific hybridization at contact zones. However, it is also possible that this process occurs more commonly between natural populations of a single species, and thus its importance in adaptive evolution may have been underestimated. In this study, we tested the consequences of genomic introgression during apparent secondary contact between phenotypically similar lineages of the riverine cichlid fish Astatotilapia calliptera. We provide population genetic evidence of a secondary contact zone in the wild, and then demonstrate using mate-choice experiments that both lineages can reproduce together successfully in laboratory conditions. Finally, we show that genomically admixed individuals display extreme phenotypes not observed in the parental lineages. Collectively, the evidence shows that secondary contact can drive the evolution of phenotypic novelty, suggesting that pulses of secondary contact may repeatedly seed genetic novelty, which when coupled with ecological opportunity could promote rapid adaptive evolution in natural circumstances.

An investigation for population maintenance mechanism in a miniature garden: genetic connectivity or independence of small islet populations of the Ryukyu five-lined skink

The Ryukyu five-lined skink (Plestiodon marginatus) is an island lizard that is even found in tiny islets with less than half a hectare of habitat area. We hypothesized that the island populations are maintained under frequent gene flow among the islands or independent of each other. To test our hypotheses, we investigated genetic structure of 21 populations from 11 land-bridge islands that were connected during the latest glacial age, and 4 isolated islands. Analyses using mitochondrial cytochrome b gene sequence (n = 67) and 10 microsatellite loci (n = 235) revealed moderate to high levels of genetic differentiation, existence of many private alleles/haplotypes in most islands, little contemporary migration, a positive correlation between genetic variability and island area, and a negative correlation between relatedness and island area. These evidences suggest a strong effect of independent genetic drift as opposed to gene flow, favoring the isolation hypothesis even in tiny islet populations. Isolation-by-distance effect was demonstrated and it became more prominent when the 4 isolated islands were excluded, suggesting that the pattern is a remnant of the land-bridge age. In a few island populations, however, the possibility of occasional overwater dispersals was partially supported and therefore could not be ruled out.

Repeated landmass reformation limits diversification in the widespread littoral zone mosquito Anopheles sundaicus sensu lato in the Indo-Oriental Region

Southeast Asia harbours abundant biodiversity, hypothesized to have been generated by Pliocene and Pleistocene climatic and environmental change. Vicariance between the island of Borneo, the remaining Indonesian archipelago and mainland Southeast Asia caused by elevated sea levels during interglacial periods has been proposed to lead to diversification in the littoral zone mosquito Anopheles (Cellia) sundaicus (Rodenwaldt) sensu lato. To test this biogeographical hypothesis, we inferred the population history and assessed gene flow of A. sundaicus s.l. sampled from 18 populations across its pan-Asian species range, using sequences from mitochondrial cytochrome c oxidase subunit 1 (CO1), the internal transcribed spacer 2 (ITS2) and the mannose phosphate isomerase (Mpi) gene. A hypothesis of ecological speciation for A. sundaicus involving divergent adaptation to brackish and freshwater larval habitats was also previously proposed, based on a deficiency of heterozygotes for Mpi allozyme alleles in sympatry. This hypothesis was not supported by Mpi sequence data, which exhibited no fixed differences between brackish and freshwater larval habitats. Mpi and CO1 supported the presence of up to eight genetically distinct population groupings. Counter to the hypothesis of three allopatric species, divergence was often no greater between Borneo, Sumatra/Java and the Southeast Asian mainland than it was between genetic groupings within these landmasses. An isolation-with-migration (IM) model indicates recurrent gene flow between the current major landmasses. Such gene flow would have been possible during glacial periods when the current landmasses merged, presenting opportunities for dispersal along expanding and contracting coastlines. Consequently, Pleistocene climatic variation has proved a homogenizing, rather than diversifying, force for A. sundaicus diversity.

Does relaxed predation drive phenotypic divergence among insular populations?

The evolution of striking phenotypes on islands is a well-known phenomenon, and there has been a long-standing debate on the patterns of body size evolution on islands. The ecological causes driving divergence in insular populations are, however, poorly understood. Reduced predator fauna is expected to lower escape propensity, increase body size and relax selection for crypsis in small-bodied, insular prey species. Here, we investigated whether escape behaviour, body size and dorsal coloration have diverged as predicted under predation release in spatially replicated islet and mainland populations of the lizard species Podarcis gaigeae. We show that islet lizards escape approaching observers at shorter distances and are larger than mainland lizards. Additionally, we found evidence for larger between-population variation in body size among the islet populations than mainland populations. Moreover, islet populations are significantly more divergent in dorsal coloration and match their respective habitats poorer than mainland lizards. These results strongly suggest that predation release on islets has driven population divergence in phenotypic and behavioural traits and that selective release has affected both trait means and variances. Relaxed predation pressure is therefore likely to be one of the major ecological factors driving body size divergence on these islands.

Mitochondrial and nuclear DNA survey of Zootoca vivipara across the eastern Italian Alps: evolutionary relationships, historical demography and conservation implications

The European common lizard Zootoca vivipara exhibits reproductive bimodality, with populations being either viviparous or oviparous. In the central-eastern Italian Alps oviparous populations (Z. v. carniolica) and viviparous populations (Z. v. vivipara) partly overlap geographically. Studying the evolutionary relationship between these taxa presents an interesting opportunity to gain insight into the evolution of this trait. We aim to: i) test whether Z. v. carniolica, which is endangered, constitutes an ESU (Evolutionary Significant Unity); ii) infer mtDNA divergence time between the Z. v. carniolica clade and all the other Z. vivipara subspecies with the aid of an external calibration point; and iii) describe the phylogeographical and demographic scenarios in the area. To do so we sequenced about 200 individuals for mitochondrial variation; 64 of them were also analysed for three nuclear genes. Furthermore, we analysed the same nuclear markers in 17 individuals from the other oviparous subspecies Z. v. louislantzi and 11 individuals of Z. v. vivipara from widespread geographical origins. The mtDNA and nDNA loci that we examined supported the monophyly of Z. v. carniolica. The mtDNA-based estimate of divergence time between Z. v. carniolica and all the other subspecies indicated a separation at 4.5 Mya (95% CI 6.1–2.6), with about 5% of sequence divergence. Considering that Z. v. carniolica harbours higher genetic diversity, while Z. v. vivipara from central-eastern Alps shows a signature of recent population and spatial expansion, we argue that Z. v. carniolica represents a distinct evolutionary unit, with a presumably long-term evolutionary history of separation. Z. v. carniolica populations, occurring at higher latitudes and altitudes than insofar supposed, live in peat bogs, a seriously threatened habitat: taking into account also its evolutionary distinctness, specific conservation measures should be considered.

Did postglacial sea-level changes initiate the evolutionary divergence of a Tasmanian endemic raptor from its mainland relative?

Populations on continental islands are often distinguishable from mainland conspecifics with respect to body size, appearance, behaviour or life history, and this is often congruent with genetic patterns. It is commonly assumed that such differences developed following the complete isolation of populations by sea-level rise following the Last Glacial Maximum (LGM). However, population divergence may predate the LGM, or marine dispersal and colonization of islands may have occurred more recently; in both cases, populations may have also diverged despite ongoing gene flow. Here, we test these alternative hypotheses for the divergence between wedge-tailed eagles from mainland Australia (Aquila audax audax) and the threatened Tasmanian subspecies (Aquila audax fleayi), based on variation at 20 microsatellite loci and mtDNA. Coalescent analyses indicate that population divergence appreciably postdates the severance of terrestrial habitat continuity and occurred without any subsequent gene flow. We infer a recent colonization of Tasmania by marine dispersal and cannot discount founder effects as the cause of differences in body size and life history. We call into question the general assumption of post-LGM marine transgression as the initiator of divergence of terrestrial lineages on continental islands and adjacent mainland, and highlight the range of alternative scenarios that should be considered.

Has the inbreeding load for a condition-dependent sexual signalling trait been purged in insular lizard populations?

Sexually selected traits are often condition-dependent and are expected to be affected by genome-wide distributed deleterious mutations and inbreeding. However, sexual selection is a powerful selective force that can counteract inbreeding through purging of deleterious mutations. Inbreeding and purging of the inbreeding load for sexually selected traits has rarely been studied across natural populations with different degrees of inbreeding. Here we investigate inbreeding effects (measured as marker-based heterozygosity) on condition-dependent sexually selected signalling trait and other morphological traits across islet- and mainland populations (n = 15) of an endemic lizard species (Podarcis gaigeae). Our data suggest inbreeding depression on a condition-dependent sexually selected signalling character among mainland subpopulations with low or intermediate levels of inbreeding, but no sign of inbreeding depression among small and isolated islet populations despite their higher overall inbreeding levels. In contrast, there was no such pattern among ten other morphological traits which are primarily naturally selected and presumably not involved in sexual signalling. These results are in line with purging of recessive deleterious alleles, or purging in combination with stochastic fixation of alleles by genetic drift, for a sexual signalling character in the islet environment, which is characterized by low population sizes and strong sexual selection. Higher clutch sizes in islet populations also raise interesting questions regarding the possibility of antagonistic pleiotropy. Purging and other non-exclusive explanations of our results are discussed.

Age of the association between Helicobacter pylori and man

When modern humans left Africa ca. 60,000 years ago (60 kya), they were already infected with Helicobacter pylori, and these bacteria have subsequently diversified in parallel with their human hosts. But how long were humans infected by H. pylori prior to the out-of-Africa event? Did this co-evolution predate the emergence of modern humans, spanning the species divide? To answer these questions, we investigated the diversity of H. pylori in Africa, where both humans and H. pylori originated. Three distinct H. pylori populations are native to Africa: hpNEAfrica in Afro-Asiatic and Nilo-Saharan speakers, hpAfrica1 in Niger-Congo speakers and hpAfrica2 in South Africa. Rather than representing a sustained co-evolution over millions of years, we find that the coalescent for all H. pylori plus its closest relative H. acinonychis dates to 88–116 kya. At that time the phylogeny split into two primary super-lineages, one of which is associated with the former hunter-gatherers in southern Africa known as the San. H. acinonychis, which infects large felines, resulted from a later host jump from the San, 43–56 kya. These dating estimates, together with striking phylogenetic and quantitative human-bacterial similarities show that H. pylori is approximately as old as are anatomically modern humans. They also suggest that H. pylori may have been acquired via a single host jump from an unknown, non-human host. We also find evidence for a second Out of Africa migration in the last 52,000 years, because hpEurope is a hybrid population between hpAsia2 and hpNEAfrica, the latter of which arose in northeast Africa 36–52 kya, after the Out of Africa migrations around 60 kya.

We previously showed that the population history of H. pylori may be used as a marker for human migrations, including the demonstration that humans carried H. pylori out of Africa 60,000 years ago during their recent global expansions. But how long were humans infected by H. pylori prior to the out-of-Africa event? Here we showed that chimpanzees in Central-East Africa do not possess Helicobacter-like bacteria, as would have been expected for pathogen-host co-evolution over millions of years. Using H. pylori gene sequences isolated from San, a group of click-speaking hunter-gatherers, and numerous other sources, we calculated that humans have been infected with H. pylori for at least 88,000–116,000 years. Phylogenetic comparisons showed similar evolutionary histories for human and H. pylori lineages and suggest that this association stemmed from a single host jump. We showed that hpAfrica2, the most divergent H. pylori population, arose in the San and that their progenitors were the source of H. acinonychis which was acquired by large felines approximately 50,000 years ago. Furthermore, our data provided clear evidence for a recent second exodus Out of Africa in the last 52,000 years which was essential for the formation of the hybrid population that currently infects Europeans.