Comparative skeletal anatomy of salt marsh and western harvest mice in relation to locomotor ecology.

The salt marsh harvest mouse (Reithrodontomys raviventris) is an endangered species, endemic to the San Francisco Bay Estuary, that co-occurs with the more broadly distributed species, the western harvest mouse (Reithrodontomys megalotis). Despite their considerable external morphological similarities, the northern subspecies of salt marsh harvest mice have relatively longer and thicker tails than do western harvest mice, which may be related to their abilities to climb emergent marsh vegetation to avoid tidal inundation. We used micro-CT to compare post-cranial skeletal anatomy between the salt marsh and western harvest mouse, to examine whether the salt marsh harvest mouse's restriction to brackish marshes is associated with skeletal adaptations for scansorial locomotion. We found that salt marsh harvest mice exhibited a deeper 3rd caudal vertebra, a more caudally located longest tail vertebra, craniocaudally longer tail vertebrae, and a longer digit III proximal phalanx than western harvest mice. These phalangeal and vertebral characteristics are known to decrease body rotations during climbing, increase contact with substrates, and decrease fall susceptibility in arboreal mammals, suggesting that the salt marsh harvest mouse may be morphologically specialized for scansorial locomotion, adaptive for its dynamic wetland environment.

Population Genetics of California Gray Foxes Clarify Origins of the Island Fox.

We used mitochondrial sequences and nuclear microsatellites to investigate population structure of gray foxes (Urocyon cinereoargenteus) and the evolutionary origins of the endemic island fox (Urocyon littoralis), which first appeared in the northern Channel Islands <13,000 years ago and in the southern Channel Islands <6000 years ago. It is unclear whether island foxes evolved directly from mainland gray foxes transported to the islands one or more times or from a now-extinct mainland population, already diverged from the gray fox. Our 345 mitochondrial sequences, combined with previous data, confirmed island foxes to be monophyletic, tracing to a most recent common ancestor approximately 85,000 years ago. Our rooted nuclear DNA tree additionally indicated genome-wide monophyly of island foxes relative to western gray foxes, although we detected admixture in northern island foxes from adjacent mainland gray foxes, consistent with some historical gene flow. Southern California gray foxes also bore a genetic signature of admixture and connectivity to a desert population, consistent with partial replacement by a late-Holocene range expansion. Using our outgroup analysis to root previous nuclear sequence-based trees indicated reciprocal monophyly of northern versus southern island foxes. Results were most consistent with island fox origins through multiple introductions from a now-extirpated mainland population.

Dietary characterization of the endangered salt marsh harvest mouse and sympatric rodents using DNA metabarcoding.

The salt marsh harvest mouse (Reithrodontomys raviventris; RERA) is an endangered species endemic to the coastal wetlands of the San Francisco Estuary, California. RERA are specialized to saline coastal wetlands, and their historical range has been severely impacted by landscape conversion and the introduction of non-native plant and rodent species. A better understanding of their diet is needed to assess habitat quality, particularly in relation to potential competitors. We investigated three questions using DNA metabarcoding with ITS2 and trnL markers: (1) Do RERA specialize on the native plant, pickleweed (Salicornia pacifica), (2) Do RERA consume non-native plants, and (3) What is the dietary niche breadth and overlap with three sympatric native and non-native rodents? RERA diet was dominated by two plants, native Salicornia and non-native salt bush (Atriplex spp.), but included 48 plant genera. RERA diet breadth was narrowest in fall, when they consumed the highest frequencies of Salicornia and Atriplex, and broadest in spring, when the frequencies of these two plants were lowest. Diet breadth was slightly lower for RERA than for co-occurring species in pairwise comparisons. All four species consumed similarly high frequencies of wetland plants, but RERA consumed fewer grasses and upland plants, suggesting that it may be less suited to fragmented habitat than sympatric rodents. Diet overlap was lowest between RERA and the native California vole (Microtis californicus). In contrast, RERA diet overlapped substantially with the native western harvest mouse (R. megalotis) and non-native house mouse (Mus musculus), suggesting potential for competition if these species become sufficiently abundant.

Next-generation phylogeography resolves post-glacial colonization patterns in a widespread carnivore, the red fox (Vulpes vulpes), in Europe.

