Demographic history shapes North American gray wolf genomic diversity and informs species' conservation.

Effective population size estimates are critical information needed for evolutionary predictions and conservation decisions. This is particularly true for species with social factors that restrict access to breeding or experience repeated fluctuations in population size across generations. We investigated the genomic estimates of effective population size along with diversity, subdivision, and inbreeding from 162,109 minimally filtered and 81,595 statistically neutral and unlinked SNPs genotyped in 437 grey wolf samples from North America collected between 1986 and 2021. We found genetic structure across North America, represented by three distinct demographic histories of western, central, and eastern regions of the continent. Further, grey wolves in the northern Rocky Mountains have lower genomic diversity than wolves of the western Great Lakes and have declined over time. Effective population size estimates revealed the historical signatures of continental efforts of predator extermination, despite a quarter century of recovery efforts. We are the first to provide molecular estimates of effective population size across distinct grey wolf populations in North America, which ranged between Ne ~ 275 and 3050 since early 1980s. We provide data that inform managers regarding the status and importance of effective population size estimates for grey wolf conservation, which are on average 5.2-9.3% of census estimates for this species. We show that while grey wolves fall above minimum effective population sizes needed to avoid extinction due to inbreeding depression in the short term, they are below sizes predicted to be necessary to avoid long-term risk of extinction.

Large variance in inbreeding within the Iberian wolf population.

The gray wolf (Canis lupus) population on the Iberian Peninsula was the largest in western and central Europe during most of the 20th century, with its size apparently never under a few hundred individuals. After partial legal protection in the 1970s in Spain, the northwest Iberian population increased to about 300-350 packs and then stabilized. In contrast to many current European wolf populations, which have been connected through gene flow, the Iberian wolf population has been isolated for decades. Here we measured changes on genomic diversity and inbreeding through the last decades in a geographic context. We find that the level of genomic diversity in Iberian wolves is low compared to other Eurasian wolf populations. Despite population expansion in the last 50 years, some modern wolves had very high inbreeding, especially in the recently recolonized and historical edge areas. These individuals contrast with others with low inbreeding within the same population. The high variance in inbreeding despite population expansion seems associated with small-scale fragmentation of the range that is revealed by the genetic similarity between modern and historical samples from close localities despite being separated by decades, remaining differentiated from other individuals that are just over 100 km away, a small distance for a species with great dispersal capacity inhabiting a continuous range. This illustrates that, despite its demographically stable condition, the population would probably benefit from favoring connectivity within the population as well as genetic exchange with other European wolf populations to avoid excessive fragmentation and local inbreeding depression.

Evolutionary history of an island endemic, the Azorean common quail.

Oceanic islands are characterized by conditions that favour diversification into endemic lineages that can be very different from their mainland counterparts. This can be the result of fast phenotypic divergence due to drift or the result of slower adaptation to local conditions. This uniqueness can obscure their evolutionary history. Here we used morphological, stable isotope, genetic and genomic data to characterize common quails (Coturnix coturnix) in the Azores archipelago and assess the divergence from neighbouring common quail populations. Historical documents suggested that these quails could have a recent origin associated with the arrival of humans in the last centuries. Our results show that Azorean quails constitute a well-differentiated lineage with small size and dark throat pigmentation that has lost the migratory ability and that diverged from mainland quail lineages more than 0.8 mya, contrary to the notion of a recent human-mediated arrival. Even though some Azorean quails carry an inversion that affects 115 Mbp of chromosome 1 and that has been associated with the loss of the migratory behaviour in other common quail populations, half of the analysed individuals do not have that inversion and still do not migrate. The long coexistence and evolution in isolation in the Azores of two chromosomal variants (with and without the inversion) is best explained by balancing selection. Thus, a unique and long evolutionary history led to the island endemic that we know today, C. c. conturbans.

Early dispersal of domestic horses into the Great Plains and northern Rockies.

The horse is central to many Indigenous cultures across the American Southwest and the Great Plains. However, when and how horses were first integrated into Indigenous lifeways remain contentious, with extant models derived largely from colonial records. We conducted an interdisciplinary study of an assemblage of historic archaeological horse remains, integrating genomic, isotopic, radiocarbon, and paleopathological evidence. Archaeological and modern North American horses show strong Iberian genetic affinities, with later influx from British sources, but no Viking proximity. Horses rapidly spread from the south into the northern Rockies and central plains by the first half of the 17th century CE, likely through Indigenous exchange networks. They were deeply integrated into Indigenous societies before the arrival of 18th-century European observers, as reflected in herd management, ceremonial practices, and culture.

