From genomics to integrative species delimitation? The case study of the Indo-Pacific Pocillopora corals

With the advent of genomics, sequencing thousands of loci from hundreds of individuals now appears feasible at reasonable costs, allowing complex phylogenies to be resolved. This is particularly relevant for cnidarians, for which insufficient data is available due to the small number of currently available markers and obscures species boundaries. Difficulties in inferring gene trees and morphological incongruences further blur the study and conservation of these organisms. Yet, can genomics alone be used to delimit species? Here, focusing on the coral genus Pocillopora, whose colonies play key roles in Indo-Pacific reef ecosystems but have challenged taxonomists for decades, we explored and discussed the usefulness of multiple criteria (genetics, morphology, biogeography and symbiosis ecology) to delimit species of this genus. Phylogenetic inferences, clustering approaches and species delimitation methods based on genome-wide single-nucleotide polymorphisms (SNP) were first used to resolve Pocillopora phylogeny and propose genomic species hypotheses from 356 colonies sampled across the Indo-Pacific (western Indian Ocean, tropical southwestern Pacific and south-east Polynesia). These species hypotheses were then compared to other lines of evidence based on genetic, morphology, biogeography and symbiont associations. Out of 21 species hypotheses delimited by genomics, 13 were strongly supported by all approaches, while six could represent either undescribed species or nominal species that have been synonymised incorrectly. Altogether, our results support (1) the obsolescence of macromorphology (i.e., overall colony and branches shape) but the relevance of micromorphology (i.e., corallite structures) to refine Pocillopora species boundaries, (2) the relevance of the mtORF (coupled with other markers in some cases) as a diagnostic marker of most species, (3) the requirement of molecular identification when species identity of colonies is absolutely necessary to interpret results, as morphology can blur species identification in the field, and (4) the need for a taxonomic revision of the genus Pocillopora. These results give new insights into the usefulness of multiple criteria for resolving Pocillopora, and more widely, scleractinian species boundaries, and will ultimately contribute to the taxonomic revision of this genus and the conservation of its species.

Integrative characterization of genetic and phenotypic differentiation in an ant species complex with strong hierarchical population structure and low dispersal abilities

Low dispersal, occurrence of asexual reproduction and geographic discontinuity increase genetic differentiation between populations, which ultimately can lead to speciation. In this work, we used a multidisciplinary framework to characterize the genetic and phenotypic differentiation between and within two cryptic ant species with restricted dispersal, Cataglyphis cursor and C. piliscapa and used behavioral experiments to test for reproductive isolation. Their distribution is segregated by the Rhône River and they have been traditionally distinguished only by hair numbers, although a statistical assessment is still lacking. We found strong genetic (microsatellites, nuclear and mitochondrial sequences), morphological (number of hairs, tibia length, male genitalia) and chemical (cuticular hydrocarbons) differentiation not only between species but also among localities within species. However, inter-specific differentiation was slightly higher than intra-specific differentiation for most markers. Overall, this pattern could either reflect reproductive isolation or could result from a longer period of geographic isolation between species than among localities within species without necessarily involving reproductive isolation. Interestingly, our behavioral experiments showed an absence of mating between species associated to a higher aggressiveness of workers towards heterospecific males. This suggests that sexual selection may, at least partially, fuel reproductive isolation. We also showed that cuticular hydrocarbons, mtDNA sequences and number of hairs provide reliable criteria allowing species discrimination. Overall, this species complex offers a case study to further investigate varying stages of a speciation continuum by estimating reproductive isolation between pairs of localities varying by their level of genetic differentiation.

From rivers to ocean basins: The role of ocean barriers and philopatry in the genetic structuring of a cosmopolitan coastal predator

The Bull Shark (Carcharhinus leucas) faces varying levels of exploitation around the world due to its coastal distribution. Information regarding population connectivity is crucial to evaluate its conservation status and local fishing impacts. In this study, we sampled 922 putative Bull Sharks from 19 locations in the first global assessment of population structure of this cosmopolitan species. Using a recently developed DNA-capture approach (DArTcap), samples were genotyped for 3400 nuclear markers. Additionally, full mitochondrial genomes of 384 Indo-Pacific samples were sequenced. Reproductive isolation was found between and across ocean basins (eastern Pacific, western Atlantic, eastern Atlantic, Indo-West Pacific) with distinct island populations in Japan and Fiji. Bull Sharks appear to maintain gene flow using shallow coastal waters as dispersal corridors, whereas large oceanic distances and historical land-bridges act as barriers. Females tend to return to the same area for reproduction, making them more susceptible to local threats and an important focus for management actions. Given these behaviors, the exploitation of Bull Sharks from insular populations, such as Japan and Fiji, may instigate local decline that cannot readily be replenished by immigration, which can in turn affect ecosystem dynamics and functions. These data also supported the development of a genetic panel to ascertain the population of origin, which will be useful in monitoring the trade of fisheries products and assessing population-level impacts of this harvest.

