Maxime Galan

Landscape connectivity influences gene flow in a roe deer population inhabiting a fragmented landscape: an individual-based approach

Changes in agricultural practices and forest fragmentation can have a dramatic effect on landscape connectivity and the dispersal of animals, potentially reducing gene flow within populations. In this study, we assessed the influence of woodland connectivity on gene flow in a traditionally forest‐dwelling species — the European roe deer — in a fragmented landscape. From a sample of 648 roe deer spatially referenced within a study area of 55 × 40 km, interindividual genetic distances were calculated from genotypes at 12 polymorphic microsatellite loci. We calculated two geographical distances between each pair of individuals: the Euclidean distance (straight line) and the ‘least cost distance’ (the trajectory that maximizes the use of wooded corridors). We tested the correlation between genetic pairwise distances and the two types of geographical pairwise distance using Mantel tests. The correlation was better using the least cost distance, which takes into account the distribution of wooded patches, especially for females (the correlation was stronger but not significant for males). These results suggest that in a fragmented woodland area roe deer dispersal is strongly linked to wooded structures and hence that gene flow within the roe deer population is influenced by the connectivity of the landscape.

Genetic structure is influenced by landscape features: empirical evidence from a roe deer population

The delimitation of population units is of primary importance in population management and conservation biology. Moreover, when coupled with landscape data, the description of population genetic structure can provide valuable knowledge about the permeability of landscape features, which is often difficult to assess by direct methods (e.g. telemetry). In this study, we investigated the genetic structuring of a roe deer population which recently recolonized a fragmented landscape. We sampled 1148 individuals from a 40 × 55‐km area containing several putative barriers to deer movements, and hence to gene flow, namely a highway, rivers and several canals. In order to assess the effect of these landscape features on genetic structure, we implemented a spatial statistical model known as geneland which analyses genetic structure, explicitly taking into account the spatial nature of the problem. Two genetic units were inferred, exhibiting a very low level of differentiation (FST = 0.008). The location of their boundaries suggested that there are no absolute barriers in this study area, but that the combination of several landscape features with low permeability can lead to population differentiation. Our analysis hence suggests that the landscape has a significant influence on the structuring of the population under study. It also illustrates the use of geneland as a powerful method to infer population structure, even in situations of young populations exhibiting low genetic differentiation.

A 454 multiplex sequencing method for rapid and reliable genotyping of highly polymorphic genes in large-scale studies

BackgroundHigh-throughput sequencing technologies offer new perspectives for biomedical, agronomical and evolutionary research. Promising progresses now concern the application of these technologies to large-scale studies of genetic variation. Such studies require the genotyping of high numbers of samples. This is theoretically possible using 454 pyrosequencing, which generates billions of base pairs of sequence data. However several challenges arise: first in the attribution of each read produced to its original sample, and second, in bioinformatic analyses to distinguish true from artifactual sequence variation. This pilot study proposes a new application for the 454 GS FLX platform, allowing the individual genotyping of thousands of samples in one run. A probabilistic model has been developed to demonstrate the reliability of this method.ResultsDNA amplicons from 1,710 rodent samples were individually barcoded using a combination of tags located in forward and reverse primers. Amplicons consisted in 222 bp fragments corresponding to DRB exon 2, a highly polymorphic gene in mammals. A total of 221,789 reads were obtained, of which 153,349 were finally assigned to original samples. Rules based on a probabilistic model and a four-step procedure, were developed to validate sequences and provide a confidence level for each genotype. The method gave promising results, with the genotyping of DRB exon 2 sequences for 1,407 samples from 24 different rodent species and the sequencing of 392 variants in one half of a 454 run. Using replicates, we estimated that the reproducibility of genotyping reached 95%.ConclusionsThis new approach is a promising alternative to classical methods involving electrophoresis-based techniques for variant separation and cloning-sequencing for sequence determination. The 454 system is less costly and time consuming and may enhance the reliability of genotypes obtained when high numbers of samples are studied. It opens up new perspectives for the study of evolutionary and functional genetics of highly polymorphic genes like major histocompatibility complex genes in vertebrates or loci regulating self-compatibility in plants. Important applications in biomedical research will include the detection of individual variation in disease susceptibility. Similarly, agronomy will benefit from this approach, through the study of genes implicated in productivity or disease susceptibility traits.

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.

Migration and recovery of the genetic diversity during the increasing density phase in cyclic vole populations

In cyclic populations, high genetic diversity is currently reported despite the periodic low numbers experienced by the populations during the low phases. Here, we report spatio‐temporal monitoring at a very fine scale of cyclic populations of the fossorial water vole (Arvicola terrestris) during the increasing density phase. This phase marks the transition from a patchy structure (demes) during low density to a continuous population in high density. We found that the genetic diversity was effectively high but also that it displayed a local increase within demes over the increasing phase. The genetic diversity remained relatively constant when considering all demes together. The increase in vole abundance was also correlated with a decrease of genetic differentiation among demes. Such results suggest that at the end of the low phase, demes are affected by genetic drift as the result of being small and geographically isolated. This leads to a loss of local genetic diversity and a spatial differentiation among demes. This situation is counterbalanced during the increasing phase by the spatial expansion of demes and the increase of the effective migration among differentiated demes. We provide evidences that in cyclic populations of the fossorial water voles, the relative influence of drift operating during low density populations and migration occurring principally while population size increases interacts closely to maintain high genetic diversity.

