Multigene phylogeny of the Mustelidae: resolving relationships, tempo and biogeographic history of a mammalian adaptive radiation

Evolutionary characteristics of the mitochondrial NADH dehydrogenase subunit 6 gene in some populations of four sympatric Mustela species (Mustelidae, Mammalia) from central Europe

BACKGROUND

Selection on or reticulate evolution of mtDNA is documented in various mammalian taxa and could lead to misleading phylogenetic conclusions if not recognized. We sequenced the MT-ND6 gene of four sympatric Mustelid species of the genus Mustela from some central European populations. We hypothesised positive selection on MT-ND6, given its functional importance and the different body sizes and life histories of the species, even though climatic differences may be unimportant for adaptation in sympatry.

METHODS AND RESULTS

MT-ND6 genes were sequenced in 187 sympatric specimens of weasels, Mustela nivalis, stoats, M. erminea, polecats, M. putorius, and steppe polecats, M. eversmannii, from eastern Austria and of fourteen allopatric polecats from eastern-central Germany. Median joining networks, neighbour joining and maximum likelihood analyses as well as Bayesian inference grouped all species according to earlier published phylogenetic models. However, polecats and steppe polecats, two very closely related species, shared the same two haplotypes. We found only negative selection within the Mustela sequences, including 131 downloaded ones covering thirteen species. Positive selection was observed on three MT-ND6 codons of other mustelid genera retrieved from GenBank.

CONCLUSIONS

Negative selection for MT-ND6 within the genus Mustela suggests absence of both environmental and species-specific effects of cellular energy metabolism despite large species-specific differences in body size. The presently found shared polymorphism in European polecats and steppe polecats may result from ancestral polymorphism before speciation and historical or recent introgressive hybridization; it may indicate mtDNA capture of steppe polecats by M. putorius in Europe.

Diversity and evolution of the MHC class II DRB gene in the Capra sibirica experienced a demographic fluctuation in China

The major histocompatibility complex (MHC) genes are the most polymorphic genes in vertebrates, and their proteins play a critical role in adaptive immunity for defense against a variety of pathogens. MHC diversity was lost in many species after experiencing a decline in size. To understand the variation and evolution of MHC genes in the Siberian ibex, Capra sibirica, which has undergone a population decline, we analyzed the variation of the second exon of MHC class II DRB genes in samples collected from five geographic localities in Xinjiang, China, that belong to three diverged mitochondrial clades. Consequently, we identified a total of 26 putative functional alleles (PFAs) with 260 bp in length from 43 individuals, and found one (for 27 individuals) to three (for 5 individuals) PFAs per individual, indicating the presence of one or two DRB loci per haploid genome. The Casi-DRB1*16 was the most frequently occurring PFA, Casi-DRB1*22 was found in only seven individuals, 14 PFAs occurred once, 7 PFAs twice, implying high frequency of rare PFAs. Interestingly, more than half (15) of the PFAs were specific to clade I, only two and three PFAs were specific to clades II and III, respectively. So, we assume that the polygamy and sexual segregation nature of this species likely contributed to the allelic diversity of DRB genes. Genetic diversity indices showed that PFAs of clade II were lower in nucleotide, amino acid, and supertype diversity compared to those of the other two clades. The pattern of allele sharing and FST values between the three clades was to some extent in agreement with the pattern observed in mitochondrial DNA divergence. In addition, recombination analyses revealed no evidence for significant signatures of recombination events. Alleles shared by clades III and the other two clades diverged 6 million years ago, and systematic neighbor grids showed Trans-species polymorphism. Together with the PAML and MEME analyses, the results indicated that the DRB gene in C. sibirica evolved under balancing and positive selection. However, by comparison, it can be clearly seen that different populations were under different selective pressures. Our results are valuable in understanding the diversity and evolution of the DRB gene in a mountain living C. sibirica and in making decisions on future long-term protection strategies.

Analysis of the Complete Mitochondrial Genome of Pteronura brasiliensis and Lontra canadensis

P. brasiliensis and L. canadensis are two otter species, which successfully occupied semi-aquatic habitats and diverged from other Mustelidae. Herein, the full-length mitochondrial genome sequences were constructed for these two otter species for the first time. Comparative mitochondrial genome, selection pressure, and phylogenetic independent contrasts (PICs) analyses were conducted to determine the structure and evolutionary characteristics of their mitochondrial genomes. Phylogenetic analyses were also conducted to confirm these two otter species' phylogenetic position. The results demonstrated that the mitochondrial genome structure of P. brasiliensis and L. canadensis were consistent across Mustelidae. However, selection pressure analyses demonstrated that the evolutionary rates of mitochondrial genome protein-coding genes (PCGs) ND1, ND4, and ND4L were higher in otters than in terrestrial Mustelidae, whereas the evolutionary rates of ND2, ND6, and COX1 were lower in otters. Additionally, PIC analysis demonstrated that the evolutionary rates of ND2, ND4, and ND4L markedly correlated with a niche type. Phylogenetic analysis showed that P. brasiliensis is situated at the base of the evolutionary tree of otters, and then L. canadensis diverged from it. This study suggests a divergent evolutionary pattern of Mustelidae mitochondrial genome PCGs, prompting the otters' adaptation to semi-aquatic habitats.

First Identification and Molecular Characterization of Trichinella britovi (Nematoda: Trichinellidae) from the Pine Marten (Martes martes Linnaeus, 1758) in Romania

Trichinella spp. are etiological zoonotic agents spread throughout the world and affect mammals, birds, and reptiles; they evolve via two cycles: domestic and sylvatic. Martes martes is a carnivorous nocturnal mammal from the family Mustelidae. In Romania, this host is widespread in all forests of the country. Martes martes has an extremely voracious appetite, feeding on fruit and also on a variety of small animals, including rodents such as mice and rats. The aim of this study was the identification and molecular characterization of Trichinella larvae isolated from the muscle tissue of Martes martes collected in different counties of Romania. The muscle samples were examined via artificial digestion, and the larvae were identified at the species level via multiplex PCR. The presence of larvae belonging to Trichinella britovi, a species frequently identified in wild carnivores in temperate zones, was observed. Although T. britovi has been already reported in several host species in Romania, this is the first time this species has been observed in a Martes martes specimen. This finding contributes to our knowledge about the host species involved in the maintenance of the Trichinella sylvatic cycle in Romania, and it confirms that this parasite is consistently present in the wild fauna of this country.

Chromosome-length genome assembly and karyotype of the endangered black-footed ferret (Mustela nigripes)

