Dietary differentiation and the evolution of population genetic structure in a highly mobile carnivore

High-throughput sequencing and fatty acid profile analyses of the Black Amur bream (Megalobrama terminalis) reveal variation in dietary niche associated with geographic segregation

Fish dietary niche is a core focus, and it reflects the diversity of resources, habitats, or environments occupied by a species. However, whether geographic segregation among different populations triggers dietary diversification and concomitant fish niche shift remains unknown. In the present study, we selected the Black Amur bream (Megalobrama terminalis) is a migratory fish species that plays an important role in the material transfer and energy cycling of river ecosystems, inhabiting southern China drainage with multiple geographic populations. Here, we utilized the combined analyses of 18S rDNA high-throughput sequencing in fish gut contents and fatty acid (FA) in muscle tissues to evaluate potential spatial patterns of habitat and resource use for M. terminalis in three rivers of southern China. Our results showed that prey items of the Xijiang (XR) population (Pearl River) exhibited the highest species diversity and richness among the three geographic populations. Moreover, diet composition of M. terminalis was affected by spatial differences associated with geographic segregation. Analyses of FA biomarkers indicated that the highest levels of C16:0, C18:3n-3, and C18:2n-6c were found in Wanquan (WS) population (Wanquan River). The XR population exhibited a distinct FA profile characterized by higher amounts of arachidonic acid (ARA) and docosahexaenoic acid (DHA). The Moyang (MY) population (Moyang River) acted as the linkage between WS and XR populations and consisted of middle levels of saturated FAs (SFAs) and polyunsaturated FAs (PUFAs). The XR population displayed a greater FA niche width compared with WS population. Furthermore, we observed a close positive relationship between the niche width and α-diversity indices of dietary resources for FA proflies. Our study provides valued information to develop different conservation strategies among different populations and improve fisheries management for M. terminalis and other endemic species in local rivers.

Evolutionary history of the extinct wolf population from France in the context of global phylogeographic changes throughout the Holocene

Phylogeographic patterns in large mammals result from natural environmental factors and anthropogenic effects, which in some cases include domestication. The grey wolf was once widely distributed across the Holarctic, but experienced phylogeographic shifts and demographic declines during the Holocene. In the 19th-20th centuries, the species became extirpated from large parts of Europe due to direct extermination and habitat loss. We reconstructed the evolutionary history of the extinct Western European wolves based on the mitogenomic composition of 78 samples from France (Neolithic-20th century) in the context of other populations of wolves and dogs worldwide. We found a close genetic similarity of French wolves from ancient, medieval and recent populations, which suggests the long-term continuity of maternal lineages. MtDNA haplotypes of the French wolves showed large diversity and fell into two main haplogroups of modern Holarctic wolves. Our worldwide phylogeographic analysis indicated that haplogroup W1, which includes wolves from Eurasia and North America, originated in Northern Siberia. Haplogroup W2, which includes only European wolves, originated in Europe ~35 kya and its frequency was reduced during the Holocene due to an expansion of haplogroup W1 from the east. Moreover, we found that dog haplogroup D, currently restricted to Europe and the Middle East, was nested within the wolf haplogroup W2. This suggests European origin of haplogroup D, probably as a result of an ancient introgression from European wolves. Our results highlight the dynamic evolutionary history of European wolves during the Holocene, with a partial lineage replacement and introgressive hybridization with local dog populations.

Maybe So, Maybe Not: Canis lepophagus at Hagerman Fossil Beds National Monument, Idaho, USA

A canid dentary is described from the Pliocene Glenns Ferry Formation at Hagerman Fossil Beds National Monument, south-central Idaho, USA. The specimen possesses traits in alliance with and measurements falling within or exceeding those of Canis lepophagus. The dentary, along with a tarsal IV (cuboid) and an exploded canine come from the base of the fossiliferous Sahara complex within the monument. Improved geochronologic control provided by new tephrochronologic mapping by the U.S. Geological Survey-National Park Service Hagerman Paleontology, Environments, and Tephrochronology Project supports an interpolated age of approximately 3.9 Ma, placing it in the early Blancan North American Land Mammal Age. It is conservatively referred to herein as Canis aff. C. lepophagus with the caveat that it is an early and robust example of that species. A smaller canid, initially assigned to Canis lepophagus and then to Canis ferox, is also known from Hagerman. Most specimens of Canis ferox, including the holotype, were recently reassigned to Eucyon ferox, but specimens from the Hagerman and Rexroad faunas were left as Canis sp. and possibly attributed to C. lepophagus. We agree that these smaller canids belong in Canis and not Eucyon but reject placing them within C. lepophagus; we refer to them here as Hagerman-Rexroad Canis. This study confirms the presence of two approximately coyote-sized canids at Hagerman and adds to the growing list of carnivorans now known from these fossil beds.