Challenging ecogeographical rules: Phenotypic variation in the Mountain Treeshrew (Tupaia montana) along tropical elevational gradients.

Bergmann's and Allen's rules were defined to describe macroecological patterns across latitudinal gradients. Bergmann observed a positive association between body size and latitude for endothermic species while Allen described shorter appendages as latitude increases. Almost two centuries later, there is still ongoing discussion about these patterns. Temperature, the common variable in these two rules, varies predictably across both latitude and elevation. Although these rules have been assessed extensively in mammals across latitude, particularly in regions with strong seasonality, studies on tropical montane mammals are scarce. We here test for these patterns and assess the variation of several other locomotory, diet-associated, body condition, and thermoregulatory traits across elevation in the Mountain Treeshrew (Tupaia montana) on tropical mountains in Borneo. Based on morphological measurements from both the field and scientific collections, we found a complex pattern: Bergmann's rule was not supported in our tropical mountain system, since skull length, body size, and weight decreased from the lowest elevations (<1000 m) to middle elevations (2000-2500 m), and then increased from middle elevations to highest elevations. Allen's rule was supported for relative tail length, which decreased with elevation, but not for ear and hindfoot length, with the former remaining constant and the latter increasing with elevation. This evidence together with changes in presumed diet-related traits (rostrum length, zygomatic breadth and upper tooth row length) along elevation suggest that selective pressures other than temperature, are playing a more important role shaping the morphological variation across the distribution of the Mountain Treeshrew. Diet, food acquisition, predation pressure, and/or intra- and inter-specific competition, are some of the potential factors driving the phenotypic variation of this study system. The lack of variation in body condition might suggest local adaptation of this species across its elevational range, perhaps due to generalist foraging strategies. Finally, a highly significant temporal effect was detected in several traits but not in others, representing the first phenotypic variation temporal trends described on treeshrews.

Massive genome inversion drives coexistence of divergent morphs in common quails.

The presence of population-specific phenotypes often reflects local adaptation or barriers to gene flow. The co-occurrence of phenotypic polymorphisms that are restricted within the range of a highly mobile species is more difficult to explain. An example of such polymorphisms is in the common quail Coturnix coturnix, a small migratory bird that moves widely during the breeding season in search of new mating opportunities, following ephemeral habitats,1,2 and whose females may lay successive clutches at different locations while migrating.3 In spite of this vagility, previous studies reported a higher frequency of heavier males with darker throat coloration in the southwest of the distribution (I. Jiménez-Blasco et al., 2015, Int. Union Game Biol., conference). We used population genomics and cytogenetics to explore the basis of this polymorphism and discovered a large inversion in the genome of the common quail. This inversion extends 115 Mbp in length and encompasses more than 7,000 genes (about 12% of the genome), producing two very different forms. Birds with the inversion are larger, have darker throat coloration and rounder wings, are inferred to have poorer flight efficiency, and are geographically restricted despite the high mobility of the species. Stable isotope analyses confirmed that birds carrying the inversion have shorter migratory distances or do not migrate. However, we found no evidence of pre- or post-zygotic isolation, indicating the two forms commonly interbreed and that the polymorphism remains locally restricted because of the effect on behavior. This illustrates a genomic mechanism underlying maintenance of geographically structured polymorphisms despite interbreeding with a lineage with high mobility.

Loss of Mitochondrial Genetic Diversity despite Population Growth: The Legacy of Past Wolf Population Declines.

Gray wolves (Canis lupus) in the Iberian Peninsula declined substantially in both range and population size in the last few centuries due to human persecution and habitat fragmentation. However, unlike many other western European populations, gray wolves never went extinct in Iberia. Since the minimum number was recorded around 1970, their numbers have significantly increased and then stabilized in recent decades. We analyzed mitochondrial genomes from 54 historical specimens of Iberian wolves from across their historical range using ancient DNA methods. We compared historical and current mitochondrial diversity in Iberian wolves at the 5' end of the control region (n = 17 and 27) and the whole mitochondrial genome excluding the control region (n = 19 and 29). Despite an increase in population size since the 1970s, genetic diversity declined. We identified 10 whole mitochondrial DNA haplotypes in 19 historical specimens, whereas only six of them were observed in 29 modern Iberian wolves. Moreover, a haplotype that was restricted to the southern part of the distribution has gone extinct. Our results illustrate a lag between demographic and genetic diversity changes, and show that after severe population declines, genetic diversity can continue to be lost in stable or even expanding populations. This suggests that such populations may be of conservation concern even after their demographic trajectory has been reversed.

Automated genotyping of microsatellite loci from feces with high throughput sequences.

