Multigenic and morphometric differentiation of ground squirrels (Spermophilus, Scuiridae, Rodentia) in Turkey, with a description of a new species

Traditional multilocus phylogeny fails to fully resolve Palearctic ground squirrels (Spermophilus) relationships but reveals a new species endemic to West Siberia

Previous efforts to reconstruct evolutionary history of Palearctic ground squirrels within the genus Spermophilus have primarily relied on a single mitochondrial marker for phylogenetic data. In this study, we present the first phylogeny with comprehensive taxon sampling of Spermophilus via a conventional multilocus approach utilizing five mitochondrial and five nuclear markers. Through application of the multispecies coalescent model, we constructed a species tree revealing four distinct clades that diverged during the Late Miocene. These clades are 1) S. alaschanicus and S. dauricus from East Asia; 2) S. musicus and S. pygmaeus from East Europe and northwestern Central Asia; 3) the subgenus Colobotis found across Central Asia and its adjacent regions and encompassing S. brevicauda, S. erythrogenys, S. fulvus, S. major, S. pallidicauda, S. ralli, S. relictus, S. selevini, and S. vorontsovi sp. nov.; and 4) a Central/Eastern Europe and Asia Minor clade comprising S. citellus, S. taurensis, S. xanthoprymnus, S. suslicus, and S. odessanus. The latter clade lacked strong support owing to uncertainty of taxonomic placement of S. odessanus and S. suslicus. Resolving relationships within the subgenus Colobotis, which radiated rapidly, remains challenging likely because of incomplete lineage sorting and introgressive hybridization. Most of modern Spermophilus species diversified during the Early-Middle Pleistocene (2.2-1.0 million years ago). We propose a revised taxonomic classification for the genus Spermophilus by recognizing 18 species including a newly identified one (S. vorontsovi sp. nov.), which is found only in a limited area in the southeast of West Siberia. Employing genome-wide single-nucleotide polymorphism genotyping, we substantiated the role of the Ob River as a major barrier ensuring robust isolation of this taxon from S. erythrogenys. Despite its inherent limitations, the traditional multilocus approach remains a valuable tool for resolving relationships and can provide important insights into otherwise poorly understood groups. It is imperative to recognize that additional efforts are needed to definitively determine phylogenetic relationships between certain species of Palearctic ground squirrels.

Evolutionary history and systematics of European blind mole rats (Rodentia: Spalacidae: Nannospalax): Multilocus phylogeny and species delimitation in a puzzling group

Species delimitation is a powerful approach to assist taxonomic decisions in challenging taxa where species boundaries are hard to establish. European taxa of the blind mole rats (genus Nannospalax) display small morphological differences and complex chromosomal evolution at a shallow evolutionary divergence level. Previous analyses led to the recognition of 25 'forms' in their distribution area. We provide a comprehensive framework to improve knowledge on the evolutionary history and revise the taxonomy of European blind mole rats based on samples from all but three of the 25 forms. We sequenced two nuclear-encoded genetic regions and the whole mitochondrial cytochrome b gene for phylogenetic tree reconstructions using concatenation and coalescence-based species-tree estimations. The phylogenetic analyses confirmed that Aegean N. insularis belongs to N. superspecies xanthodon, and that it represents the second known species of this superspecies in Europe. Mainland taxa reached Europe from Asia Minor in two colonisation events corresponding to two superspecies-level taxa: N. superspecies monticola (taxon established herewith) reached Europe c. 2.1 million years ago (Mya) and was followed by N. superspecies leucodon (re-defined herewith) c. 1.5 Mya. Species delimitation allowed the clarification of the taxonomic contents of the above superspecies. N. superspecies monticola contains three species geographically confined to the western periphery of the distribution of blind mole rats, whereas N. superspecies leucodon is more speciose with six species and several additional subspecies. The observed geographic pattern hints at a robust peripatric speciation process and rapid chromosomal evolution. The present treatment is thus regarded as the minimum taxonomic content of each lineage, which can be further refined based on other sources of information such as karyological traits, crossbreeding experiments, etc. The species delimitation models also allowed the recognition of a hitherto unnamed blind mole rat taxon from Albania, described here as a new subspecies.

Geographic Genetic Structure of Alectoris chukar in Türkiye: Post-LGM-Induced Hybridization and Human-Mediated Contaminations

