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.

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.

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.

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.

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.

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.

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.

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.

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'.

Mesoamerican tree squirrels evolution (Rodentia: Sciuridae): a molecular phylogenetic analysis

The tribe Sciurini comprehends the genera Sciurus, Syntheosiurus, Microsciurus, Tamiasciurus and Rheinthrosciurus. The phylogenetic relationships within Sciurus have been only partially done, and the relationship between Mesoamerican species remains unsolved. The phylogenetic relationships of the Mesoamerican tree squirrels were examined using molecular data. Sequence data publicly available (12S, 16S, CYTB mitochondrial genes and IRBP nuclear gene) and cytochrome B gene sequences of four previously not sampled Mesoamerican Sciurus species were analyzed under a Bayesian multispecies coalescence model. Phylogenetic analysis of the multilocus data set showed the neotropical tree squirrels as a monophyletic clade. The genus Sciurus was paraphyletic due to the inclusion of Microsciurus species (M. alfari and M. flaviventer). The South American species S. aestuans and S. stramineus showed a sister taxa relationship. Single locus analysis based on the most compact and complete data set (i.e. CYTB gene sequences), supported the monophyly of the South American species and recovered a Mesoamerican clade including S. aureogaster, S. granatensis and S. variegatoides. These results corroborated previous findings based on cladistic analysis of cranial and post-cranial characters. Our data support a close relationship between Mesoamerican Sciurus species and a sister relationship with South American species, and corroborates previous findings in relation to the polyphyly of Microsciurus and Syntheosciurus paraphyly.

Differentiation in cranial variables among six species of Hylopetes (Sciurinae : Pteromyini)

There is some discrepancy in the classification of different species of Hylopetes, particularly regarding systematic status of H. electilis and H. phayrei and their relationship to other species. In the present study, for the first time we have brought together six of the nine Hylopetes species and performed statistical analysis of 14 measurable cranial variables, analyzing in total 89 specimens, including H. electilis, H. alboniger, H. phayrei, H. lepidus, H. spadiceus, and H. nigripes. Both univariate and multivariate analysis results indicate that H. electilis can not only be obviously distinguished from H. phayrei, but also clearly differs from the other four Hylopetes species. These results sustain the contention that H. electilis is neither a synonym nor subspecies of H. phayrei,but should be considered a distinct and valid species. Subsequently, a straightforward discussion on the biogeography of Hylopetes in southeastern Asia gives further insight into the differentiation and variety of species belonging to this genus. There is some discrepancy in the classification of different species of Hylopetes, particularly regarding systematic status of H. electilis and H. phayrei and their relationship to other species. In the present study, for the first time we have brought together six of the nine Hylopetes speciesand performed statistical analysis of 14 measurable cranial variables, analyzing in total 89 specimens, including H. electilis, H. alboniger, H. phayrei, H. lepidus, H. spadiceus, and H. nigripes. Both univariate and multivariate analysis results indicate that H. electilis can not only be obviously distinguished from H. phayrei, but also clearly differs from the other four Hylopetes species. These results sustain the contention that H. electilis is neither a synonym nor subspecies of H. phayrei,but should be considered a distinct and valid species. Subsequently, a straightforward discussion on the biogeography of Hylopetes in southeastern Asia gives further insight into the differentiation and variety of species belonging to this genus.

Conservatism and adaptability during squirrel radiation: what is mandible shape telling us?

Both functional adaptation and phylogeny shape the morphology of taxa within clades. Herein we explore these two factors in an integrated way by analyzing shape and size variation in the mandible of extant squirrels using landmark-based geometric morphometrics in combination with a comparative phylogenetic analysis. Dietary specialization and locomotion were found to be reliable predictors of mandible shape, with the prediction by locomotion probably reflecting the underlying diet. In addition a weak but significant allometric effect could be demonstrated. Our results found a strong phylogenetic signal in the family as a whole as well as in the main clades, which is in agreement with the general notion of squirrels being a conservative group. This fact does not preclude functional explanations for mandible shape, but rather indicates that ancient adaptations kept a prominent role, with most genera having diverged little from their ancestral clade morphologies. Nevertheless, certain groups have evolved conspicuous adaptations that allow them to specialize on unique dietary resources. Such adaptations mostly occurred in the Callosciurinae and probably reflect their radiation into the numerous ecological niches of the tropical and subtropical forests of Southeastern Asia. Our dietary reconstruction for the oldest known fossil squirrels (Eocene, 36 million years ago) show a specialization on nuts and seeds, implying that the development from protrogomorphous to sciuromorphous skulls was not necessarily related to a change in diet.

