Survival and divergence in a small group: The extraordinary genomic history of the endangered Apennine brown bear stragglers

Large scale genome reconstructions illuminate Wolbachia evolution

Wolbachia is an iconic example of a successful intracellular bacterium. Despite its importance as a manipulator of invertebrate biology, its evolutionary dynamics have been poorly studied from a genomic viewpoint. To expand the number of Wolbachia genomes, we screen over 30,000 publicly available shotgun DNA sequencing samples from 500 hosts. By assembling over 1000 Wolbachia genomes, we provide a substantial increase in host representation. Our phylogenies based on both core-genome and gene content provide a robust reference for future studies, support new strains in model organisms, and reveal recent horizontal transfers amongst distantly related hosts. We find various instances of gene function gains and losses in different super-groups and in cytoplasmic incompatibility inducing strains. Our Wolbachia-host co-phylogenies indicate that horizontal transmission is widespread at the host intraspecific level and that there is no support for a general Wolbachia-mitochondrial synchronous divergence.

By greatly expanding the number of assembled genomes for Wolbachia (a group of intracellular bacteria) and constructing robust phylogenies, this study finds strong rate heterogeneity among Wolbachiapopulations and no support for synchronous divergence between Wolbachia and host mitochondria.

Distinguishing among complex evolutionary models using unphased whole-genome data through random forest approximate Bayesian computation

Inferring past demographic histories is crucial in population genetics, and the amount of complete genomes now available should in principle facilitate this inference. In practice, however, the available inferential methods suffer from severe limitations. Although hundreds complete genomes can be simultaneously analysed, complex demographic processes can easily exceed computational constraints, and the procedures to evaluate the reliability of the estimates contribute to increase the computational effort. Here we present an approximate Bayesian computation framework based on the random forest algorithm (ABC-RF), to infer complex past population processes using complete genomes. To this aim, we propose to summarize the data by the full genomic distribution of the four mutually exclusive categories of segregating sites (FDSS), a statistic fast to compute from unphased genome data and that does not require the ancestral state of alleles to be known. We constructed an efficient ABC pipeline and tested how accurately it allows one to recognize the true model among models of increasing complexity, using simulated data and taking into account different sampling strategies in terms of number of individuals analysed, number and size of the genetic loci considered. We also compared the FDSS with the unfolded and folded site frequency spectrum (SFS), and for these statistics we highlighted the experimental conditions maximizing the inferential power of the ABC-RF procedure. We finally analysed real data sets, testing models on the dispersal of anatomically modern humans out of Africa and exploring the evolutionary relationships of the three species of Orangutan inhabiting Borneo and Sumatra.

Dental Calculus as a Tool to Study the Evolution of the Mammalian Oral Microbiome

Dental calculus, the calcified form of the mammalian oral microbial plaque biofilm, is a rich source of oral microbiome, host, and dietary biomolecules and is well preserved in museum and archaeological specimens. Despite its wide presence in mammals, to date, dental calculus has primarily been used to study primate microbiome evolution. We establish dental calculus as a valuable tool for the study of nonhuman host microbiome evolution, by using shotgun metagenomics to characterize the taxonomic and functional composition of the oral microbiome in species as diverse as gorillas, bears, and reindeer. We detect oral pathogens in individuals with evidence of oral disease, assemble near-complete bacterial genomes from historical specimens, characterize antibiotic resistance genes, reconstruct components of the host diet, and recover host genetic profiles. Our work demonstrates that metagenomic analyses of dental calculus can be performed on a diverse range of mammalian species, which will allow the study of oral microbiome and pathogen evolution from a comparative perspective. As dental calculus is readily preserved through time, it can also facilitate the quantification of the impact of anthropogenic changes on wildlife and the environment.

Adaptive and maladaptive genetic diversity in small populations: Insights from the Brook Charr (Salvelinus fontinalis) case study