Ecological and conservation genetic studies often use noninvasive sampling, especially with elusive or endangered species. Because microsatellites are generally short in length, they can be amplified from low quality samples such as feces. Microsatellites are highly polymorphic so few markers are enough for reliable individual identification, kinship determination, or population characterization. However, the genotyping process from feces is expensive and time consuming. Given next-generation sequencing (NGS) and recent software developments, automated microsatellite genotyping from NGS data may now be possible. These software packages infer the genotypes directly from sequence reads, increasing throughput. Here we evaluate the performance of four software packages to genotype microsatellite loci from Iberian wolf (Canis lupus) feces using NGS. We initially combined 46 markers in a single multiplex reaction for the first time, of which 19 were included in the final analyses. Megasat was the software that provided genotypes with fewer errors. Coverage over 100X provided little additional information, but a relatively high number of PCR replicates were necessary to obtain a high quality genotype from highly unoptimized, multiplexed reactions (10 replicates for 18 of the 19 loci analyzed here). This could be reduced through optimization. The use of new bioinformatic tools and next-generation sequencing data to genotype these highly informative markers may increase throughput at a reasonable cost and with a smaller amount of laboratory work. Thus, high throughput sequencing approaches could facilitate the use of microsatellites with fecal DNA to address ecological and conservation questions.

A sliver of the past: The decimation of the genetic diversity of the Mexican wolf.

The endangered Mexican wolf (Canis lupus baileyi) is known to carry exceedingly low levels of genetic diversity. This could be (i) the result of long-term evolutionary patterns as they exist at the southernmost limit of the species distribution at a relatively reduced effective size, or (ii) due to rapid population decline caused by human persecution over the last century. If the former, purifying selection is expected to have minimized the impact of inbreeding. If the latter, rapid and recent declines in genetic diversity may have resulted in severe fitness consequences. To differentiate these hypotheses, we conducted comparative whole-genome analyses of five historical Mexican wolves (1907-1917) and 18 contemporary Mexican and grey wolves from North America and Eurasia. Based on whole-genome data, historical and modern Mexican wolves together form a discrete unit. Moreover, we found that modern Mexican wolves have reduced genetic diversity and increased inbreeding relative to the historical population, which was widespread across the southwestern United States and not restricted to Mexico as previously assumed. Finally, although Mexican wolves have evolved in sympatry with coyotes (C. latrans), we observed lower introgression between historical Mexican wolves and coyotes than with modern Mexican wolves, despite similarities in body size. Taken together, our data show that recent population declines probably caused the reduced level of genetic diversity, but not the observed differentiation of the Mexican wolves from other North American wolves.

Phylogenomics and evolutionary history of Oreobates (Anura: Craugastoridae) Neotropical frogs along elevational gradients.

Mountain ranges offer opportunities for understanding how species evolved and diversified across different environmental conditions. Neotropical frogs of the genus Oreobates (Anura: Craugastoridae) are adapted to highland and lowland habitats along the Andes, but many aspects of their evolution remain unknown. We studied their evolutionary history using ~18,000 exons enriched by targeted sequence-capture. Since capture success was very variable across samples, we evaluated to what degree differing data filtering produced robust inferences. The inferred evolutionary framework evidenced phylogenetic discordances among lowland species that can be explained by taxonomic misidentification or admixture of ancestral lineages. Highland species showed smaller effective populations than lowland frogs, probably due to greater habitat fragmentation in montane environments. Stronger genetic drift likely decreased the power of purifying selection and led to an increased proportion of nonsynonymous mutations in highland populations that could play an important role in their adaptation. Overall, our work sheds light on the evolutionary history and diversification of this group of Neotropical frogs along elevational gradients in the Andes as well as on their patterns of intraspecific diversity.

Markers for genetic change.

BACKGROUND AND AIMS: Wildlife conservation has focused primarily on species for the last decades. Recently, popular perception and laws have begun to recognize the central importance of genetic diversity in the conservation of biodiversity. How to incorporate genetic diversity in ongoing monitoring and management of wildlife is still an open question. METHODS: We tested a panel of multiplexed, high-throughput sequenced introns in the small mammal communities of two UNESCO World Heritage Sites on different continents to assess their viability for large-scale monitoring of genetic variability in a spectrum of diverse species. To enhance applicability across other systems, the bioinformatic pipeline for primer design was outlined. RESULTS: The number of loci amplified and amplification evenness decreased as phylogenetic distance increased from the reference taxa, yet several loci were still variable across multiple mammal orders. CONCLUSIONS: Genetic variability found is informative for population genetic analyses and for addressing phylogeographic and phylogenetic questions, illustrated by small mammal examples here.