Diversity, seasonality, and context of mammalian roadkills in the southern Great Plains

Thousands of mammals are killed annually from vehicle collisions, making the issue an important one for conservation biologists and environmental managers. We recorded all readily identifiable kills on or immediately adjacent to roads in the southern Great Plains from March 2004-March 2007. We also recorded distance traveled, whether a road was paved or divided, the number of lanes, and prevailing habitat. Surveys were opportunistic and were conducted by car during conditions of good visibility. Over our 239 surveys and > 16,500 km traveled, we recorded 1412 roadkills from 18 different mammal species (size ranged from Sciurus squirrels to the white-tailed deer, Odocolieus virginianus). The overall kill rate was 8.50/100 km. Four species were prone to collisions: the Virginia opossum (Didelphis virginiana), nine-banded armadillo (Dasypus novemcinctus), striped skunk (Mephitis mephitis), and northern raccoon (Procyon lotor). Together they accounted for approximately 85% (1198) of all roadkills. Mortality rate differed significantly between 2- and 4-lane roads (8.39 versus 7.79/100 km). Kill rates were significantly higher on paved versus unpaved roads (8.60 versus 3.65/100 km), but did not depend on whether a road was divided. Roadkills were higher in spring than in fall (1.5 x), winter (1.4x), or summer (1.3x). The spring peak (in kills/100 km) was driven chiefly by the armadillo (2.76 in spring/summer versus 0.73 in autumn/winter) and opossum (2.65 versus 1.47). By contrast, seasonality was dampened by a late winter/early spring peak in skunk mortalities, for which 41% occurred in the 6-week period of mid-February through March. The raccoon did not exhibit a strong seasonal pattern. Our data are consistent with dispersal patterns of these species. Our results underscore the high rate of highway mortality in the southern plains, as well as differences in seasonality and road type that contribute to mortality. Conservation and management efforts should focus on creating underpasses or using other means to reduce roadkill rates.

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

This study reports on the molecular phylogenetics of ground squirrels, genus Spermophilus, in Turkey using cytochrome b (1140bp), part of the D-loop and flanking tRNAs (572bp), X chromosome (867-1051bp) and Y chromosome (983-989bp) DNA sequences. Individuals also were characterized by karyotype and with geometric morphometric analyses of mandibles and skulls. Two hundred fourteen individuals from 91 localities were studied. All the data support the recognition of a new species in SW Anatolia: the Taurus ground squirrel Spermophilus taurensis sp. nov. The new species has a small distribution in the Taurus Mountains in an area that is a hotspot for biodiversity. Molecular clock analysis suggests that the new species diverged from the European ground squirrel, Spermophilus citellus, about 2.5 million years ago and that the ancestor of these two species diverged from the widespread Anatolian ground squirrel, Spermophilus xanthoprymnus, about 5 million years ago. Morphometric differentiation in skull and mandible shape among the three species is incomplete, but statistically significant. S. xanthoprymnus is subdivided into five cytochrome b phylogroups and we use these data to infer the location of glacial refugia where the species lived during the last glacial maximum. This study illustrates the potential of combined molecular and morphometric studies to uncover new Anatolian species and to reconstruct their phylogeographic history. The new species is important for squirrel taxonomy and for understanding Eurasian mammal evolution.

Phylogeny and biogeography of the Petaurista philippensis complex (Rodentia: Sciuridae), inter- and intraspecific relationships inferred from molecular and morphometric analysis