Türkiye is considered an important evolutionary area for Chukar partridge (Alectoris chukar), since it is both a potential ancestral area and a diversification center for the species. Using 2 mitochondrial (Cty-b and D-loop) and 13 polymorphic microsatellite markers, we investigated the geographic genetic structure of A. chukar populations to determine how past climatic fluctuations and human activities have shaped the gene pool of this species in Türkiye. Our results indicate, firstly, that only A. chukar of the genus Alectoris is present in Türkiye (Anatolia and Thrace), with no natural or artificial gene flow from congenerics. Secondly, the geographic genetic structure of the species in Türkiye has been shaped by topographic heterogeneity, Pleistocene climatic fluctuations, and artificial transport by humans. Third, there appears to be three genetic clusters: Thracian, Eastern, and Western. Fourth, the post-LGM demographic expansion of the Eastern and Western populations has formed a hybrid zone in Central Anatolia (~8 kyBP). Fifth, the rate of China clade-B contamination in Türkiye is about 8% in mtDNA and about 12% in nuDNA, with the Southeastern Anatolian population having the highest contamination. Sixth, the Thracian population was the most genetically distinct, with the lowest genetic diversity and highest level of inbreeding and no China clad-B contamination. These results can contribute to the conservation regarding A. chukar populations, especially the Thracian population.

Complete mitogenomes of Turkish tree squirrels, Sciurus anomalus and S. vulgaris, (Sciuridae: Rodentia: Mammalia) and their phylogenetic status within the tribe Sciurini

The genus Sciurus, a member of the family Sciuridae, is widely distributed in the Holarctic region. To better understand mitogenomic characteristics and to reveal internal phylogenetic relationships of the genus, 20 complete mitogenomes of Turkish tree squirrels were successfully sequenced for the first time, including 19 for S. anomalus (from 16,505 bp to 16,510 bp) and one for S. vulgaris (16,511 bp). The mitogenomes of two species were AT-biased. All tRNAs for two species displayed a typical clover-leaf structure, except for tRNA^Ser(AGY). The tRNA Serine1 (S1)-GCT structure lacked the dihydrouridine (DHU) loop and stem. Based on mitogenomic dataset for phylogeny of Sciurinae, phylogenetic analyses (Bayesian Inference and Maximum Likelihood) did not support monophyly of Sciurus and proposed that S. anomalus, the most basal taxa in the Sciurini tribe, had at least five mitogenome lineages, which were also supported by network analysis. The dissimilarities among the five linegaes of S. anomalus ranged from 0.0042 (0.42%) to 0.0062 (0.62%) using K2P sequence pairwise distances. In addition to this mitogenomic analysis result, phylogenetic analyses using the CYTB + D-loop dataset proposed the existence of at least nine lineages for S. anomalus, which was different than those of the previous studies. The current study proposed that the use of mitogenomic data for reconstructing the phylogeny of Turkey' Sciurus holds an important value for revealing evolutionary relationships.

The story of a rock-star: multilocus phylogeny and species delimitation in the starred or roughtail rock agama, Laudakia stellio (Reptilia: Agamidae)

Situated at the junction of three continents, Europe, Asia and Africa, the Eastern Mediterranean is an ideal region to study the effects of palaeogeography, ecology and long human presence on animal evolution. Laudakia stellio (Squamata: Agamidae) is found across this region and offers an excellent opportunity for such studies. The high morphological variation across their range suggests that these lizards might represent a species complex. This is the first study exploring their evolutionary history, using molecular markers and individuals from all described subspecies. We employed the latest phylogenetic and species-delimitation methods to identify all distinct evolutionary lineages, their genetic variation and divergence times. The phenotypical diversity of L. stellio matches its genetic differentiation: almost all subspecies correspond to well-supported retrieved subclades and additional distinct lineages representing intermediate morphs have been retrieved. ‘Laudakia stellio’ represents three distinct evolutionary entities that diverged during the Plio-Pleistocene transition, which we propose as distinct species. One includes Greek and Turkish populations, as well as cryptic Anatolian lineages. The second comprises all other Near East populations and the third is endemic to Cyprus. Our results indicate a role of humans in shaping present distribution patterns, and highlight the importance of the Aegean, Anatolia and the Levant as glacial refugia and diversity hotspots.

The roles of possible geographic barriers and geological events on the phylogeographic structure of the Eastern broad toothed field mouse (Apodemus mystacinus)

