Present-day genetic introgression from domestic pigs into European wild boar has been suggested in various studies. However, no hybrids have been identified beyond doubt mainly because available methods were unable to quantify the extent of introgression and rule out natural processes. Genetic introgression from domestic pigs may have far-reaching ecological consequences by altering traits like the reproduction rate or immunology of wild boar. In this study, we demonstrate a novel approach to investigate genetic introgression in a Northwest (NW) European wild boar data set using a genome-wide single nucleotide polymorphism (SNP) assay developed for domestic pigs. We quantified the extent of introgression using allele frequency spectrum analysis, in silico hybridization simulations and genome distribution patterns of introgressed SNPs. Levels of recent introgression in the study area were expected to be low, as pig farming practices are prevailingly intensive and indoors. However, evidence was found for geographically widespread presence of domestic pig SNPs in 10% of analysed wild boar. This was supported by the identification of two different pig mitochondrial DNA haplotypes in three of the identified hybrid wild boar, suggesting that introgression had occurred from multiple sources (pig breeds). In silico hybridization simulations showed that the level of introgression in the identified hybrid wild boar is equivalent to first-generation hybrids until fifth-generation backcrosses with wild boar. The distribution pattern of introgressed SNPs supported these assignments in four of nine hybrids. The other five hybrids are considered advanced-generation hybrids, resulting from interbreeding among hybrid individuals. Three of nine hybrids were genetically associated with a different wild boar population than the one in which they were sampled. This discrepancy suggests that genetic introgression has occurred through the escape or release of an already hybridized farmed wild boar stock. We conclude that genetic introgression from domestic pigs into NW European wild boar populations is more recent and more common than expected and that genome-wide SNP analysis is a promising tool to quantify recent hybridization in free-living populations.

Complex microbial communities inhabit vertebrate digestive systems but thorough understanding of the ecological dynamics and functions of host-associated microbiota within natural habitats is limited. We investigate the role of environmental conditions in shaping gut and skin microbiota under natural conditions by performing a field survey and reciprocal transfer experiments with salamander larvae inhabiting two distinct habitats (ponds and streams). We show that gut and skin microbiota are habitat-specific, demonstrating environmental factors mediate community structure. Reciprocal transfer reveals that gut microbiota, but not skin microbiota, responds differentially to environmental change. Stream-to-pond larvae shift their gut microbiota to that of pond-to-pond larvae, whereas pond-to-stream larvae change to a community structure distinct from both habitat controls. Predicted functions, however, match that of larvae from the destination habitats in both cases. Thus, microbial function can be matched without taxonomic coherence and gut microbiota appears to exhibit metagenomic plasticity.

BackgroundPopulation genetic studies focus on natural dispersal and isolation by landscape barriers as the main drivers of genetic population structure. However, anthropogenic factors such as reintroductions, translocations and wild x domestic hybridization may also have strong effects on genetic population structure. In this study we genotyped 351 Single Nucleotide Polymorphism markers evenly spread across the genome in 645 wild boar (Sus scrofa) from Northwest Europe to evaluate determinants of genetic population structure.ResultsWe show that wild boar genetic population structure is influenced by historical reintroductions and by genetic introgression from domestic pigs. Six genetically distinct and geographically coherent wild boar clusters were identified in the Netherlands and Western Germany. The Dutch Veluwe cluster is known to be reintroduced, and three adjacent Dutch and German clusters are suspected to be a result of reintroduction, based on clustering results, low levels of heterozygosity and relatively high genetic distances to nearby populations. Recent wild x domestic hybrids were found geographically widespread across clusters and at low frequencies (average 3.9%). The relationship between pairwise kinship coefficients and geographic distance showed male-biased dispersal at the population genetic level.ConclusionsOur results demonstrate that wildlife and landscape management by humans are shaping the genetic diversity of an iconic wildlife species. Historical reintroductions, translocation and recent restocking activities with farmed wild boar have all influenced wild boar genetic population structure. The current trend of wild boar population growth and range expansion has recently led to a number of contact zones between clusters, and further admixture between the different wild boar clusters is to be expected.

The evolution of island populations in natural systems is driven by local adaptation and genetic drift. However, evolutionary pathways may be altered by humans in several ways. The wild boar (WB) (Sus scrofa) is an iconic game species occurring in several islands, where it has been strongly managed since prehistoric times. We examined genomic diversity at 49 803 single-nucleotide polymorphisms in 99 Sardinian WBs and compared them with 196 wild specimens from mainland Europe and 105 domestic pigs (DP; 11 breeds). High levels of genetic variation were observed in Sardinia (80.9% of the total number of polymorphisms), which can be only in part associated to recent genetic introgression. Both Principal Component Analysis and Bayesian clustering approach revealed that the Sardinian WB population is highly differentiated from the other European populations (FST=0.126–0.138), and from DP (FST=0.169). Such evidences were mostly unaffected by an uneven sample size, although clustering results in reference populations changed when the number of individuals was standardized. Runs of homozygosity (ROHs) pattern and distribution in Sardinian WB are consistent with a past expansion following a bottleneck (small ROHs) and recent population substructuring (highly homozygous individuals). The observed effect of a non-random selection of Sardinian individuals on diversity, FST and ROH estimates, stressed the importance of sampling design in the study of structured or introgressed populations. Our results support the heterogeneity and distinctiveness of the Sardinian population and prompt further investigations on its origins and conservation status.

The utilization of similar habitats by different species provides an ideal opportunity to identify genes underlying adaptation and acclimatization. Here, we analysed the gene expression of two closely related salamander species: Salamandra salamandra in Central Europe and Salamandra infraimmaculata in the Near East. These species inhabit similar habitat types: 'temporary ponds' and 'permanent streams' during larval development. We developed two species-specific gene expression microarrays, each targeting over 12 000 transcripts, including an overlapping subset of 8331 orthologues. Gene expression was examined for systematic differences between temporary ponds and permanent streams in larvae from both salamander species to establish gene sets and functions associated with these two habitat types. Only 20 orthologues were associated with a habitat in both species, but these orthologues did not show parallel expression patterns across species more than expected by chance. Functional annotation of a set of 106 genes with the highest effect size for a habitat suggested four putative gene function categories associated with a habitat in both species: cell proliferation, neural development, oxygen responses and muscle capacity. Among these high effect size genes was a single orthologue (14-3-3 protein zeta/YWHAZ) that was downregulated in temporary ponds in both species. The emergence of four gene function categories combined with a lack of parallel expression of orthologues (except 14-3-3 protein zeta) suggests that parallel habitat adaptation or acclimatization by larvae from S. salamandra and S. infraimmaculata to temporary ponds and permanent streams is mainly realized by different genes with a converging functionality.