Renaud Vitalis

rehh: an R package to detect footprints of selection in genome-wide SNP data from haplotype structure

UNLABELLED With the development of next-generation sequencing and genotyping approaches, large single nucleotide polymorphism haplotype datasets are becoming available in a growing number of both model and non-model species. Identifying genomic regions with unexpectedly high local haplotype homozygosity relatively to neutral expectation represents a powerful strategy to ascertain candidate genes responding to natural or artificial selection. To facilitate genome-wide scans of selection based on the analysis of long-range haplotypes, we developed the R package rehh. It provides a versatile tool to detect the footprints of recent or ongoing selection with several graphical functions that help visual interpretation of the results. AVAILABILITY AND IMPLEMENTATION Stable version is available from CRAN: http://cran.r-project.org/. Development version is available from the R-forge repository: http://r-forge.r-project.org/projects/rehh. Both versions can be installed directly from R. Function documentation and example data files are provided within the package and a tutorial is available as Supplementary Material. rehh is distributed under the GNU General Public Licence (GPL ≥ 2).

In defence of model-based inference in phylogeography

Recent papers have promoted the view that model‐based methods in general, and those based on Approximate Bayesian Computation (ABC) in particular, are flawed in a number of ways, and are therefore inappropriate for the analysis of phylogeographic data. These papers further argue that Nested Clade Phylogeographic Analysis (NCPA) offers the best approach in statistical phylogeography. In order to remove the confusion and misconceptions introduced by these papers, we justify and explain the reasoning behind model‐based inference. We argue that ABC is a statistically valid approach, alongside other computational statistical techniques that have been successfully used to infer parameters and compare models in population genetics. We also examine the NCPA method and highlight numerous deficiencies, either when used with single or multiple loci. We further show that the ages of clades are carelessly used to infer ages of demographic events, that these ages are estimated under a simple model of panmixia and population stationarity but are then used under different and unspecified models to test hypotheses, a usage the invalidates these testing procedures. We conclude by encouraging researchers to study and use model‐based inference in population genetics.

Inferring Population Decline and Expansion From Microsatellite Data: A Simulation-Based Evaluation of the Msvar Method

Reconstructing the demographic history of populations is a central issue in evolutionary biology. Using likelihood-based methods coupled with Monte Carlo simulations, it is now possible to reconstruct past changes in population size from genetic data. Using simulated data sets under various demographic scenarios, we evaluate the statistical performance of Msvar, a full-likelihood Bayesian method that infers past demographic change from microsatellite data. Our simulation tests show that Msvar is very efficient at detecting population declines and expansions, provided the event is neither too weak nor too recent. We further show that Msvar outperforms two moment-based methods (the M-ratio test and Bottleneck) for detecting population size changes, whatever the time and the severity of the event. The same trend emerges from a compilation of empirical studies. The latest version of Msvar provides estimates of the current and the ancestral population size and the time since the population started changing in size. We show that, in the absence of prior knowledge, Msvar provides little information on the mutation rate, which results in biased estimates and/or wide credibility intervals for each of the demographic parameters. However, scaling the population size parameters with the mutation rate and scaling the time with current population size, as coalescent theory requires, significantly improves the quality of the estimates for contraction but not for expansion scenarios. Finally, our results suggest that Msvar is robust to moderate departures from a strict stepwise mutation model.

Inbreeding Depression Is Purged in the Invasive Insect Harmonia axyridis

Bottlenecks in population size reduce genetic diversity and increase inbreeding, which can lead to inbreeding depression. It is thus puzzling how introduced species, which typically pass through bottlenecks, become such successful invaders. However, under certain theoretical conditions, bottlenecks of intermediate size can actually purge the alleles that cause inbreeding depression. Although this process has been confirmed in model laboratory systems, it has yet to be observed in natural invasive populations. We evaluate whether such purging could facilitate biological invasions by using the world-wide invasion of the ladybird (or ladybug) Harmonia axyridis. We first show that invasive populations endured a bottleneck of intermediate intensity. We then demonstrate that replicate introduced populations experience almost none of the inbreeding depression suffered by native populations. Thus, rather than posing a barrier to invasion as often assumed, bottlenecks, by purging deleterious alleles, can enable the evolution of invaders that maintain high fitness even when inbred.

Relationships Between Transposable Elements Based Upon the Integrase-Transposase Domains: Is There a Common Ancestor?

The integrase domain of RNA-mediated elements (class I) and the transposase domain of DNA-mediated transposable elements (class II) were compared. A number of elements contain the DDE signature, which plays an important role in their integration. The possible relationships betweenmariner-Tc1 andIS elements, retrotransposons, and retroviruses were analyzed from an alignment of this region. Themariner-Tc1 superfamily, and LTR retrotransposons and retroviruses were found to be monophyletic groups. However, theIS elements of bacteria were found in several groups. These results were used to propose an evolutionary history that suggests a common ancestor for some integrases and transposases.

LARGE-SCALE CANDIDATE GENE SCAN REVEALS THE ROLE OF CHEMORECEPTOR GENES IN HOST PLANT SPECIALIZATION AND SPECIATION IN THE PEA APHID