The Eastern broad toothed field mouse, Apodemus mystacinus, is a rodent species distributed in Turkey, the Middle East, and a few Aegean Islands. The aim of this study is to analyse the phylogeographic structure of A. mystacinus and possible causes of its differentiation, on the basis of mitochondrial and nuclear DNA sequences using a large number of new samples from Turkey. In this context, partial mitochondrial sequences of cytochrome b (Cytb), control region (D-loop) and a nuclear interphotoreceptor retinoid-binding protein (IRBP) gene were used to reveal the geographical differentiation among A. mystacinus populations and the validity of its subspecies. The estimated divergence times revealed that the first separation of A. mystacinus into three distinct groups (subspecies of A. mystacinus: A. m. mystacinus, A. m. smyrnensis, and A. m. euxinus) begun 0.641 Mya. The possible physical barriers in Anatolia such as high mountains and rivers could interrupt the gene flow between A. mystacinus populations. The results of the present study indicated that A. mystacinus might have used the high rocky areas along the Anatolian Diagonal as a dispersal way. Moreover, mitochondrial data in this study suggested for the first time that A. m. rhodius is synonymous with the nominative subspecies A. m. mystacinus.

Mitogenome Phylogeny Including Data from Additional Subspecies Provides New Insights into the Historical Biogeography of the Eurasian lynx Lynx lynx

Previous molecular studies of the wide-ranging Eurasian lynx Lynx lynx focused mainly on its northern Palearctic populations, with the consequence that the reconstruction of this species’ evolutionary history did not include genetic variation present in its southern Palearctic distribution. We sampled a previously not considered Asian subspecies (L. l. dinniki), added published data from another Asian subspecies (L. l. isabellinus), and reassessed the Eurasian lynx mtDNA phylogeny along with previously published data from northern Palearctic populations. Our mitogenome-based analyses revealed the existence of three major clades (A: Central Asia, B: SE Europe/SW Asia, C: Europe and Northern Asia) and at least five lineages, with diversification in Lynx lynx commencing at least 28kyr earlier than hitherto estimated. The subspecies L. l. isabellinus harbors the most basal matriline, consistent with the origin of Lynx lynx in this subspecies’ current range. L. l. dinniki harbors the second most basal matriline, which is related to, and may be the source of, the mtDNA diversity of the critically endangered Balkan lynx L. l. balcanicus. Our results suggest that the Anatolian peninsula was a glacial refugium for Eurasian lynx, with previously unconsidered implications for the colonization of Europe by this species.

From refugia to contact: Pine processionary moth hybrid zone in a complex biogeographic setting

Contact zones occur at the crossroad between specific dispersal routes and are facilitated by biogeographic discontinuities. Here, we focused on two Lepidoptera sister species that come in contact near the Turkish Straits System (TSS). We aimed to infer their phylogeographic histories in the Eastern Mediterranean and finely analyze their co-occurrence and hybridization patterns in this biogeographic context. We used molecular mitochondrial and nuclear markers to study 224 individuals from 42 localities. We used discordances between markers and complementary assignment methods to identify and map hybrids and parental individuals. We confirmed the parapatric distribution of Thaumetopoea pityocampa (Lepidoptera: Notodontidae) in the west and Thaumetopoea wilkinsoni in the east and identified a narrow contact zone. We identified several glacial refugia of T. wilkinsoni in southern Turkey with a strong east-west differentiation in this species. Unexpectedly, T. pityocampa crossed the TSS and occur in northern Aegean Turkey and some eastern Greek islands. We found robust evidence of introgression between the two species in a restricted zone in northwestern Turkey, but we did not identify any F1 individuals. The identified hybrid zone was mostly bimodal. The distributions and genetic patterns of the studied species were strongly influenced both by the Quaternary climatic oscillations and the complex geological history of the Aegean region. T. pityocampa and T. wilkinsoni survived the last glacial maximum in disjoint refugia and met in western Turkey at the edge of the recolonization routes. Expanding population of T. wilkinsoni constrained T. pityocampa to the western Turkish shore. Additionally, we found evidence of recurrent introgression by T. wilkinsoni males in several T. pityocampa populations. Our results suggest that some prezygotic isolation mechanisms, such as differences in timing of the adult emergences, might be a driver of the isolation between the sister species.