Population genetic structure of wolves in the northwestern Dinaric-Balkan region

The Balkan Peninsula and the Dinaric Mountains possess extraordinary biodiversity and support one of the largest and most diverse wolf (Canis lupus) populations in Europe. Results obtained with diverse genetic markers show west-east substructure, also seen in various other species, despite the absence of obvious barriers to movement. However, the spatial extent of the genetic clusters remains unresolved, and our aim was to combine fine-scale sampling with population and spatial genetic analyses to improve resolution of wolf genetic clusters. We analyzed 16 autosomal microsatellites from 255 wolves sampled in Slovenia, Croatia, Bosnia and Herzegovina (BIH), and Serbia and documented three genetic clusters. These comprised (1) Slovenia and the regions of Gorski kotar and Lika in Croatia, (2) the region of Dalmatia in southern Croatia and BIH, and (3) Serbia. When we mapped the clusters geographically, we observed west-east genetic structure across the study area, together with some specific structure in BIH-Dalmatia. We observed that cluster 1 had a smaller effective population size, consistent with earlier reports of population recovery since the 1980s. Our results provide foundation for future genomic studies that would further resolve the observed west-east population structure and its evolutionary history in wolves and other taxa in the region and identify focal areas for habitat conservation. They also have immediate importance for conservation planning for the wolves in one of the most important parts of the species' European range.

The Indispensable Dog

Dogs' remarkable success in living in a human-dominated world rests on a set of adaptations to cohabitation with humans. In this paper, I review the nature of these adaptations. They include changes in reproductive and foraging behavior from their ancestor species, wolves, which can be understood as adaptations to the change from hunting live prey to feeding on human food residues. Dogs also show several changes in social behavior which are more controversial and even somewhat paradoxical. Contrary to theories of canine domestication which view dogs as less aggressive and more cooperative than wolves, several studies show that dogs' social interactions with conspecifics are more hierarchical and competitive than are wolves'. As scavengers rather than hunters, dogs do not need to cooperate with conspecifics the way that wolves do. But how then can we understand dogs' willingness to cooperate with humans? I propose an integrated account of dogs' social behavior that does not assume that dogs need to recognize the species-identity of the individuals with whom they interact. Because of the overlap in formal signals of dominance and submission between dog and human and people's complete control over the resources dogs need, I propose that people occupy a status of "super-dominance" over dogs. This conception suggests several new lines of research which could shed light on the human-dog relationship to the benefit of both partners.

The dynamics of prey selection by the trap-building predator Gasteracantha hasselti

Prey selection by generalist predators can be highly dynamic depending on the prey community structure. However, the dynamics of prey selection at the stage of prey entrapping are rarely investigated in trap-building predators, probably because their traps have been previously considered to intercept mobile prey proportionally to its availability in environment. Here we investigated the dynamics of prey selection by the orb-weaving spider Gasteracantha hasselti (Araneidae) depending on the composition of the available prey in tropical lowland forests located in north-eastern Thailand. We found that Gasteracantha captured a wide variety of prey but selected, on average, mostly Coleoptera and Diptera. The selectivity of Gasteracantha’s webs for Coleoptera was constant across the changes in overall prey availability and prey composition. The web selectivity for Hemiptera decreased rapidly with increasing relative densities of Hemiptera in the environment. The selectivity for Diptera and Hymenoptera increased and decreased, respectively, with their absolute densities in the environment. The relative selectivity of Gasteracantha’s traps for a particular prey type was driven by the presence and density of the highly selected prey rather than overall prey density. The results show that the selectivity of Gasteracantha’s traps for prey had both fixed and dynamic components and the dynamic component was determined by the relative as well as absolute densities of the particular prey types in the environment.

Cryptic population structure reveals low dispersal in Iberian wolves