Carnivores tend to exhibit a lack of (or less pronounced) genetic structure at continental scales in both a geographic and temporal sense and this can confound the identification of post-glacial colonization patterns in this group. In this study we used genome-wide data (using genotyping by sequencing [GBS]) to reconstruct the phylogeographic history of a widespread carnivore, the red fox (Vulpes vulpes), by investigating broad-scale patterns of genomic variation, differentiation and admixture amongst contemporary populations in Europe. Using 15,003 single nucleotide polymorphisms (SNPs) from 524 individuals allowed us to identify the importance of refugial regions for the red fox in terms of endemism (e.g., Iberia). In addition, we tested multiple post-glacial recolonization scenarios of previously glaciated regions during the Last Glacial Maximum using an Approximate Bayesian Computation (ABC) approach that were unresolved from previous studies. This allowed us to identify the role of admixture from multiple source population post-Younger Dryas in the case of Scandinavia and ancient land-bridges in the colonization of the British Isles. A natural colonization of Ireland was deemed more likely than an ancient human-mediated introduction as has previously been proposed and potentially points to a larger mammalian community on the island in the early post-glacial period. Using genome-wide data has allowed us to tease apart broad-scale patterns of structure and diversity in a widespread carnivore in Europe that was not evident from using more limited marker sets and provides a foundation for next-generation phylogeographic studies in other non-model species.

Pleistocene origins, western ghost lineages, and the emerging phylogeographic history of the red wolf and coyote.

The red wolf (Canis rufus) of the eastern US was driven to near-extinction by colonial-era persecution and habitat conversion, which facilitated coyote (C. latrans) range expansion and widespread hybridization with red wolves. The observation of some grey wolf (C. lupus) ancestry within red wolves sparked controversy over whether it was historically a subspecies of grey wolf with its predominant "coyote-like" ancestry obtained from post-colonial coyote hybridization (2-species hypothesis) versus a distinct species closely related to the coyote that hybridized with grey wolf (3-species hypothesis). We analysed mitogenomes sourced from before the 20th century bottleneck and coyote invasion, along with hundreds of modern amplicons, which led us to reject the 2-species model and to investigate a broader phylogeographic 3-species model suggested by the fossil record. Our findings broadly support this model, in which red wolves ranged the width of the American continent prior to arrival of the grey wolf to the mid-continent 60-80 ka; red wolves subsequently disappeared from the mid-continent, relegated to California and the eastern forests, which ushered in emergence of the coyote in their place (50-30 ka); by the early Holocene (12-10 ka), coyotes had expanded into California, where they admixed with and phenotypically replaced western red wolves in a process analogous to the 20th century coyote invasion of the eastern forests. Findings indicate that the red wolf pre-dated not only European colonization but human, and possibly coyote, presence in North America. These findings highlight the urgency of expanding conservation efforts for the red wolf.

Phylogeography of Cedros and Tiburón Island Mule Deer in North America's Desert Southwest.

Though mule deer (Odocoileus hemionus) persist in robust populations throughout most of their North American distribution, habitat loss, unregulated hunting, and other factors have reduced their historical range in México. Two of the 6 putative subspecies inhabiting México's deserts and Baja California peninsula are of conservation concern, occupying islands in the Pacific Ocean (Odocoileus hemionus cerrosensis on Cedros Island: endangered) and Sea of Cortés (Odocoileus hemionus sheldoni on Tiburón Island: threatened). Focusing on the desert southwest (n = 448), we sampled Tiburón (n = 22) and Cedros (n = 15) Island mule deer using contemporary samples and natural history museum specimens to complete a phylogeographic evaluation of the species complex, and assess the phylogeography of these insular subspecies. Both insular subspecies formed endemic haplotype lineages, consistent with island biogeographic theory. Bayesian skyline plots were consistent with Holocene demographic expansion. Cedros Island deer were genetically most similar to adjacent mainland Baja California deer, but exhibited a suite of unique haplotypes and reduced genetic variation. Tiburón Island deer haplotypes unexpectedly nested within a mainland lineage found in distant New Mexico, rather than the adjacent mainland Sonoran lineage. Such findings suggest the importance of postglacial climate fluctuations and biotic community turnover in the phylogeographic history of mule deer in the desert southwest. Our genetic data corroborates cultural, archaeological, and phenotypic evidence supporting Cedros and Tiburón deer endemicity and subspecies status. Reduced genetic variation, divergence from mainland populations, and demographic trends on both islands indicate that conservation, monitoring, and management are critical to ensure persistence of these endemic insular subspecies.