With modified DNA extraction and purification protocols, the complete cytochrome b gene sequences (1140 bp) were determined from degraded museum specimens. Molecular analysis and morphological examination of cranial characteristics of the giant flying squirrels of Petaurista philippensis complex (P. grandis, P. hainana, and P. yunanensis) and other Petaurista species yielded new insights into long-standing controversies in the Petaurista systematics. Patterns of genetic variations and morphological differences observed in this study indicate that P. hainana, P. albiventer, and P. yunanensis can be recognized as distinct species, and P. grandis and P. petaurista are conspecific populations. Phylogenetic relationships reconstructed by using parsimony, likelihood, and Bayesian methods reveal that, with P. leucogenys as the basal branch, all Petaurista groups formed two distinct clades. Petaurista philippensis, P. hainana, P. yunanensis, and P. albiventer are clustered in the same clade, while P. grandis shows a close relationship to P. petaurista. Deduced divergence times based on Bayesian analysis and the transversional substitution at the third codon suggest that the retreating of glaciers and upheavals or movements of tectonic plates in the Pliocene-Pleistocene were the major factors responsible for the present geographical distributions of Petaurista groups.

Intraspecific Phylogeography of Red Squirrels (Tamiasciurus hudsonicus) in the Central Rocky Mountain Region of North America

We used variation in a portion of the mitochondrial DNA control region to examine phylogeography of Tamiasciurus hudsonicus, a boreal-adapted small mammal in the central Rocky Mountain region. AMOVA revealed that 65.66% of genetic diversity was attributable to variation within populations, 16.93% to variation among populations on different mountain ranges, and 17.41% to variation among populations within mountain ranges. Nested clade analysis revealed two major clades that likely diverged in allopatry during the Pleistocene: a southern clade from southern Colorado and a northern clade comprising northern Colorado, Wyoming, eastern Utah, and eastern Idaho. Historically restricted gene flow as a result of geographic barriers was indicated between populations on opposite sides of the Green River and Wyoming Basin and among populations in eastern Wyoming. In some instances genetic structure indicated isolation by distance.

The effects of Cenozoic global change on squirrel phylogeny

By modifying habitats and creating bridges and barriers between landmasses, climate change and tectonic events are believed to have important consequences for diversification of terrestrial organisms. Such consequences should be most evident in phylogenetic histories of groups that are ancient, widespread, and diverse. The squirrel family (Sciuridae) is one of very few mammalian families endemic to Eurasia, Africa, and North and South America and is ideal for examining these issues. Through phylogenetic and molecular-clock analyses, we infer that arrival and diversification of squirrels in Africa, on Sunda Shelf islands, across Beringea, and across the Panamanian isthmus coincide in timing and location with multiple well-documented sea-level, tectonic, and paleontological events. These precise correspondences point to an important role for global change in the diversification of a major group of mammals.

Extensive mtDNA variation within the yellow-pine chipmunk, Tamias amoenus (Rodentia: Sciuridae), and phylogeographic inferences for northwest North America

The yellow-pine chipmunk, Tamias amoenus, is common in xerophytic forests throughout much of northwest North America. We analyzed cytochrome b sequence variation from 155 individuals representing 57 localities across the distribution of T. amoenus including 10 additional species of Tamias. Maximum likelihood and parsimony tree estimation methods were used in conjunction with nested clade analysis to infer both deep and population-level processes. Our results indicate that two currently recognized subspecies of T. amoenus (T. a. canicaudus and T. a. cratericus) are not nested within other samples of T. amoenus. Maximum uncorrected levels of intraspecific sequence divergence within remaining samples of T. amoenus are >7%. Substantial geographic variation is characterized by 12 well-supported clades that correspond to distinct mountain ranges, but do not necessarily follow existing subspecific taxonomy. Significant association between geography and genealogy was detected within many of these clades and can be attributed to different population-level processes including past fragmentation, recent range expansion, and isolation by distance.

Choline acetyltransferase is expressed by non-starburst amacrine cells in the ground squirrel retina