Investigating the relative importance of neutral versus selective processes governing the accumulation of genetic variants is a key goal in both evolutionary and conservation biology. This is particularly true in the context of small populations, where genetic drift can counteract the effect of selection. Using Brook Charr (Salvelinus fontinalis) from Québec, Canada, as a case study, we investigated the importance of demographic versus selective processes governing the accumulation of both adaptive and maladaptive mutations in closed versus open and connected populations to assess gene flow effect. This was achieved by using 14,779 high-quality filtered SNPs genotyped among 1,416 fish representing 50 populations from three life history types: lacustrine (closed populations), riverine and anadromous (connected populations). Using the PROVEAN algorithm, we observed a considerable accumulation of putative deleterious mutations across populations. The absence of correlation between the occurrence of putatively beneficial or deleterious mutations and local recombination rate supports the hypothesis that genetic drift might be the main driver of the accumulation of such variants. However, despite a lower genetic diversity observed in lacustrine than in riverine or anadromous populations, lacustrine populations do not exhibit more deleterious mutations than the two other history types, suggesting that the negative effect of genetic drift in lacustrine populations may be mitigated by that of relaxed purifying selection. Moreover, we also identified genomic regions associated with anadromy, as well as an overrepresentation of transposable elements associated with variation in environmental variables, thus supporting the importance of transposable elements in adaptation.

Analyses of key genes involved in Arctic adaptation in polar bears suggest selection on both standing variation and de novo mutations played an important role

BACKGROUND

Polar bears are uniquely adapted to an Arctic existence. Since their relatively recent divergence from their closest living relative, brown bears, less than 500,000 years ago, the species has evolved an array of novel traits suited to its Arctic lifestyle. Previous studies sought to uncover the genomic underpinnings of these unique characteristics, and disclosed the genes showing the strongest signal of positive selection in the polar bear lineage. Here, we survey a comprehensive dataset of 109 polar bear and 33 brown bear genomes to investigate the genomic variants within these top genes present in each species. Specifically, we investigate whether fixed homozygous variants in polar bears derived from selection on standing variation in the ancestral gene pool or on de novo mutation in the polar bear lineage.

RESULTS

We find that a large number of sites fixed in polar bears are biallelic in brown bears, suggesting selection on standing variation. Moreover, we uncover sites in which polar bears are fixed for a derived allele while brown bears are fixed for the ancestral allele, which we suggest may be a signal of de novo mutation in the polar bear lineage.

CONCLUSIONS

Our findings suggest that, among other mechanisms, natural selection acting on changes in genes derived from a combination of variation already in the ancestral gene pool, and from de novo missense mutations in the polar bear lineage, may have enabled the rapid adaptation of polar bears to their new Arctic environment.

Drift, selection and adaptive variation in small populations of a threatened rattlesnake

An important goal of conservation genetics is to determine if the viability of small populations is reduced by a loss of adaptive variation due to genetic drift. Here, we assessed the impact of drift and selection on direct measures of adaptive variation (toxin loci encoding venom proteins) in the eastern massasauga rattlesnake (Sistrurus catenatus), a threatened reptile that exists in small isolated populations. We estimated levels of individual polymorphism in 46 toxin loci and 1,467 control loci across 12 populations of this species, and compared the results with patterns of selection on the same loci following speciation of S. catenatus and its closest relative, the western massasauga (S. tergeminus). Multiple lines of evidence suggest that both drift and selection have had observable impacts on standing adaptive variation. In support of drift effects, we found little evidence for selection on toxin variation within populations and a significant positive relationship between current levels of adaptive variation and long- and short-term estimates of effective population size. However, we also observed levels of directional selection on toxin loci among populations that are broadly similar to patterns predicted from interspecific selection analyses that pre-date the effects of recent drift, and that functional variation in these loci persists despite small short-term effective sizes. This suggests that much of the adaptive variation present in populations may represent an example of "drift debt," a nonequilibrium state where present-day levels of variation overestimate the amount of functional genetic diversity present in future populations.

Survival in the Tropics despite isolation, inbreeding and asexual reproduction: insights from the genome of the world's southernmost poplar (Populus ilicifolia)

Species are becoming extinct at unprecedented rates as a consequence of human activity. Hence it is important to understand the evolutionary dynamics of species with already small population sizes. Populus ilicifolia is a vulnerable poplar species that is isolated from other poplar species and is uniquely adapted to the Tropics. It has a very limited size, reproduces partly clonally and is therefore an excellent case study for conservation genomics. We present here the first annotated draft genome of P. ilicifolia, characterize genome-wide patterns of polymorphisms and compare those to other poplar species with larger natural ranges. P. ilicifolia experienced a more prolonged and severe decline of effective population size (N_e ) and signs of genetic erosion than any other poplar species with which it was compared. At present, the species has the lowest genome-wide genetic diversity, the highest abundance of long runs of homozygosity, high inbreeding levels as well as a high overall accumulation of deleterious variants. However, more effective purging of severely deleterious variants and adaptation to the Tropics may have contributed to its survival. Hence, in spite of its limited genetic variation, it is certainly worth pursuing the conservation efforts of this unique species.