Density-related changes in selection pattern for major histocompatibility complex genes in fluctuating populations of voles.

Host–pathogen interactions are of particular interest in studies of the interplay between population dynamics and natural selection. The major histocompatibility complex (MHC) genes of demographically fluctuating species are highly suitable markers for such studies, because they are involved in initiating the immune response against pathogens and display a high level of adaptive genetic variation. We investigated whether two MHC class II genes (DQA1, DRB) were subjected to contemporary selection during increases in the density of fossorial water vole (Arvicola terrestris) populations, by comparing the neutral genetic structure of seven populations with that estimated from MHC genes. Tests for heterozygosity excess indicated that DQA1 was subject to intense balancing selection. No such selection operated on neutral markers. This pattern of selection became more marked with increasing abundance. In the low‐abundance phase, when populations were geographically isolated, both overall differentiation and isolation‐by‐distance were more marked for MHC genes than for neutral markers. Model‐based simulations identified DQA1 as an outlier (i.e. under selection) in a single population, suggesting the action of local selection in fragmented populations. The differences between MHC and neutral markers gradually disappeared with increasing effective migration between sites. In the high‐abundance year, DQA1 displayed significantly lower levels of overall differentiation than the neutral markers. This gene therefore displayed stronger homogenization than observed under drift and migration alone. The observed signs of selection were much weaker for DRB. Spatial and temporal fluctuations in parasite pressure and locus‐specific selection are probably the most plausible mechanisms underlying the observed changes in selection pattern during the demographic cycle.

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.

ChickRH6: a chicken whole-genome radiation hybrid panel

As a first step towards the development of radiation hybrid maps, we have produced a radiation hybrid panel in the chicken by fusing female embryonic diploid fibroblasts irradiated at 6 000 rads with HPRT-deficient hamster Wg3hCl2 cells. Due to the low retention frequency of the chicken fragments, a high number of clones was produced from which the best ones were selected. Thus, 452 fusion clones were tested for retention frequencies with a panel of 46 markers. Based on these results, 103 clones with a mean marker retention of 23.8% were selected for large scale culture to produce DNA in sufficient quantities for the genotyping of numerous markers. Retention frequency was tested again with the same 46 markers and the 90 best clones, with a final mean retention frequency of 21.9%, were selected for the final panel. This panel will be a valuable resource for fine mapping of markers and genes in the chicken, and will also help in building BAC contigs.

Genetic structure of the cyclic fossorial water vole (Arvicola terrestris): landscape and demographic influences.

Genetic structure can be strongly affected by landscape features and variation through time and space of demographic parameters such as population size and migration rate. The fossorial water vole (Arvicola terrestris) is a cyclic species characterized by large demographic fluctuations over short periods of time. The outbreaks do not occur everywhere at the same time but spread as a wave at a regional scale. This leads to a pattern of large areas (i.e. some hundreds of km2), each with different vole abundances, at any given time. Here, we describe the abundance and genetic structures in populations of the fossorial water vole. We use the data to try to understand how landscape and demographic features act to shape the genetic structure. The spatial variability of vole abundance was assessed from surface indices, collected in spring 2002 (April) in eastern central France. Genetic variability was analysed using eight microsatellite loci at 23 localities sampled between October 2001 and April 2002. We found some congruence between abundance and genetic structures. At a regional scale, the genetic disruptions were associated with both sharp relief and transition between an area of low abundance and another of high abundance. At a local scale, we observed a variation of the isolation‐by‐distance pattern according to the abundance level of vole populations. From these results we suggest that the dispersal pattern in cyclic rodent populations varies throughout the demographic cycle.

Next-generation sequencing for rodent barcoding: species identification from fresh, degraded and environmental samples.

Rodentia is the most diverse order among mammals, with more than 2,000 species currently described. Most of the time, species assignation is so difficult based on morphological data solely that identifying rodents at the specific level corresponds to a real challenge. In this study, we compared the applicability of 100 bp mini-barcodes from cytochrome b and cytochrome c oxidase 1 genes to enable rodent species identification. Based on GenBank sequence datasets of 115 rodent species, a 136 bp fragment of cytochrome b was selected as the most discriminatory mini-barcode, and rodent universal primers surrounding this fragment were designed. The efficacy of this new molecular tool was assessed on 946 samples including rodent tissues, feces, museum samples and feces/pellets from predators known to ingest rodents. Utilizing next-generation sequencing technologies able to sequence mixes of DNA, 1,140 amplicons were tagged, multiplexed and sequenced together in one single 454 GS-FLX run. Our method was initially validated on a reference sample set including 265 clearly identified rodent tissues, corresponding to 103 different species. Following validation, 85.6% of 555 rodent samples from Europe, Asia and Africa whose species identity was unknown were able to be identified using the BLASTN program and GenBank reference sequences. In addition, our method proved effective even on degraded rodent DNA samples: 91.8% and 75.9% of samples from feces and museum specimens respectively were correctly identified. Finally, we succeeded in determining the diet of 66.7% of the investigated carnivores from their feces and 81.8% of owls from their pellets. Non-rodent species were also identified, suggesting that our method is sensitive enough to investigate complete predator diets. This study demonstrates how this molecular identification method combined with high-throughput sequencing can open new realms of possibilities in achieving fast, accurate and inexpensive species identification.