Hybridization between wild species and their domestic counterparts may represent a major threat to natural populations. However, high genetic similarity between the hybridizing taxa makes the detection of hybrids a difficult task and may hinder attempts to assess the impact of hybridization in conservation biology. In this work, we used a combination of 42 autosomal microsatellites together with Y-chromosome microsatellite-defined haplotypes and mtDNA sequences to investigate the occurrence and dynamics of wolf-dog hybridization in the Iberian Peninsula. To do this, we applied a variety of Bayesian analyses and a parallel set of simulation studies to evaluate (i) the differences between Iberian wolves and dogs, (ii) the frequency and geographical distribution of hybridization and (iii) the directionality of hybridization. First, we show that Iberian wolves and dogs form two well-differentiated genetic entities, suggesting that introgressive hybridization is not a widespread phenomenon shaping both gene pools. Second, we found evidence for the existence of hybridization that is apparently restricted to more peripheral and recently expanded wolf populations. Third, we describe compelling evidence suggesting that the dynamics of hybridization in wolf populations is mediated by crosses between male dogs and female wolves. More importantly, the observation of a population showing the occurrence of a continuum of hybrid classes forming mixed packs may indicate that we have underestimated hybridization. If future studies confirm this pattern, then an intriguing avenue of research is to investigate how introgression from free-ranging domestic dogs is enabling wolf populations to adapt to the highly humanized habitats of southern Europe while still maintaining their genetic differentiation.

Wolves and dogs provide a paradigmatic example of the ecological and conservation implications of hybridization events between wild and domesticated forms. However, our understanding of such implications has been traditionally hampered by both high genetic similarity and the difficulties in obtaining tissue samples (TS), which limit our ability to assess ongoing hybridization events. To assess the occurrence and extension of hybridization in a pack of wolf-dog hybrids in northwestern Iberia, we compared the power of 52 nuclear markers implemented on TS with a subset of 13 ancestry informative markers (AIMs) typed in noninvasive samples (NIS). We demonstrate that the 13 AIMs are as accurate as the 52 markers that were chosen without regard to the power to differentiate between wolves and dogs, also having the advantage of being rapidly screened on NIS. The efficiency of AIMs significantly outperformed ten random sets of similar size and an additional commercial set of 18 markers. Bayesian clustering analysis implemented on AIMs and NIS identified nine hybrids, two wolves and two dogs. Four hybrids were unambiguously assigned to F1xWolf backcrosses. Our approach (AIMs + NIS) overcomes previous difficulties related to sample availability and informative power of markers, allowing a quick identification of wolf-dog hybrids in the first phases of hybridization episodes. This provides managers with a reliable tool to evaluate hybridization and estimate the success of their actions. This approach may be easily adapted for other pairs of wild/domesticated species, thus improving our understanding of the introgression of domestication genes into natural populations.

Over thousands of years humans changed the genetic and phenotypic composition of several organisms and in the process transformed wild species into domesticated forms. From this close association, domestic animals emerged as important models in biomedical and fundamental research, in addition to their intrinsic economical and cultural value. The domestic rabbit is no exception but few studies have investigated the impact of domestication on its genetic variability. In order to study patterns of genetic structure in domestic rabbits and to quantify the genetic diversity lost with the domestication process, we genotyped 45 microsatellites for 471 individuals belonging to 16 breeds and 13 wild localities. We found that both the initial domestication and the subsequent process of breed formation, when averaged across breeds, culminated in losses of ~20% of genetic diversity present in the ancestral wild population and domestic rabbits as a whole, respectively. Despite the short time elapsed since breed diversification we uncovered a well-defined structure in domestic rabbits where the FST between breeds was 22%. However, we failed to detect deeper levels of structure, probably consequence of a recent and single geographic origin of domestication together with a non-bifurcating process of breed formation, which were often derived from crosses between two or more breeds. Finally, we found evidence for intrabreed stratification that is associated with demographic and selective causes such as formation of strains, colour morphs within the same breed, or country/breeder of origin. These additional layers of population structure within breeds should be taken into account in future mapping studies.

Although whole-genome sequencing is becoming more accessible and feasible for nonmodel organisms, microsatellites have remained the markers of choice for various population and conservation genetic studies. However, the criteria for choosing microsatellites are still controversial due to ascertainment bias that may be introduced into the genetic inference. An empirical study of red deer (Cervus elaphus) populations, in which cross-specific and species-specific microsatellites developed through pyrosequencing of enriched libraries, was performed for this study. Two different strategies were used to select the species-specific panels: randomly vs. highly polymorphic markers. The results suggest that reliable and accurate estimations of genetic diversity can be obtained using random microsatellites distributed throughout the genome. In addition, the results reinforce previous evidence that selecting the most polymorphic markers leads to an ascertainment bias in estimates of genetic diversity, when compared with randomly selected microsatellites. Analyses of population differentiation and clustering seem less influenced by the approach of microsatellite selection, whereas assigning individuals to populations might be affected by a random selection of a small number of microsatellites. Individual multilocus heterozygosity measures produced various discordant results, which in turn had impacts on the heterozygosity-fitness correlation test. Finally, we argue that picking the appropriate microsatellite set should primarily take into account the ecological and evolutionary questions studied. Selecting the most polymorphic markers will generally overestimate genetic diversity parameters, leading to misinterpretations of the real genetic diversity, which is particularly important in managed and threatened populations.

