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.

Photographic capture–recapture is a valuable tool for obtaining demographic information on wildlife populations due to its noninvasive nature and cost‐effectiveness. Recently, several computer‐aided photo‐matching algorithms have been developed to more efficiently match images of unique individuals in databases with thousands of images. However, the identification accuracy of these algorithms can severely bias estimates of vital rates and population size. Therefore, it is important to understand the performance and limitations of state‐of‐the‐art photo‐matching algorithms prior to implementation in capture–recapture studies involving possibly thousands of images. Here, we compared the performance of four photo‐matching algorithms; Wild‐ID, I3S Pattern+, APHIS, and AmphIdent using multiple amphibian databases of varying image quality. We measured the performance of each algorithm and evaluated the performance in relation to database size and the number of matching images in the database. We found that algorithm performance differed greatly by algorithm and image database, with recognition rates ranging from 100% to 22.6% when limiting the review to the 10 highest ranking images. We found that recognition rate degraded marginally with increased database size and could be improved considerably with a higher number of matching images in the database. In our study, the pixel‐based algorithm of AmphIdent exhibited superior recognition rates compared to the other approaches. We recommend carefully evaluating algorithm performance prior to using it to match a complete database. By choosing a suitable matching algorithm, databases of sizes that are unfeasible to match “by eye” can be easily translated to accurate individual capture histories necessary for robust demographic estimates.

The postembryonic development of amphibians has been characterized as divided into three predominant periods, hereafter named primary developmental stages: premetamorphosis (PreM), prometamorphosis (ProM), metamorphic climax (Meta), and completion of metamorphosis (PostM), largely based on examination of anuran development. Here, we categorized the postembryonic development of larvae of a poisonous fire salamander (Salamandra salamandra) by integrating morphology and gene expression (transcriptomic) data. Morphological analysis revealed three distinct clusters suggestive of PreM, ProM, and Meta, which were confirmed in parallel by microarray-derived gene expression analysis. In total, 3,510 probes targeted transcripts differentially expressed between the clusters we identified. Genes upregulated in PreM related to organogenesis, and those upregulated in Meta underlie structural proteins and related to development of anatomical structures and pigmentation. Biosynthesis pathways of pigments (pteridines and melanin) were upregulated during late ProM and Meta. Gas chromatographic analysis of alkaloids indicated the onset of steroidal alkaloid biosynthesis at ProM. When comparing gene expression in the fire salamander to that in other amphibians-three anurans, Xenopus laevis, X. tropicalis, and Michrohyla fissipes, and one caudate, Ambystoma mexicanum- we identified genes with conserved expression patterns involved in basic metamorphic processes such as skin restructuring and tail fin resorption. Our results support that primary stages of postembryonic development in caudates are homologous to those of anurans, and offer a baseline for the study of the evolution of developmental modes.

Aim We compared the power of different nuclear markers to investigate genetic structure of southern Balkan wild boar. We distinguished between historic events, such as isolation in different refugia during glacial periods, from recent demographic processes, such as naturally occurring expansions.Location Southern Balkans/Greece. MethodsWe sampled 555 wild boars from 20 different locations in southern Balkans/Greece. All individuals were analysed with 10 microsatellites and a subgroup of 91 with 49,508 single nucleotide polymorphisms (SNPs). Patterns of genetic structure and demographic processes were assessed with Bayesian clustering, linkage disequilibrium and past effective population size estimation analysis.Results Both microsatellite and SNP data analyses detected genetic structure caused by historic events and support the existence of three groups in the studied area. A hybrid zone between two of the groups was also detected. We also showed that genome-wide SNP data analysis can identify recent events in bottlenecked populations. Main conclusionsWe inferred the three groups diverged~50,000-10,000yr bp when populations contracted to different refugia. Our findings strengthened the evidence that the southern Balkan area was a glacial refugium including further local smaller refugia. Genome-wide genotyping inferred a recent population expansion that can mimic a 'refugium within refugium' scenario. It seems that microsatellite data tend to overestimate genetic structure when genetic drift happens in bottlenecked populations over a short distance. Therefore, genome-wide SNPs are more powerful at inferring phylogeography in natural populations, resolving inconsistencies from mitochondrial and microsatellite data sets.