Genetic studies of rare and endangered species often focus on defining and preserving genetically distinct populations, especially those having unique adaptations [1, 2]. Much less attention is directed at understanding the landscape of deleterious variation, an insidious consequence of geographic isolation and the inefficiency of natural selection to eliminate harmful variants in small populations [3-5]. With population sizes of many vertebrates decreasing and isolation increasing through habitat fragmentation and loss, understanding the extent and nature of deleterious variation in small populations is essential for predicting and enhancing population persistence. The Channel Island fox (Urocyon littoralis) is a dwarfed species that inhabits six of California's Channel Islands and is derived from the mainland gray fox (U. cinereoargenteus). These isolated island populations have persisted for thousands of years at extremely small population sizes [6, 7] and, consequently, are a model for testing ideas about the accumulation of deleterious variation in small populations under natural conditions. Analysis of complete genome sequence data from island foxes shows a dramatic decrease in genome-wide variation and a sharp increase in the homozygosity of deleterious variants. The San Nicolas Island population has a near absence of variation, demonstrating a unique genetic flatlining that is punctuated by heterozygosity hotspots, enriched for olfactory receptor genes and other genes with high levels of ancestral variation. These findings question the generality of the small-population paradigm that maintains substantial genetic variation is necessary for short- and long-term persistence.

The golden jackal of Africa (Canis aureus) has long been considered a conspecific of jackals distributed throughout Eurasia, with the nearest source populations in the Middle East. However, two recent reports found that mitochondrial haplotypes of some African golden jackals aligned more closely to gray wolves (Canis lupus), which is surprising given the absence of gray wolves in Africa and the phenotypic divergence between the two species. Moreover, these results imply the existence of a previously unrecognized phylogenetically distinct species despite a long history of taxonomic work on African canids. To test the distinct-species hypothesis and understand the evolutionary history that would account for this puzzling result, we analyzed extensive genomic data including mitochondrial genome sequences, sequences from 20 autosomal loci (17 introns and 3 exon segments), microsatellite loci, X- and Y-linked zinc-finger protein gene (ZFX and ZFY) sequences, and whole-genome nuclear sequences in African and Eurasian golden jackals and gray wolves. Our results provide consistent and robust evidence that populations of golden jackals from Africa and Eurasia represent distinct monophyletic lineages separated for more than one million years, sufficient to merit formal recognition as different species: C. anthus (African golden wolf) and C. aureus (Eurasian golden jackal). Using morphologic data, we demonstrate a striking morphologic similarity between East African and Eurasian golden jackals, suggesting parallelism, which may have misled taxonomists and likely reflects uniquely intense interspecific competition in the East African carnivore guild. Our study shows how ecology can confound taxonomy if interspecific competition constrains size diversification.

Summary 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 (U. cinereoargenteus), have inhabited California’s Channel Islands for >9,000 years [2-4]. Previous genomic analyses revealed 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 following 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-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 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.

Genome admixture in two endemic North American wolf species.

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.

19The observation that small, isolated populations often suffer reduced fitness as a result of 20 inbreeding depression has guided conservation theory and practice for decades. However, 21investigating the genome-wide dynamics associated with inbreeding depression in natural 22 populations is only now feasible with relatively inexpensive sequencing technology and 23 annotated reference genomes. To characterize the genome-wide effects of intense inbreeding and 24 isolation, we sequenced complete genomes from an iconic inbred population, the gray wolves 25 (Canis lupus) of Isle Royale. Through comparison with other wolf genomes from a variety of 26 demographic histories, we found that Isle Royale wolf genomes contain extensive runs of 27 homozygosity, but neither the overall level of heterozygosity nor the number of deleterious 28 variants per genome were reliable predictors of inbreeding depression. These findings are 29 consistent with the hypothesis that severe inbreeding depression results from increased 30 homozygosity of strongly deleterious recessive mutations, which are more prevalent in 31 historically large source populations. Our results have particular relevance in light of the recently 32

In cases of severe wildlife population decline, a key question is whether recovery efforts will be impeded by genetic factors, such as inbreeding depression. Decades of excess mortality from gillnet fishing have driven Mexico’s vaquita porpoise ( Phocoena sinus ) to ~10 remaining individuals. We analyzed whole-genome sequences from 20 vaquitas and integrated genomic and demographic information into stochastic, individual-based simulations to quantify the species’ recovery potential. Our analysis suggests that the vaquita’s historical rarity has resulted in a low burden of segregating deleterious variation, reducing the risk of inbreeding depression. Similarly, genome-informed simulations suggest that the vaquita can recover if bycatch mortality is immediately halted. This study provides hope for vaquitas and other naturally rare endangered species and highlights the utility of genomics in predicting extinction risk.

The vaquita is the most critically endangered marine mammal, with fewer than 19 remaining in the wild. First described in 1958, the vaquita has been in rapid decline for more than 20 years resulting from inadvertent deaths due to the increasing use of large‐mesh gillnets. To understand the evolutionary and demographic history of the vaquita, we used combined long‐read sequencing and long‐range scaffolding methods with long‐ and short‐read RNA sequencing to generate a near error‐free annotated reference genome assembly from cell lines derived from a female individual. The genome assembly consists of 99.92% of the assembled sequence contained in 21 nearly gapless chromosome‐length autosome scaffolds and the X‐chromosome scaffold, with a scaffold N50 of 115 Mb. Genome‐wide heterozygosity is the lowest (0.01%) of any mammalian species analysed to date, but heterozygosity is evenly distributed across the chromosomes, consistent with long‐term small population size at genetic equilibrium, rather than low diversity resulting from a recent population bottleneck or inbreeding. Historical demography of the vaquita indicates long‐term population stability at less than 5,000 (Ne) for over 200,000 years. Together, these analyses indicate that the vaquita genome has had ample opportunity to purge highly deleterious alleles and potentially maintain diversity necessary for population health.