Like a rolling stone: Colonization and migration dynamics of the gray reef shark (Carcharhinus amblyrhynchos)

Designing appropriate management plans requires knowledge of both the dispersal ability and what has shaped the current distribution of the species under consideration. Here, we investigated the evolutionary history of the endangered gray reef shark (Carcharhinus amblyrhynchos) across its range by sequencing thousands of RADseq loci in 173 individuals in the Indo-Pacific (IP). We first bring evidence of the occurrence of a range expansion (RE) originating close to the Indo-Australian Archipelago (IAA) where two stepping-stone waves (east and westward) colonized almost the entire IP. Coalescent modeling additionally highlighted a homogenous connectivity (Nm ~ 10 per generation) throughout the range, and isolation by distance model suggested the absence of barriers to dispersal despite the affinity of C. amblyrhynchos to coral reefs. This coincides with long-distance swims previously recorded, suggesting that the strong genetic structure at the IP scale (F ST ~ 0.56 between its ends) is the consequence of its broad current distribution and organization in a large number of demes. Our results strongly suggest that management plans for the gray reef shark should be designed on a range-wide rather than a local scale due to its continuous genetic structure. We further contrasted these results with those obtained previously for the sympatric but strictly lagoon-associated Carcharhinus melanopterus, known for its restricted dispersal ability. Carcharhinus melanopterus exhibits a similar RE dynamic but is characterized by a stronger genetic structure and a nonhomogeneous connectivity largely dependent on local coral reefs availability. This sheds new light on shark evolution, emphasizing the roles of IAA as source of biodiversity and of life-history traits in shaping the extent of genetic structure and diversity.

Ecological and biogeographic features shaped the complex evolutionary history of an iconic apex predator (Galeocerdo cuvier)

BACKGROUND

The tiger shark (Galeocerdo cuvier) is a large iconic marine predator inhabiting worldwide tropical and subtropical waters. So far, only mitochondrial markers and microsatellites studies have investigated its worldwide historical demography with inconclusive outcomes. Here, we assessed for the first time the genomic variability of tiger shark based on RAD-seq data for 50 individuals from five sampling sites in the Indo-Pacific (IP) and one in the Atlantic Ocean (AO) to decipher the extent of the species' global connectivity and its demographic history.

RESULTS

Clustering algorithms (PCA and NMF), F_ST and an approximate Bayesian computation framework revealed the presence of two clusters corresponding to the two oceanic basins. By modelling the two-dimensional site frequency spectrum, we tested alternative isolation/migration scenarios between these two identified populations. We found the highest support for a divergence time between the two ocean basins of ~ 193,000 years before present (B.P) and an ongoing but limited asymmetric migration ~ 176 times larger from the IP to the AO (Nm ~ 3.9) than vice versa (Nm ~ 0.02).

CONCLUSIONS

The two oceanic regions are isolated by a strong barrier to dispersal more permeable from the IP to the AO through the Agulhas leakage. We finally emphasized contrasting recent demographic histories for the two regions, with the IP characterized by a recent bottleneck around 2000 years B.P. and the AO by an expansion starting 6000 years B.P. The large differentiation between the two oceanic regions and the absence of population structure within each ocean basin highlight the need for two large management units and call for future conservation programs at the oceanic rather than local scale, particularly in the Indo-Pacific where the population is declining.

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.

One panel to rule them all: DArTcap genotyping for population structure, historical demography, and kinship analyses, and its application to a threatened shark

With recent advances in sequencing technology, genomic data are changing how important conservation management decisions are made. Applications such as Close-Kin Mark-Recapture demand large amounts of data to estimate population size and structure, and their full potential can only be realised through ongoing improvements in genotyping strategies. Here we introduce DArTcap, a cost-efficient method that combines DArTseq and sequence capture, and illustrate its use in a high resolution population analysis of Glyphis garricki, a rare, poorly known and threatened euryhaline shark. Clustering analyses and spatial distribution of kin pairs from four different regions across northern Australia and one in Papua New Guinea, representing its entire known range, revealed that each region hosts at least one distinct population. Further structuring is likely within Van Diemen Gulf, the region that included the most rivers sampled, suggesting additional population structuring would be found if other rivers were sampled. Coalescent analyses and spatially explicit modelling suggest that G. garricki experienced a recent range expansion during the opening of the Gulf of Carpentaria following the conclusion of the Last Glacial Maximum. The low migration rates between neighbouring populations of a species that is found only in restricted coastal and riverine habitats show the importance of managing each population separately, including careful monitoring of local and remote anthropogenic activities that may affect their environments. Overall we demonstrated how a carefully chosen SNP panel combined with DArTcap can provide highly accurate kinship inference and also support population structure and historical demography analyses, therefore maximising cost-effectiveness.