The black-footed ferret (Mustela nigripes) narrowly avoided extinction to become an oft-cited example of the benefits of intensive management, research, and collaboration to save a species through ex situ conservation breeding and reintroduction into its former range. However, the species remains at risk due to possible inbreeding, disease susceptibility, and multiple fertility challenges. Here, we report the de novo genome assembly of a male black-footed ferret generated through a combination of linked-read sequencing, optical mapping, and Hi-C proximity ligation. In addition, we report the karyotype for this species, which was used to anchor and assign chromosome numbers to the chromosome-length scaffolds. The draft assembly was ~2.5 Gb in length, with 95.6% of it anchored to 19 chromosome-length scaffolds, corresponding to the 2n = 38 chromosomes revealed by the karyotype. The assembly has contig and scaffold N50 values of 148.8 kbp and 145.4 Mbp, respectively, and is up to 96% complete based on BUSCO analyses. Annotation of the assembly, including evidence from RNA-seq data, identified 21,406 protein-coding genes and a repeat content of 37.35%. Phylogenomic analyses indicated that the black-footed ferret diverged from the European polecat/domestic ferret lineage 1.6 million yr ago. This assembly will enable research on the conservation genomics of black-footed ferrets and thereby aid in the further restoration of this endangered species.

How was genetic diversity transferred with translocations from ex situ to in situ? A case study of the European mink translocation to Hiiumaa Island in Estonia

The European mink (Mustela lutreola) is one of the most threatened small carnivores, listed as Critically Endangered on the IUCN Red List. Tallinn Zoo started a conservation breeding operation in 1980, which in 1992 was shaped into the European mink EEP Program to maintain a demographically and genetically healthy population in captivity. Since 2000, mink have been translocated on a yearly basis from the breeding facility in the zoo to Hiiumaa Island (Estonia) until the formation of the wild island population in 2016. Maintaining a healthy genetic structure in a captive population was a priority, so genetically the least valuable animals, according to calculations made by a population management program, were used for translocation. This study aims to assess the amount of genetic diversity passed from an ex situ population on to the island population. Comparisons of the genetic diversity were made by mitochondrial, microsatellite and nuclear markers. In addition, our results were combined with the pedigree data from the European mink EEP Studbook to further evaluate the flow of genetic diversity from the founder population to the established wild population. According to the findings, the island population's allelic richness was comparable to that of the founder population, and no evidence was found that its genetic structure had diverged from that of the founder population. It seems that the formation of the island population has been a gradual process of no more than the last seven yearly translocations.

Functional and phylogenetic interpretation of the forelimb myology of two South American carnivorans, the ring-tailed coati (Nasua nasua) and crab-eating raccoon (Procyon cancrivorus)

A comparative analysis of the forelimb myology of two neotropical procyonids (Nasua nasua and Procyon cancrivorus) was performed to assess how observed differences in their myological configuration would be related to their diverse ecological behaviors and phylogeny. Although both species are associated with the arboreal substrate, N. nasua is a more agile climber that usually digs; whereas P. cancrivorus spends most of its time on the ground foraging, climbing on the trees as shelter and is a good swimmer. Here, myological descriptions, muscle maps, phylogenetic optimizations, and muscle mass data of the forelimb of these two procyonids are presented. The main functional muscular groups are discussed in a comparative framework with other carnivorans that present a wide ecological diversity. Also, muscular characters were mapped onto a phylogeny to explore their evolution and to obtain ancestral state reconstructions. Results indicate clear myological differences among the two neotropical procyonids associated with their ecological preferences. One of the most remarkable anatomical differences is the arrangement and relative mass of the extrinsic musculature, mainly the musculus rhomboideus and the delto-pectoral complexes. In Nasua nasua, these suggested a greater stability in their shoulder girdle for climbing and digging and probably would provide stronger neck and head movements when they use them for foraging on the ground. Conversely, P. cancrivorus has a different extrinsic muscular configuration, which would allow an increment on the stride length and faster movements of the forelimb associated with more frequent terrestrial gaits. Also, significant differences are observed in the distal musculature, associated with strong movements of forepaws when climbing and digging in N. nasua; whereas, P. cancrivorus configuration suggested precise forearm and digits movements, related to manipulation of food items when they are catching prey or feeding. Most of the codified features of P. cancrivorus would reflect retention of plesiomorphies acquired in the common ancestor of caniforms or arctoids, whereas N. nasua shows derived traits, particularly in the proximal forelimb region. The present work increases the information available on the myology of these particular taxa and extant generalized arctoid models in general. The analyses presented here will be useful both for other comparative myological studies (morpho-functional and phylogenetic) and for muscular reconstruction in extinct procyonids, as well as other carnivorans.

Molecular Identification of Sarcocystis rileyi and Sarcocystis sp. (Closely Related to Sarcocystis wenzeli) in Intestines of Mustelids from Lithuania

The genus Sarcocystis is a group of numerous protozoan parasites having a two-host life cycle. Based on laboratory experiments and/or phylogenetic analysis results it was shown that seven Sarcocystis spp. producing sarcocsyts in bird tissues are transmitted via predatory placental mammals. To date the role of small mammals of the family Mustelidae in the distribution of avian Sarcocystis spp. have not been studied. During the current investigation, intestinal mucosa scrapings of 115 mustelids belonging to five species were tested for S. albifronsi, S. anasi, S. rileyi, and S. wenzeli infecting anseriforms and chickens. Microscopically, free sporocysts, sporulating oocysts, and loose oocysts were found in 61 samples (53.0%). Using nested PCR targeting the ITS1 region and sequencing, S. rileyi was confirmed in eight American minks, two European polecats and single European badger. Sarcocystis sp. was identified in one American mink and one European pine marten. Based on the partial ITS1 region this parasite showed that 100% identity to pathogenic Sarcocystis sp. caused a fatal infection in backyard chickens from Brazil. Phylogenetically, the Sarcocystis sp. identified in our study was most closely related to S. wenzeli parasitising domestic fowl (Gallus domesticus).

A chromosome-level genome assembly reveals genomic characteristics of the American mink (Neogale vison)

Availability of a contiguous chromosome-level genome assembly is the foundational step to develop genome-based studies in American mink (Neogale vison). The main objective of this study was to provide a high quality chromosome-level genome assembly for American mink. An initial draft of the genome assembly was generated using 2,884,047 PacBio long reads. Integration of Hi-C data into the initial draft led to an assembly with 183 scaffolds and scaffold N50 of 220 Mb. This gap-free genome assembly of American mink (ASM_NN_V1) had a length of 2.68 Gb in which about 98.6% of the whole genome was covered by 15 chromosomes. In total, 25,377 genes were predicted across the American mink genome using the NCBI Eukaryotic Genome Annotation Pipeline. In addition, gene orthology, demographic history, synteny blocks, and phylogenetic relationships were studied in connection with the genomes of other related Carnivora. Furthermore, population-based statistics of 100 sequenced mink were presented using the newly assembled genome. Remarkable improvements were observed in genome contiguity, the number of scaffolds, and annotation compared to the first draft of mink genome assembly (NNQGG.v01). This high-quality genome assembly will support the development of efficient breeding strategies as well as conservation programs for American mink.