Understanding the drivers of speciation is critical to interpreting patterns of biodiversity. The identification of the genetic changes underlying adaptation and reproductive isolation is necessary to link barriers to gene flow to the causal origins of divergence. Here, we present a novel approach to the genetics of speciation, which should complement the commonly used approaches of quantitative trait locus mapping and genome‐wide scans for selection. We present a large‐scale candidate gene approach by means of sequence capture, applied to identifying the genetic changes underlying reproductive isolation in the pea aphid, a model system for the study of ecological speciation. Targeted resequencing enabled us to scale up the candidate gene approach, specifically testing for the role of chemosensory gene families in host plant specialization. Screening for the signature of divergence under selection at 172 candidate and noncandidate loci, we revealed a handful of loci that show high levels of differentiation among host races, which almost all correspond to odorant and gustatory receptor genes. This study offers the first indication that some chemoreceptor genes, often tightly linked together in the genome, could play a key role in local adaptation and reproductive isolation in the pea aphid and potentially other phytophagous insects. Our approach opens a new route toward the functional genomics of ecological speciation.

Sex-Specific Genetic Structure and Social Organization in Central Asia: Insights from a Multi-Locus Study

In the last two decades, mitochondrial DNA (mtDNA) and the non-recombining portion of the Y chromosome (NRY) have been extensively used in order to measure the maternally and paternally inherited genetic structure of human populations, and to infer sex-specific demography and history. Most studies converge towards the notion that among populations, women are genetically less structured than men. This has been mainly explained by a higher migration rate of women, due to patrilocality, a tendency for men to stay in their birthplace while women move to their husbands house. Yet, since population differentiation depends upon the product of the effective number of individuals within each deme and the migration rate among demes, differences in male and female effective numbers and sex-biased dispersal have confounding effects on the comparison of genetic structure as measured by uniparentally inherited markers. In this study, we develop a new multi-locus approach to analyze jointly autosomal and X-linked markers in order to aid the understanding of sex-specific contributions to population differentiation. We show that in patrilineal herder groups of Central Asia, in contrast to bilineal agriculturalists, the effective number of women is higher than that of men. We interpret this result, which could not be obtained by the analysis of mtDNA and NRY alone, as the consequence of the social organization of patrilineal populations, in which genetically related men (but not women) tend to cluster together. This study suggests that differences in sex-specific migration rates may not be the only cause of contrasting male and female differentiation in humans, and that differences in effective numbers do matter.

When genes go to sleep: the population genetic consequences of seed dormancy and monocarpic perenniality.

In many annual plant populations, seeds may be dormant for several seasons before they germinate. Here, we investigate the consequences of both conditional (dispersed seeds cannot enter a dormant stage) and unconditional seed dormancy on the amount and the distribution of neutral genetic diversity within and among populations. We present joint demographic and population genetics models for single and subdivided populations and derive the effective size and population differentiation at both local and metapopulation scales. We suggest that a Wahlund effect is unlikely to result from age structure alone. Furthermore, the differentiation among populations is decreased by the presence of seed banks. We also extend these models to describe monocarpic (semelparous) perennial life cycle, where the nonreproductive stages are vegetative rosettes instead of dormant seeds. The main difference between the models relies in the way the density‐dependent regulation is acting. The effective size of monocarpic perennial species may be less than the census number of individuals, and among‐population differentiation is always larger than in annual species. We discuss our results in the light of recent population genetics surveys of annual plants with seed banks.

Genomic pattern of adaptive divergence in Arabidopsis halleri, a model species for tolerance to heavy metal

Pollution by heavy metals is one of the strongest environmental constraints in human‐altered environments that only a handful of species can cope with. Identifying the genes conferring to those species the ability to grow in polluted areas is a first step towards a global understanding of the evolutionary processes involved and will eventually improve phytoremediation practices. We used a genome‐scan approach to detect loci under divergent selection among four populations of Arabidopsis halleri growing on either polluted or nonpolluted habitats. Based on a high density of amplified fragment length polymorphism (AFLP) markers (820 AFLP markers, i.e. ~1 marker per 0.3 Mb), evidence for selection was found for some markers in every sampled population. Four loci departed from neutrality in both metallicolous populations and thus constitute high‐quality candidates for general adaptation to pollution. Interestingly, some candidates differed between the two metallicolous populations, suggesting the possibility that different loci may be involved in adaptation in the different metallicolous populations.

Chromosomal Inversions, Natural Selection and Adaptation in the Malaria Vector Anopheles funestus

Chromosomal polymorphisms, such as inversions, are presumably involved in the rapid adaptation of populations to local environmental conditions. Reduced recombination between alternative arrangements in heterozygotes may protect sets of locally adapted genes, promoting ecological divergence and potentially leading to reproductive isolation and speciation. Through a comparative analysis of chromosomal inversions and microsatellite marker polymorphisms, we hereby present biological evidence that strengthens this view in the mosquito Anopheles funestus s.s, one of the most important and widespread malaria vectors in Africa. Specimens were collected across a wide range of geographical, ecological, and climatic conditions in Cameroon. We observed a sharp contrast between population structure measured at neutral microsatellite markers and at chromosomal inversions. Microsatellite data detected only a weak signal for population structuring among geographical zones (F(ST) < 0.013, P < 0.01). By contrast, strong differentiation among ecological zones was revealed by chromosomal inversions (F(ST) > 0.190, P < 0.01). Using standardized estimates of F(ST), we show that inversions behave at odds with neutral expectations strongly suggesting a role of environmental selection in shaping their distribution. We further demonstrate through canonical correspondence analysis that heterogeneity in eco-geographical variables measured at specimen sampling sites explained 89% of chromosomal variance in A. funestus. These results are in agreement with a role of chromosomal inversions in ecotypic adaptation in this species. We argue that this widespread mosquito represents an interesting model system for the study of chromosomal speciation mechanisms and should provide ample opportunity for comparative studies on the evolution of reproductive isolation and speciation in major human malaria vectors.