We have used immunostaining techniques to reveal a new type of amacrine cell that is immunoreactive for choline acetyltransferase (ChAT), the acetylcholine synthesizing enzyme, in the Ground Squirrel (Spermophilus beecheyi) retina. Cryostat sections and double immunostained wholemount preparations were examined by confocal microscopy. This new ChAT type III cell is distinct in morphology and neurotransmitter content from the well know 'starburst' amacrine cells (types I and II) that are so well represented in the ground squirrel retina [J. Comp. Neurol. 365 (1996) 173-216]. The type III cell colocalizes glycine with the acetylcholine and does not appear to be GABAergic or exhibit calcium-binding proteins like the well-known starburst type. As well, type III cells do not occur as a mirror-symmetric pair with normally placed and displaced varieties. The type III cell is probably a small field amacrine type branching broadly in upper sublamina b of the inner plexiform layer, and is most likely A6 of the Ground Squirrel retina [J. Comp. Neurol. 365 (1996) 173-216]. Type III cells are ideally placed in the architecture of the Ground Squirrel retina to influence ON directionally selective ganglion cell types.

Parapoxvirus causes a deleterious disease in red squirrels associated with UK population declines

The disease implications of novel pathogens need to be considered when investigating the ecological impact of species translocations on native fauna. Traditional explanations based on competition or predation may often not be the whole story. Evidence suggests that an emerging infectious disease, caused by a parapoxvirus, may be a significant component of the impact that the introduced grey squirrel has had on UK red squirrel populations. Here we validate the potential role of parapoxvirus by proving that the virus is highly pathogenic in the red squirrel while having no detectable effect on grey squirrel health.

Phylogeography of the red-tailed chipmunk (Tamias ruficaudus), a northern Rocky Mountain endemic

The northern Rocky Mountains have experienced a complex history of geological events and environmental fluctuation, including Pleistocene glaciation. To provide an initial assessment of the genetic impact of this history on the regional biota we estimated phylogenetic relationships within Tamias ruficaudus, a regional endemic, from cytochrome b sequence variation using parsimony, maximum likelihood, and nested clade analysis. Analyses of sequence variation in 187 individuals from 43 localities across the distribution of T. ruficaudus indicate a history of vicariance events and range fluctuation consistent with successive periods of extensive Pleistocene glaciation in the northern Rocky Mountains. Intraspecific divergence levels (c. 4.7% uncorrected) and phylogenetic structure are consistent with a genealogical vicariance initiated prior to the Late Pleistocene, whereas nested clade analyses indicate more recent population history structured by both fragmentation and range expansion. A comparison of sequence variation with bacular morphology indicates that the two genetically and morphologically differentiated entities exhibit a zone of differential character introgression. Sequence data support a multiple refugia hypothesis and provide a phylogeographical case study for the ongoing synthesis of regional biogeography for northern Rocky Mountain endemics.

Molecular phylogeny of the chipmunks inferred from Mitochondrial cytochrome b and cytochrome oxidase II gene sequences

There are currently 25 recognized species of the chipmunk genus Tamias. In this study we sequenced the complete mitochondrial cytochrome b (cyt b) gene of 23 Tamias species. We analyzed the cyt b sequence and then analyzed a combined data set of cyt b along with a previous data set of cytochrome oxidase subunit II (COII) sequence. Maximum-likelihood was used to further test the fit of models of evolution to the cyt b data. Other sciurid cyt b sequence was added to examine the evolution of Tamias in the context of other sciurids. Relationships among Tamias species are discussed, particularly the possibility of a current sorting event among taxa of the southwestern United States and the extreme divergences among the three subgenera (Neotamias, Eutamias, and Tamias).

[Genetic connections and differentiation of ground squirrels Marmotinae Pocock, 1923 (Rodentia, Sciuridae) from Palearctics]