Molecular Dietary Analysis of Three Sympatric Mustelidae in Northeast China

Diet analysis is essential to fully understand the biology of a species and its function within the ecosystem, as well as being key in identifying food web interactions and the population dynamics of predators and prey. The understanding of the diet of small to mid-sized carnivores remains generally lacking or uninformative due to the inability for taxonomic resolution based on morphology. The yellow-throated marten (Martes flavigula), Eurasian river otter (Lutra lutra), and Siberian weasel (Mustela sibirica) are three important Mustelidae species in ecosystems of northeast China. Based on fecal DNA and a next-generation sequencing (NGS) approach, we analyzed the vertebrate prey of these three sympatric Mustelidae. Prey included 7 mammalian taxa, 10 fishes, 2 birds, and 2 amphibians, with 85% of the taxa assigned to the species level. In total, twenty-one vertebrate prey taxa were identified from seven yellow-throated martens, eight Eurasian river otters, and two Siberian weasels. Concerning identified dietary species, 10 taxa were consumed by yellow-throated martens, 14 by Eurasian river otters, and 4 by Siberian weasels. Some prey species were identified in more than one species. Amphibians and fishes were the most dominant Eurasian river otter prey categories, whereas Eurasian badger (Meles leucurus), birds, and rodents were the main yellow-throated marten prey; amphibians and rodents were largely contained in Siberian weasel prey. Among prey items, Dybowski's frog (Rana dybowskii) and Korean field mouse (Apodemus peninsulae) were identified in all three Mustelidae species but our analyses suggest potential diet preferences among Mustelidae species. Future studies should focus on understanding the trophic relationships of these three Mustelidae species, providing valuable information for their conservation planning.

Role of feeding specialization in taste receptor loss: insights from sweet and umami receptor evolution in Carnivora

Controversy and misunderstanding surround the role of feeding specialization in taste receptor loss in vertebrates. We refined and tested the hypothesis that this loss is caused by feeding specializations. Specifically, feeding specializations were proposed to trigger time-dependent process of taste receptor loss through deprivation of benefit of using the receptor's gustatory function. We propose that this process may be accelerated by abiotic environmental conditions or decelerated/stopped because of extragustatory functions of the receptor's protein(s). As test case we used evolution of the sweet (TAS1R2+TAS1R3) and umami (TAS1R1+TAS1R3) receptors in Carnivora (dogs, cats, and kin). We predicted these receptors' absence/presence using data on presence/absence of inactivating mutations in these receptors' genes and data from behavioral sweet/umami preference tests. We identified 20 evolutionary events of sweet (11) or umami (9) receptor loss. These events affected species with feeding specializations predicted to favor sweet/umami receptor loss (27 and 22 species, respectively). All species with feeding habits predicted to favor sweet/umami receptor retention (11 and 24, respectively) were found to retain that receptor. Six species retained the sweet (5) or umami (1) receptor despite feeding specialization predicted to favor loss of that receptor, which can be explained by the time dependence of sweet/umami receptor loss process and the possible decelerating effect of TAS1R extragustatory functions so that the sweet/umami receptor process is ongoing in these species. Our findings support the idea that feeding specialization leads to taste receptor loss and is the main if not only triggering factor for evolutionary loss of taste receptors.

The impact of environmental factors on the evolution of brain size in carnivorans

The reasons why some animals have developed larger brains has long been a subject of debate. Yet, it remains unclear which selective pressures may favour the encephalization and how it may act during evolution at different taxonomic scales. Here we studied the patterns and tempo of brain evolution within the order Carnivora and present large-scale comparative analysis of the effect of ecological, environmental, social, and physiological variables on relative brain size in a sample of 174 extant carnivoran species. We found a complex pattern of brain size change between carnivoran families with differences in both the rate and diversity of encephalization. Our findings suggest that during carnivorans’ evolution, a trade-off have occurred between the cognitive advantages of acquiring a relatively large brain allowing to adapt to specific environments, and the metabolic costs of the brain which may constitute a disadvantage when facing the need to colonize new environments.

The brain size of carnivores has evolved to balance a trade-off between increased cognitive function and increased metabolic cost.

The impact of environmental factors on the evolution of brain size in carnivorans

The reasons why some animals have developed larger brains has long been a subject of debate. Yet, it remains unclear which selective pressures may favour the encephalization and how it may act during evolution at different taxonomic scales. Here we studied the patterns and tempo of brain evolution within the order Carnivora and present large-scale comparative analysis of the effect of ecological, environmental, social, and physiological variables on relative brain size in a sample of 174 extant carnivoran species. We found a complex pattern of brain size change between carnivoran families with differences in both the rate and diversity of encephalization. Our findings suggest that during carnivorans' evolution, a trade-off have occurred between the cognitive advantages of acquiring a relatively large brain allowing to adapt to specific environments, and the metabolic costs of the brain which may constitute a disadvantage when facing the need to colonize new environments.

Comparative Genome Analysis Reveals the Genomic Basis of Semi-Aquatic Adaptation in American Mink (Neovison vison)

Simple Summary

In nature, every animal lives in its suitable environment to support life and reproduction. Due to survival pressure, the ancient ancestors of land animals changed from living in water to living on land through a long time of evolution. But as more and more animals live on land, the pressure to survive increases, and some animals continue to evolve and re-enter the water. Species evolved from water to land, fins became limbs, and re-entering organisms evolved webbing, all in an effort to better adapt to their environment. The American mink’s life is extremely dependent on the water environment, but its external changes are less pronounced than those of other water-dependent animals. Since the limited external changes are not sufficient to indicate that the American mink is a semi-aquatic mammal, we can explore the evidence of aquatic adaptation at the molecular level. Through comparative genomic analysis, we obtained that American mink has certain adaptive evolution to aquatic environment in olfactory, coagulation, immunity and other aspects. The results of this study have important reference significance for exploring biological evolution and species conservation.

Abstract