Coping with Pleistocene climatic fluctuations: Demographic responses in remote endemic reef fishes

Elucidating demographic history during the settlement of ecological communities is crucial for properly inferring the mechanisms that shape patterns of species diversity and their persistence through time. Here, we used genomic data and coalescent-based approaches to elucidate for the first time the demographic dynamics associated with the settlement by endemic reef fish fauna of one of the most remote peripheral islands of the Pacific Ocean, Rapa Nui (Easter Island). We compared the demographic history of nine endemic species in order to explore their demographic responses to Pleistocene climatic fluctuations. We found that species endemic to Rapa Nui share a common demographic history, as signatures of population expansions were retrieved for almost all of the species studied here, and synchronous demographic expansions initiated during the last glacial period were recovered for more than half of the studied species. These results suggest that eustatic fluctuations associated with Milankovitch cycles have played a central role in species demographic histories and in the final stage of the community assembly of many Rapa Nui reef fishes. Specifically, sea level lowstands resulted in the maximum reef habitat extension for Rapa Nui endemic species; we discuss the potential role of seamounts in allowing endemic species to cope with Pleistocene climatic fluctuations, and we highlight the importance of local historical processes over regional ones. Overall, our results shed light on the mechanisms by which endemism arises and is maintained in peripheral reef fish fauna.

Urbanization without isolation: the absence of genetic structure among cities and forests in the tiny acorn ant Temnothorax nylanderi

Urban alteration of neutral and adaptive evolutionary processes is still underexplored. Using a genome-wide SNP dataset, we investigated (i) urban-induced modifications of population demography, genetic diversity and population structure and (ii) signature of divergent selection between urban and forest populations in the ant species, Temnothorax nylanderi. Our results did not reveal an impact of urbanization on neutral processes since we observed: (i) analogous genetic diversity among paired urban/forest sites and two control populations; (ii) weak population genetic structure explained neither by habitat (urban versus forest) nor by geography; (iii) a remarkably similar demographic history across populations with an ancestral growth followed by a recent decline, regardless of their current habitat or geographical location. The micro-geographical home range of ants may explain their resilience to urbanization. Finally, we detected 19 candidate loci discriminating urban/forest populations and associated with core cellular components, molecular function or biological process. Two of these loci were associated with a gene ontology term that was previously found to belong to a module of co-expressed genes related to caste phenotype. These results call for transcriptomics analyses to identify genes associated with ant social traits and to infer their potential role in urban adaptation.

Out of Africa: demographic and colonization history of the Algerian mouse (Mus spretus Lataste)

North Africa is now recognized as a major area for the emergence and dispersal of anatomically modern humans from at least 315 kya. The Mediterranean Basin is thus particularly suited to study the role of climate versus human-mediated changes on the evolutionary history of species. The Algerian mouse (Mus spretus Lataste) is an endemic species from this basin, with its distribution restricted to North Africa (from Libya to Morocco), Iberian Peninsula and South of France. A rich paleontological record of M. spretus exists in North Africa, suggesting hypotheses concerning colonization pathways, and the demographic and morphologic history of this species. Here we combined genetic (3 mitochondrial DNA loci and 18 microsatellites) and climatic niche modeling data to infer the evolutionary history of the Algerian mouse. We collected 646 new individuals in 51 localities. Our results are consistent with an anthropogenic translocation of the Algerian mouse from North Africa to the Iberian Peninsula via Neolithic navigators, probably from the Tingitane Peninsula. Once arrived in Spain, suitable climatic conditions would then have favored the dispersion of the Algerian mice to France. The morphological differentiation observed between Spanish, French and North African populations could be explained by a founder effect and possibly local adaptation. This article helps to better understand the role of climate versus human-mediated changes on the evolutionary history of mammal species in the Mediterranean Basin.

Survival and divergence in a small group: The extraordinary genomic history of the endangered Apennine brown bear stragglers

