Pristionchus pacificus, recently established as a model organism in evolutionary biology, is a cosmopolitan nematode that has a necromenic association with scarab beetles. The diverse array of host beetle species and habitat types occupied by P. pacificus make it a good model for investigating local adaptation to novel environments. Presence of P. pacificus on La Réunion Island, a young volcanic island with a dynamic geological history and a wide variety of ecozones, facilitates such investigation in an island biogeographic setting. Microsatellite data from 20 markers and 223 strains and mitochondrial sequence data from 272 strains reveal rich genetic diversity among La Réunion P. pacificus isolates, shaped by differentially timed introductions from diverse sources and in association with different beetle species. Distinctions between volcanic zones and between arid western and wet eastern climatic zones have likely limited westward dispersal of recently colonized lineages and maintained a genetic distinction between eastern and western clades. The highly selfing lifestyle of P. pacificus contributes to the strong fine-scale population structure detected, with each beetle host harbouring strongly differentiated assemblages of strains. Periodic out-crossing generates admixture between genetically diverse lineages, creating a diverse array of allelic combinations likely to increase the evolutionary potential of the species and facilitate adaptation to local environments and beetle hosts.

Historical trajectories of tree species during the late Quaternary have been well reconstructed through genetic and palaeobotanical studies. However, many congeneric tree species are interfertile, and the timing and contribution of introgression to species divergence during their evolutionary history remains largely unknown. We quantified past and current gene flow events between four morphologically divergent oak species (Quercus petraea, Q. robur, Q. pyrenaica, Q. pubescens), by two independent inference methods: diffusion approximation to the joint frequency spectrum (∂a∂i) and approximate Bayesian computation (ABC). For each pair of species, alternative scenarios of speciation allowing gene flow over different timescales were evaluated. Analyses of 3524 single nucleotide polymorphisms (SNPs) randomly distributed in the genome, showed that these species evolved in complete isolation for most of their history, but recently came into secondary contact, probably facilitated by the most recent period of postglacial warming. We demonstrated that: there was sufficient genetic differentiation before secondary contact for the accumulation of barriers to gene flow; and current European white oak genomes are a mosaic of genes that have crossed species boundaries and genes impermeable to gene flow.

An Illumina Infinium SNP genotyping array was constructed for European white oaks. Six individuals of Quercus petraea and Q. robur were considered for SNP discovery using both previously obtained Sanger sequences across 676 gene regions (1371 in vitro SNPs) and Roche 454 technology sequences from 5112 contigs (6542 putative in silico SNPs). The 7913 SNPs were genotyped across the six parental individuals, full-sib progenies (one within each species and two interspecific crosses between Q. petraea and Q. robur) and three natural populations from south-western France that included two additional interfertile white oak species (Q. pubescens and Q. pyrenaica). The genotyping success rate in mapping populations was 80.4% overall and 72.4% for polymorphic SNPs. In natural populations, these figures were lower (54.8% and 51.9%, respectively). Illumina genotype clusters with compression (shift of clusters on the normalized x-axis) were detected in ~25% of the successfully genotyped SNPs and may be due to the presence of paralogues. Compressed clusters were significantly more frequent for SNPs showing a priori incorrect Illumina genotypes, suggesting that they should be considered with caution or discarded. Altogether, these results show a high experimental error rate for the Infinium array (between 15% and 20% of SNPs potentially unreliable and 10% when excluding all compressed clusters), and recommendations are proposed when applying this type of high-throughput technique. Finally, results on diversity levels and shared polymorphisms across targeted white oaks and more distant species of the Quercus genus are discussed, and perspectives for future comparative studies are proposed.