Pleistocene climate fluctuations drove demographic history of African golden wolves (Canis lupaster).

Pleistocene climate change impacted entire ecosystems throughout the world. In the northern hemisphere, the distribution of Arctic species expanded during glacial periods, while more temperate and mesic species contracted into climatic refugia, where isolation drove genetic divergence. Cycles of local cooling and warming in the Sahara region of northern Africa caused repeated contractions and expansions of savannah-like environments which connected mesic species isolated in refugia during interglacial times, possibly driving population expansions and contractions; divergence and geneflow in the associated fauna. Here, we use whole genome sequences of African golden wolves (Canis lupaster), a generalist mesopredator with a wide distribution in northern Africa to estimate their demographic history and past episodes of geneflow. We detect a correlation between divergence times and cycles of increased aridity-associated Pleistocene glacial cycles. A complex demographic history with responses to local climate change in different lineages was found, including a relict lineage north of the High Atlas Mountains of Morocco that has been isolated for more than 18,000 years, possibly a distinct ecotype.

Little genetic structure in a Bornean endemic small mammal across a steep ecological gradient.

Janzen's influential "mountain passes are higher in the tropics" hypothesis predicts restricted gene flow and genetic isolation among populations spanning elevational gradients in the tropics. Few studies have tested this prediction, and studies that focus on population genetic structure in Southeast Asia are particularly underrepresented in the literature. Here, we test the hypothesis that mountain treeshrews (Tupaia montana) exhibit limited dispersal across their broad elevational range which spans ~2,300 m on two peaks in Kinabalu National Park (KNP) in Borneo: Mt Tambuyukon (MT) and Mt Kinabalu (MK). We sampled 83 individuals across elevations on both peaks and performed population genomics analyses on mitogenomes and single nucleotide polymorphisms from 4,106 ultraconserved element loci. We detected weak genetic structure and infer gene flow both across elevations and between peaks. We found higher genetic differentiation on MT than MK despite its lower elevation and associated environmental variation. This implies that, contrary to our hypothesis, genetic structure in this system is not primarily shaped by elevation. We propose that this pattern may instead be the result of historical processes and limited upslope gene flow on MT. Importantly, our results serve as a foundational estimate of genetic diversity and population structure from which to track potential future effects of climate change on mountain treeshrews in KNP, an important conservation stronghold for the mountain treeshrew and other montane species.

Spatiotemporal analyses suggest the role of glacial history and the ice-free corridor in shaping American badger population genetic variation.

Recurring glacial cycles through the Quaternary period drastically altered the size and distribution of natural populations of North American flora and fauna. The "southerly refugia model" has been the longstanding framework for testing the effects of glaciation on contemporary genetic patterns; however, insights from ancient DNA have contributed to the reconstruction of more complex histories for some species. The American badger, Taxidea taxus, provides an interesting species for exploring the genetic legacy of glacial history, having been hypothesized to have postglacially emerged from a single, southerly refugium to recolonize northern latitudes. However, previous studies have lacked genetic sampling from areas where distinct glacial refugia have been hypothesized, including the Pacific Northwest and American Far North (Yukon, Alaska). In order to further investigate the phylogeographic history of American badgers, we collected mitochondrial DNA sequence data from ancient subfossil material collected within the historical range (Alaska, Yukon) and combined them with new and previously published data from across the species' contemporary distribution (n = 1,207). We reconstructed a mostly unresolved phylogenetic tree and star-like haplotype network indicative of emergence from a largely panmictic glacial refugium and recent population expansion, the latter further punctuated by significantly negative Tajima's D and Fu's Fs values. Although directionality of migration cannot be unequivocally inferred, the moderate to high levels of genetic variation exhibited by American badgers, alongside the low frequency of haplotypes with indels in the Midwest, suggest a potential recolonization into central North America after the hypothesized ice-free corridor reopened ~13,000 years ago. Overall, the expanded reconstruction of phylogeographic history of American badgers offers a broader understanding of contemporary range-wide patterns and identifies unique genetic units that can likely be used to inform conservation of at-risk populations at the northern periphery.

Mitogenomes Reveal Multiple Colonization of Mountains by Rattus in Sundaland.

