Julien Foucaud

Clonal reproduction by males and females in the little fire ant

Sexual reproduction can lead to major conflicts between sexes and within genomes. Here we report an extreme case of such conflicts in the little fire ant Wasmannia auropunctata. We found that sterile workers are produced by normal sexual reproduction, whereas daughter queens are invariably clonally produced. Because males usually develop from unfertilized maternal eggs in ants and other haplodiploid species, they normally achieve direct fitness only through diploid female offspring. Hence, although the clonal production of queens increases the queens relatedness to reproductive daughters, it potentially reduces male reproductive success to zero. In an apparent response to this conflict between sexes, genetic analyses reveal that males reproduce clonally, most likely by eliminating the maternal half of the genome in diploid eggs. As a result, all sons have nuclear genomes identical to those of their father. The obligate clonal production of males and queens from individuals of the same sex effectively results in a complete separation of the male and female gene pools. These findings show that the haplodiploid sex-determination system provides grounds for the evolution of extraordinary genetic systems and new types of sexual conflict.

The effect of RAD allele dropout on the estimation of genetic variation within and between populations

Inexpensive short‐read sequencing technologies applied to reduced representation genomes is revolutionizing genetic research, especially population genetics analysis, by allowing the genotyping of massive numbers of single‐nucleotide polymorphisms (SNP) for large numbers of individuals and populations. Restriction site–associated DNA (RAD) sequencing is a recent technique based on the characterization of genomic regions flanking restriction sites. One of its potential drawbacks is the presence of polymorphism within the restriction site, which makes it impossible to observe the associated SNP allele (i.e. allele dropout, ADO). To investigate the effect of ADO on genetic variation estimated from RAD markers, we first mathematically derived measures of the effect of ADO on allele frequencies as a function of different parameters within a single population. We then used RAD data sets simulated using a coalescence model to investigate the magnitude of biases induced by ADO on the estimation of expected heterozygosity and FST under a simple demographic model of divergence between two populations. We found that ADO tends to overestimate genetic variation both within and between populations. Assuming a mutation rate per nucleotide between 10−9 and 10−8, this bias remained low for most studied combinations of divergence time and effective population size, except for large effective population sizes. Averaging FST values over multiple SNPs, for example, by sliding window analysis, did not correct ADO biases. We briefly discuss possible solutions to filter the most problematic cases of ADO using read coverage to detect markers with a large excess of null alleles.

Estimation of population allele frequencies from next-generation sequencing data: pool-versus individual-based genotyping

Molecular markers produced by next‐generation sequencing (NGS) technologies are revolutionizing genetic research. However, the costs of analysing large numbers of individual genomes remain prohibitive for most population genetics studies. Here, we present results based on mathematical derivations showing that, under many realistic experimental designs, NGS of DNA pools from diploid individuals allows to estimate the allele frequencies at single nucleotide polymorphisms (SNPs) with at least the same accuracy as individual‐based analyses, for considerably lower library construction and sequencing efforts. These findings remain true when taking into account the possibility of substantially unequal contributions of each individual to the final pool of sequence reads. We propose the intuitive notion of effective pool size to account for unequal pooling and derive a Bayesian hierarchical model to estimate this parameter directly from the data. We provide a user‐friendly application assessing the accuracy of allele frequency estimation from both pool‐ and individual‐based NGS population data under various sampling, sequencing depth and experimental error designs. We illustrate our findings with theoretical examples and real data sets corresponding to SNP loci obtained using restriction site–associated DNA (RAD) sequencing in pool‐ and individual‐based experiments carried out on the same population of the pine processionary moth (Thaumetopoea pityocampa). NGS of DNA pools might not be optimal for all types of studies but provides a cost‐effective approach for estimating allele frequencies for very large numbers of SNPs. It thus allows comparison of genome‐wide patterns of genetic variation for large numbers of individuals in multiple populations.

Empirical assessment of RAD sequencing for interspecific phylogeny.

Next-generation sequencing opened up new possibilities in phylogenetics; however, choosing an appropriate method of sample preparation remains challenging. Here, we demonstrate that restriction-site-associated DNA sequencing (RAD-seq) generates useful data for phylogenomics. Analysis of our RAD library using current bioinformatic and phylogenetic tools produced 400× more sites than our Sanger approach (2,262,825 nt/species), fully resolving relationships between 18 species of ground beetles (divergences up to 17 My). This suggests that RAD-seq is promising to infer phylogeny of eukaryotic species, though potential biases need to be evaluated and new methodologies developed to take full advantage of such data.

Worldwide invasion by the little fire ant: routes of introduction and eco-evolutionary pathways

Biological invasions are generally thought to occur after human aided migration to a new range. However, human activities prior to migration may also play a role. We studied here the evolutionary genetics of introduced populations of the invasive ant Wasmannia auropunctata at a worldwide scale. Using microsatellite markers, we reconstructed the main routes of introduction of the species. We found three main routes of introduction, each of them strongly associated to human history and trading routes. We also demonstrate the overwhelming occurrence of male and female clonality in introduced populations of W. auropunctata, and suggest that this particular reproduction system is under selection in human‐modified habitats. Together with previous researches focused on native populations, our results suggest that invasive clonal populations may have evolved within human modified habitats in the native range, and spread further from there. The evolutionarily most parsimonious scenario for the emergence of invasive populations of the little fire ant might thus be a two‐step process. The W. auropunctata case illustrates the central role of humans in biological change, not only due to changes in migration patterns, but also in selective pressures over species.