About 100 km east of Rome, in the central Apennine Mountains, a critically endangered population of ∼50 brown bears live in complete isolation. Mating outside this population is prevented by several 100 km of bear-free territories. We exploited this natural experiment to better understand the gene and genomic consequences of surviving at extremely small population size. We found that brown bear populations in Europe lost connectivity since Neolithic times, when farming communities expanded and forest burning was used for land clearance. In central Italy, this resulted in a 40-fold population decline. The overall genomic impact of this decline included the complete loss of variation in the mitochondrial genome and along long stretches of the nuclear genome. Several private and deleterious amino acid changes were fixed by random drift; predicted effects include energy deficit, muscle weakness, anomalies in cranial and skeletal development, and reduced aggressiveness. Despite this extreme loss of diversity, Apennine bear genomes show nonrandom peaks of high variation, possibly maintained by balancing selection, at genomic regions significantly enriched for genes associated with immune and olfactory systems. Challenging the paradigm of increased extinction risk in small populations, we suggest that random fixation of deleterious alleles (i) can be an important driver of divergence in isolation, (ii) can be tolerated when balancing selection prevents random loss of variation at important genes, and (iii) is followed by or results directly in favorable behavioral changes.

Population expansions dominate demographic histories of endemic and widespread Pacific reef fishes

Despite the unique nature of endemic species, their origin and population history remain poorly studied. We investigated the population history of 28 coral reef fish species, close related, from the Gambier and Marquesas Islands, from five families, with range size varying from widespread to small-range endemic. We analyzed both mitochondrial and nuclear sequence data using neutrality test and Bayesian analysis (EBSP and ABC). We found evidence for demographic expansions for most species (24 of 28), irrespective of range size, reproduction strategy or archipelago. The timing of the expansions varied greatly among species, from 8,000 to 2,000,000 years ago. The typical hypothesis for reef fish that links population expansions to the Last Glacial Maximum fit for 14 of the 24 demographic expansions. We propose two evolutionary processes that could lead to expansions older than the LGM: (a) we are retrieving the signature of an old colonization process for widespread, large-range endemic and paleoendemic species or (b) speciation; the expansion reflects the birth of the species for neoendemic species. We show for the first time that the demographic histories of endemic and widespread reef fish are not distinctly different and suggest that a number of processes drive endemism.

On the role played by the carrying capacity and the ancestral population size during a range expansion

Most species are structured and various population genetics models have been proposed to investigate their history. For mathematical tractability, most of these models make the simplifying assumption of equilibrium. Here we focus on the properties of a nonequilibrium spatial explicit model, range expansions (REs). Despite their abundance, many details of their genetic consequences need yet to be fully investigated. The model we studied is characterized by four main parameters: the effective population size of each deme (N), the migration rate per generation per deme (m), the time of the expansion (T_exp) and the effective size of the deme from which the expansion started (N_anc). By means of extensive coalescent simulations, we focused on two aspects of range expansions for fixed Nm: (1) the separate influence of N and m and (2) the role of N_anc. We compared our results with an equilibrium stepping stone model and found two main features typical of REs: an excess of rare variants for larger N and a complex interaction between N, T_exp and N_anc in shaping the degree of population differentiation (which depends only on Nm in the stepping stone model). Finally, we developed an approximate Bayesian computation approach to jointly estimate N and m and to infer N_anc. When applied to pseudo-observed data sets, we could correctly recover both N and m (but not N_anc), provided a large number of demes were sampled. These findings highlight how it will be possible to estimate the dispersal rate in nonequilibrium metapopulations by using population genetics approaches.

Inferring Population Size History from Large Samples of Genome-Wide Molecular Data - An Approximate Bayesian Computation Approach

Inferring the ancestral dynamics of effective population size is a long-standing question in population genetics, which can now be tackled much more accurately thanks to the massive genomic data available in many species. Several promising methods that take advantage of whole-genome sequences have been recently developed in this context. However, they can only be applied to rather small samples, which limits their ability to estimate recent population size history. Besides, they can be very sensitive to sequencing or phasing errors. Here we introduce a new approximate Bayesian computation approach named PopSizeABC that allows estimating the evolution of the effective population size through time, using a large sample of complete genomes. This sample is summarized using the folded allele frequency spectrum and the average zygotic linkage disequilibrium at different bins of physical distance, two classes of statistics that are widely used in population genetics and can be easily computed from unphased and unpolarized SNP data. Our approach provides accurate estimations of past population sizes, from the very first generations before present back to the expected time to the most recent common ancestor of the sample, as shown by simulations under a wide range of demographic scenarios. When applied to samples of 15 or 25 complete genomes in four cattle breeds (Angus, Fleckvieh, Holstein and Jersey), PopSizeABC revealed a series of population declines, related to historical events such as domestication or modern breed creation. We further highlight that our approach is robust to sequencing errors, provided summary statistics are computed from SNPs with common alleles.