Highly mobile mammalian carnivores are expected to have the capability to maintain high levels of gene flow across large geographic scales. Nonetheless, surprising levels of genetic structure have been found in many such populations. We combined genetic and spatial behavioural information from wolves (Canis lupus) in the Iberian Peninsula (Western Europe) during the last two decades to present a particular case of low dispersal levels in a large carnivore population persisting in human-dominated landscapes. We found an exceptionally reticulated pattern of cryptic population structure emerging at two hierarchical levels, in which four or eleven meaningful genetic clusters can be recognized, respectively. These clusters were characterized by moderate-high levels of differentiation (average pairwise F_ST = 0.09–0.19), low levels of admixture and varying degrees of genetic diversity. The number of dispersers identified among the 11 clusters was very low (<4% out of 218 wolves). Spatial information of tracked wolves further confirmed the geographical genetic patterns (only 2 out of 85 collared wolves overlapped with more than one genetic cluster). The high levels of genetic structure in this population may be determined by the recent demographic history of this population, among other factors. The identification of meaningful genetic clusters has implications for the delineation of conservation units and, consequently, on the conservation and management actions for Iberian wolves.

The host age related occurrence of Alaria alata in wild canids in Latvia

The trematode Alaria alata has a complex life cycle in which carnivore mammals may become infected by feeding on mesocercariae-infected, second intermediate hosts or paratenic hosts. Afterwards, young flukes migrate through various organs of the definitive host, including the diaphragm and the lungs, before reaching the small intestine. The aim of the present study was to establish the prevalence and infection intensity of A. alata in carnivore mammals, the host age and related occurrence of different A. alata development stages in particular animals in Latvia. Overall, samples from 992 animals (539 red foxes, 411 raccoon dogs and 42 grey wolves) were examined using sedimentation and counting techniques to observe the presence of A. alata metacercariae and adult trematodes in each particular host. The present study, based on data collected over five hunting seasons (2010/2011-2014/2015), is a comprehensive study on the circulation of A. alata in wild canids in Latvia. A. alata was present in all species examined and the overall prevalence reached 86.2%. A significant positive correlation was observed between the burden of A. alata metacercariae and adult trematodes in raccoon dogs (r_s = 0.27, p = 0.001) and red foxes (r_s = 0.17, p = 0.0002). A significantly higher (p = 0.003) A. alata metacercariae mean intensity was observed in juvenile raccoon dogs compared to the one observed in adult raccoon dogs and red foxes from both age groups. The presence of both A. alata development stages in a particular animal concurrently and at different intensity levels may indicate primary infection and re-cured infection of the host.

Genetic diversity, genetic structure and diet of ancient and contemporary red deer (Cervus elaphus L.) from north-eastern France

In north-eastern France, red deer (Cervus elaphus L.) populations were rebuilt from a few hundred individuals, which have subsisted in remote valleys of the Vosges mountains, and to a lesser extent from individuals escaped from private enclosures; at present times, this species occupies large areas, mainly in the Vosges Mountains. In this study, we examined the population dynamics of red deer in the Vosges Mountains using ancient and contemporary mitochondrial DNA (mtDNA) from 140 samples (23 ancient + 117 modern) spanning the last 7'000 years. In addition, we reconstructed the feeding habits and the habitat of red deer since the beginning of agriculture applying isotopic analyses in order to establish a basis for current environmental management strategies. We show that past and present red deer in the Vosges Mountains belong to mtDNA haplogroup A, suggesting that they originated from the Iberian refugium after the last glacial maximum (LGM). Palaeogenetic analysis of ancient bone material revealed the presence of two distinct haplotypes with different temporal distributions. Individuals belonging to the two haplotype groups apparently occupied two different habitats over at least 7'000 years. AM6 correlates with an ecological type that feeds in densely forested mountain landscapes, while AM235 correlates with feeding in lowland landscapes, composed of a mixture of meadows and riverine, herb-rich woodlands. Our results suggest that red deer of north-eastern France was able to adapt, over the long term, to these different habitat types, possibly due to efficient ethological barriers. Modern haplotype patterns support the historical record that red deer has been exposed to strong anthropogenic influences as a major game species.

Wolf population genetics in Europe: a systematic review, meta-analysis and suggestions for conservation and management