Genetic distances between eight species of sousliks (Spermophilus) and five species of marmots (Marmota) were estimated on the basis of 39 biochemical loci. All taxa were shown to be genetically discrete. The genetic differentiation was minimal (Pfd = 11.3) between parapatric species of Palearctic sousliks of the suslicus pigmaeus group and Marmota species, intermediate (Pfd = 34.7) between allopatric sousliks species, and maximal (Pfd = 56.7) between representatives of different genera. The following trends were revealed in the geographic differentiation of the genus Spermophilus: (1) genetic similarity was associated with the geographic distance; (2) the eastern and western Palearctic phyla were markedly different genetically; (3) the eastern Palearctic forms exhibited higher differentiation than the western ones. The revealed speciation pattern is consistent with the general trend of temporal differentiation in Palearctic phyla and confirms the periodic speciation mode in the Palearctics.

Scaling of the mandible in squirrels

We compared the shape of the mandible among New World tree squirrels and selected outgroup taxa using linear measurements and areas defined by the median axis and conventional anatomical landmarks. We modified the median axis technique to define novel measurements, which proved complementary to those obtained from conventional landmarks. Allometric analyses showed that the scaling of the mandible among the New World tree squirrels is generally isometric (as has been observed in other groups of mammals), but diverges from isometry in a tendency in smaller animals for the masseteric ridge to be displaced anteriorly, the condylar process and posterior portion of the ascending ramus to be relatively elongated, and the coronoid process to be shortened. Allometric analyses also revealed the ways and extent that outgrowth taxa deviated from the scaling pattern observed for the New World tree squirrels. A flying squirrel (subfamily Pteromyinae), a moderate-sized callosciurine squirrel, and three species of pygmy tree squirrels from Asia and Africa show mandibular proportions very similar to those predicted for New World tree squirrels of corresponding size. Ground squirrels (tribe Marmotini) and successively more distant relatives such as Aplodontia, two myomorph rodents, and a rabbit show greater differences from the New World tree squirrels in their mandibular proportions. Combining the use of median-axis and conventional measurements makes it possible to examine changing relationships between locations of anatomically homologous landmarks and the geometry of the form.

Phylogenetic relationships among Japanese species of the family Sciuridae (Mammalia, Rodentia), inferred from nucleotide sequences of mitochondrial 12S ribosomal RNA genes

In order to investigate phylogenetic relationships of the family Sciuridae living in Japan, we sequenced partial regions (379 bases) of mitochondrial 12S rRNA genes in six species of Japanese and other Asian squirrels. Phylogenetic trees constructed by sequence data indicated that two genera of flying squirrels (Petaurista and Pteromys) were clustered in a group distinct from non-flying squirrels, suggesting a possible monophyletic relationships of these flying squirrels. The evolutionary distance between the Japanese squirrel (Sciurus lis) from Honshu island and the Eurasian red squirrel (Sciurus vulgaris) from Hokkaido island was comparable to intraspecific distances of the remaining species examined.

Phylogeny of six Sciurus aberti subspecies based on nucleotide sequences of cytochrome b

The tassel-eared squirrel, Sciurus aberti, is separated into six subspecies which occupy restricted and comparable habitats in ponderosa pine forests in the south-western United States and Mexico. These forests and squirrel populations are currently isolated by large arid areas and, as such, S. aberti appears to offer an example of incipient speciation. We sequenced the complete mitochondrial cytochrome b gene to construct a molecular phylogeny for S. aberti and to determine whether subspecific genetic structure and geographic patterns are correlative. Twenty alleles were identified among 612 squirrels throughout the species' range. Nucleotide divergence between alleles ranged from 0.009 to 0.0233, whereas average sequence divergence between S. aberti and an outgroup species, Sciurus niger, was 0.1823. Neighbor-joining and maximum parsimony analyses identified three major genetic assemblages composed of the following subspecies groups: (1) barberi and durangi; (2) aberti and kaibabensis; and (3) ferreus. The chuscensis samples were the only population with divergent sequences; one sequence was identical to an aberti sequence and a second unique sequence clustered with the ferreus sequences. The presence of divergent sequences in the chuscensis population, coupled with its central geographic position between aberti and ferreus, suggests a relatively recent influx of aberti mtDNA. Estimates of the times separating sequences in subspecies within different groups ranged from 0.94 to 1.52 x 10(6) years, based on a rate estimate of 7.15 x 10(-9) substitutions/year/site. The limited divergence observed between (1) aberti and kaibabensis as well as (2) barberi and durangi suggests relatively recent separations of subspecies within each assemblage. In fact, populations defined morphologically and geographically as ferreus exhibited greater sequence divergence than the aforementioned groups, identifying ferreus as the subspecies with the greatest genetic substructuring. The levels of cytochrome b divergence observed for the three distinct groups argues against a significant role for late Pleistocene glaciation in dispersal of this particular species; however, the proposed intermixing of aberti and chuscensis populations may well have been associated with such glacial events. Nucleotide diversity within subspecies ranked chuscensis >> aberti > barberi approximately kaibabensis approximately ferreus subpopulations; the relatively high level of diversity of chuscensis samples likely results from the apparent introgression of an aberti haplotype. The comparative levels of diversity in the aberti, barberi, kaibabensis, and ferreus sample populations do not correlate with respective habitat size (and presumably population size), suggesting that relatively recent forces, e.g., glaciation and inconsistent timber harvests, may have influenced diversity in these populations without apparent alterations in population size.