Molecular data sampled from extant individuals contains considerable information about their demographic history. In particular, one classical question in population genetics is to reconstruct past population size changes from such data. Relating these changes to various climatic, geological or anthropogenic events allows characterizing the main factors driving genetic diversity and can have major outcomes for conservation. Until recently, mostly very simple histories, including one or two population size changes, could be estimated from genetic data. This has changed with the sequencing of entire genomes in many species, and several methods allow now inferring complex histories consisting of several tens of population size changes. However, analyzing entire genomes, while accounting for recombination, remains a statistical and numerical challenge. These methods, therefore, can only be applied to small samples with a few diploid genomes. We overcome this limitation by using an approximate estimation approach, where observed genomes are summarized using a small number of statistics related to allele frequencies and linkage disequilibrium. In contrast to previous approaches, we show that our method allows us to reconstruct also the most recent part (the last 100 generations) of the population size history. As an illustration, we apply it to large samples of whole-genome sequences in four cattle breeds.

Population genetic history of Aristeus antennatus (Crustacea: Decapoda) in the Western and Central Mediterranean Sea

Aristeus antennatus is an ecologically and economically important deep-water species in the Mediterranean Sea. In this study we investigated the genetic variability of A. antennatus sampled from 10 sampling stations in the Western and Central Mediterranean. By comparing our new samples with available data from the Western area, we aim to identify potential genetic stocks of A. antennatus and to reconstruct its historical demography in the Mediterranean. We analyzed two regions of mitochondrial DNA in 319 individuals, namely COI and 16S. We found two main results: i) the genetic diversity values consistent with previous data within the Mediterranean and the absence of barriers to gene flow within the Mediterranean Sea; ii) a constant long-term effective population size in almost all demes but a strong signature of population expansion in the pooled sample about 50,000 years B.P./ago. We propose two explanation for our results. The first is based on the ecology of A. antennatus. We suggest the existence of a complex meta-population structured into two layers: a deeper-dwelling stock, not affected by fishing, which preserves the pattern of historical demography; and genetically homogeneous demes inhabiting the fishing grounds. The larval dispersal, adult migration and continuous movements of individuals from “virgin” deeper grounds not affected by fishing to upper fishing areas support an effective ‘rescue effect’ contributing to the recovery of the exploited stocks and explain their genetic homogeneity throughout the Mediterranean Sea. The second is based on the reproduction model of this shrimp: the high variance in offspring production calls for a careful interpretation of the data observed under classical population genetics and Kingman’s coalescent. In both cases, management policies for A. antennatus will therefore require careful evaluation of the meta-population dynamics of all stocks in the Mediterranean. In the future, it will be particularly relevant to sample the deepest ones directly.

Evolutionary history and global spread of the Mycobacterium tuberculosis Beijing lineage

Mycobacterium tuberculosis strains of the Beijing lineage are globally distributed and are associated with the massive spread of multidrug-resistant (MDR) tuberculosis in Eurasia. Here we reconstructed the biogeographical structure and evolutionary history of this lineage by genetic analysis of 4,987 isolates from 99 countries and whole-genome sequencing of 110 representative isolates. We show that this lineage initially originated in the Far East, from where it radiated worldwide in several waves. We detected successive increases in population size for this pathogen over the last 200 years, practically coinciding with the Industrial Revolution, the First World War and HIV epidemics. Two MDR clones of this lineage started to spread throughout central Asia and Russia concomitantly with the collapse of the public health system in the former Soviet Union. Mutations identified in genes putatively under positive selection and associated with virulence might have favored the expansion of the most successful branches of the lineage.

Genetic consequences of habitat fragmentation during a range expansion

We investigate the effect of habitat fragmentation on the genetic diversity of a species experiencing a range expansion. These two evolutionary processes have not been studied yet, at the same time, owing to the difficulties of deriving analytic results for non-equilibrium models. Here we provide a description of their interaction by using extensive spatial and temporal coalescent simulations and we suggest guidelines for a proper genetic sampling to detect fragmentation. To model habitat fragmentation, we simulated a two-dimensional lattice of demes partitioned into groups (patches) by adding barriers to dispersal. After letting a population expand on this grid, we sampled lineages from the lattice at several scales and studied their coalescent history. We find that in order to detect fragmentation, one needs to extensively sample at a local level rather than at a landscape level. This is because the gene genealogy of a scattered sample is less sensitive to the presence of genetic barriers. Considering the effect of temporal changes of fragmentation intensities, we find that at least 10, but often >100, generations are needed to affect local genetic diversity and population structure. This result explains why recent habitat fragmentation does not always lead to detectable signatures in the genetic structure of populations. Finally, as expected, long-distance dispersal increases local genetic diversity and decreases levels of population differentiation, efficiently counteracting the effects of fragmentation.