Highly mobile mammalian carnivores are expected to have the capability to maintain high levels of gene flow across large geographic scales. Nonetheless, surprising levels of genetic structure have been found in many such populations. We combined genetic and spatial behavioural information from wolves (Canis lupus) in the Iberian Peninsula (Western Europe) during the last two decades to present a particular case of low dispersal levels in a large carnivore population persisting in human-dominated landscapes. We found an exceptionally reticulated pattern of cryptic population structure emerging at two hierarchical levels, in which four or eleven meaningful genetic clusters can be recognized, respectively. These clusters were characterized by moderate-high levels of differentiation (average pairwise F = 0.09-0.19), low levels of admixture and varying degrees of genetic diversity. The number of dispersers identified among the 11 clusters was very low (<4% out of 218 wolves). Spatial information of tracked wolves further confirmed the geographical genetic patterns (only 2 out of 85 collared wolves overlapped with more than one genetic cluster). The high levels of genetic structure in this population may be determined by the recent demographic history of this population, among other factors. The identification of meaningful genetic clusters has implications for the delineation of conservation units and, consequently, on the conservation and management actions for Iberian wolves.

We report new primers for the amplification and sequencing of 11 nuclear markers in squamate reptiles and anuran amphibians (five in squamates, six in anurans). Ten out of the 11 loci are introns (three of which are linked) that were amplified using an exon-primed, introncrossing (EPIC) PCR strategy, whereas an eleventh locus spans part of a protein-coding gene. Squamate and anuran primers were initially developed for Lacerta schreiberi (Squamata: Lacertidae) and Pelodytes spp. (Anura: Pelodytidae), respectively. Cross-species amplification of the squamate markers was evaluated in four genera representing two additional families, whereas for anurans three genera corresponding to three additional families were tested. Three out of the five loci were successfully sequenced in all squamate taxa tested. Cross-amplification of the six anuran markers had lower, but still significant, success. We predict these markers will be of great utility for both population genetics and phylogenetic studies.

Faba bean (Vicia faba L.) is a facultative cross-pollinating legume crop with a great importance for food and feed due to its high protein content as well as the important role in soil fertility and nitrogen fixation. In this work we evaluated genetic diversity and population structure of faba bean accessions from the Western Mediterranean basin and wild related species. For that purpose we screened 53 V. faba, 2 V. johannis and 7 V. narbonensis accessions from Portugal, Spain and Morocco with 28 faba bean Single Sequence Repeats (SSR). SSR genotyping showed that the number of alleles detected per locus for the polymorphic markers ranged between 2 and 10, with Polymorphic Information Content (PIC) values between 0.662 and 0.071, and heterozygosity (HO) between 0–0.467. Heterozygosity and inbreeding coefficient levels indicate a higher level of inbreeding in wild related species than in cultivated Vicia. The analysis of molecular variance (AMOVA) showed a superior genetic diversity within accessions than between accessions even from distant regions. These results are in accordance to population structure analysis showing that individuals from the same accession can be genetically more similar to individuals from far away accessions, than from individuals from the same accession. In all three levels of analysis (whole panel of cultivated and wild accessions, cultivated faba bean accessions and Portuguese accessions) no population structure was observed based on geography or climatic factors. Differences between V. narbonensis and V. johannis are undetectable although these wild taxa are clearly distinct from V. faba accessions. Thus, a limited gene flow occurred between cultivated accessions and wild relatives. Contrastingly, the lack of population structure seems to indicate a high degree of gene flow between V. faba accessions, possibly explained by the partially allogamous habit in association with frequent seed exchange/introduction.

Large mammals are seriously threatened in North Africa, with emblematic cases of extinction reported during the twentieth century. The Dorcas gazelle (Gazella dorcas) is an endangered species whose populations drastically declined in the last few decades. In this work we applied both invasive and non-invasive molecular methods to document for the first time patterns of genetic diversity and population structure of G. dorcas in its northwestern range, using 13 microsatellite loci and a 716 bp fragment of mitochondrial DNA in seven populations (one wild, four semi-captive and two captive). Genetic diversity at both mtDNA and nuclear markers showed a clear dichotomy among G. dorcas populations. While the wild population and the captive population maintained in Almeria (Spain) exhibited appreciable levels of diversity, all five captive and semi-captive populations across Morocco were fixed for a single mtDNA haplotype and showed low values of nuclear diversity. The allele frequency spectrum analysis of these five populations revealed profiles expected under a bottleneck scenario, but statistical tests performed to investigate this situation were not significant. Genetic differentiation measured by summary statistics (F ST and D est ) and population structure revealed by Bayesian clustering analysis suggest that the Sidi Chiker Reserve, in the northern plains of Morocco, may harbour the last individuals belonging to G. d. massaesyla, while individuals from the El Kheng Reserve exhibited a moderate degree of differentiation and could not be unambiguously associated with one of the two remaining putative subspecies. These data should be taken into consideration in order to implement a conservation action plan for G. dorcas in Morocco.