Although the American mink is extremely dependent on water and has evolved a range of aquatic characteristics, its structural adaptation to water is still less obvious than that of other typical semi-aquatic mammals, such as otters. Therefore, many scholars consider it not to be a semi-aquatic mammal. In order to make the point that minks are semi-aquatic mammals more convincing, we provide evidence at the micro (genome)-level. In particular, we used the genomes of the American mink and 13 mammalian species to reconstruct their evolutionary history, identified genes that affect aquatic adaptation, and examined the evolution of aquatic adaptation. By analyzing unique gene families, the expansion and contraction of gene families, and positive selection genes, we found that the American mink genome has evolved specifically for aquatic adaptation. In particular, we found that the main adaptive characteristics of the American mink include the external structural characteristics of bone and hair development, as well as the internal physiological characteristics of immunity, olfaction, coagulation, lipid metabolism, energy metabolism, and nitrogen metabolism. We also observed that the genomic characteristics of the American mink are similar to those of other aquatic and semi-aquatic mammals. This not only provides solid genomic evidence for the idea that minks are semi-aquatic mammals, but also leads to a clearer understanding of semi-aquatic species. At the same time, this study also provides a reference for the protection and utilization of the American mink.

Extensive genome introgression between domestic ferret and European polecat during population recovery in Great Britain

The European polecat (Mustela putorius) is a mammalian predator which occurs across much of Europe east to the Ural Mountains. In Great Britain, following years of persecution the range of the European polecat contracted and by the early 1900s was restricted to unmanaged forests of central Wales. The European polecat has recently undergone a population increase due to legal protection and its range now overlaps that of feral domestic ferrets (Mustela putorius furo). During this range expansion, European polecats hybridized with feral domestic ferrets producing viable offspring. Here, we carry out population-level whole-genome sequencing on 8 domestic ferrets, 19 British European polecats, and 15 European polecats from the European mainland. We used a range of population genomics methods to examine the data, including phylogenetics, phylogenetic graphs, model-based clustering, phylogenetic invariants, ABBA-BABA tests, topology weighting, and Fst. We found high degrees of genome introgression in British polecats outside their previous stronghold, even in those individuals phenotyped as “pure” polecats. These polecats ranged from presumed F1 hybrids (gamma = 0.53) to individuals that were much less introgressed (gamma = 0.2). We quantify this introgression and find introgressed genes containing Fst outliers associated with cognitive function and sight.

The exceptionally high diversity of small carnivorans from the Late Miocene hominid locality of Hammerschmiede (Bavaria, Germany)

The present study deals with new material of carnivorans (Mustelidae, Mephitidae, Ailuridae, Potamotheriinae and Viverridae) from the basal Tortonian (Late Miocene, late Astaracian) hominid-bearing locality of Hammerschmiede (Bavaria, Germany). The small carnivoran fauna includes 20 species belonging to nine different subfamilies (Guloninae, Lutrinae, Mellivorinae, Potamotheriinae, Leptarctinae, Mephitinae, Simocyoninae, Genettinae and Viverrinae). The identified forms include: “Martes” sansaniensis, “Martes” cf. munki, “Martes” sp., Circamustela hartmanni n. sp., Laphyctis mustelinus, Guloninae indet., Eomellivora moralesi, Vishnuonyx neptuni, Paralutra jaegeri, Lartetictis cf. dubia, Trocharion albanense, Palaeomeles pachecoi, Proputorius sansaniensis, Proputorius pusillus, Alopecocyon goeriachensis, Simocyoninae indet., Potamotherium sp., Semigenetta sansaniensis, Semigenetta grandis and Viverrictis modica. The new species Circamustela hartmanni n. sp. is differentiated from the other members of the genus by its small size and the morphology of its dental cusps in the upper and lower carnassials. This is one of the highest reported taxonomic diversities for fossil small carnivorans in the Miocene of Europe, including also first and last occurrences for several genera and species. Additionally, the assemblage comprises some rare taxa such as Palaeomeles pachecoi and Eomellivora moralesi. An ecomorphological comparison of the discovered taxa reveals possible cases of competition and niche partitioning.

Phylogenomics of the world's otters

Comparative whole-genome analyses hold great power to illuminate commonalities and differences in the evolution of related species that share similar ecologies. The mustelid subfamily Lutrinae includes 13 currently recognized extant species of otters,^1-5 a semiaquatic group whose evolutionary history is incompletely understood. We assembled a dataset comprising 24 genomes from all living otter species, 14 of which were newly sequenced. We used this dataset to infer phylogenetic relationships and divergence times, to characterize patterns of genome-wide genealogical discordance, and to investigate demographic history and current genomic diversity. We found that genera Lutra, Aonyx, Amblonyx, and Lutrogale form a coherent clade that should be synonymized under Lutra, simplifying the taxonomic structure of the subfamily. The poorly known tropical African Aonyx congicus and the more widespread Aonyx capensis were found to be reciprocally monophyletic (having diverged 440,000 years ago), supporting the validity of the former as a distinct species. We observed variable changes in effective population sizes over time among otters within and among continents, although several species showed similar trends of expansions and declines during the last 100,000 years. This has led to different levels of genomic diversity assessed by overall heterozygosity, genome-wide SNV density, and run of homozygosity burden. Interestingly, there were cases in which diversity metrics were consistent with the current threat status (mostly based on census size), highlighting the potential of genomic data for conservation assessment. Overall, our results shed light on otter evolutionary history and provide a framework for further in-depth comparative genomic studies targeting this group.

Characterizing the oral and distal gut microbiota of the threatened southern sea otter (Enhydra lutris nereis) to enhance conservation practice

The southern sea otter (Enhydra lutris nereis) is a threatened sub-species in coastal ecosystems. To understand better the role of diet, monitor health, and enhance management of this and other marine mammal species, we characterized the oral (gingival) and distal gut (rectal and fecal) microbiota of 158 wild southern sea otters living off the coast of central California, USA, and 12 captive sea otters, some of which were included in a diet shift experiment. We found that the sea otter fecal microbiota was distinct from that of three other otter species, and that captivity does not significantly alter the community structure of the sea otter gingival or distal gut microbiota. Metagenomic analysis unexpectedly revealed that the majority of sea otter fecal DNA is derived from prey, rather than from indigenous bacteria or host cells as with most other mammals. We speculate that a reduced bacterial biomass in the sea otter gut reflects rapid gut transit time and a particular strategy for foraging and energy harvest. This study establishes a reference for the healthy sea otter microbiota, highlights how a marine lifestyle may shape the mammalian microbiota, and may inform future health assessments and conservation management of sea otter populations.

Multiple types of genomic variation contribute to adaptive traits in the mustelid subfamily Guloninae