Phylogenetic relationship and variation of alarm call traits of populations of red-cheeked ground squirrels (Spermophilus erythrogenys sensu lato) suggest taxonomic delineation

Distribution area and taxonomic borders within the species complex Spermophilus erythrogenys sensu lato remain questionable. Early evidence suggests that red-cheeked ground squirrels of Southeast Kazakhstan are remarkably different in terms of the acoustic structure of their alarm calls from the red-cheeked ground squirrels of the Kurgan region in Russia. In this study, we analyzed the differences in the acoustic structure of the alarm call and mitochondrial DNA (complete control region, 1005-1006 bp and complete cytochrome b gene, 1140 bp) in 3 populations of red-cheeked ground squirrels (Tara, Altyn-Emel and Balkhash), all located within areas isolated by geographical barriers in Southeast Kazakhstan. We found that the alarm call variables were similar between the 3 study populations and differed by the maximum fundamental frequency (8.46 ± 0.75 kHz) from the values (5.62 ± 0.06 kHz) reported for the red-cheeked ground squirrels from the Kurgan region of Russia. Variation in mtDNA control region was only 3% and variation in cytochrome b gene was only 2.5%. Phylogenetic trees based on cytochrome b gene polymorphism of 44 individuals from the study area and adjacent territories indicated 3 clades with high (98-100%) bootstrap support: "intermedius," "brevicauda" and "iliensis"). We conclude that the 3 study populations in Southeast Kazakhstan belong to the clade intermedius and suggest a taxonomical revision of the species complex Spermophilus erythrogenys sensu lato, including analyses of nuclear DNA and alarm calls for populations of the brevicauda and iliensis clades.

Eyes on the size: accuracy of visual length estimates of white sharks, Carcharodon carcharias

Visual estimates have been used extensively to determine the length of large organisms that are logistically challenging to measure. However, there has been little effort to quantify the accuracy or validity of this technique despite inaccurate size estimates leading to incorrect population assessments and misinformed management strategies. Here, we compared visually estimated total length measurements of white sharks, Carcharodon carcharias, during cage-diving operations with measurements obtained from stereo-video cameras and assessed the accuracy of those estimates in relation to suspected biases (shark size, and observer experience and gender) using generalized linear mixed-models and linear regressions. Observer experience on board cage-diving vessels had the greatest effect on the accuracy of visual length estimates, with scientists being more accurate (mean accuracy ± standard error: 23.0 ± 16.5 cm) than crew (39.9 ± 33.8 cm) and passengers (49.4 ± 38.5 cm). Observer gender and shark size had no impact on the overall accuracy of visual length estimates, but passengers overestimated sharks less than 3 m and underestimated sharks greater than 3 m. Our findings show that experience measuring animals is the most substantial driver of accurate visual length estimates regardless of the amount of exposure to the species being measured. Scientists were most accurate, even though crew observe white sharks more frequently. Our results show that visual length estimates are not impacted by shark size and are a valid measurement tool for many aspects of C. carcharias research, provided they come from people who have previously been involved in measuring animals, i.e. scientists.

The rodent fauna from the Neolithic human settlement of Tepecik-Çiftlik (Niğde, Turkey)

The primary aim of this study was to establish for the very first time the results on the rodent fauna gathered from the ongoing Neolithic excavations at the Tepecik-Çiftlik site in southern Cappadocia (Niğde, Turkey). So far, the fauna being study is represented by species of arvicolines (Microtus cf. arvalis and Arvicola cf. amphibius), murines (Mus cf. musculus), sousliks (Spermophilus xanthoprymnus), hamsters (Mesocricetus brandti) and mole-rats (Spalax xanthodon), and their remains are described in detail. Paleoenvironmental assumption based on common preferences of the fauna elements indicates, at first sight, a rather dry steppe environment with sparse plant cover or perennial short grasses. On the other hand, the presence of Arvicola cf. amphibius points out streams and marsh-like vegetation cover with bodies of water . The predation or burn marks observed on some specimens and the presence of subterraneous rodents raise questions concerning their taphonomy.

Quantifying the three-dimensional facial morphology of the laboratory rat with a focus on the vibrissae