Tropical mountains are cradles of biodiversity and endemism. Sundaland, tropical Southeast Asia, hosts 3 species of Rattus endemic to elevations above 2000 m with an apparent convergence in external morphology: Rattus korinchi and R. hoogerwerfi from Sumatra, and R. baluensis from Borneo. A fourth one, R. tiomanicus, is restricted to lowland elevations across the whole region. The origins of these endemics are little known due to the absence of a robust phylogenetic framework. We use complete mitochondrial genomes from the 3 high altitude Rattus, and several related species to determine their relationships, date divergences, reconstruct their history of colonization, and test for selection on the mitochondrial DNA. We show that mountain colonization happened independently in Borneo (<390 Kya) and Sumatra (~1.38 Mya), likely from lowland lineages. The origin of the Bornean endemic R. baluensis is very recent and its genetic diversity is nested within the diversity of R. tiomanicus. We found weak evidence of positive selection in the high-elevation lineages and attributed the greater nonsynonymous mutations on these branches (specially R. baluensis) to lesser purifying selection having acted on the terminal branches in the phylogeny.

Endemism and diversity of small mammals along two neighboring Bornean mountains.

Mountains offer replicated units with large biotic and abiotic gradients in a reduced spatial scale. This transforms them into well-suited scenarios to evaluate biogeographic theories. Mountain biogeography is a hot topic of research and many theories have been proposed to describe the changes in biodiversity with elevation. Geometric constraints, which predict the highest diversity to occur in mid-elevations, have been a focal part of this discussion. Despite this, there is no general theory to explain these patterns, probably because of the interaction among different predictors with the local effects of historical factors. We characterize the diversity of small non-volant mammals across the elevational gradient on Mount (Mt.) Kinabalu (4,095 m) and Mt. Tambuyukon (2,579 m), two neighboring mountains in Borneo, Malaysia. We documented a decrease in species richness with elevation which deviates from expectations of the geometric constraints and suggests that spatial factors (e.g., larger diversity in larger areas) are important. The lowland small mammal community was replaced in higher elevations (from above ~1,900 m) with montane communities consisting mainly of high elevation Borneo endemics. The positive correlation we find between elevation and endemism is concordant with a hypothesis that predicts higher endemism with topographical isolation. This supports lineage history and geographic history could be important drivers of species diversity in this region.

On the path to extinction: Inbreeding and admixture in a declining grey wolf population.

Allee effects reduce the viability of small populations in many different ways, which act synergistically to lead populations towards extinction vortexes. The Sierra Morena wolf population, isolated in the south of the Iberian Peninsula and composed of just one or few packs for decades, represents a good example of how diverse threats act additively in very small populations. We sequenced the genome of one of the last wolves identified (and road-killed) in Sierra Morena and that of another wolf in the Iberian Wolf Captive Breeding Program and compared them with other wolf and dog genomes from around the world (including two previously published genome sequences from northern Iberian wolves). The results showed relatively low overall genetic diversity in Iberian wolves, but diverse population histories including past introgression of dog genes. The Sierra Morena wolf had an extraordinarily high level of inbreeding and long runs of homozygosity, resulting from the long isolation. In addition, about one-third of the genome was of dog origin. Despite the introgression of dog genes, heterozygosity remained low because of continued inbreeding after several hybridization events. The results thus illustrate the case of a small and isolated wolf population where the low population density may have favoured hybridization and introgression of dog alleles, but continued inbreeding may have resulted in large chromosomal fragments of wolf origin completely disappearing from the population, and being replaced by chromosomal fragments of dog origin. The latest population surveys suggest that this population may have gone extinct.

Evolutionary History of Saber-Toothed Cats Based on Ancient Mitogenomics.

Saber-toothed cats (Machairodontinae) are among the most widely recognized representatives of the now largely extinct Pleistocene megafauna. However, many aspects of their ecology, evolution, and extinction remain uncertain. Although ancient-DNA studies have led to huge advances in our knowledge of these aspects of many other megafauna species (e.g., mammoths and cave bears), relatively few ancient-DNA studies have focused on saber-toothed cats [1-3], and they have been restricted to short fragments of mitochondrial DNA. Here we investigate the evolutionary history of two lineages of saber-toothed cats (Smilodon and Homotherium) in relation to living carnivores and find that the Machairodontinae form a well-supported clade that is distinct from all living felids. We present partial mitochondrial genomes from one S. populator sample and three Homotherium sp. samples, including the only Late Pleistocene Homotherium sample from Eurasia [4]. We confirm the identification of the unique Late Pleistocene European fossil through ancient-DNA analyses, thus strengthening the evidence that Homotherium occurred in Europe over 200,000 years later than previously believed. This in turn forces a re-evaluation of its demography and extinction dynamics. Within the Machairodontinae, we find a deep divergence between Smilodon and Homotherium (∼18 million years) but limited diversity between the American and European Homotherium specimens. The genetic data support the hypothesis that all Late Pleistocene (or post-Villafrancian) Homotherium should be considered a single species, H. latidens, which was previously proposed based on morphological data [5, 6].