Species of the mustelid subfamily Guloninae inhabit diverse habitats on multiple continents, and occupy a variety of ecological niches. They differ in feeding ecologies, reproductive strategies and morphological adaptations. To identify candidate loci associated with adaptations to their respective environments, we generated a de novo assembly of the tayra (Eira barbara), the earliest diverging species in the subfamily, and compared this with the genomes available for the wolverine (Gulo gulo) and the sable (Martes zibellina). Our comparative genomic analyses included searching for signs of positive selection, examining changes in gene family sizes, as well as searching for species-specific structural variants (SVs). Among candidate loci associated with phenotypic traits, we observed many related to diet, body condition and reproduction. For example, for the tayra, which has an atypical gulonine reproductive strategy of aseasonal breeding, we observe species-specific changes in many pregnancy-related genes. For the wolverine, a circumpolar hypercarnivore that must cope with seasonal food scarcity, we observed many changes in genes associated with diet and body condition. All types of genomic variation examined (single nucleotide polymorphisms, gene family expansions, structural variants) contributed substantially to the identification of candidate loci. This strongly argues for consideration of variation other than single nucleotide polymorphisms in comparative genomics studies aiming to identify loci of adaptive significance.

Cranial variation in the Altai weasel Mustela altaica (Carnivora: Mustelidae) and its possible taxonomic implications

The Altai weasel, Mustela altaicaPallas 1811 (Carnivora: Mustelidae), is widely distributed across open landscapes of central and eastern Asia. The geographical variation in skull morphology and taxonomic composition of this small mustelid remain poorly known. Based on extensive sampling in museum collections, we analyzed morphometric variation in 23 cranial characters of 232 M. altaica specimens from across its entire range. The multivariate analysis (nonmetric multidimensional scaling) revealed a geographical pattern of variation in cranial size and shape. Based on this result, two morphological groups have been comprehensively characterized. Weasels of the first group differ from those of the second in having larger skulls, wider neurocrania, elongated toothrows, and larger carnassial teeth, as well as in sexual size dimorphism. The first group is confined to the southwestern part of the species range (Tien Shan Mountains, Pamir, the Himalayas, and southern Tibet). The second group is widespread from the foothills of Tien Shan in the west to Khingan and Manchuria in the east. The distribution of the two groups of M. altaica is consistent with the “Out-of-Tibet” hypothesis. The first group presumably can be attributed to M. a. temonHodgson, 1857, whereas the second group can be treated as nominotypical M. a. altaicaPallas, 1811.

The locomotion of extinct secondarily aquatic tetrapods

The colonisation of freshwater and marine ecosystems by land vertebrates has repeatedly occurred in amphibians, reptiles, birds and mammals over the course of 300 million years. Functional interpretations of the fossil record are crucial to understanding the forces shaping these evolutionary transitions. Secondarily aquatic tetrapods have acquired a suite of anatomical, physiological and behavioural adaptations to locomotion in water. However, much of this information is lost for extinct clades, with fossil evidence often restricted to osteological data and a few extraordinary specimens with soft tissue preservation. Traditionally, functional morphology in fossil secondarily aquatic tetrapods was investigated through comparative anatomy and correlation with living functional analogues. However, in the last two decades, biomechanics in palaeobiology has experienced a remarkable methodological shift. Anatomy-based approaches are increasingly rigorous, informed by quantitative techniques for analysing shape. Moreover, the incorporation of physics-based methods has enabled objective tests of functional hypotheses, revealing the importance of hydrodynamic forces as drivers of evolutionary innovation and adaptation. Here, we present an overview of the latest research on the locomotion of extinct secondarily aquatic tetrapods, with a focus on amniotes, highlighting the state-of-the-art experimental approaches used in this field. We discuss the suitability of these techniques for exploring different aspects of locomotory adaptation, analysing their advantages and limitations and laying out recommendations for their application, with the aim to inform future experimental strategies. Furthermore, we outline some unexplored research avenues that have been successfully deployed in other areas of palaeobiomechanical research, such as the use of dynamic models in feeding mechanics and terrestrial locomotion, thus providing a new methodological synthesis for the field of locomotory biomechanics in extinct secondarily aquatic vertebrates. Advances in imaging technology and three-dimensional modelling software, new developments in robotics, and increased availability and awareness of numerical methods like computational fluid dynamics make this an exciting time for analysing form and function in ancient vertebrates.

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.

Museomics dissects the genetic basis for adaptive seasonal colouration in the least weasel

Dissecting the link between genetic variation and adaptive phenotypes provides outstanding opportunities to understand fundamental evolutionary processes. Here, we use a museomics approach to investigate the genetic basis and evolution of winter coat colouration morphs in least weasels (Mustela nivalis), a repeated adaptation for camouflage in mammals with seasonal pelage colour moults across regions with varying winter snow. Whole-genome sequence data was obtained from biological collections and mapped onto a newly assembled reference genome for the species. Sampling represented two replicate transition zones between nivalis and vulgaris colouration morphs in Europe, which typically develop white or brown winter coats, respectively. Population analyses showed that the morph distribution across transition zones is not a by-product of historical structure. Association scans linked a 200 kb genomic region to colouration morph, which was validated by genotyping museum specimens from inter-morph experimental crosses. Genotyping the wild populations narrowed down the association to pigmentation gene MC1R and pinpointed a candidate amino acid change co-segregating with colouration morph. This polymorphism replaces an ancestral leucine residue by lysine at the start of the first extracellular loop of the protein in the vulgaris morph. A selective sweep signature overlapped the association region in vulgaris, suggesting that past adaptation favoured winter-brown morphs and can anchor future adaptive responses to decreasing winter snow. Using biological collections as valuable resources to study natural adaptations, our study showed a new evolutionary route generating winter colour variation in mammals and that seasonal camouflage can be modulated by changes at single key genes.

Complete mitochondrial genome sequencing of Lutra lutra (Linnaeus, 1758) (Carnivora: Mustelidae) and its phylogenetic status in Mustelidae

We report the complete mitochondrial genome sequence of the endangered Eurasian otter, Lutra lutra. The complete mitochondrial genome is 16,537 bp in length and contains 13 protein-coding genes, 22 transfer RNA, two ribosomal RNA, and one control region. The mitogenome is A + T rich, with a composition of 32.2% A, 27.5% C, 14.5% G, and 25.8% T. Phylogenetic analysis based on 13 protein-coding mitochondrial genes of Mustelidae supports the conventional systematic treatment with eight subfamilies. Lutra lutra is closely related to Lutra sumatrana and the subfamily Lutrinae was closely grouped with the Ictonychinae. This study provides genetic and taxonomic information for future studies of Eurasian otters and the Mustelidae.