The grey wolf (Canis lupus) is an iconic large carnivore that has increasingly been recognized as an apex predator with intrinsic value and a keystone species. However, wolves have also long represented a primary source of human-carnivore conflict, which has led to long-term persecution of wolves, resulting in a significant decrease in their numbers, genetic diversity and gene flow between populations. For more effective protection and management of wolf populations in Europe, robust scientific evidence is crucial. This review serves as an analytical summary of the main findings from wolf population genetic studies in Europe, covering major studies from the 'pre-genomic era' and the first insights of the 'genomics era'. We analyse, summarize and discuss findings derived from analyses of three compartments of the mammalian genome with different inheritance modes: maternal (mitochondrial DNA), paternal (Y chromosome) and biparental [autosomal microsatellites and single nucleotide polymorphisms (SNPs)]. To describe large-scale trends and patterns of genetic variation in European wolf populations, we conducted a meta-analysis based on the results of previous microsatellite studies and also included new data, covering all 19 European countries for which wolf genetic information is available: Norway, Sweden, Finland, Estonia, Latvia, Lithuania, Poland, Czech Republic, Slovakia, Germany, Belarus, Russia, Italy, Croatia, Bulgaria, Bosnia and Herzegovina, Greece, Spain and Portugal. We compared different indices of genetic diversity in wolf populations and found a significant spatial trend in heterozygosity across Europe from south-west (lowest genetic diversity) to north-east (highest). The range of spatial autocorrelation calculated on the basis of three characteristics of genetic diversity was 650-850 km, suggesting that the genetic diversity of a given wolf population can be influenced by populations up to 850 km away. As an important outcome of this synthesis, we discuss the most pressing issues threatening wolf populations in Europe, highlight important gaps in current knowledge, suggest solutions to overcome these limitations, and provide recommendations for science-based wolf conservation and management at regional and Europe-wide scales.

A European Concern? Genetic Structure and Expansion of Golden Jackals (Canis aureus) in Europe and the Caucasus

In the first continent-wide study of the golden jackal (Canis aureus), we characterised its population genetic structure and attempted to identify the origin of European populations. This provided a unique insight into genetic characteristics of a native carnivore population with rapid large-scale expansion. We analysed 15 microsatellite markers and a 406 base-pair fragment of the mitochondrial control region. Bayesian-based and principal components methods were applied to evaluate whether the geographical grouping of samples corresponded with genetic groups. Our analysis revealed low levels of genetic diversity, reflecting the unique history of the golden jackal among Europe’s native carnivores. The results suggest ongoing gene flow between south-eastern Europe and the Caucasus, with both contributing to the Baltic population, which appeared only recently. The population from the Peloponnese Peninsula in southern Greece forms a common genetic cluster with samples from south-eastern Europe (ΔK approach in STRUCTURE, Principal Components Analysis [PCA]), although the results based on BAPS and the estimated likelihood in STRUCTURE indicate that Peloponnesian jackals may represent a distinct population. Moreover, analyses of population structure also suggest either genetic distinctiveness of the island population from Samos near the coast of Asia Minor (BAPS, most STRUCTURE, PCA), or possibly its connection with the Caucasus population (one analysis in STRUCTURE). We speculate from our results that ancient Mediterranean jackal populations have persisted to the present day, and have merged with jackals colonising from Asia. These data also suggest that new populations of the golden jackal may be founded by long-distance dispersal, and thus should not be treated as an invasive alien species, i.e. an organism that is “non-native to an ecosystem, and which may cause economic or environmental harm or adversely affect human health”. These insights into the genetic structure and ancestry of Baltic jackals have important implications for management and conservation of jackals in Europe. The golden jackal is listed as an Annex V species in the EU Habitats Directive and as such, considering also the results presented here, should be legally protected in all EU member states.

Ecological opportunities and specializations shaped genetic divergence in a highly mobile marine top predator

Environmental conditions can shape genetic and morphological divergence. Release of new habitats during historical environmental changes was a major driver of evolutionary diversification. Here, forces shaping population structure and ecotype differentiation ('pelagic' and 'coastal') of bottlenose dolphins in the North-east Atlantic were investigated using complementary evolutionary and ecological approaches. Inference of population demographic history using approximate Bayesian computation indicated that coastal populations were likely founded by the Atlantic pelagic population after the Last Glacial Maxima probably as a result of newly available coastal ecological niches. Pelagic dolphins from the Atlantic and the Mediterranean Sea likely diverged during a period of high productivity in the Mediterranean Sea. Genetic differentiation between coastal and pelagic ecotypes may be maintained by niche specializations, as indicated by stable isotope and stomach content analyses, and social behaviour. The two ecotypes were only weakly morphologically segregated in contrast to other parts of the World Ocean. This may be linked to weak contrasts between coastal and pelagic habitats and/or a relatively recent divergence. We suggest that ecological opportunity to specialize is a major driver of genetic and morphological divergence. Combining genetic, ecological and morphological approaches is essential to understanding the population structure of mobile and cryptic species.

Genome-wide analyses suggest parallel selection for universal traits may eclipse local environmental selection in a highly mobile carnivore