Historical Population Size Change and Differentiation of Relict Populations of the Endangered Giant Kangaroo Rat.

From a conservation management perspective it is important to understand how genetic diversity is partitioned across a species' range, including 1) identification of evolutionarily distinct units versus those recently isolated through anthropogenic activities and 2) the relative genetic contributions among components of fragmented (meta)populations. To address these questions, we investigated the phylogeography and metapopulation structure among relict populations of the endangered giant kangaroo rat (Dipodomys ingens) in the highly altered San Joaquin Desert Ecosystem. This keystone species underwent a ~97% range reduction over the past century, resulting in a current range that is highly fragmented, with 2 dominant northern and southern populations occurring 150 km apart. We sequenced >800 bp of mitochondrial DNA and genotyped 17 nuclear microsatellites in >275 D. ingens to assess the evolutionary relationship of these populations as well as the genetic structure within the northern metapopulation. A Bayesian Skyline Plot indicated that the species experienced a demographic expansion toward the end of the Pleistocene, with a recent population decline. Northern and southern D. ingens split 1857-13 443 years ago, prior to the massive conversion of the San Joaquin Valley to irrigated agriculture. We recommend that the northern and southern populations of D. ingens be re-classified as distinct population segments under the United States Endangered Species Act. We also observed population structure and asymmetrical migration within northern D. ingens where the Tumey Hills acted as a source contributing gene flow to all peripheral populations. This emphasized the importance of this location in the conservation of the metapopulation as a whole.

Nuclear Genetic Analysis of the Red Fox Across its Trans-Pacific Range.

The red fox (Vulpes vulpes) occurs on multiple continents in diverse habitats, making it an informative system for evolutionary genomic research. However, its phylogeography remains unclear. Previously, mitochondrial DNA and small numbers of nuclear loci provided discordant views. Both markers indicated deep divergence (~ 0.5 million years [MY]) between Eurasian and southern North American populations but differed in the apparent continental affinity of Alaskan red foxes, implying some degree of gene exchange during secondary contact (~0.1 MY). We assayed >173000 nuclear genomic sites in 52 red foxes, along with 2 Rueppell's foxes (Vulpes rueppellii) and a gray wolf (Canis lupus) using the Illumina CanineHD BeadChip. We obtained 5107 single nucleotide polymorphisms (SNPs) in the foxes. Consistent with the Afro-Eurasian origins of red foxes, genetic diversity was higher in Eurasian than North American samples. Phylogenetic trees indicated that Alaskan and southern North American red foxes formed a monophyletic group nested within the Eurasian clade. However, admixture models suggested Alaskan red foxes contained up to 40% Eurasian ancestry. We hypothesize that North American red foxes either hybridized with Eurasian foxes in Beringia at the start of the last glaciation or merged with a Beringian population after the last glaciation. Future work is needed to test between these scenarios and assess speciation.

Mitochondrial Analysis of the Most Basal Canid Reveals Deep Divergence between Eastern and Western North American Gray Foxes (Urocyon spp.) and Ancient Roots in Pleistocene California.

Pleistocene aridification in central North America caused many temperate forest-associated vertebrates to split into eastern and western lineages. Such divisions can be cryptic when Holocene expansions have closed the gaps between once-disjunct ranges or when local morphological variation obscures deeper regional divergences. We investigated such cryptic divergence in the gray fox (Urocyon cinereoargenteus), the most basal extant canid in the world. We also investigated the phylogeography of this species and its diminutive relative, the island fox (U. littoralis), in California. The California Floristic Province was a significant source of Pleistocene diversification for a wide range of taxa and, we hypothesized, for the gray fox as well. Alternatively, gray foxes in California potentially reflected a recent Holocene expansion from further south. We sequenced mitochondrial DNA from 169 gray foxes from the southeastern and southwestern United States and 11 island foxes from three of the Channel Islands. We estimated a 1.3% sequence divergence in the cytochrome b gene between eastern and western foxes and used coalescent simulations to date the divergence to approximately 500,000 years before present (YBP), which is comparable to that between recognized sister species within the Canidae. Gray fox samples collected from throughout California exhibited high haplotype diversity, phylogeographic structure, and genetic signatures of a late-Holocene population decline. Bayesian skyline analysis also indicated an earlier population increase dating to the early Wisconsin glaciation (~70,000 YBP) and a root height extending back to the previous interglacial (~100,000 YBP). Together these findings support California's role as a long-term Pleistocene refugium for western Urocyon. Lastly, based both on our results and re-interpretation of those of another study, we conclude that island foxes of the Channel Islands trace their origins to at least 3 distinct female founders from the mainland rather than to a single matriline, as previously suggested.