Airborne vocal communication in adult neotropical otters (Lontra longicaudis)

Most aquatic mammals have complex social and communication systems. Interestingly, little is known about otters' vocal communication compared to other aquatic mammals. Here, for the first time, we acoustically describe vocalizations of the neotropical otter (Lontra longicaudis), a solitary and endangered New World otter species. We recorded vocalizations and behavioral contexts from six captive neotropical otters at Projeto Lontra, Santa Catarina Island, Brazil. Analysis of acoustic parameters were used to classify the vocalizations according to structure and context. We describe six call types with highly tonal as well as chaotic vocalizations with fundamental frequencies ranging from 90 to 2500 Hz. Additionally, we identified sex differences in the usage of calls. Results suggest that the neotropical river otter has a rich vocal repertoire, similar in complexity to other solitary otter species, but less complex than that of the social giant otter. Despite differences in sociality, phylogeny and ecology, L. longicaudis seems to possess vocalizations homologous to those found in other otters (e.g. hah and chirp), suggesting phylogenetic inertia in otter communicative repertoire. Otters thus offer an interesting but neglected group to explore the evolution of communication systems.

Pliocene-Early Pleistocene Geological Events Structure Pacific Martens (Martes caurina)

The complex topography, climate, and geological history of Western North America have shaped contemporary patterns of biodiversity and species distributions in the region. Pacific martens (Martes caurina) are distributed along the northern Pacific Coast of North America with disjunct populations found throughout the Northwestern Forested Mountains and Marine West Coast Forest ecoregions of the West Coast. Martes in this region have been classified into subspecies; however, the subspecific designation has been extensively debated. In this study, we use genomic data to delineate conservation units of Pacific marten in the Sierra-Cascade-Coastal montane belt in the western United States. We analyzed the mitochondrial genome for 94 individuals to evaluate the spatial distribution and divergence times of major lineages. We further genotyped 401 individuals at 13 microsatellite loci to investigate major patterns of population structure. Both nuclear and mitochondrial DNA suggest substantial genetic substructure concordant with historical subspecies designations. Our results revealed that the region contains 2 distinct mitochondrial lineages: a Cascades/Sierra lineage that diverged from the Cascades/coastal lineage 2.23 (1.48-3.14 mya), consistent with orogeny of the Cascade Mountain chain. Interestingly, Pacific Martes share phylogeographic patterns similar with other sympatric taxa, suggesting that the complex geological history has shaped the biota of this region. The information is critical for conservation and management efforts, and further investigation of adaptive diversity is warranted following appropriate revision of conservation management designations.

Social organization of otters in relation to their ecology

Otter species are known to fluctuate intraspecifically from a solitary lifestyle to group-living arrangements. By examining what is known about habitat use and foraging style in otters of 13 different species, based on 93 studied sites, we assessed (1) the relationship between social habits and preferred habitats, (2) the relationship between species and prey preferences, and (3) the effect of predator avoidance on their social organization in order to assess the socio-ecological factors influencing otters. Females remain the core of their social stability. We show the major influence of habitats and feeding strategies (i.e. socio-ecology) of otters. The different species of solitary otters most often inhabit linear environments, such as freshwater ecosystems or wave-exposed marine coasts, and their habitat is often subject to disturbances that fragment their functional continuity. Social otters are more often found in extensive habitats with high plant cover, regular food resources and in areas with large predators compared to solitary species. The maintenance of regular resources and the fact that the main trophic resources are replenished rapidly might be determining factors driving sociality. Group-living and bachelor congregations among otters can also respond to pressure from large predators. This suggests that foraging, habitat use and the presence of large predators may be the drivers of sociality in otters. We conclude that most otters have a greater social potential than previously assumed, which is confirmed by their various vocalizations recently described.

First mitochondrial genome of the marbled polecat Vormela peregusna (Carnivora, Mustelidae)

The marbled polecat, Vormela peregusna, is one of the least studied species in the Mustelidae family, especially with regard to phylogeography and genetic diversity. In this study, we determined the mitochondrial genome sequence of V. peregusna and investigated its position within the Mustelidae phylogeny. The generated mitogenome is 15,982 bp in length; it consists of 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and a control region.

Differing limb functions and their potential influence upon the diversification of the mustelid hindlimb skeleton

Though form-function relationships of the mammalian locomotor system have been investigated for over a century, recent models of trait evolution have hitherto been seldom used to identify likely evolutionary processes underlying the locomotor system’s morphological diversity. Using mustelids, an ecologically diverse carnivoran lineage, I investigated whether variation in hindlimb skeletal morphology functionally coincides with climbing, digging, swimming and generalized locomotor habits by using 15 linear traits of the femur, tibia, fibula, calcaneum and metatarsal III across 44 species in a principal component analysis. I subsequently fit different models of Brownian motion and adaptive trait diversification individually to each trait. Climbing, digging and swimming mustelids occupy distinct regions of phenotypic space characterized by differences in bone robustness. Models of adaptive and neutral evolution are, respectively, the best fits for long bone lengths and muscle in-levers, suggesting that different kinds of traits may be associated with different evolutionary processes. However, simulations based upon models of best fit reveal low statistical power to rank the models. Though differences in mustelid hindlimb skeletal morphology appear to coincide with locomotor habits, further study, with sampling expanded beyond the Mustelidae, is necessary to better understand to what degree adaptive evolution shapes morphological diversity of the locomotor system.

Fifty Years of Research on European Mink Mustela lutreola L., 1761 Genetics: Where Are We Now in Studies on One of the Most Endangered Mammals?

The purpose of this review is to present the current state of knowledge about the genetics of European mink Mustela lutreola L., 1761, which is one of the most endangered mammalian species in the world. This article provides a comprehensive description of the studies undertaken over the last 50 years in terms of cytogenetics, molecular genetics, genomics (including mitogenomics), population genetics of wild populations and captive stocks, phylogenetics, phylogeography, and applied genetics (including identification by genetic methods, molecular ecology, and conservation genetics). An extensive and up-to-date review and critical analysis of the available specialist literature on the topic is provided, with special reference to conservation genetics. Unresolved issues are also described, such as the standard karyotype, systematic position, and whole-genome sequencing, and hotly debated issues are addressed, like the origin of the Southwestern population of the European mink and management approaches of the most distinct populations of the species. Finally, the most urgent directions of future research, based on the research questions arising from completed studies and the implementation of conservation measures to save and restore M. lutreola populations, are outlined. The importance of the popularization of research topics related to European mink genetics among scientists is highlighted.

Are subspecies (of Eira barbara) real?

