Whole-genome analysis of Mustela erminea finds that pulsed hybridization impacts evolution at high latitudes

Tracing evolutionary trajectories in the presence of gene flow in South American temperate lizards (Squamata: Liolaemus kingii group)

Evolutionary processes behind lineage divergence often involve multidimensional differentiation. However, in the context of recent divergences, the signals exhibited by each dimension may not converge. In such scenarios, incomplete lineage sorting, gene flow, and scarce phenotypic differentiation are pervasive. Here, we integrated genomic (RAD loci of 90 individuals), phenotypic (linear and geometric traits of 823 and 411 individuals, respectively), spatial, and climatic data to reconstruct the evolutionary history of a speciation continuum of liolaemid lizards (Liolaemus kingii group). Specifically, we (a) inferred the population structure of the group and contrasted it with the phenotypic variability; (b) assessed the role of postdivergence gene flow in shaping phylogeographic and phenotypic patterns; and (c) explored ecogeographic drivers of diversification across time and space. We inferred eight genomic clusters exhibiting leaky genetic borders coincident with geographic transitions. We also found evidence of postdivergence gene flow resulting in transgressive phenotypic evolution in one species. Predicted ancestral niches unveiled suitable areas in southern and eastern Patagonia during glacial and interglacial periods. Our study underscores integrating different data and model-based approaches to determine the underlying causes of diversification, a challenge faced in the study of recently diverged groups. We also highlight Liolaemus as a model system for phylogeographic and broader evolutionary studies.

Evaluating hybrid speciation and swamping in wild carnivores with a decision-tree approach

Hybridization is an important evolutionary force with a principal role in the origin of new species, known as hybrid speciation. However, ongoing hybridization can create hybrid swamping, in which parental genomes are completely lost. This can become a biodiversity threat if it involves species that have adapted to certain environmental conditions and occur nowhere else. Because conservation scientists commonly have a negative attitude toward hybrids, it is important to improve understanding of the influence of interspecific gene flow on the persistence of species. We reviewed the literature on species hybridization to build a list of all known cases in the order Carnivora. To examine the relative impact, we also noted level of introgression, whether fertile offspring were produced, and whether there was mention of negative or positive evolutionary effects (hybrid speciation and swamping). To evaluate the conservation implications of hybrids, we developed a decision-making tree with which to determine which actions should be taken to manage hybrid species. We found 53 hybrids involving 68 unique taxa, which is roughly 23% of all carnivore species. They mainly involved monophyletic (83%) and sympatric species (75%). For 2 species, the outcome of the assessment was to eliminate or restrict the hybrids: Ethiopian wolf (Canis simensis) and Scottish wildcat (Felis silvestris silvestris). Both species hybridize with their domestic conspecifics. For all other cases, we suggest hybrids be protected in the same manner as native species. We found no evidence of genomic extinction in Carnivora. To the contrary, some species appear to be of hybrid origin, such as the Asiatic black bear (Ursus thibetanus) and African golden wolf (Canis lupaster). Other positive outcomes of hybridization are novel genetic diversity, adaptation to extreme environments, and increased reproductive fitness. These outcomes are particularly valuable for counterbalancing genetic drift and enabling adaptive introgression in a human-dominated world.

Population genomics of the muskox' resilience in the near absence of genetic variation

Genomic studies of species threatened by extinction are providing crucial information about evolutionary mechanisms and genetic consequences of population declines and bottlenecks. However, to understand how species avoid the extinction vortex, insights can be drawn by studying species that thrive despite past declines. Here, we studied the population genomics of the muskox (Ovibos moschatus), an Ice Age relict that was at the brink of extinction for thousands of years at the end of the Pleistocene yet appears to be thriving today. We analysed 108 whole genomes, including present-day individuals representing the current native range of both muskox subspecies, the white-faced and the barren-ground muskox (O. moschatus wardi and O. moschatus moschatus) and a ~21,000-year-old ancient individual from Siberia. We found that the muskox' demographic history was profoundly shaped by past climate changes and post-glacial re-colonizations. In particular, the white-faced muskox has the lowest genome-wide heterozygosity recorded in an ungulate. Yet, there is no evidence of inbreeding depression in native muskox populations. We hypothesize that this can be explained by the effect of long-term gradual population declines that allowed for purging of strongly deleterious mutations. This study provides insights into how species with a history of population bottlenecks, small population sizes and low genetic diversity survive against all odds.

Analyzing the Evolution and Host Adaptation of the Rabies Virus from the Perspective of Codon Usage Bias