Ecological and environmental heterogeneity can produce genetic differentiation in highly mobile species. Accordingly, local adaptation may be expected across comparatively short distances in the presence of marked environmental gradients. Within the European continent, wolves (Canis lupus) exhibit distinct north–south population differentiation. We investigated more than 67‐K single nucleotide polymorphism (SNP) loci for signatures of local adaptation in 59 unrelated wolves from four previously identified population clusters (northcentral Europe n = 32, Carpathian Mountains n = 7, Dinaric‐Balkan n = 9, Ukrainian Steppe n = 11). Our analyses combined identification of outlier loci with findings from genome‐wide association study of individual genomic profiles and 12 environmental variables. We identified 353 candidate SNP loci. We examined the SNP position and neighboring megabase (1 Mb, one million bases) regions in the dog (C. lupus familiaris) genome for genes potentially under selection, including homologue genes in other vertebrates. These regions included functional genes for, for example, temperature regulation that may indicate local adaptation and genes controlling for functions universally important for wolves, including olfaction, hearing, vision, and cognitive functions. We also observed strong outliers not associated with any of the investigated variables, which could suggest selective pressures associated with other unmeasured environmental variables and/or demographic factors. These patterns are further supported by the examination of spatial distributions of the SNPs associated with universally important traits, which typically show marked differences in allele frequencies among population clusters. Accordingly, parallel selection for features important to all wolves may eclipse local environmental selection and implies long‐term separation among population clusters.

Population genetic structure of gray wolves (Canis lupus) in a marine archipelago suggests island-mainland differentiation consistent with dietary niche

BACKGROUND

Emerging evidence suggests that ecological heterogeneity across space can influence the genetic structure of populations, including that of long-distance dispersers such as large carnivores. On the central coast of British Columbia, Canada, wolf (Canis lupus L., 1758) dietary niche and parasite prevalence data indicate strong ecological divergence between marine-oriented wolves inhabiting islands and individuals on the coastal mainland that interact primarily with terrestrial prey. Local holders of traditional ecological knowledge, who distinguish between mainland and island wolf forms, also informed our hypothesis that genetic differentiation might occur between wolves from these adjacent environments.

RESULTS

We used microsatellite genetic markers to examine data obtained from wolf faecal samples. Our results from 116 individuals suggest the presence of a genetic cline between mainland and island wolves. This pattern occurs despite field observations that individuals easily traverse the 30 km wide study area and swim up to 13 km among landmasses in the region.

CONCLUSIONS

Natal habitat-biased dispersal (i.e., the preference for dispersal into familiar ecological environments) might contribute to genetic differentiation. Accordingly, this working hypothesis presents an exciting avenue for future research where marine resources or other components of ecological heterogeneity are present.

Genetic variability of the grey wolf Canis lupus in the Caucasus in comparison with Europe and the Middle East: distinct or intermediary population?

Despite continuous historical distribution of the grey wolf (Canis lupus) throughout Eurasia, the species displays considerable morphological differentiation that resulted in delimitation of a number of subspecies. However, these morphological discontinuities are not always consistent with patterns of genetic differentiation. Here we assess genetic distinctiveness of grey wolves from the Caucasus (a region at the border between Europe and West Asia) that have been classified as a distinct subspecies C. l. cubanensis. We analysed their genetic variability based on mtDNA control region, microsatellite loci and genome-wide SNP genotypes (obtained for a subset of the samples), and found similar or higher levels of genetic diversity at all these types of loci as compared with other Eurasian populations. Although we found no evidence for a recent genetic bottleneck, genome-wide linkage disequilibrium patterns suggest a long-term demographic decline in the Caucasian population--a trend consistent with other Eurasian populations. Caucasian wolves share mtDNA haplotypes with both Eastern European and West Asian wolves, suggesting past or ongoing gene flow. Microsatellite data also suggest gene flow between the Caucasus and Eastern Europe. We found evidence for moderate admixture between the Caucasian wolves and domestic dogs, at a level comparable with other Eurasian populations. Taken together, our results show that Caucasian wolves are not genetically isolated from other Eurasian populations, share with them the same demographic trends, and are affected by similar conservation problems.

Large-scale genetic structuring of a widely distributed carnivore--the Eurasian lynx (Lynx lynx)