Genetic footprints reveal geographic patterns of expansion in Fennoscandian red foxes.

Population expansions of boreal species are among the most substantial ecological consequences of climate change, potentially transforming both structure and processes of northern ecosystems. Despite their importance, little is known about expansion dynamics of boreal species. Red foxes (Vulpes vulpes) are forecasted to become a keystone species in northern Europe, a process stemming from population expansions that began in the 19th century. To identify the relative roles of geographic and demographic factors and the sources of northern European red fox population expansion, we genotyped 21 microsatellite loci in modern and historical (1835-1941) Fennoscandian red foxes. Using Bayesian clustering and Bayesian inference of migration rates, we identified high connectivity and asymmetric migration rates across the region, consistent with source-sink dynamics, whereby more recently colonized sampling regions received immigrants from multiple sources. There were no clear clines in allele frequency or genetic diversity as would be expected from a unidirectional range expansion from south to north. Instead, migration inferences, demographic models and comparison to historical red fox genotypes suggested that the population expansion of the red fox is a consequence of dispersal from multiple sources, as well as in situ demographic growth. Together, these findings provide a rare glimpse into the anatomy of a boreal range expansion and enable informed predictions about future changes in boreal communities.

Range-wide multilocus phylogeography of the red fox reveals ancient continental divergence, minimal genomic exchange and distinct demographic histories.

Widely distributed taxa provide an opportunity to compare biogeographic responses to climatic fluctuations on multiple continents and to investigate speciation. We conducted the most geographically and genomically comprehensive study to date of the red fox (Vulpes vulpes), the world's most widely distributed wild terrestrial carnivore. Analyses of 697 bp of mitochondrial sequence in ~1000 individuals suggested an ancient Middle Eastern origin for all extant red foxes and a 400 kya (SD = 139 kya) origin of the primary North American (Nearctic) clade. Demographic analyses indicated a major expansion in Eurasia during the last glaciation (~50 kya), coinciding with a previously described secondary transfer of a single matriline (Holarctic) to North America. In contrast, North American matrilines (including the transferred portion of Holarctic clade) exhibited no signatures of expansion until the end of the Pleistocene (~12 kya). Analyses of 11 autosomal loci from a subset of foxes supported the colonization time frame suggested by mtDNA (and the fossil record) but, in contrast, reflected no detectable secondary transfer, resulting in the most fundamental genomic division of red foxes at the Bering Strait. Endemic continental Y-chromosome clades further supported this pattern. Thus, intercontinental genomic exchange was overall very limited, consistent with long-term reproductive isolation since the initial colonization of North America. Based on continental divergence times in other carnivoran species pairs, our findings support a model of peripatric speciation and are consistent with the previous classification of the North American red fox as a distinct species, V. fulva.

Temporal genetic variation of the red fox, Vulpes vulpes, across western Europe and the British Isles.

Quaternary climatic fluctuations have had profound effects on the phylogeographic structure of many species. Classically, species were thought to have become isolated in peninsular refugia, but there is limited evidence that large, non-polar species survived outside traditional refugial areas. We examined the phylogeographic structure of the red fox (Vulpes vulpes), a species that shows high ecological adaptability in the western Palaearctic region. We compared mitochondrial DNA sequences (cytochrome b and control region) from 399 modern and 31 ancient individuals from across Europe. Our objective was to test whether red foxes colonised the British Isles from mainland Europe in the late Pleistocene, or whether there is evidence that they persisted in the region through the Last Glacial Maximum. We found red foxes to show a high degree of phylogeographic structuring across Europe and, consistent with palaeontological and ancient DNA evidence, confirmed via phylogenetic indicators that red foxes were persistent in areas outside peninsular refugia during the last ice age. Bayesian analyses and tests of neutrality indicated population expansion. We conclude that there is evidence that red foxes from the British Isles derived from central European populations that became isolated after the closure of the landbridge with Europe.

On the origin of a domesticated species: Identifying the parent population of Russian silver foxes (Vulpes vulpes).