Investigating sex-specific dynamics using uniparental markers: West New Guinea as a case study

Mitochondrial DNA (mtDNA) and Y chromosome (NRY) genetic markers have been often contrasted to investigate sex-specific dynamics. Traditionally, isolation by distance, intrapopulation genetic diversity and population differentiation are estimated from both markers and compared. Two possible sources of bias are often neglected. First, kilometric distances are frequently used as predictor of the connectivity between groups, hiding the role played by environmental features at a microgeographic scale. Second, the comparison of intrapopulation diversity and population differentiation between mtDNA and NRY is hampered by their different mutational mechanisms and rates. Here, we show how to account for these biases by analyzing from a different perspective a published dataset of eight West New Guinea (WNG) populations for which mtDNA control region sequences and seven linked NRY microsatellites had been typed. First, we modeled the connectivity among sampled populations by computing the number of days required to travel between groups. Then, we investigated the differences between the two sexes accounting for the molecular characteristics of the markers examined to obtain estimates on the product of the effective population size and the migration rate among demes (Nm). We achieved this goal by studying the shape of the gene genealogy at several sampling levels and using spatial explicit simulations. Both the direction and the rate of migration differ between male and females, with an Nm estimated to be >6 times higher in the latter under many evolutionary scenarios. We finally highlight the importance of applying metapopulation models when analyzing the genetic diversity of a species. We have applied the prediction of the sampling theory in a meta-population and we have corroborated our finding using spatial explicit simulations. Both approaches are fundamentally meant to deal with structured populations: we strongly believe in the importance of tacking structure into account when inferring the demographic history of a species.

The integration of multiple independent data reveals an unusual response to Pleistocene climatic changes in the hard tick Ixodes ricinus

In the last few years, improved analytical tools and the integration of genetic data with multiple sources of information have shown that temperate species exhibited more complex responses to ice ages than previously thought. In this study, we investigated how Pleistocene climatic changes affected the current distribution and genetic diversity of European populations of the tick Ixodes ricinus, an ectoparasite with high ecological plasticity. We first used mitochondrial and nuclear genetic markers to investigate the phylogeographic structure of the species and its Pleistocene history using coalescent-based methods; then we used species distribution modelling to infer the climatic niche of the species at last glacial maximum; finally, we reviewed the literature on the I. ricinus hosts to identify the locations of their glacial refugia. Our results support the scenario that during the last glacial phase, I. ricinus never experienced a prolonged allopatric divergence in separate glacial refugia, but persisted with interconnected populations across Southern and Central Europe. The generalist behaviour in host choice of I. ricinus would have played a major role in maintaining connections between its populations. Although most of the hosts persisted in separate refugia, from the point of view of I. ricinus, they represented a continuity of 'bridges' among populations. Our study highlights the importance of species-specific ecology in affecting responses to Pleistocene glacial-interglacial cycles. Together with other cases in Europe and elsewhere, it contributes to setting new hypotheses on how species with wide ecological plasticity coped with Pleistocene climatic changes.

Geographical delimitation of a partial selective sweep in African Drosophila melanogaster

Positive selection leaves characteristic footprints on DNA variation but detecting such patterns is challenging as the age, the intensity and the mode of selection as well as demography and evolutionary parameters (mutation and recombination rates) all play roles and these are difficult to disentangle. We recorded nucleotide variation in a sample of isogenic chromosomes from a western African population of Drosophila melanogaster at a locus (Fbp2) for which a partial selective sweep had previously been reported. We compared this locus to four other genes from the same chromosomes and from a European and an East African population. Then, we assessed Fbp2 variation in a sample of 370 chromosomes covering a comprehensive geographic sampling of 16 African localities. The signature of selection was tested while accounting for the demographic history of the populations. We found a significant signal of selection in two West African localities including Ivory Coast. Variation at Fpb2 would thus represent a case of an ongoing selective sweep in the range of this species. A weaker, nonsignificant, signal of selection was, however, apparent in some other populations, thus leaving open several possibilities: (i) the selective sweep originated in Ivory Coast and has spread to the rest of the continent; (ii) several African populations report the signature of a selective event having occurred in an ancestral population; (iii) this genome region is subject to independent selective events in African populations; and (iv) A neutral scenario with population subdivision and local bottleneck cannot be fully excluded to explain the molecular patterns observed in some populations.

ABC as a flexible framework to estimate demography over space and time: some cons, many pros