Accumulation of Deleterious Mutations in Landlocked Threespine Stickleback Populations

Colonization of new habitats often reduces population sizes and may result in the accumulation of deleterious mutations by genetic drift. Compared with the genomic basis for adaptation to new environments, genome-wide analysis of deleterious mutations in isolated populations remains limited. In the present study, we investigated the accumulation of deleterious mutations in five endangered freshwater populations of threespine stickleback (Gasterosteus aculeatus) in the central part of the mainland of Japan. Using whole-genome resequencing data, we first conducted phylogenomic analysis and confirmed at least two independent freshwater colonization events in the central mainland from ancestral marine ecotypes. Next, analyses of single nucleotide polymorphisms showed a substantial reduction of heterozygosity in freshwater populations compared with marine populations. Reduction in heterozygosity was more apparent at the center of each chromosome than the peripheries and on X chromosomes compared with autosomes. Third, bioinformatic analysis of deleterious mutations showed increased accumulation of putatively deleterious mutations in the landlocked freshwater populations compared with marine populations. For the majority of populations examined, the frequencies of putatively deleterious mutations were higher on X chromosomes than on autosomes. The interpopulation comparison indicated that the majority of putatively deleterious mutations may have accumulated independently. Thus, whole-genome resequencing of endangered populations can help to estimate the accumulation of deleterious mutations and inform us of which populations are the most severely endangered. Furthermore, analysis of variation among chromosomes can give insights into whether any particular chromosomes are likely to accumulate deleterious mutations.

Genomic patterns in the widespread Eurasian lynx shaped by Late Quaternary climatic fluctuations and anthropogenic impacts

Disentangling the contribution of long-term evolutionary processes and recent anthropogenic impacts to current genetic patterns of wildlife species is key to assessing genetic risks and designing conservation strategies. Here, we used 80 whole nuclear genomes and 96 mitogenomes from populations of the Eurasian lynx covering a range of conservation statuses, climatic zones and subspecies across Eurasia to infer the demographic history, reconstruct genetic patterns, and discuss the influence of long-term isolation and/or more recent human-driven changes. Our results show that Eurasian lynx populations shared a common history until 100,000 years ago, when Asian and European populations started to diverge and both entered a period of continuous and widespread decline, with western populations, except Kirov, maintaining lower effective sizes than eastern populations. Population declines and increased isolation in more recent times probably drove the genetic differentiation between geographically and ecologically close westernmost European populations. By contrast, and despite the wide range of habitats covered, populations are quite homogeneous genetically across the Asian range, showing a pattern of isolation by distance and providing little genetic support for the several proposed subspecies. Mitogenomic and nuclear divergences and population declines starting during the Late Pleistocene can be mostly attributed to climatic fluctuations and early human influence, but the widespread and sustained decline since the Holocene is more probably the consequence of anthropogenic impacts which intensified in recent centuries, especially in western Europe. Genetic erosion in isolated European populations and lack of evidence for long-term isolation argue for the restoration of lost population connectivity.

Paths for colonization or exodus? New insights from the brown bear (Ursus arctos) population of the Cantabrian Mountains

Over the centuries, the geographical distribution of brown bear (Ursus arctos) across the Iberian Peninsula has been decreasing, with the species currently confined to North Iberia. The Cantabrian brown bear population is one of the smallest in Europe and is structured into two subpopulations, positioned along an east-west axis. Given the current critically endangered status of this population, it is essential to have a clear picture of its within-population genetic patterns and processes. We use a set of three molecular markers (mitochondrial DNA, autosomal microsatellites and sex markers) to clarify the genetic origins and assess the migration patterns and gene flow of the Cantabrian brown bear population. Our results reveal the presence of two different mitochondrial (matrilineal) haplotypes in the Cantabrian population, both belonging to European brown bear clade 1a. The two haplotypes are geographically structured between Eastern (haplotype CanE) and Western Cantabrian (haplotype CanW) subpopulations, which is consistent with the genetic structure previously identified using nuclear markers. Additionally, we show that CanE is closer to the historical Pyrenean (Pyr) haplotype than to CanW. Despite strong structuring at the levels of mtDNA and nuclear loci, there is evidence of bidirectional gene flow and admixture among subpopulations. Gene flow is asymmetrical and significantly more intense from the Eastern to the Western Cantabrian subpopulation. In fact, we only detected first generation male migrants from the Eastern to the Western Cantabrian subpopulation. These results suggest more intense migration from the smaller and more vulnerable Eastern Cantabrian subpopulation towards the larger and more stable Western Cantabrian subpopulation. These new insights are relevant for assessments of on-going conservation measures, namely the role of dispersal corridors and enhanced connectivity. We discuss the importance of complementary conservation measures, such as human-wildlife conflict mitigation and habitat improvement, for the conservation of a viable Cantabrian brown bear population.