The subspecies concept is one of the most controversial in Linnean taxonomy. In the past, subspecies were described without a clear conceptual framework, triggering confusion and motivating criticism of the very concept of a subspecies. At present, subspecies are conceived as aggregates of populations that are geographically isolated, are composed of interfertile individuals, and are morphologically diagnosable. The tayra, Eira barbara, was described in 1758 and has had a stable taxonomic history at the species level. However, below the species level, 16 subspecies have been named, with from two to seven subspecies recognized as valid by different authors. None of the subspecies were, however, described within a clear conceptual framework. Using the modern concept of a subspecies, I analyzed subspecies of E. barbara recognized by recent authors. I gathered morphometric data from 155 specimens in mammal collections, georeferenced each specimen, and recorded membership to subspecies assigned by different references and by its location. I gathered climate and geographic data for each location. I analyzed data using Principal Components Analysis (PCA) and analysis of variance (ANOVA). Specimens exhibited sexual dimorphism in size but not in skull shape. I used regression analysis to test for associations between skull shape and size and climate data. Geographic analyses documented that subspecies are not allopatric, violating one of the main properties of the subspecies concept. ANOVA showed significant differences in skull morphology between some pairs of recognized subspecies but not others. However, none of the subspecies segregated in the PCA. Thus, the recognized subspecies could not be diagnosed from morphological data, violating another property of the subspecies concept. Size varied greatly between the sexes using different schemes for recognized subspecies. Climate variables explained between 4% and 6% of size variation for males and females. Skull shape proved not to be geographically variable.

Adaptations for amphibious vision in sea otters (Enhydra lutris): structural and functional observations

Sea otters (Enhydra lutris) are amphibious mammals that maintain equal in-air and underwater visual acuity. However, their lens-based underwater accommodative mechanism presumably requires a small pupil that may limit sensitivity across light levels. In this study, we consider adaptations for amphibious living by assessing the tapetum lucidum, retina, and pupil dynamics in sea otters. The sea otter tapetum lucidum resembles that of terrestrial carnivores in thickness and fundic coverage. A heavily rod-dominated retina appears qualitatively similar to the ferret and domestic cat, and a thick outer nuclear layer relative to a thinner inner nuclear layer is consistent with nocturnal vertebrates and other amphibious carnivores. Pupil size range in two living sea otters is smaller relative to other amphibious marine carnivores (pinnipeds) when accounting for test conditions. The pupillary light response seems slower than other aquatic and terrestrial species tested in comparable brightness, although direct comparisons require further assessment. Our results suggest that sea otters have retained features for low-light vision but rapid adjustments and acute underwater vision may be constrained across varying light levels by a combination of pupil shape, absolute eye size, and the presumed coupling between anterior lens curvature and pupil size during accommodation.

New insights into the giant mustelids (Mammalia, Carnivora, Mustelidae) from Langebaanweg fossil site (West Coast Fossil Park, South Africa, early Pliocene)

Giant mustelids are a paraphyletic group of mustelids found in the Neogene of Eurasia, Africa and North America. Most are known largely from dental remains, with their postcranial skeleton mostly unknown. Here, we describe new craniodental and postcranial remains of the large lutrine Sivaonyx hendeyi and the leopard-size gulonine Plesiogulo aff. monspessulanus from the early Pliocene site Langebaanweg, South Africa. The new material of the endemic S. hendeyi, includes upper incisors and premolars, and fragmentary humerus, ulna and a complete astragalus. Its postcrania shares more traits with the living Aonyx capensis than the late Miocene Sivaonyx beyi from Chad. Sivaonyx hendeyi could therefore be tentatively interpreted as a relatively more aquatic taxon than the Chadian species, comparable to A. capensis. The new specimens of Plesiogulo comprise two edentulous maxillae, including one of a juvenile individual with incomplete decidual dentition, and a fragmentary forelimb of an adult individual. The new dental measurements point to this form being amongst the largest specimens of the genus. Both P3-4 differs from the very large species Plesiogulo botori from late Miocene of Kenya and Ethiopia. This confirms the existence of two distinct large species of Plesiogulo in Africa during the Mio/Pliocene, P. botori in the Late Miocene of Eastern Africa (6.1–5.5 Ma) and Plesiogulo aff. monspessulanus at the beginning of the Pliocene in southern Africa (5.2 Ma). Lastly, we report for the first time the presence of both Sivaonyx and Plesiogulo in MPPM and LQSM at Langebaanweg, suggesting that the differences observed from the locality may be produced by sedimentation or sampling biases instead of temporal replacement within the carnivoran guild.

First Draft Genome of the Sable, Martes zibellina

Members of genus Martes provide early warning signals about forest ecosystem health and are designated as a Management Indicator Species. As one of the most widespread members in Martes, the sable (Martes zibellina) is a circumboreal small predator found throughout all taiga zoogeographical zones of Eurasia and shows distinct population differentiation and morphological variations. To support further studies on striking local adaptation and population evolution, we present the first sable genome, assembled de novo from an individual originating in the Great Khingan Mountains (China). The assembled genome is 2.42 Gb, consisting of 15,814 scaffolds with a scaffold N50 of 5.20 Mb. Searches for complete Mammalia BUSCO (Benchmarking Universal Single-Copy Ortholog) gene groups found that 95.15% of the curated single-copy orthologs were assembled as complete, suggesting a high level of completeness of the genome. We totally predicted 19,413 protein-coding genes, and 0.82 Gb of repeat sequences was annotated. We also detected 1,257 olfactory receptor genes and found more functional olfactory receptor genes in sable than in other Mustelidae species, which provide a possible genetic explanation for the acute sense of smell of the sable for searching the preys under deep snow. Phylogenetic analyses revealed that the ferret (Mustela putorius furo) and sea otter (Enhydra lutris) form a clade that is sister to the sable, which was dated ∼16.4 Ma. Overall, our study provided the first reference genome for research in a broad range of areas including local adaptations, population evolution, conservation, and management for sable.

Functional Morphology and Morphological Diversification of Hind Limb Cross-Sectional Traits in Mustelid Mammals