The morphology of an animal's face will have large effects on the sensory information it can acquire. Here we quantify the arrangement of cranial sensory structures of the rat, with special emphasis on the mystacial vibrissae (whiskers). Nearly all mammals have vibrissae, which are generally arranged in rows and columns across the face. The vibrissae serve a wide variety of important behavioral functions, including navigation, climbing, wake following, anemotaxis, and social interactions. To date, however, there are few studies that compare the morphology of vibrissal arrays across species, or that describe the arrangement of the vibrissae relative to other facial sensory structures. The few studies that do exist have exploited the whiskers' grid-like arrangement to quantify array morphology in terms of row and column identity. However, relying on whisker identity poses a challenge for comparative research because different species have different numbers and arrangements of whiskers. The present work introduces an approach to quantify vibrissal array morphology regardless of the number of rows and columns, and to quantify the array's location relative to other sensory structures. We use the three-dimensional locations of the whisker basepoints as fundamental parameters to generate equations describing the length, curvature, and orientation of each whisker. Results show that in the rat, whisker length varies exponentially across the array, and that a hard limit on intrinsic curvature constrains the whisker height-to-length ratio. Whiskers are oriented to "fan out" approximately equally in dorsal-ventral and rostral-caudal directions. Quantifying positions of the other sensory structures relative to the whisker basepoints shows remarkable alignment to the somatosensory cortical homunculus, an alignment that would not occur for other choices of coordinate systems (e.g., centered on the midpoint of the eyes). We anticipate that the quantification of facial sensory structures, including the vibrissae, will ultimately enable cross-species comparisons of multi-modal sensing volumes.

Determination of genetic variations in the genus Dryomys Thomas, 1906 (Rodentia: Gliridae) distributed in Turkey using NADH dehydrogenase 1 (ND1) gene

Genetic diversity and phylogeny of Dryomys nitedula and Dryomys laniger from Turkey was described in the present study by using mitochondrial DNA NADH dehydrogenase 1 gene (ND1). Genetic variation in ND1 gene was determined by two model-based phylogenetic analyses and a network analysis revealed 27 haplotypes of D. nitedula constructing four main lineages (Thrace, Anatolia, North-eastern Anatolia and Şavşat) that have non-overlapping geographic distributions and no shared haplotypes, but on the other hand, three haplotypes were detected in four samples of D. laniger from Turkey. It was determined that nucleotide diversity was low but haplotype diversity was high in D. laniger, whereas, D. nitedula has both high level of haplotype and nucleotide diversity. Characterization of Thrace lineage of D. nitedula with low nucleotide diversity and determination of the total nucleotide diversity of Anatolian lineages (Anatolia + North-eastern Anatolia+Şavşat) to be approximately four times higher than that of Thrace lineage indicated that Anatolia may have served as a refuge for D. nitedula. Divergence times and high level of nucleotide differences between D. nitedula lineages showed that diversification of the lineages may have occurred before and during ice ages in Turkey, thought to be a refuge for post-glacial colonization and biodiversity resource of Europe. Additionally, estimated divergence times and calculated genetic distances yielded compatible results with the previous paleontological and genomic data for the diversification time of two species in the genus.

Genetic characterization of the Turkish gray hamster (Cricetulus migratorius) based on mitochondrial cytochrome b and 12S rRNA sequences

Although genetic diversity and phylogenetic status of the gray hamster (Cricetulus migratorius) have been investigated from different regions in previous studies, genetic data on this species from Turkey are still lacking, since previous data have been based on a limited number of gray hamsters sampled across the Anatolian part of Turkey. The aim of this study was to determine the genetic diversity of the Anatolian population and to reveal the phylogenetic relationships among the Anatolian population and conspecific populations of the gray hamster. The complete and partial fragments of mitochondrial Cyt b and 12S rRNA from the 20 Turkish samples were amplified and sequenced. Ten 12S rRNA (901 bp) and 15 Cyt b (1140 bp) haplotypes found in this work were not previously reported. Based on Bayesian, Maximum Likelihood, Neighbour-Joining and Median-Joining network analyses by using mitochondrial data under the name Cricetulus, the results of phylogenetic and network analyses indicated that there was a deep separation among the distinct lineages within the genus Cricetulus. When considering the position of the Turkish haplotypes in median joining network, the Anatolian part of Turkey may have hosted a source population of the gray hamster for expansion to adjacent regions in the past period. Additionally, the Anatolian population of gray hamster had relatively high haplotype diversity and the present study propounded the importance of data obtained from the Anatolian population of gray hamster to reveal the phylogenetic relationships among conspecific populations of the gray hamster.