Rabies virus (RABV) is a highly pathogenic virus that causes a fatal disease in humans and other mammals, but the mechanism of its evolution, spread, and spillover remains unknown. In this study, we analyzed the codon usage pattern of 2,018 RABV full-length genome sequences from 79 countries collected between 1931 and 2021 to provide an insight into its molecular evolution and unravel its unknown host-adapted pattern. We found that RABV exhibited a weak codon usage bias, with a preference for the codons ending in A (28.10 ± 0.01) or U (26.43 ± 0.02). Moreover, natural selection plays a major role in shaping the codon usage bias of the RABV. Notably, nearly half of the 18 codons in the virus were best matched to the hosts' most abundant isoacceptor tRNAs, which might account for the wide range of RABV hosts. Furthermore, significant differences were observed in the codon usage patterns of RABV for different host species, suggesting that codon usage bias may be influenced by host-specific factors. In conclusion, our study reveals codon usage patterns of RABV that may help in the development of control strategies and effective vaccines and therapies against this deadly virus.

Surrounded by Kindred: Spermophilus major Hybridization with Other Spermophilus Species in Space and Time

Among the numerous described cases of hybridization in mammals, the most intriguing are (a) cases of introgressive hybridization deeply affecting the evolutionary history of species, and (b) models involving not a pair of species but a multi-species complex. Therefore, the hybridization history of the russet ground squirrel Spermophilus major, whose range has repeatedly changed due to climatic fluctuations and now borders the ranges of four related species, is of great interest. The main aims of this study were to determine the direction and intensity of gene introgression, the spatial depth of the infiltration of extraneous genes into the S. major range, and to refine the hypothesis of the hybridogenic replacement of mitochondrial genomes in the studied group. Using phylogenetic analysis of the variability of mitochondrial (CR, cytb) and nuclear (SmcY, BGN, PRKCI, c-myc, i6p53) markers, we determined the contribution of neighboring species to the S. major genome. We showed that 36% of S. major individuals had extraneous alleles. All peripheral species that were in contact with S. major contributed towards its genetic variability. We also proposed a hypothesis for the sequence and localization of serial hybridization events. Our assessment of the S. major genome implications of introgression highlights the importance of implementing conservation measures to protect this species.

Ancient bears provide insights into Pleistocene ice age refugia in Southeast Alaska

During the Late Pleistocene, major parts of North America were periodically covered by ice sheets. However, there are still questions about whether ice-free refugia were present in the Alexander Archipelago along the Southeast (SE) Alaska coast during the last glacial maximum (LGM). Numerous subfossils have been recovered from caves in SE Alaska, including American black (Ursus americanus) and brown (U. arctos) bears, which today are found in the Alexander Archipelago but are genetically distinct from mainland bear populations. Hence, these bear species offer an ideal system to investigate long-term occupation, potential refugial survival and lineage turnover. Here, we present genetic analyses based on 99 new complete mitochondrial genomes from ancient and modern brown and black bears spanning the last ~45,000 years. Black bears form two SE Alaskan subclades, one preglacial and another postglacial, that diverged >100,000 years ago. All postglacial ancient brown bears are closely related to modern brown bears in the archipelago, while a single preglacial brown bear is found in a distantly related clade. A hiatus in the bear subfossil record around the LGM and the deep split of their pre- and postglacial subclades fail to support a hypothesis of continuous occupancy in SE Alaska throughout the LGM for either species. Our results are consistent with an absence of refugia along the SE Alaska coast, but indicate that vegetation quickly expanded after deglaciation, allowing bears to recolonize the area after a short-lived LGM peak.

A comparison of microsatellites and genome-wide SNPs for the detection of admixture brings the first molecular evidence for hybridization between Mustela eversmanii and M. putorius (Mustelidae, Carnivora)

Introgressive hybridization can pose a serious threat to endangered species which have an overlapping distribution such as in the case of two polecat species, Mustela eversmanii and M. putorius, in Europe. The population size of steppe polecat is known to continuously shrink, whereas its sister species, the European polecat, is still somehow widespread. In this study, we perform an analysis using microsatellite (SSR) and genomic (SNP) data sets to identify natural hybrids between polecats. Four populations were genotyped for eight polymorphic SSR loci, and thousands of unlinked SNPs were generated using a reduced-representation sequencing approach, RADseq, to characterize the genetic make-up of allopatric populations and to identify hybrids in the sympatric area. We applied standard population genetic analyses to characterize the populations based on their SSR allelic frequency. Only a single sample out of 48 sympatric samples showed exact intermediacy that we identified as an F1 hybrid. Additionally, one specimen was indicated in the genomic data sets as backcrossed. Other backcrosses, indicated by SSRs, were not validated by SNPs, which highlights the higher efficacy of the genomic method to identify backcrossed individuals. The low frequency of hybridization suggests that the difference in habitat preference of the two species may act as a barrier to admixture. Therefore, it is apparently unlikely that polecat populations are threatened by significant introgression. The two species showed a clear genetic differentiation using both techniques. We found higher genetic diversity values in the sympatric steppe polecat population than in the other studies on polecat populations. Although M. putorius is a hunted species in most countries, genetic diversity values indicate worse conditions in Europe than in the protected sibling species M. eversmanii. Suspending hunting and providing protected status of the former seems to be reasonable and timely.