A recent origin of marmots

The complete mitochondrial cytochrome b genes from four ground squirrels, a marmot, and a tree squirrel were sequenced and compared. A level of divergence three times lower than expected based on fossil evidence was found between the marmot and the ground squirrels. The observed divergence was similar to that found within the ground squirrels and is consistent with a recent divergence of marmots (Marmota) from true ground squirrels (Spermophilus).

Variation in renal structure and urine concentrating capacity among ground squirrels of the Spermophilus townsendii complex (Rodentia: Sciuridae)

Renal structure and urine concentrating capacity varied significantly among nine populations of desert ground squirrels representing four closely related taxa. Variation in renal structure was weakly associated with taxonomic grouping. Urine concentrating capacity was not correlated with standard measures of habitat were aridity, but it was correlated with soil salinity. Maximum urine osmolalities for one population from the most saline habitat were among the highest documented for sciurid rodents. Results do not support the hypothesis that hibernating sciurids have poor renal efficiency.

Hybridization in the Mexican and 13-lined ground squirrels, Spermophilus mexicanus and Spermophilus tridecemlineatus

Evidence of hybridization between the ground squirrels, Spermophilus tridecemlineatus and S. mexicanus is presented on the basis of chromosomal and protein data. The hybrids produced appear to be completely interfertile. Based on the evidence of hybridization and the recent reestablishment of contact between these two species, they are considered to be semispecies.

The chromosomes of Spermophilus townsendi (Rodentia: Sciuridae) and report of a new subspecies

Chromosomes were analyzed from five subspecies of the ground squirrel, <i>Spermophilus townsendi</i>. The diploid chromosome number was 46 in <i>S</i>. <i>t</i>. <i>vigilis </i>and <i>S</i>. <i>t</i>. <i>canus, </i>38 in <i>S</i>. <i>t</i>. <i>idahoensis </i>and <i>S</i>. <i>t</i>. <i>mollis </i>and 36 in <i>S</i>. <i>t</i>. <i>townsendi</i>. A chromosomally homogenous population of ground squirrels, separated from <i>S</i>. <i>t</i>. <i>townsendi </i>by the Yakima River and having 38 instead of 36 chromosomes, was recognized as a new subspecies, <i>S</i>. <i>townsendi nancyae</i>. Chromosomal evidence linked with zoogeography suggest that the subspecies of <i>S</i>. <i>townsendi </i>radiated from an ancestral population (2n = 46), probably originating in Oregon, with eventual spread north across the Columbia River <i>(S</i>. <i>t</i>. <i>townsendi </i>and <i>nancyae</i>) and southeast into Idaho, Utah and Nevada <i>(S</i>. <i>t</i>. <i>mollis </i>and <i>idahoensis</i>). Robertsonian centric fusion rather than dissociation of metacentric chromosomes receives strong zoogeographic support as the mechanism responsible for karyotype evolution in <i>S</i>. <i>townsendi.</i>