Consensify: A Method for Generating Pseudohaploid Genome Sequences from Palaeogenomic Datasets with Reduced Error Rates

A standard practise in palaeogenome analysis is the conversion of mapped short read data into pseudohaploid sequences, frequently by selecting a single high-quality nucleotide at random from the stack of mapped reads. This controls for biases due to differential sequencing coverage, but it does not control for differential rates and types of sequencing error, which are frequently large and variable in datasets obtained from ancient samples. These errors have the potential to distort phylogenetic and population clustering analyses, and to mislead tests of admixture using D statistics. We introduce Consensify, a method for generating pseudohaploid sequences, which controls for biases resulting from differential sequencing coverage while greatly reducing error rates. The error correction is derived directly from the data itself, without the requirement for additional genomic resources or simplifying assumptions such as contemporaneous sampling. For phylogenetic and population clustering analysis, we find that Consensify is less affected by artefacts than methods based on single read sampling. For D statistics, Consensify is more resistant to false positives and appears to be less affected by biases resulting from different laboratory protocols than other frequently used methods. Although Consensify is developed with palaeogenomic data in mind, it is applicable for any low to medium coverage short read datasets. We predict that Consensify will be a useful tool for future studies of palaeogenomes.

Recent evolutionary history predicts population but not ecosystem-level patterns

In the face of rapid anthropogenic environmental change, it is increasingly important to understand how ecological and evolutionary interactions affect the persistence of natural populations. Augmented gene flow has emerged as a potentially effective management strategy to counteract negative consequences of genetic drift and inbreeding depression in small and isolated populations. However, questions remain about the long-term impacts of augmented gene flow and whether changes in individual and population fitness are reflected in ecosystem structure, potentiating eco-evolutionary feedbacks. In this study, we used Trinidadian guppies (Poecilia reticulata) in experimental outdoor mesocosms to assess how populations with different recent evolutionary histories responded to a scenario of severe population size reduction followed by expansion in a high-quality environment. We also investigated how variation in evolutionary history of the focal species affected ecosystem dynamics. We found that evolutionary history (i.e., gene flow vs. no gene flow) consistently predicted variation in individual growth. In addition, gene flow led to faster population growth in populations from one of the two drainages, but did not have measurable impacts on the ecosystem variables we measured: zooplankton density, algal growth, and decomposition rates. Our results suggest that benefits of gene flow may be long-term and environment-dependent. Although small in replication and duration, our study highlights the importance of eco-evolutionary interactions in determining population persistence and sets the stage for future work in this area.

Small population size and low genomic diversity have no effect on fitness in experimental translocations of a wild fish

Little empirical work in nature has quantified how wild populations with varying effective population sizes and genetic diversity perform when exposed to a gradient of ecologically important environmental conditions. To achieve this, juvenile brook trout from 12 isolated populations or closed metapopulations that differ substantially in population size and genetic diversity were transplanted to previously fishless ponds spanning a wide gradient of ecologically important variables. We evaluated the effect of genome-wide variation, effective population size (N_e), pond habitat, and initial body size on two fitness correlates (survival and growth). Genetic variables had no effect on either fitness correlate, which was determined primarily by habitat (pond temperature, depth, and pH) and initial body size. These results suggest that some vertebrate populations with low genomic diversity, low N_e, and long-term isolation can represent important sources of variation and are capable of maintaining fitness in, and ultimately persisting and adapting to, changing environments. Our results also reinforce the paramount importance of improving available habitat and slowing habitat degradation for species conservation.

Conservation through the lens of (mal)adaptation: Concepts and meta-analysis