Phylogenomic systematics of the spotted skunks (Carnivora, Mephitidae, Spilogale): Additional species diversity and Pleistocene climate change as a major driver of diversification

Four species of spotted skunks (Carnivora, Mephitidae, Spilogale) are currently recognized: Spilogale angustifrons, S. gracilis, S. putorius, and S. pygmaea. Understanding species boundaries within this group is critical for effective conservation given that regional populations or subspecies (e.g., S. p. interrupta) have experienced significant population declines. Further, there may be currently unrecognized diversity within this genus as some taxa (e.g., S. angustifrons) and geographic regions (e.g., Central America) never have been assessed using DNA sequence data. We analyzed species limits and diversification patterns in spotted skunks using multilocus nuclear (ultraconserved elements) and mitochondrial (whole mitogenomes and single gene analysis) data sets from broad geographic sampling representing all currently recognized species and subspecies. We found a high degree of genetic divergence among Spilogale that reflects seven distinct species and eight unique mitochondrial lineages. Initial divergence between S. pygmaea and all other Spilogale occurred in the Early Pliocene (∼ 5.0 million years ago). Subsequent diversification of the remaining Spilogale into an "eastern" and a "western" lineage occurred during the Early Pleistocene (∼1.5 million years ago). These two lineages experienced temporally coincident patterns of diversification at ∼0.66 and ∼0.35 million years ago into two and ultimately three distinct evolutionary units, respectively. Diversification was confined almost entirely within the Pleistocene during a timeframe characterized by alternating glacial-interglacial cycles, with the origin of this diversity occurring in northeastern Mexico and the southwestern United States of America. Mitochondrial-nuclear discordance was recovered across three lineages in geographic regions consistent with secondary contact, including a distinct mitochondrial lineage confined to the Sonoran Desert. Our results have direct consequences for conservation of threatened populations, or species, as well as for our understanding of the evolution of delayed implantation in this enigmatic group of small carnivores.

Dynamic landscapes in northwestern North America structured populations of wolverines (Gulo gulo)

Cyclic climatic and glacial fluctuations of the Late Quaternary produced a dynamic biogeographic history for high latitudes. To refine our understanding of this history in northwestern North America, we explored geographic structure in a wide-ranging carnivore, the wolverine (Gulo gulo). We examined genetic variation in populations across mainland Alaska, coastal Southeast Alaska, and mainland western Canada using nuclear microsatellite genotypes and sequence data from the mitochondrial DNA (mtDNA) control region and Cytochrome b (Cytb) gene. Data from maternally inherited mtDNA reflect stable populations in Northwest Alaska, suggesting the region harbored wolverine populations since at least the Last Glacial Maximum (LGM; 21 Kya), consistent with their persistence in the fossil record of Beringia. Populations in Southeast Alaska are characterized by minimal divergence, with no genetic signature of long-term refugial persistence (consistent with the lack of pre-Holocene fossil records there). The Kenai Peninsula population exhibits mixed signatures depending on marker type: mtDNA data indicate stability (i.e., historical persistence) and include a private haplotype, whereas biparentally inherited microsatellites exhibit relatively low variation and a lack of private alleles consistent with a more recent Holocene colonization of the peninsula. Our genetic work is largely consistent with the early 20th century taxonomic hypothesis that wolverines on the Kenai Peninsula belong to a distinct subspecies. Our finding of significant genetic differentiation of wolverines inhabiting the Kenai Peninsula, coupled with the peninsula’s burgeoning human population and the wolverine’s known sensitivity to anthropogenic impacts, provides valuable foundational data that can be used to inform conservation and management prescriptions for wolverines inhabiting these landscapes.

Harmonizing hybridization dissonance in conservation

A dramatic increase in the hybridization between historically allopatric species has been induced by human activities. However, the notion of hybridization seems to lack consistency in two respects. On the one hand, it is inconsistent with the biological species concept, which does not allow for interbreeding between species, and on the other hand, it is considered either as an evolutionary process leading to the emergence of new biodiversity or as a cause of biodiversity loss, with conservation implications. In the first case, we argue that conservation biology should avoid the discussion around the species concept and delimit priorities of conservation units based on the impact on biodiversity if taxa are lost. In the second case, we show that this is not a paradox but an intrinsic property of hybridization, which should be considered in conservation programmes. We propose a novel view of conservation guidelines, in which human-induced hybridization may also be a tool to enhance the likelihood of adaptation to changing environmental conditions or to increase the genetic diversity of taxa affected by inbreeding depression. The conservation guidelines presented here represent a guide for the development of programmes aimed at protecting biodiversity as a dynamic evolutionary system.