The analysis of genetic variation to estimate demographic and historical parameters and to quantitatively compare alternative scenarios recently gained a powerful and flexible approach: the Approximate Bayesian Computation (ABC). The likelihood functions does not need to be theoretically specified, but posterior distributions can be approximated by simulation even assuming very complex population models including both natural and human-induced processes. Prior information can be easily incorporated and the quality of the results can be analysed with rather limited additional effort. ABC is not a statistical analysis per se, but rather a statistical framework and any specific application is a sort of hybrid between a simulation and a data-analysis study. Complete software packages performing the necessary steps under a set of models and for specific genetic markers are already available, but the flexibility of the method is better exploited combining different programs. Many questions relevant in ecology can be addressed using ABC, but adequate amount of time should be dedicated to decide among alternative options and to evaluate the results. In this paper we will describe and critically comment on the different steps of an ABC analysis, analyse some of the published applications of ABC and provide user guidelines.

Population dynamic of the extinct European aurochs: genetic evidence of a north-south differentiation pattern and no evidence of post-glacial expansion

Background

The aurochs (Bos primigenius) was a large bovine that ranged over almost the entirety of the Eurasian continent and North Africa. It is the wild ancestor of the modern cattle (Bos taurus), and went extinct in 1627 probably as a consequence of human hunting and the progressive reduction of its habitat. To investigate in detail the genetic history of this species and to compare the population dynamics in different European areas, we analysed Bos primigenius remains from various sites across Italy.

Results

Fourteen samples provided ancient DNA fragments from the mitochondrial hypervariable region. Our data, jointly analysed with previously published sequences, support the view that Italian aurochsen were genetically similar to modern bovine breeds, but very different from northern/central European aurochsen. Bayesian analyses and coalescent simulations indicate that the genetic variation pattern in both Italian and northern/central European aurochsen is compatible with demographic stability after the last glaciation. We provide evidence that signatures of population expansion can erroneously arise in stable aurochsen populations when the different ages of the samples are not taken into account.

Conclusions

Distinct groups of aurochsen probably inhabited Italy and northern/central Europe after the last glaciation, respectively. On the contrary, Italian and Fertile Crescent aurochsen likely shared several mtDNA sequences, now common in modern breeds. We argue that a certain level of genetic homogeneity characterized aurochs populations in Southern Europe and the Middle East, and also that post-glacial recolonization of northern and central Europe advanced, without major demographic expansions, from eastern, and not southern, refugia.

The genetic impact of translocations and habitat fragmentation in chamois (Rupicapra) spp

The chamois is a useful species with which to investigate the combined genetic impact of habitat fragmentation, over hunting, and translocations. Genetic variation within and between chamois (genus Rupicapra) populations was analyzed in 259 individuals from 16 sampling sites located in Italy, Spain, Slovakia, and the Czech Republic. Two mitochondrial DNA markers (control region and cytochrome b) and 11 nuclear microsatellites were typed. The principal results of this study can be summarized as follows: 1) high and significant differentiation between almost all chamois populations is observed even on a microgeographical scale, probably caused by the patchy distribution of this species, sharp geographical barriers to gene flow, and drift effects related to recent bottlenecks; 2) historical translocation events have left a clear genetic signature, including interspecific hybridization in some Alpine localities; 3) the Apennine subspecies of chamois, Rupicapra pyrenaica ornata, shows a high and similar level of divergence (about 1.5 My) from the Pyrenean (Rupicapra pyrenaica pyrenaica) and the Alpine (Rupicapra rupicapra) chamois; therefore, the specific status of these taxa should be revised. These results confirm the potential of population genetic analyses to dissect and interpret complex patterns of diversity in order to define factors important to conservation and management.

Disentangling the effects of recombination, selection, and demography on the genetic variation at a major histocompatibility complex class II gene in the alpine chamois

The major histocompatibility complex (MHC) harbours some of the most polymorphic loci in vertebrate genomes. MHC genes are thought to be subject to some form of balancing selection, most likely pathogen-mediated selection. Hence, MHC genes are excellent candidates for exploring adaptive processes. In this study, we investigated the genetic variation at exon 2 of the DRB class II MHC locus in 191 alpine chamois (Rupicapra rupicapra) from 10 populations in the eastern Alps of Italy. In particular, we were interested in distinguishing and estimating the relative impact of selective and demographic factors, while taking into account the confounding effect of recombination. The extremely high d(n)/d(s) ratio and the presence of trans-species polymorphisms suggest that a strong long-term balancing selection effect has been operating at this locus throughout the evolutionary history of this species. We analysed patterns of genetic variation within and between populations, and the mitochondrial D-loop polymorphism patterns were analysed to provide a baseline indicator of the effects of demographic processes. These analyses showed that (i) the chamois experienced a demographic decline in the last 5000-30 000 years, most likely related to the postglacial elevation in temperature; (ii) this demographic process can explain the results of neutrality tests applied to MHC variation within populations, but cannot justify the much weaker divergence between populations implied by MHC as opposed to mitochondrial DNA; (iii) similar sets of divergent alleles are probably maintained with similar frequencies by balancing selection in different populations, and this mechanism is also operating in small isolated populations, which are strongly affected by drift.