Over the last decades the phylogeography and genetic structure of a multitude of species inhabiting Europe and North America have been described. The flora and fauna of the vast landmasses of north-eastern Eurasia are still largely unexplored in this respect. The Eurasian lynx is a large felid that is relatively abundant over much of the Russian sub-continent and the adjoining countries. Analyzing 148 museum specimens collected throughout its range over the last 150 years we have described the large-scale genetic structuring in this highly mobile species. We have investigated the spatial genetic patterns using mitochondrial DNA sequences (D-loop and cytochrome b) and 11 microsatellite loci, and describe three phylogenetic clades and a clear structuring along an east-west gradient. The most likely scenario is that the contemporary Eurasian lynx populations originated in central Asia and that parts of Europe were inhabited by lynx during the Pleistocene. After the Last Glacial Maximum (LGM) range expansions lead to colonization of north-western Siberia and Scandinavia from the Caucasus and north-eastern Siberia from a refugium further east. No evidence of a Berinigan refugium could be detected in our data. We observed restricted gene flow and suggest that future studies of the Eurasian lynx explore to what extent the contemporary population structure may be explained by ecological variables.

North-South differentiation and a region of high diversity in European wolves (Canis lupus)

European wolves (Canis lupus) show population genetic structure in the absence of geographic barriers, and across relatively short distances for this highly mobile species. Additional information on the location of and divergence between population clusters is required, particularly because wolves are currently recolonizing parts of Europe. We evaluated genetic structure in 177 wolves from 11 countries using over 67K single nucleotide polymorphism (SNP) loci. The results supported previous findings of an isolated Italian population with lower genetic diversity than that observed across other areas of Europe. Wolves from the remaining countries were primarily structured in a north-south axis, with Croatia, Bulgaria, and Greece (Dinaric-Balkan) differentiated from northcentral wolves that included individuals from Finland, Latvia, Belarus, Poland and Russia. Carpathian Mountain wolves in central Europe had genotypes intermediate between those identified in northcentral Europe and the Dinaric-Balkan cluster. Overall, individual genotypes from northcentral Europe suggested high levels of admixture. We observed high diversity within Belarus, with wolves from western and northern Belarus representing the two most differentiated groups within northcentral Europe. Our results support the presence of at least three major clusters (Italy, Carpathians, Dinaric-Balkan) in southern and central Europe. Individuals from Croatia also appeared differentiated from wolves in Greece and Bulgaria. Expansion from glacial refugia, adaptation to local environments, and human-related factors such as landscape fragmentation and frequent killing of wolves in some areas may have contributed to the observed patterns. Our findings can help inform conservation management of these apex predators and the ecosystems of which they are part.

A range-wide synthesis and timeline for phylogeographic events in the red fox (Vulpes vulpes)

BACKGROUND

Many boreo-temperate mammals have a Pleistocene fossil record throughout Eurasia and North America, but only few have a contemporary distribution that spans this large area. Examples of Holarctic-distributed carnivores are the brown bear, grey wolf, and red fox, all three ecological generalists with large dispersal capacity and a high adaptive flexibility. While the two former have been examined extensively across their ranges, no phylogeographic study of the red fox has been conducted across its entire Holarctic range. Moreover, no study included samples from central Asia, leaving a large sampling gap in the middle of the Eurasian landmass.

RESULTS

Here we provide the first mitochondrial DNA sequence data of red foxes from central Asia (Siberia), and new sequences from several European populations. In a range-wide synthesis of 729 red fox mitochondrial control region sequences, including 677 previously published and 52 newly obtained sequences, this manuscript describes the pattern and timing of major phylogeographic events in red foxes, using a Bayesian coalescence approach with multiple fossil tip and root calibration points. In a 335 bp alignment we found in total 175 unique haplotypes. All newly sequenced individuals belonged to the previously described Holarctic lineage. Our analyses confirmed the presence of three Nearctic- and two Japan-restricted lineages that were formed since the Mid/Late Pleistocene.

CONCLUSIONS

The phylogeographic history of red foxes is highly similar to that previously described for grey wolves and brown bears, indicating that climatic fluctuations and habitat changes since the Pleistocene had similar effects on these highly mobile generalist species. All three species originally diversified in Eurasia and later colonized North America and Japan. North American lineages persisted through the last glacial maximum south of the ice sheets, meeting more recent colonizers from Beringia during postglacial expansion into the northern Nearctic. Both brown bears and red foxes colonized Japan's northern island Hokkaido at least three times, all lineages being most closely related to different mainland lineages. Red foxes, grey wolves, and brown bears thus represent an interesting case where species that occupy similar ecological niches also exhibit similar phylogeographic histories.