Evolutionary approaches are gaining popularity in conservation science, with diverse strategies applied in efforts to support adaptive population outcomes. Yet conservation strategies differ in the type of adaptive outcomes they promote as conservation goals. For instance, strategies based on genetic or demographic rescue implicitly target adaptive population states whereas strategies utilizing transgenerational plasticity or evolutionary rescue implicitly target adaptive processes. These two goals are somewhat polar: adaptive state strategies optimize current population fitness, which should reduce phenotypic and/or genetic variance, reducing adaptability in changing or uncertain environments; adaptive process strategies increase genetic variance, causing maladaptation in the short term, but increase adaptability over the long term. Maladaptation refers to suboptimal population fitness, adaptation refers to optimal population fitness, and (mal)adaptation refers to the continuum of fitness variation from maladaptation to adaptation. Here, we present a conceptual classification for conservation that implicitly considers (mal)adaptation in the short-term and long-term outcomes of conservation strategies. We describe cases of how (mal)adaptation is implicated in traditional conservation strategies, as well as strategies that have potential as a conservation tool but are relatively underutilized. We use a meta-analysis of a small number of available studies to evaluate whether the different conservation strategies employed are better suited toward increasing population fitness across multiple generations. We found weakly increasing adaptation over time for transgenerational plasticity, genetic rescue, and evolutionary rescue. Demographic rescue was generally maladaptive, both immediately after conservation intervention and after several generations. Interspecific hybridization was adaptive only in the F1 generation, but then rapidly leads to maladaptation. Management decisions that are made to support the process of adaptation must adequately account for (mal)adaptation as a potential outcome and even as a tool to bolster adaptive capacity to changing conditions.

Population correlates of rapid captive-induced maladaptation in a wild fish

Understanding the extent to which captivity generates maladaptation in wild species can inform species recovery programs and elucidate wild population responses to novel environmental change. Although rarely quantified, effective population size (N e ) and genetic diversity should influence the magnitude of plastic and genetic changes manifested in captivity that reduce wild fitness. Sexually dimorphic traits might also mediate consequences of captivity. To evaluate these relationships, we generated >600 full- and half-sibling families from nine wild brook trout populations, reared them for one generation under common, captive environmental conditions and contrasted several fitness-related traits in wild versus captive lines. We found substantial variation in lifetime success (lifetime survival and reproductive success) and life history traits among wild populations after just one captive generation (fourteen- and threefold ranges across populations, respectively). Populations with lower heterozygosity showed lower captive lifetime success, suggesting that captivity generates maladaptation within one generation. Greater male-biased mortality in captivity occurred in populations having disproportionately higher growth rates in males than females. Wild population N e and allelic diversity had little or no influence on captive trait expression and lifetime success. Our results have four conservation implications: (i) Trait values and lifetime success were highly variable across populations following one generation of captivity. (ii) Maladaptation induced by captive breeding might be particularly intense for the very populations practitioners are most interested in conserving, such as those with low heterozygosity. (iii) Maladaptive sex differences in captivity might be associated with population-dependent growth costs of reproduction. (iv) Heterozygosity can be a good indicator of short-term, intraspecific responses to novel environmental change.

Genetic diversity of CHC22 clathrin impacts its function in glucose metabolism

CHC22 clathrin plays a key role in intracellular membrane traffic of the insulin-responsive glucose transporter GLUT4 in humans. We performed population genetic and phylogenetic analyses of the CHC22-encoding CLTCL1 gene, revealing independent gene loss in at least two vertebrate lineages, after arising from gene duplication. All vertebrates retained the paralogous CLTC gene encoding CHC17 clathrin, which mediates endocytosis. For vertebrates retaining CLTCL1, strong evidence for purifying selection supports CHC22 functionality. All human populations maintained two high frequency CLTCL1 allelic variants, encoding either methionine or valine at position 1316. Functional studies indicated that CHC22-V1316, which is more frequent in farming populations than in hunter-gatherers, has different cellular dynamics than M1316-CHC22 and is less effective at controlling GLUT4 membrane traffic, altering its insulin-regulated response. These analyses suggest that ancestral human dietary change influenced selection of allotypes that affect CHC22's role in metabolism and have potential to differentially influence the human insulin response.

Genomic signatures of extensive inbreeding in Isle Royale wolves, a population on the threshold of extinction