Ecologically heterogeneous populations of the invasive ant Wasmannia auropunctata within its native and introduced ranges

Abstract 1. The biology of most invasive species in their native geographical areas remains largely unknown. Such studies are, however, crucial in shedding light on the ecological and evolutionary processes underlying biological invasions.

RARE SEXUAL REPRODUCTION EVENTS IN THE CLONAL REPRODUCTION SYSTEM OF INTRODUCED POPULATIONS OF THE LITTLE FIRE ANT

Abstract A unique reproductive system has previously been described in Wasmannia auropunctata, a widespread invasive ant species, where males are produced clonally, female queens are parthenogens, and female workers are produced sexually. However, these findings were mostly based on samples originating from only a limited part of the native range of the species in South America. We used microsatellite markers to uncover the reproductive modes displayed by a large number of nests collected in various invasive W. auropunctata populations introduced 40 years ago into New Caledonia, where the species now forms a single 450‐km‐long supercolony. Although the main reproduction system in New Caledonia remained clonality for both male and female reproductives, we found evidence of rare sexual reproduction events that led to the production of both new queen and male clonal lineages. All clonal lineages observed in New Caledonia potentially derived from sexual reproduction, recombination, and mutation events from a single female and a single male genotype. Hence, the male and female gene pools are not strictly separated in New Caledonia and the two sexes do not follow independent evolutionary trajectories. Our results also suggest genetic determination for both parthenogenesis and caste. We discuss the evolutionary implications of the emergence of sex in the clonal reproduction system of introduced populations of W. auropunctata.

Where do adaptive shifts occur during invasion? A multidisciplinary approach to unravelling cold adaptation in a tropical ant species invading the Mediterranean area

Evolution may improve the invasiveness of populations, but it often remains unclear whether key adaptation events occur after introduction into the recipient habitat (i.e. post-introduction adaptation scenario), or before introduction within the native range (i.e. prior-adaptation scenario) or at a primary site of invasion (i.e. bridgehead scenario). We used a multidisciplinary approach to determine which of these three scenarios underlies the invasion of the tropical ant Wasmannia auropunctata in a Mediterranean region (i.e. Israel). Species distribution models (SDM), phylogeographical analyses at a broad geographical scale and laboratory experiments on appropriate native and invasive populations indicated that Israeli populations followed an invasion scenario in which adaptation to cold occurred at the southern limit of the native range before dispersal to Israel. We discuss the usefulness of combining SDM, genetic and experimental approaches for unambiguous determination of eco-evolutionary invasion scenarios.

Reproductive system, social organization, human disturbance and ecological dominance in native populations of the little fire ant, Wasmannia auropunctata.

The invasive ant species Wasmannia auropunctata displays both ecologically dominant and non‐dominant populations within its native range. Three factors could theoretically explain the ecological dominance of some native populations of W. auropunctata: (i) its clonal reproductive system, through demographic and/or adaptive advantages; (ii) its unicolonial social organization, through lower intraspecific and efficient interspecific competition; (iii) the human disturbance of its native range, through the modification of biotic and abiotic environmental conditions. We used microsatellite markers and behavioural tests to uncover the reproductive modes and social organization of dominant and non‐dominant native populations in natural and human‐modified habitats. Microsatellite and mtDNA data indicated that dominant and non‐dominant native populations (supercolonies as determined by aggression tests) of W. auropunctata did not belong to different evolutionary units. We found that the reproductive system and the social organization are neither necessary nor sufficient to explain W. auropunctata ecological dominance. Dominance rather seems to be set off by unknown ecological factors altered by human activities, as all dominant populations were recorded in human‐modified habitats. The clonal reproductive system found in some populations of W. auropunctata may however indirectly contribute to its ecological dominance by allowing the species to expand its environmental niche, through the fixation over time of specific combinations of divergent male and female genotypes. Unicoloniality may rather promote the range expansion of already dominant populations than actually trigger ecological dominance. The W. auropunctata model illustrates the strong impact of human disturbance on species’ ecological features and the adaptive potential of clonal reproductive systems.

Deciphering the routes of invasion of Drosophila suzukii by means of ABC random forest

Abstract Deciphering invasion routes from molecular data is crucial to understanding biological invasions, including identifying bottlenecks in population size and admixture among distinct populations. Here, we unravel the invasion routes of the invasive pest Drosophila suzukii using a multi-locus microsatellite dataset (25 loci on 23 worldwide sampling locations). To do this, we use approximate Bayesian computation (ABC), which has improved the reconstruction of invasion routes, but can be computationally expensive. We use our study to illustrate the use of a new, more efficient, ABC method, ABC random forest (ABC-RF) and compare it to a standard ABC method (ABC-LDA). We find that Japan emerges as the most probable source of the earliest recorded invasion into Hawaii. Southeast China and Hawaii together are the most probable sources of populations in western North America, which then in turn served as sources for those in eastern North America. European populations are genetically more homogeneous than North American populations, and their most probable source is northeast China, with evidence of limited gene flow from the eastern US as well. All introduced populations passed through bottlenecks, and analyses reveal five distinct admixture events. These findings can inform hypotheses concerning how this species evolved between different and independent source and invasive populations. Methodological comparisons indicate that ABC-RF and ABC-LDA show concordant results if ABC-LDA is based on a large number of simulated datasets but that ABC-RF out-performs ABC-LDA when using a comparable and more manageable number of simulated datasets, especially when analyzing complex introduction scenarios.