The foxes at Novosibirsk, Russia, are the only population of domesticated foxes in the world. These domesticated foxes originated from farm-bred silver foxes (Vulpes vulpes), whose genetic source is unknown. In this study we examined the origin of the domesticated strain of foxes and two other farm-bred fox populations (aggressive and unselected) maintained in Novosibirsk. To identify the phylogenetic origin of these populations we sequenced two regions of mtDNA, cytochrome b and D-loop, from 24 Novosibirsk foxes (8 foxes from each population) and compared them with corresponding sequences of native red foxes from Europe, Asia, Alaska and Western Canada, Eastern Canada, and the Western Mountains of the USA. We identified seven cytochrome b - D-loop haplotypes in Novosibirsk populations, four of which were previously observed in Eastern North America. The three remaining haplotypes differed by one or two base change from the most common haplotype in Eastern Canada. Φ(ST) analysis showed significant differentiation between Novosibirsk populations and red fox populations from all geographic regions except Eastern Canada. No haplotypes of Eurasian origin were identified in the Novosibirsk populations. These results are consistent with historical records indicating that the original breeding stock of farm-bred foxes originated from Prince Edward Island, Canada. Mitochondrial DNA data together with historical records indicate two stages in the selection of domesticated foxes: the first includes captive breeding for ~50 years with unconscious selection for behaviour; the second corresponds to over 50 further years of intensive selection for tame behaviour.

A restricted hybrid zone between native and introduced red fox (Vulpes vulpes) populations suggests reproductive barriers and competitive exclusion.

Introduced species can threaten native taxa in multiple ways, including competition and hybridization, which can reduce fitness, alter ecological niches or swamp native genomes. Encroachment and hybridization by introduced species also provide opportunities to study the dynamics of invasiveness and hybridization during early stages following contact. We used 33 microsatellites, 51 single nucleotide polymorphisms and a mtDNA marker to characterize the extent and spatial pattern of encroachment and hybridization between a native, endemic subspecies of red fox (Vulpes vulpes patwin) and an introduced red fox population composed of highly admixed, phylogenetically divergent stock, resulting from a century of domestication. Both nuclear and mtDNA markers indicated that hybridization was primarily restricted to a narrow zone where the two populations came into contact. Although a few introgressed genotypes were detected in the interior of the native range, we found no immigrant foxes or F(1) or F(2) hybrids there, suggesting native foxes excluded introduced individuals. We speculate that the observed interbreeding at the periphery was facilitated by low densities. In total, 98% of mtDNA haplotypes in the native range were native and 96% of the nuclear ancestry was estimated to be native. Although the introduced range had expanded fivefold over the past four decades, native and non-native haplotypes from museum samples collected in and near the native range three decades earlier showed a similar geographic distribution as today, suggesting that the native range and hybrid zone were relatively stable. We hypothesize that the monogamous mating system of red foxes and other wild canids may enhance their resistance to hybridization because of greater fitness consequences associated with mate discrimination.

Phylogeography of the North American red fox: vicariance in Pleistocene forest refugia.

Fossil, archaeological, and morphometric data suggest that indigenous red foxes in North America were derived from vicariance in two disjunct refugia during the last glaciation: one in Beringia and one in the contiguous USA. To test this hypothesis, we conducted a phylogeographical analysis of the North American red fox within its presettlement range. We sequenced portions of the mitochondrial cytochrome b (354 bp) gene and D-loop (342 bp) from 220 historical red fox specimens. Phylogenetic analysis of the cytochrome b gene produced two clades that diverged c. 400,000 years before present (bp): a Holarctic and a Nearctic clade. D-loop analyses of the Nearctic clade indicated three distinct subclades (> or = 99% Bayesian posterior probability); two that were more recently derived (rho estimate c. 20,000 bp) and were restricted to the southwestern mountains and the eastern portion of North America, and one that was older (rho estimate c. 45,000 bp) and more widespread in North America. Populations that migrated north from the southern refugium following deglaciation were derived from the colonization of North America during or prior to the Illinoian glaciation (300,000-130,000 bp), whereas populations that migrated south from the northern refugium represent a more recent colonization event during the Wisconsin glaciation (100,000-10,000 bp). Our findings indicate that Nearctic clade red foxes are phylogenetically distinct from their Holarctic counterparts, and reflect long-term isolation in two disjunct forest refugia during the Pleistocene. The montane lineage, which includes endangered populations, may be ecologically and evolutionarily distinct.