Phylogenetic and morphological significance of an overlooked flying squirrel (Pteromyini, Rodentia) from the eastern Himalayas with the description of a new genus.

The flying squirrels (Pteromyini, Rodentia) are the most diverse and widely distributed group of gliding mammals. Taxonomic boundaries and relationships within flying squirrels remain an area of active research in mammalogy. The discovery of new specimens of Pteromys ( Hylopetes) leonardi Thomas, 1921 previously considered a synonym of Hylopetes alboniger, in Yunnan Province, China allowed a morphological and genetic reassessment of the status of this taxon. Phylogenetic reconstruction was implemented using sequences of two mitochondrial (12S ribosomal DNA and 16S ribosomal DNA) and one nuclear (interphotoreceptor retinoid-binding protein) gene fragments. Morphological assessments involved examinations of features preserved on skins, skulls, and penises of museum specimens, supplemented with principal component analysis of craniometric data. Together these assessments revealed that this taxon should be recognized not only as a distinct species, and should also be placed within a new genus, described here as Priapomys.

Museomics of tree squirrels: a dense taxon sampling of mitogenomes reveals hidden diversity, phenotypic convergence, and the need of a taxonomic overhaul.

BACKGROUND: Tree squirrels (Sciuridae, Sciurini), in particular the highly diverse Neotropical lineages, are amongst the most rapidly diversifying branches of the mammal tree of life but also some of the least known. Negligence of this group by systematists is likely a product of the difficulties in assessing morphological informative traits and of the scarcity or unavailability of fresh tissue samples for DNA sequencing. The highly discrepant taxonomic arrangements are a consequence of the lack of phylogenies and the exclusive phenotypic-based classifications, which can be misleading in a group with conservative morphology. Here we used high-throughput sequencing and an unprecedented sampling of museum specimens to provide the first comprehensive phylogeny of tree squirrels, with a special emphasis on Neotropical taxa. RESULTS: We obtained complete or partial mitochondrial genomes from 232 historical and modern samples, representing 40 of the 43 currently recognized species of Sciurini. Our phylogenetic analyses-performed with datasets differing on levels of missing data and taxa under distinct analytical methods-strongly support the monophyly of Sciurini and consistently recovered 12 major clades within the tribe. We found evidence that the diversity of Neotropical tree squirrels is underestimated, with at least six lineages that represent taxa to be named or revalidated. Ancestral state reconstructions of number of upper premolars and number of mammae indicated that alternative conditions of both characters must have evolved multiple times throughout the evolutionary history of tree squirrels. CONCLUSIONS: Complete mitogenomes were obtained from museum specimens as old as 120 years, reinforcing the potential of historical samples for phylogenetic inferences of elusive lineages of the tree of life. None of the taxonomic arrangements ever proposed for tree squirrels fully corresponded to our phylogenetic reconstruction, with only a few of the currently recognized genera recovered as monophyletic. By investigating the evolution of two morphological traits widely employed in the taxonomy of the group, we revealed that their homoplastic nature can help explain the incongruence between phylogenetic results and the classification schemes presented so far. Based on our phylogenetic results we suggest a tentative supraspecific taxonomic arrangement for Sciurini, employing 13 generic names used in previous taxonomic classifications.