Locomotor habits in mammals are strongly tied to limb bones’ lengths, diameters, and proportions. By comparison, fewer studies have examined how limb bone cross-sectional traits relate to locomotor habit. Here, we tested whether climbing, digging, and swimming locomotor habits reflect biomechanically meaningful differences in three cross-sectional traits rendered dimensionless— cross-sectional area (CSA), second moments of area (SMA), and section modulus (MOD)—using femora, tibiae, and fibulae of 28 species of mustelid. CSA and SMA represent resistance to axial compression and bending, respectively, whereas MOD represents structural strength. Given the need to counteract buoyancy in aquatic environments and soil’s high density, we predicted that natatorial and fossorial mustelids have higher values of cross-sectional traits. For all three traits, we found that natatorial mustelids have the highest values, followed by fossorial mustelids, with both of these groups significantly differing from scansorial mustelids. However, phylogenetic relatedness strongly influences diversity in cross-sectional morphology, as locomotor habit strongly correlates with phylogeny. Testing whether hind limb bone cross-sectional traits have evolved adaptively, we fit Ornstein–Uhlenbeck (OU) and Brownian motion (BM) models of trait diversification to cross-sectional traits. The cross-sectional traits of the femur, tibia, and fibula appear to have, respectively, diversified under a multi-rate BM model, a single rate BM model, and a multi-optima OU model. In light of recent studies on mustelid body size and elongation, our findings suggest that the mustelid body plan—and perhaps that of other mammals—is likely the sum of a suite of traits evolving under different models of trait diversification.

Adaptive Changes in Hemoglobin Function in High-Altitude Tibetan Canids Were Derived via Gene Conversion and Introgression

A key question in evolutionary biology concerns the relative importance of different sources of adaptive genetic variation, such as de novo mutations, standing variation, and introgressive hybridization. A corollary question concerns how allelic variants derived from these different sources may influence the molecular basis of phenotypic adaptation. Here, we use a protein-engineering approach to examine the phenotypic effect of putatively adaptive hemoglobin (Hb) mutations in the high-altitude Tibetan wolf that were selectively introgressed into the Tibetan mastiff, a high-altitude dog breed that is renowned for its hypoxia tolerance. Experiments revealed that the introgressed coding variants confer an increased Hb–O₂ affinity in conjunction with an enhanced Bohr effect. We also document that affinity-enhancing mutations in the β-globin gene of Tibetan wolf were originally derived via interparalog gene conversion from a tandemly linked β-globin pseudogene. Thus, affinity-enhancing mutations were introduced into the β-globin gene of Tibetan wolf via one form of intragenomic lateral transfer (ectopic gene conversion) and were subsequently introduced into the Tibetan mastiff genome via a second form of lateral transfer (introgression). Site-directed mutagenesis experiments revealed that the increased Hb–O₂ affinity requires a specific two-site combination of amino acid replacements, suggesting that the molecular underpinnings of Hb adaptation in Tibetan mastiff (involving mutations that arose in a nonexpressed gene and which originally fixed in Tibetan wolf) may be qualitatively distinct from functionally similar changes in protein function that could have evolved via sequential fixation of de novo mutations during the breed’s relatively short duration of residency at high altitude.

The Anatolian glacial refugium and human-mediated colonization: a phylogeographical study of the stone marten (Martes foina) in Turkey

The Anatolian Peninsula, comprising most of modern Turkey, is situated at the junction of Europe, the Middle East and Asia. Together with its complex geomorphological and climatic history, this has given rise to a rich fauna and flora, which exhibits a wide range of historical biogeographical patterns. The stone marten (Martes foina) is a small carnivore that is found across the temperate Palaearctic region including Anatolia, where it is often associated with habitats modified by humans, but few genetic data exist for this species. We sequenced a 1840-bp region of the mitochondrial genome from 97 martens sampled across the peninsula and intron 7 of the nuclear β-fibrinogen gene from 53 of these. Two mitochondrial lineages were recovered, with overlapping eastern and western distributions, but there was no geographical structure for the autosomal marker. Coalescent analyses indicated that both of the lineages originated during the Last Glacial Maximum, one of them within an eastern Anatolian refugium and the other in a western Anatolian or Balkan refugium. The western lineage colonized most of Europe in the Holocene, while the eastern lineage may be endemic to Anatolia, from where it colonized the Iberian Peninsula via human translocation. The presence of at least one refugial stone marten population highlights the importance of Anatolia to the preservation of genetic variation and biodiversity.

Aquatic Adaptation and Depleted Diversity: A Deep Dive into the Genomes of the Sea Otter and Giant Otter

Despite its recent invasion into the marine realm, the sea otter (Enhydra lutris) has evolved a suite of adaptations for life in cold coastal waters, including limb modifications and dense insulating fur. This uniquely dense coat led to the near-extinction of sea otters during the 18th-20th century fur trade and an extreme population bottleneck. We used the de novo genome of the southern sea otter (E. l. nereis) to reconstruct its evolutionary history, identify genes influencing aquatic adaptation, and detect signals of population bottlenecks. We compared the genome of the southern sea otter with the tropical freshwater-living giant otter (Pteronura brasiliensis) to assess common and divergent genomic trends between otter species, and with the closely related northern sea otter (E. l. kenyoni) to uncover population-level trends. We found signals of positive selection in genes related to aquatic adaptations, particularly limb development and polygenic selection on genes related to hair follicle development. We found extensive pseudogenization of olfactory receptor genes in both the sea otter and giant otter lineages, consistent with patterns of sensory gene loss in other aquatic mammals. At the population level, the southern sea otter and the northern sea otter showed extremely low genomic diversity, signals of recent inbreeding, and demographic histories marked by population declines. These declines may predate the fur trade and appear to have resulted in an increase in putatively deleterious variants that could impact the future recovery of the sea otter.

The earliest mustelid in North America

Until now, the pre-Miocene fossil record of mustelids in North America has been restricted to specimens attributable to oligobunine taxa and isolated remains tentatively allocated to the genus Plesictis. In the present study, we report on a nearly complete cranium and a referred dentary of a new genus and species of mustelid. The specimens were recovered from the Turtle Cove and Kimberly Members of the John Day Formation, Oregon, USA.

These excellently preserved specimens more confidently confirm the presence of mustelids in the Early and Late Oligocene (Early and Late Arikareean) of North America. Like the holotype specimen of ‘Plesictis’ julieni, the new species lacks an alisphenoid canal and a postprotocrista on the M1 (synapomorphies of Mustelidae), but retains a dorsally deep suprameatal fossa (a feature occasionally suggested to be unique to Procyonidae). Phylogenetic analyses, applying parsimony and Bayesian inference to combined molecular (five genes totalling 5490 bp) and morphological data, recover this new species of mustelid as sister-species to ‘Plesictis’ julieni. The results of these analyses reveal that the new genus is a close relative of other species of Plesictis and several taxa traditionally allied with Oligobuninae, thereby rendering Oligobuninae paraphyletic. We further discuss the significance of the relatively small size of this new mustelid as it relates to predictions based on increased aridification of the palaeoclimate and the expansion of open habitats in the Oligocene.