The impact of the Austronesian expansion: evidence from mtDNA and Y chromosome diversity in the Admiralty Islands of Melanesia

The genetic ancestry of Polynesians can be traced to both Asia and Melanesia, which presumably reflects admixture occurring between incoming Austronesians and resident non-Austronesians in Melanesia before the subsequent occupation of the greater Pacific; however, the genetic impact of the Austronesian expansion to Melanesia remains largely unknown. We therefore studied the diversity of nonrecombining Y chromosomal (NRY) and mitochondrial (mt) DNA in the Admiralty Islands, located north of mainland Papua New Guinea, and updated our previous data from Asia, Melanesia, and Polynesia with new NRY markers. The Admiralties are occupied today solely by Austronesian-speaking groups, but their human settlement history goes back 20,000 years prior to the arrival of Austronesians about 3,400 years ago. On the Admiralties, we found substantial mtDNA and NRY variation of both Austronesian and non-Austronesian origins, with higher frequencies of Asian mtDNA and Melanesian NRY haplogroups, similar to previous findings in Polynesia and perhaps as a consequence of Austronesian matrilocality. Thus, the Austronesian language replacement on the Admiralties (and elsewhere in Island Melanesia and coastal New Guinea) was accompanied by an incomplete genetic replacement that is more associated with mtDNA than with NRY diversity. These results provide further support for the "Slow Boat" model of Polynesian origins, according to which Polynesian ancestors originated from East Asia but genetically mixed with Melanesians before colonizing the Pacific. We also observed that non-Austronesian groups of coastal New Guinea and Island Melanesia had significantly higher frequencies of Asian mtDNA haplogroups than of Asian NRY haplogroups, suggesting sex-biased admixture perhaps as a consequence of non-Austronesian patrilocality. We additionally found that the predominant NRY haplogroup of Asian origin in the Admiralties (O-M110) likely originated in Taiwan, thus providing the first direct Y chromosome evidence for a Taiwanese origin of the Austronesian expansion. Furthermore, we identified a NRY haplogroup (K-P79, also found on the Admiralties) in Polynesians that most likely arose in the Bismarck Archipelago, providing the first direct link between northern Island Melanesia and Polynesia. These results significantly advance our understanding of the impact of the Austronesian expansion and human history in the Pacific region.

Patterns of Y-Chromosome Diversity Intersect with the Trans-New Guinea Hypothesis

The island of New Guinea received part of the first human expansion out of Africa (>40,000 years ago), but its human genetic history remains poorly understood. In this study, we examined Y-chromosome diversity in 162 samples from the Bird's Head region of northwest New Guinea (NWNG) and compared the results with previously obtained data from other parts of the island. NWNG harbors a high level of cultural and linguistic diversity and is inhabited by non-Austronesian (i.e., Papuan)-speaking groups as well as harboring most of West New Guinea's (WNG) Austronesian-speaking groups. However, 97.5% of its Y-chromosomes belong to 5 haplogroups that originated in Melanesia; hence, the Y-chromosome diversity of NWNG (and, according to available data, of New Guinea as a whole) essentially reflects a local history. The remaining 2.5% belong to 2 haplogroups (O-M119 and O-M122) of East Asian origin, which were brought to New Guinea by Austronesian-speaking migrants around 3,500 years ago. Thus, the Austronesian expansion had only a small impact on shaping Y-chromosome diversity in NWNG, although the linguistic impact of this expansion to this region was much higher. In contrast, the expansion of Trans-New Guinea (TNG) speakers (non-Austronesian) starting about 6,000-10,000 years ago from the central highlands of what is now Papua New Guinea, presumably in combination with the expansion of agriculture, played a more important role in determining the Y-chromosome diversity of New Guinea. In particular, we identified 2 haplogroups (M-P34 and K-M254) as suggestive markers for the TNG expansion, whereas 2 other haplogroups (C-M38 and K-M9) most likely reflect the earlier local Y-chromosome diversity. We propose that sex-biased differences in the social structure and cultural heritage of the people involved in the Austronesian and the TNG expansions played an important role (among other factors) in shaping the New Guinean Y-chromosome landscape.