Pleistocene range shifts, refugia and the origin of widespread species in western Palaearctic water beetles

Quaternary glacial cycles drove major shifts in both the extent and location of the geographical ranges of many organisms. During glacial maxima, large areas of central and northern Europe were inhospitable to temperate species, and these areas are generally assumed to have been recolonized during interglacials by range expansions from Mediterranean refugia. An alternative is that this recolonization was from non-Mediterranean refugia, in central Europe or western Asia, but data on the origin of widespread central and north European species remain fragmentary, especially for insects. We studied three widely distributed lineages of freshwater beetles (the Platambus maculatus complex, the Hydraena gracilis complex, and the genus Oreodytes), all restricted to running waters and including both narrowly distributed southern endemics and widespread European species, some with distributions spanning the Palearctic. Our main goal was to determine the role of the Pleistocene glaciations in shaping the diversification and current distribution of these lineages. We sequenced four mitochondrial and two nuclear genes in populations drawn from across the ranges of these taxa, and used Bayesian probabilities and Maximum Likelihood to reconstruct their phylogenetic relationships, age and geographical origin. Our results suggest that all extant species in these groups are of Pleistocene origin. In the H. gracilis complex, the widespread European H. gracilis has experienced a rapid, recent range expansion from northern Anatolia, to occupy almost the whole of Europe. However, in the other two groups widespread central and northern European taxa appear to originate from central Asia, rather than the Mediterranean. These widespread species of eastern origin typically have peripherally isolated forms in the southern Mediterranean peninsulas, which may be remnants of earlier expansion-diversification cycles or result from incipient isolation of populations during the most recent Holocene expansion. The accumulation of narrow endemics of such lineages in the Mediterranean may result from successive cycles of range expansion, with subsequent speciation (and local extinction in glaciated areas) through multiple Pleistocene climatic cycles.

Implications of hybridization, NUMTs, and overlooked diversity for DNA Barcoding of Eurasian ground squirrels

The utility of DNA Barcoding for species identification and discovery has catalyzed a concerted effort to build the global reference library; however, many animal groups of economical or conservational importance remain poorly represented. This study aims to contribute DNA barcode records for all ground squirrel species (Xerinae, Sciuridae, Rodentia) inhabiting Eurasia and to test efficiency of this approach for species discrimination. Cytochrome c oxidase subunit 1 (COI) gene sequences were obtained for 97 individuals representing 16 ground squirrel species of which 12 were correctly identified. Taxonomic allocation of some specimens within four species was complicated by geographically restricted mtDNA introgression. Exclusion of individuals with introgressed mtDNA allowed reaching a 91.6% identification success rate. Significant COI divergence (3.5-4.4%) was observed within the most widespread ground squirrel species (Spermophilus erythrogenys, S. pygmaeus, S. suslicus, Urocitellus undulatus), suggesting the presence of cryptic species. A single putative NUMT (nuclear mitochondrial pseudogene) sequence was recovered during molecular analysis; mitochondrial COI from this sample was amplified following re-extraction of DNA. Our data show high discrimination ability of 100 bp COI fragments for Eurasian ground squirrels (84.3%) with no incorrect assessments, underscoring the potential utility of the existing reference librariy for the development of diagnostic 'mini-barcodes'.

Congruent morphological and genetic differentiation as a signature of range expansion in a fragmented landscape

Phenotypic differentiation is often interpreted as a result of local adaptation of individuals to their environment. Here, we investigated the skull morphological differentiation in 11 populations of the white-footed mouse (Peromyscus leucopus). These populations were sampled in an agricultural landscape in the Montérégie region (Québec, Canada), at the northern edge of the distribution of the white-footed mouse. We found a strong pattern of phenotypic differentiation matching the genetic structure across these populations. Landscape fragmentation and the presence of geographic barriers, in particular north-south oriented rivers, contribute to this differentiation and modulate the pattern of rapid ongoing northward range expansion of the white-footed mouse in response to climate warming. We conclude that while large rivers and postglacial recolonization routes have shaped the current pattern of distribution and differentiation of white-footed mouse populations, further local differentiation is occurring, at the scale of the landscape. We posit that the northern expansion of the white-footed mouse is achieved through successive independent founder events in a fragmented landscape at the northern range edge of the species. The phenotypic differentiation we observe is thus a result of a number of mechanisms operating at different spatial and temporal scales.