Ecology and Infection Dynamics of Multi-Host Amdoparvoviral and Protoparvoviral Carnivore Pathogens

Amdoparvovirus and Protoparvovirus are monophyletic viral genera that infect carnivores. We performed surveillance for and sequence analyses of parvoviruses in mustelids in insular British Columbia to investigate parvoviral maintenance and cross-species transmission among wildlife. Overall, 19.1% (49/256) of the tested animals were parvovirus-positive. Aleutian mink disease virus (AMDV) was more prevalent in mink (41.6%, 32/77) than martens (3.1%, 4/130), feline panleukopenia virus (FPV) was more prevalent in otters (27.3%, 6/22) than mink (5.2%, 4/77) or martens (2.3%, 3/130), and canine parvovirus 2 (CPV-2) was found in one mink, one otter, and zero ermines (N = 27). Viruses were endemic and bottleneck events, founder effects, and genetic drift generated regional lineages. We identified two local closely related AMDV lineages, one CPV-2 lineage, and five FPV lineages. Highly similar viruses were identified in different hosts, demonstrating cross-species transmission. The likelihood for cross-species transmission differed among viruses and some species likely represented dead-end spillover hosts. We suggest that there are principal maintenance hosts (otters for FPV, raccoons for CPV-2/FPV, mink for AMDV) that enable viral persistence and serve as sources for other susceptible species. In this multi-host system, viral and host factors affect viral persistence and distribution, shaping parvoviral ecology and evolution, with implications for insular carnivore conservation.

Conservation Genomics in a Changing Arctic

Although logistically challenging to study, the Arctic is a bellwether for global change and is becoming a model for questions pertinent to the persistence of biodiversity. Disruption of Arctic ecosystems is accelerating, with impacts ranging from mixing of biotic communities to individual behavioral responses. Understanding these changes is crucial for conservation and sustainable economic development. Genomic approaches are providing transformative insights into biotic responses to environmental change, but have seen limited application in the Arctic due to a series of limitations. To meet the promise of genome analyses, we urge rigorous development of biorepositories from high latitudes to provide essential libraries to improve the conservation, monitoring, and management of Arctic ecosystems through genomic approaches.

Fluctuating fortunes: genomes and habitat reconstructions reveal global climate-mediated changes in bats' genetic diversity

Over the last approximately 2.6 Myr, Earth's climate has been dominated by cyclical ice ages that have profoundly affected species' population sizes, but the impact of impending anthropogenic climate change on species' extinction potential remains a worrying problem. We investigated 11 bat species from different taxonomic, ecological and geographical backgrounds using combined information from palaeoclimatic habitat reconstructions and genomes to analyse biotic impacts of historic climate change. We discover tightly correlated fluctuations between species' historic distribution and effective population size, identify frugivores as particularly susceptible to global warming, pinpoint large insectivores as having overall low effective population size and flag the onset of the Holocene (approx. 10-12 000 years ago) as the period with the generally lowest effective population sizes across the last approximately 1 Myr. Our study shows that combining genomic and palaeoclimatological approaches reveals effects of climatic shifts on genetic diversity and may help predict impacts of future climate change.

Living on the edge: Exploring the role of coastal refugia in the Alexander Archipelago of Alaska

Although islands are of long-standing interest to biologists, only a handful of studies have investigated the role of climatic history in shaping evolutionary diversification in high-latitude archipelagos. In this study of the Alexander Archipelago (AA) of Southeast Alaska, we address the impact of glacial cycles on geographic genetic structure for three mammals co-distributed along the North Pacific Coast. We examined variation in mitochondrial and nuclear loci for long-tailed voles (Microtus longicaudus), northwestern deermice (Peromyscus keeni), and dusky shrews (Sorex monticola), and then tested hypotheses derived from Species Distribution Models, reconstructions of paleoshorelines, and island area and isolation. In all three species, we identified paleoendemic clades that likely originated in coastal refugia, a finding consistent with other paleoendemic lineages identified in the region such as ermine. Although there is spatial concordance at the regional level for endemism, finer scale spatial and temporal patterns are less clearly defined. Demographic expansion across the region for these distinctive clades is also evident and highlights the dynamic history of Late Quaternary contraction and expansion that characterizes high-latitude species.