The observation that small isolated populations often suffer reduced fitness from inbreeding depression has guided conservation theory and practice for decades. However, investigating the genome-wide dynamics associated with inbreeding depression in natural populations is only now feasible with relatively inexpensive sequencing technology and annotated reference genomes. To characterize the genome-wide effects of intense inbreeding and isolation, we performed whole-genome sequencing and morphological analysis of an iconic inbred population, the gray wolves (Canis lupus) of Isle Royale. Through population genetic simulations and comparison with wolf genomes from a variety of demographic histories, we find evidence that severe inbreeding depression in this population is due to increased homozygosity of strongly deleterious recessive mutations. Our results have particular relevance in light of the recent translocation of wolves from the mainland to Isle Royale, as well as broader implications for management of genetic variation in the fragmented landscape of the modern world.

The intrepid urban coyote: a comparison of bold and exploratory behavior in coyotes from urban and rural environments

Coyotes (Canis latrans) are highly adaptable, medium-sized carnivores that now inhabit nearly every large city in the United States and Canada. To help understand how coyotes have adapted to living in urban environments, we compared two ecologically and evolutionarily important behavioral traits (i.e., bold-shy and exploration-avoidance behavior) in two contrasting environments (i.e., rural and urban). Boldness is an individual’s reaction to a risky situation and exploration is an individual’s willingness to explore novel situations. Our results from both tests indicate that urban coyotes are bolder and more exploratory than rural coyotes and that within both populations there are individuals that vary across both spectrums. Bolder behavior in urban coyotes emerged over several decades and we speculate on possible processes (e.g., learning and selection) and site differences that could be playing a role in this behavioral adaptation. We hypothesize that an important factor is how people treat coyotes; in the rural area coyotes were regularly persecuted whereas in the urban area coyotes were rarely persecuted and sometimes positively rewarded to be in close proximity of people. Negative consequences of this behavioral adaptation are coyotes that become bold enough to occasionally prey on pets or attack humans.

Genetic differentiation and asymmetric gene flow among Carpathian brown bear (Ursus arctos) populations-Implications for conservation of transboundary populations

The abundance and distribution of large carnivores in Europe have been historically reduced. Their recovery requires multilevel coordination, especially regarding transboundary populations. Here, we apply nuclear and mitochondrial genetic markers to test for admixture level and its impact on population genetic structure of contemporary brown bears (Ursus arctos) from the Eastern, Southern, and Western Carpathians. Carpathian Mountains (Europe). Nearly 400 noninvasive brown bear DNA samples from the Western (Poland) and Eastern Carpathians (Bieszczady Mountains in Poland, Slovakia, Ukraine) were collected. Together with DNA isolates from Slovakia and Romania, they were analyzed using the set of eight microsatellite loci and two mtDNA regions (control region and cytochrome b). A set of 113 individuals with complete genotypes was used to investigate genetic differentiation across national boundaries, genetic structuring within and between populations, and movement between populations. Transboundary brown bear subpopulations (Slovakia and Poland) did not show significant internal genetic structure, and thus were treated as cohesive units. All brown bears from the Western Carpathians carried mitochondrial haplotypes from the Eastern lineage, while the Western lineage prevailed in the brown bears from the Bieszczady Mountains. Despite similar levels of microsatellite variability, we documented significant differentiation among the studied populations for nuclear markers and mtDNA. We also detected male-biased and asymmetrical movement into the Bieszczady Mountains population from the Western Carpathians. Our findings suggest initial colonization of the Western Carpathians by brown bears possessing mtDNA from the Eastern lineage. Genetic structuring among populations at microsatellite loci could be a result of human-mediated alterations. Detected asymmetric gene flow suggests ongoing expansion from more abundant populations into the Bieszczady Mountains and thus supports a metapopulation model. The knowledge concerning this complex pattern can be implemented in a joint Carpathian brown bear management plan that should allow population mixing by dispersing males.

Purging of Strongly Deleterious Mutations Explains Long-Term Persistence and Absence of Inbreeding Depression in Island Foxes