A Novel Retrovirus (Gunnison's Prairie Dog Retrovirus) Associated With Thymic Lymphoma in Gunnison's Prairie Dogs in Colorado, USA.

As part of research and wildlife disease surveillance efforts, we performed necropsy examinations of 125 free-ranging (n = 114) and captive (n = 11) prairie dogs in Colorado from 2009 to 2017. From these cases, we identified three cases of thymic lymphoma in free-ranging Gunnison's prairie dogs (Cynomys gunnisoni), and we identified a novel retroviral sequence associated with these tumors. The viral sequence is 7700 nucleotides in length and exhibits a genetic organization that is consistent with the characteristics of a type D betaretrovirus. The proposed name of this virus is Gunnison's prairie dog retrovirus (GPDRV). We screened all 125 prairie dogs for the presence of GPDRV using PCR with envelope-specific primers and DNA extracted from spleen samples. Samples were from Gunnison's prairie dogs (n = 59), black-tailed prairie dogs (Cynomys ludovicianus) (n = 40), and white-tailed prairie dogs (Cynomys leucurus) (n = 26). We identified GPDRV in a total of 7/125 (5.6%) samples including all three of the prairie dogs with thymic lymphoma, as well as spleen from an additional four Gunnison's prairie dogs with no tumors recognized at necropsy. None of the GPDRV-negative Gunnison's prairie dogs had thymic lymphomas. We also identified a related, apparently endogenous retroviral sequence in all prairie dog samples. These results suggest that GPDRV infection may lead to development of thymic lymphoma in Gunnison's prairie dogs.

Search for polyoma-, herpes-, and bornaviruses in squirrels of the family Sciuridae.

BACKGROUND: Squirrels (family Sciuridae) are globally distributed members of the order Rodentia with wildlife occurrence in indigenous and non-indigenous regions (as invasive species) and frequent presence in zoological gardens and other holdings. Multiple species introductions, strong inter-species competition as well as the recent discovery of a novel zoonotic bornavirus resulted in increased research interest on squirrel pathogens. Therefore we aimed to test a variety of squirrel species for representatives of three virus families. METHODS: Several species of the squirrel subfamilies Sciurinae, Callosciurinae and Xerinae were tested for the presence of polyomaviruses (PyVs; family Polyomaviridae) and herpesviruses (HVs; family Herpesviridae), using generic nested polymerase chain reaction (PCR) with specificity for the PyV VP1 gene and the HV DNA polymerase (DPOL) gene, respectively. Selected animals were tested for the presence of bornaviruses (family Bornaviridae), using both a broad-range orthobornavirus- and a variegated squirrel bornavirus 1 (VSBV-1)-specific reverse transcription-quantitative PCR (RT-qPCR). RESULTS: In addition to previously detected bornavirus RNA-positive squirrels no more animals tested positive in this study, but four novel PyVs, four novel betaherpesviruses (BHVs) and six novel gammaherpesviruses (GHVs) were identified. For three PyVs, complete genomes could be amplified with long-distance PCR (LD-PCR). Splice sites of the PyV genomes were predicted in silico for large T antigen, small T antigen, and VP2 coding sequences, and experimentally confirmed in Vero and NIH/3T3 cells. Attempts to extend the HV DPOL sequences in upstream direction resulted in contiguous sequences of around 3.3 kilobase pairs for one BHV and two GHVs. Phylogenetic analysis allocated the novel squirrel PyVs to the genera Alpha- and Betapolyomavirus, the BHVs to the genus Muromegalovirus, and the GHVs to the genera Rhadinovirus and Macavirus. CONCLUSIONS: This is the first report on molecular identification and sequence characterization of PyVs and HVs and the detection of bornavirus coinfections with PyVs or HVs in two squirrel species. Multiple detection of PyVs and HVs in certain squirrel species exclusively indicate their potential host association to a single squirrel species. The novel PyVs and HVs might serve for a better understanding of virus evolution in invading host species in the future.

Pharmacokinetic Profiles of Gabapentin after Oral and Subcutaneous Administration in Black-tailed Prairie Dogs (Cynomys ludovicianus).

In veterinary and human medicine, gabapentin (a chemical analog of γ-aminobutyric acid) is commonly prescribed to treat postoperative and chronic neuropathic pain. This study explored the pharmacokinetics of oral and subcutaneous administration of gabapentin at high (80 mg/kg) and low (30 mg/kg) doses as a potential analgesic in black-tailed prairie dogs (Cynomys ludovicianus; n = 24). The doses (30 and 80 mg/kg) and half maximal effective concentration (1.4 to 16.7 ng/mL) for this study were extrapolated from pharmacokinetic efficacy studies in rats, rabbits, and cats. Gabapentin in plasma was measured by using an immunoassay, and data were evaluated using noncompartmental analysis. The peak plasma concentrations (mean ±1 SD) were 42.6 ±14.8 and 115.5 ±15.2 ng/mL, respectively, after 30 and 80 mg/kg SC and 14.5 ±3.5 and 20.7 ±6.1 ng/mL after the low and high oral dosages, respectively. All peak plasma concentrations of gabapentin occurred within 5 h of administration. Disappearance half-lives for the low and high oral doses were 7.4 ± 6.0 h and 5.0 ± 0.8 h, respectively. The results of this study demonstrate that oral administration of gabapentin at low (30 mg/kg) doses likely would achieve and maintain plasma concentrations at half maximum effective concentration for 12 h, making it a viable option for an every 12-h treatment.

Gliding Dragons and Flying Squirrels: Diversifying versus Stabilizing Selection on Morphology following the Evolution of an Innovation.

Evolutionary innovations and ecological competition are factors often cited as drivers of adaptive diversification. Yet many innovations result in stabilizing rather than diversifying selection on morphology, and morphological disparity among coexisting species can reflect competitive exclusion (species sorting) rather than sympatric adaptive divergence (character displacement). We studied the innovation of gliding in dragons (Agamidae) and squirrels (Sciuridae) and its effect on subsequent body size diversification. We found that gliding either had no impact (squirrels) or resulted in strong stabilizing selection on body size (dragons). Despite this constraining effect in dragons, sympatric gliders exhibit greater size disparity compared with allopatric gliders, a pattern consistent with, although not exclusively explained by, ecological competition changing the adaptive landscape of body size evolution to induce character displacement. These results show that innovations do not necessarily instigate further differentiation among species, as is so often assumed, and suggest that competition can be a powerful force generating morphological divergence among coexisting species, even in the face of strong stabilizing selection.

Exclusion of interspecific competition reduces scatter-hoarding of Siberian chipmunk Tamias sibiricus: A field study.

Although food availability and the abundance of seed predators have been postulated to affect seed dispersal, it is not clear how seed-eating animals modify their scatter-hoarding strategies in response to different levels of interspecific competition. We placed paired germinated and ungerminated acorns of Quercus mongolica on 30-cm high platforms to exclude potential interspecific competition of the predominant larder hoarders Apodemus peninsulae and Myodes rufocanus, to investigate seed dispersal by a predominant scatter-hoarder, Tamias sibiricus, in the field in north-eastern China. Our results showed that T. sibiricus ate more acorns in situ in the absence of interspecific competition. In the presence of interspecific competition of A. peninsulae and C. rufocanus, however, more acorns were scatter-hoarded by T. sibiricus. Regardless of interspecific competition, germination of acorns showed no significant effects on seed dispersal patterns, inconsistent with the "seed perishability hypothesis" that animals avoid hoarding seeds with high perishability. Exclusion of interspecific competition, though relatively increasing the per capita seed abundance, appears to reduce seed dispersal, scatter-hoarding and seedling establishment. Therefore, we propose that moderate interspecific competition rather than competition exclusion may benefit seed scatter-hoarding and seedling establishment.

No evidence for phylosymbiosis in western chipmunk species.

Phylosymbiosis refers to a congruent pattern between the similarity of microbiomes of different species and the branching pattern of the host phylogeny. Phylosymbiosis has been detected in a variety of vertebrate and invertebrate hosts, but has only been assessed in geographically isolated populations. We tested for phylosymbiosis in eight (sub)species of western chipmunks with overlapping ranges and ecological niches; we used a nuclear (Acrosin) and a mitochondrial (CYTB) phylogenetic marker because there are many instances of mitochondrial introgression in chipmunks. We predicted that similarity among microbiomes increases with: (1) increasing host mitochondrial relatedness, (2) increasing host nuclear genome relatedness and (3) decreasing geographic distance among hosts. We did not find statistical evidence supporting phylosymbiosis in western chipmunks. Furthermore, in contrast to studies of other mammalian microbiomes, similarity of chipmunk microbiomes is not predominantly determined by host species. Sampling site explained most variation in microbiome composition, indicating an important role of local environment in shaping microbiomes. Fecal microbiomes of chipmunks were dominated by Bacteroidetes (72.2%), followed by Firmicutes (24.5%), which is one of the highest abundances of Bacteroidetes detected in wild mammals. Future work will need to elucidate the effects of habitat, ecology and host genomics on chipmunk microbiomes.

Temporal genomic contrasts reveal rapid evolutionary responses in an alpine mammal during recent climate change.

Many species have experienced dramatic changes in their abundance and distribution during recent climate change, but it is often unclear whether such ecological responses are accompanied by evolutionary change. We used targeted exon sequencing of 294 museum specimens (160 historic, 134 modern) to generate independent temporal genomic contrasts spanning a century of climate change (1911-2012) for two co-distributed chipmunk species: an endemic alpine specialist (Tamias alpinus) undergoing severe range contraction and a stable mid-elevation species (T. speciosus). Using a novel analytical approach, we reconstructed the demographic histories of these populations and tested for evidence of recent positive directional selection. Only the retracting species showed substantial population genetic fragmentation through time and this was coupled with positive selection and substantial shifts in allele frequencies at a gene, Alox15, involved in regulation of inflammation and response to hypoxia. However, these rapid population and gene-level responses were not detected in an analogous temporal contrast from another area where T. alpinus has also undergone severe range contraction. Collectively, these results highlight that evolutionary responses may be variable and context dependent across populations, even when they show seemingly synchronous ecological shifts. Our results demonstrate that temporal genomic contrasts can be used to detect very recent evolutionary responses within and among contemporary populations, even in the face of complex demographic changes. Given the wealth of specimens archived in natural history museums, comparative analyses of temporal population genomic data have the potential to improve our understanding of recent and ongoing evolutionary responses to rapidly changing environments.

Alternative treatments of genetic distances for species delimitation in Callosciurinae and Sciurinae (Rodentia: Sciuridae).

There are many DNA taxonomy methods available, but single-locus data coupled with distance-based analyses still dominate species identification in practice. Genetic distance thresholds are often used for assigning genetic diversity into entities corresponding to species but are criticized for non-universality or arbitrary. This study investigated an alternative approach for determining genetic thresholds that used iteratively sister-clade clustering for sister species boundary detection within a phylogenetic framework. This method was separately applied to two close subfamilies of Callosciurinae and Sciurinae in Rodentia using the cyt b gene. Our results showed that genetic thresholds for Callosciurinae and Sciurinae were 4.0% and 6.0%, respectively, indicating that the optimized thresholds could be lineage-specific. The use of these thresholds for taxa partitioning yielded a very similar putative species number as the prevailing ABGD method, and increased species diversity by 74.2% and 20.7% in Callosciurinae and Sciurinae, respectively. This suggested that additional cryptic species were present that warrant further investigation. We further tested the performance limitation of our method by simulating the impacts of tree construction and sampling limitation. The results showed that it preformed equally well for different trees but failed to work when inter- or intraspecific sampling is insufficient. These findings support the feasibility of this approach as an alternative tool for species delineation when only single-locus information is available for large datasets.

Impacts of Inference Method and Data set Filtering on Phylogenomic Resolution in a Rapid Radiation of Ground Squirrels (Xerinae: Marmotini).

Phylogenomic data sets are illuminating many areas of the Tree of Life. However, the large size of these data sets alone may be insufficient to resolve problematic nodes in the most rapid evolutionary radiations, because inferences in zones of extraordinarily low phylogenetic signal can be sensitive to the model and method of inference, as well as the information content of loci employed. We used a data set of $>$3950 ultraconserved element (UCE) loci from a classic mammalian radiation, ground-dwelling squirrels of the tribe Marmotini (Sciuridae: Xerinae), to assess sensitivity of phylogenetic estimates to varying per-locus information content across four different inference methods (RAxML, ASTRAL, NJst, and SVDquartets). Persistent discordance was found in topology and bootstrap support between concatenation- and coalescent-based inferences; among methods within the coalescent framework; and within all methods in response to different filtering scenarios. Contrary to some recent empirical UCE-based studies, filtering by information content did not promote complete among-method concordance. Nevertheless, filtering did improve concordance relative to randomly selected locus sets, largely via improved consistency of two-step summary methods (particularly NJst) under conditions of higher average per-locus variation (and thus increasing gene tree precision). The benefits of phylogenomic data set filtering are variable among classes of inference methods and across different evolutionary scenarios, reiterating the complexities of resolving rapid radiations, even with robust taxon and character sampling.

Oldest skeleton of a fossil flying squirrel casts new light on the phylogeny of the group.

Flying squirrels are the only group of gliding mammals with a remarkable diversity and wide geographical range. However, their evolutionary story is not well known. Thus far, identification of extinct flying squirrels has been exclusively based on dental features, which, contrary to certain postcranial characters, are not unique to them. Therefore, fossils attributed to this clade may indeed belong to other squirrel groups. Here we report the oldest fossil skeleton of a flying squirrel (11.6 Ma) that displays the gliding-related diagnostic features shared by extant forms and allows for a recalibration of the divergence time between tree and flying squirrels. Our phylogenetic analyses combining morphological and molecular data generally support older dates than previous molecular estimates (~23 Ma), being congruent with the inclusion of some of the earliest fossils (~36 Ma) into this clade. They also show that flying squirrels experienced little morphological change for almost 12 million years.

Convergent evolution in the Euarchontoglires.

Convergence-the independent evolution of similar phenotypes in distantly related clades-is a widespread and much-studied phenomenon. An often-cited, but hitherto untested, case of morphological convergence is that between the aye-aye and squirrels. The aye-aye (Daubentonia madagascariensis) is a highly unusual lemuriform primate that has evolved a dentition similar to that of rodents: it possesses large, ever-growing incisors which it uses to strip the bark from trees in order to feed on wood-boring beetle larvae. Indeed, such is the similarity that some of the earliest classifications of the aye-aye placed it in the squirrel genus Sciurus Here, we aimed to test the degree of convergence between the skulls and lower jaws of squirrels and the aye-aye. Three-dimensional landmarks were recorded from the crania and mandibles of 46 taxa representing the majority of families in the Euarchontoglires. Results were plotted as phylomorphospaces and convergence measures were calculated. The convergence between squirrels and the aye-aye was shown to be statistically significant for both the cranium and mandible, although the mandibles seem to converge more closely in shape. The convergence may indicate strong functional drivers of morphology in these taxa, i.e. the use of the incisors to produce high bite forces during feeding. Overall, we have shown that this classic case of convergence stands up to quantitative analysis.

The Taxonomic Status of Spermophilus in the Plague Area of Dingbian County, Shaanxi Province, China.

This study was conducted to define the taxonomic status of Spermophilus in the plague area of Dingbian County in Shaanxi Province, China, through the two-factor variance analysis of morphological characteristics, DNA barcoding, and chromosome karyotype analysis. The Spermophilus samples collected from Dingbian and Zhengxiang Baiqi Counties exhibited significant differences in their morphological measurements. All Spermophilus samples form two distinct branches in neighbor-joining (NJ) tree. One branch included the Spermophilus samples collected from Inner Mongolia, and the other branch included samples collected from the plague foci of Shaanxi Province and the Ningxia Region. The Spermophilus samples collected from Dingbian County had a chromosome number of 2n = 38 in 84.40% of all their cells. The Spermophilus species collected from the plague area of Dingbian County was categorized as Spermophilus alashanicus (S.alashamicus). The findings reported in this study are epidemiologically significant for monitoring plague in this region of west-central China.

Impacts of late Quaternary environmental change on the long-tailed ground squirrel (Urocitellus undulatus) in Mongolia.

Impacts of Quaternary environmental changes on mammal faunas of central Asia remain poorly understood due to a lack of geographically comprehensive phylogeographic sampling for most species. To help address this knowledge gap, we conducted the most extensive molecular analysis to date of the long-tailed ground squirrel (Urocitellus undulatus Pallas 1778) in Mongolia, a country that comprises the southern core of this species' range. Drawing on material from recent collaborative field expeditions, we genotyped 128 individuals at 2 mitochondrial genes (cytochrome b and cytochrome oxidase I; 1 797 bp total). Phylogenetic inference supports the existence of two deeply divergent infraspecific lineages (corresponding to subspecies U. u. undulatus and U. u. eversmanni), a result in agreement with previous molecular investigations but discordant with patterns of range-wide craniometric and external phenotypic variation. In the widespread westerneversmanni lineage, we recovered geographically-associated clades from the: (a) Khangai, (b) Mongolian Altai, and (c) Govi Altai mountain ranges. Phylogeographic structure in U. u. eversmanni is consistent with an isolation-by-distance model; however, genetic distances are significantly lower than among subspecies, and intra-clade relationships are largely unresolved. The latter patterns, as well as the relatively higher nucleotide polymorphism of populations from the Great Lakes Depression of northwestern Mongolia, suggest a history of range shifts into these lowland areas in response to Pleistocene glaciation and environmental change, followed by upslope movements and mitochondrial lineage sorting with Holocene aridification. Our study illuminates possible historical mechanisms responsible for U. undulatus genetic structure and contributes to a framework for ongoing exploration of mammalian response to past and present climate change in central Asia.

Distribution and prevalence of vector-borne diseases in California chipmunks (Tamias spp.).

California, with 13 chipmunk (Tamias) species, has more than any other state or country, occupying habitats ranging from chaparral to the high peaks of the Sierra Nevada. Chipmunks host zoonotic pathogens including Yersinia pestis, Anaplasma phagocytophilum, relapsing fever (RF) Borrelia spp., Borrelia burgdorferi, and spotted fever group (SFG) Rickettsia species. Chipmunk species are often not differentiated by public health workers, yet different species utilize different ecological niches and may have intrinsically different capacities for maintaining vector-borne pathogens and infecting vectors. We surveyed over 700 individuals from nine species of chipmunks throughout California for exposure to and infection by Y. pestis, A. phagocytophilum, RF Borrelia spp., Borrelia burgdorferi, and SFG Rickettsia species. DNA of all five pathogens was found and all chipmunks except Merriam's chipmunk (T. merriami) were PCR-positive for at least one of the pathogens. Anaplasma phagocytophilum was most common (40.0%, 2/5) in Sonoma chipmunks (T. sonomae) from Marin county and B. burgdorferi most common (37.5%, 27/72) in redwood chipmunks (T. ochrogenys) from Mendocino county. RF Borrelia spp. was detected in 2% (6/297) of redwood chipmunks in Mendocino county and 10% (1/10) of both least (T. minimus) and lodgepole (T. speciosus) chipmunks in the western Sierra. Exposure to SFG Rickettsia spp. was found in the Northern Coastal region (Del Norte, Humboldt and Mendocino counties) and in the northern and western Sierra in several species of chipmunks. Y. pestis infection was found only in the western Sierra-in a yellow-pine (T. amoenus) and a long-eared (T. quadrimaculatus) chipmunk. Though more data are needed to thoroughly understand the roles that different chipmunk species play in disease transmission, our findings suggest that some chipmunk species may be more important to the maintenance of vector-borne diseases than others within each geographic area.

DNA analysis of a 30,000-year-old Urocitellus glacialis from northeastern Siberia reveals phylogenetic relationships between ancient and present-day arctic ground squirrels.

In contrast to the abundant fossil record of arctic ground squirrels, Urocitellus parryii, from eastern Beringia, only a limited number of fossils is known from its western part. In 1946, unnamed GULAG prisoners discovered a nest with three mummified carcasses of arctic ground squirrels in the permafrost sediments of the El'ga river, Yakutia, Russia, that were later attributed to a new species, Citellus (Urocitellus) glacialis Vinogr. To verify this assignment and to explore phylogenetic relationships between ancient and present-day arctic ground squirrels, we performed 14C dating and ancient DNA analyses of one of the El'ga mummies and four contemporaneous fossils from Duvanny Yar, northeastern Yakutia. Phylogenetic reconstructions, based on complete cytochrome b gene sequences of five Late Pleistocene arctic ground squirrels and those of modern U. parryii from 21 locations across western Beringia, provided no support for earlier proposals that ancient arctic ground squirrels from Siberia constitute a distinct species. In fact, we observed genetic continuity of the glacialis mitochondrial DNA lineage in modern U. parryii of the Kamchatka peninsula. When viewed in a broader geographic perspective, our findings provide new insights into the genetic history of U. parryii in Late Pleistocene Beringia.

Comparative Phylogenomic Assessment of Mitochondrial Introgression among Several Species of Chipmunks (Tamias).

Many species are not completely reproductively isolated, resulting in hybridization and genetic introgression. Organellar genomes, such as those derived from mitochondria (mtDNA) and chloroplasts, introgress frequently in natural systems; however, the forces shaping patterns of introgression are not always clear. Here, we investigate extensive mtDNA introgression in western chipmunks, focusing on species in the Tamias quadrivittatus group from the central and southern Rocky Mountains. Specifically, we investigate the role of selection in driving patterns of introgression. We sequenced 51 mtDNA genomes from six species and combine these sequences with other published genomic data to yield annotated mitochondrial reference genomes for nine species of chipmunks. Genomic characterization was performed using a series of molecular evolutionary and phylogenetic analyses to test protein-coding genes for positive selection. We fit a series of maximum likelihood models using a model-averaging approach, assessed deviations from neutral expectations, and performed additional tests to search for codons under the influence of selection. We found no evidence for positive selection among these genomes, suggesting that selection has not been the driving force of introgression in these species. Thus, extensive mtDNA introgression among several species of chipmunks likely reflects genetic drift of introgressed alleles in historically fluctuating populations.

Directional selection effects on patterns of phenotypic (co)variation in wild populations.

Phenotypic (co)variation is a prerequisite for evolutionary change, and understanding how (co)variation evolves is of crucial importance to the biological sciences. Theoretical models predict that under directional selection, phenotypic (co)variation should evolve in step with the underlying adaptive landscape, increasing the degree of correlation among co-selected traits as well as the amount of genetic variance in the direction of selection. Whether either of these outcomes occurs in natural populations is an open question and thus an important gap in evolutionary theory. Here, we documented changes in the phenotypic (co)variation structure in two separate natural populations in each of two chipmunk species (Tamias alpinus and T. speciosus) undergoing directional selection. In populations where selection was strongest (those of T. alpinus), we observed changes, at least for one population, in phenotypic (co)variation that matched theoretical expectations, namely an increase of both phenotypic integration and (co)variance in the direction of selection and a re-alignment of the major axis of variation with the selection gradient.

Rapid divergence and gene flow at high latitudes shape the history of Holarctic ground squirrels (Urocitellus).

Across the animal tree of life, the prevalence and evolutionary role(s) of hybridization remain incompletely understood. Rapidly radiating clades can serve as important systems for investigating these issues; however, such groups are often characterized by additional, widespread sources of gene tree discordance (e.g., incomplete lineage sorting). In this paper, we employed a multilocus dataset, Bayesian gene tree inference, and multiple species tree reconstruction methods to infer phylogeny of Holarctic ground squirrels (Urocitellus). We tested phylogenetic hypotheses based on previous morphological, cytological and single-locus datasets, and began to parse the causes of pervasive gene tree discordance that was observed. There is widespread incomplete lineage sorting in Urocitellus, consistent with rapid diversification embedded within the larger radiation of marmotine ground squirrels. We also recovered strong support for 2 instances of mitonuclear discord due to ancient hybridization among members of the high-latitude parryii-richardsonii-elegans clade. These results add to a growing number of documented hybridization events in ground squirrels, suggesting their radiation is a fertile system for understanding the interplay of diversification and hybridization in animal evolution.