The recovery and persistence of rare and endangered species are often threatened by genetic factors, such as the accumulation of deleterious mutations, loss of adaptive potential, and inbreeding depression [1]. Island foxes (Urocyon littoralis), the dwarfed descendants of mainland gray foxes (Urocyon cinereoargenteus), have inhabited California's Channel Islands for >9,000 years [2-4]. Previous genomic analyses revealed that island foxes have exceptionally low levels of diversity and elevated levels of putatively deleterious variation [5]. Nonetheless, all six populations have persisted for thousands of generations, and several populations rebounded rapidly after recent severe bottlenecks [6, 7]. Here, we combine morphological and genomic data with population-genetic simulations to determine the mechanism underlying the enigmatic persistence of these foxes. First, through analysis of genomes from 1929 to 2009, we show that island foxes have remained at small population sizes with low diversity for many generations. Second, we present morphological data indicating an absence of inbreeding depression in island foxes, confirming that they are not afflicted with congenital defects common to other small and inbred populations. Lastly, our population-genetic simulations suggest that long-term small population size results in a reduced burden of strongly deleterious recessive alleles, providing a mechanism for the absence of inbreeding depression in island foxes. Importantly, the island fox illustrates a scenario in which genetic restoration through human-assisted gene flow could be a counterproductive or even harmful conservation strategy. Our study sheds light on the puzzle of island fox persistence, a unique success story that provides a model for the preservation of small populations.

Impact of supplementation on deleterious mutation distribution in an exploited salmonid

Deleterious mutations have important implications for the evolutionary trajectories of populations. While several studies recently investigated the dynamics of deleterious mutations in wild populations, no study has yet explored the fate of deleterious mutations in a context of populations managed by supplementation. Here, based on a dataset of nine wild and 15 supplemented Lake Trout populations genotyped at 4,982 single nucleotide polymorphisms (SNP)s by means of genotype by sequencing (GBS), we explored the effect of supplementation on the frequency of putatively deleterious variants. Three main findings are consequential for the management of fish populations. First, an increase in neutral genetic diversity in stocked populations compared with unstocked ones was observed. Second, putatively deleterious mutations were more likely to be found in unstocked than in stocked populations, suggesting a lower efficiency to purge deleterious mutations in unstocked lakes. Third, a population currently used as a major source for supplementation is characterized by several fixed putatively deleterious alleles. Therefore, other source populations with lower abundance of putatively deleterious mutations should be favored as sources of supplementation. We discuss management implications of our results, especially pertaining to the joint identification of neutral and deleterious mutations that could help refining the choice of source and sink populations for supplementation in order to maximize their evolutionary potential and to limit their mutation load.

Wolf outside, dog inside? The genomic make-up of the Czechoslovakian Wolfdog

BACKGROUND

Genomic methods can provide extraordinary tools to explore the genetic background of wild species and domestic breeds, optimize breeding practices, monitor and limit the spread of recessive diseases, and discourage illegal crossings. In this study we analysed a panel of 170k Single Nucleotide Polymorphisms with a combination of multivariate, Bayesian and outlier gene approaches to examine the genome-wide diversity and inbreeding levels in a recent wolf x dog cross-breed, the Czechoslovakian Wolfdog, which is becoming increasingly popular across Europe.

RESULTS

Pairwise FST values, multivariate and assignment procedures indicated that the Czechoslovakian Wolfdog was significantly differentiated from all the other analysed breeds and also well-distinguished from both parental populations (Carpathian wolves and German Shepherds). Coherently with the low number of founders involved in the breed selection, the individual inbreeding levels calculated from homozygosity regions were relatively high and comparable with those derived from the pedigree data. In contrast, the coefficient of relatedness between individuals estimated from the pedigrees often underestimated the identity-by-descent scores determined using genetic profiles. The timing of the admixture and the effective population size trends estimated from the LD patterns reflected the documented history of the breed. Ancestry reconstruction methods identified more than 300 genes with excess of wolf ancestry compared to random expectations, mainly related to key morphological features, and more than 2000 genes with excess of dog ancestry, playing important roles in lipid metabolism, in the regulation of circadian rhythms, in learning and memory processes, and in sociability, such as the COMT gene, which has been described as a candidate gene for the latter trait in dogs.

CONCLUSIONS

In this study we successfully applied genome-wide procedures to reconstruct the history of the Czechoslovakian Wolfdog, assess individual wolf ancestry proportions and, thanks to the availability of a well-annotated reference genome, identify possible candidate genes for wolf-like and dog-like phenotypic traits typical of this breed, including commonly inherited disorders. Moreover, through the identification of ancestry-informative markers, these genomic approaches could provide tools for forensic applications to unmask illegal crossings with wolves and uncontrolled trades of recent and undeclared wolfdog hybrids.