Basal Clades and Molecular Systematics of Heteromyid Rodents

Springing into action: Comparing escape responses between bipedal and quadrupedal rodents

Predation is a fundamental selective pressure on animal morphology, as morphology is directly linked with physical performance and evasion. Bipedal heteromyid rodents, which are characterized by unique morphological traits such as enlarged hindlimbs, appear to be more successful than sympatric quadrupedal rodents at escaping predators such as snakes and owls, but no studies have directly compared the escape performance of bipedal and quadrupedal rodents. We used simulated predator attacks to compare the evasive jumping ability of bipedal kangaroo rats (Dipodomys) to that of three quadrupedal rodent groups-pocket mice (Chaetodipus), woodrats (Neotoma), and ground squirrels (Otospermophilus). Jumping performance of pocket mice was remarkably similar to that of kangaroo rats, which may be driven by their shared anatomical features (such as enlarged hindlimb muscles) and facilitated by their relatively small body size. Woodrats and ground squirrels, in contrast, almost never jumped as a startle response, and they took longer to perform evasive escape maneuvers than the heteromyid species (kangaroo rats and pocket mice). Among the heteromyids, take-off velocity was the only jump performance metric that differed significantly between species. These results support the idea that bipedal body plans facilitate vertical leaping in larger-bodied rodents as a means of predator escape and that vertical leaping likely translates to better evasion success.

Tectonic extension and paleoelevation influence mammalian diversity dynamics in the Basin and Range Province of western North America

Landscape properties have a profound influence on the diversity and distribution of biota, with present-day biodiversity hot spots occurring in topographically complex regions globally. Complex topography is created by tectonic processes and further shaped by interactions between climate and land-surface processes. These processes enrich diversity at the regional scale by promoting speciation and accommodating increased species richness along strong environmental gradients. Synthesis of the mammalian fossil record and a geophysical model of topographic evolution of the Basin and Range Province in western North America enable us to directly quantify relationships between mammal diversity and landscape dynamics over the past 30 million years. We analyze the covariation between tectonic history (extensional strain rates, paleotopography, and ruggedness), global temperature, and diversity dynamics. Mammal species richness and turnover exhibit stronger responses to rates of change in landscape properties than to the specific properties themselves, with peaks in diversity coinciding with high tectonic strain rates and large changes in elevation across spatial scales.

Updated list of the mammals of Costa Rica, with notes on recent taxonomic changes

Although Costa Rica occupies a mere 0.03% of the Earths land area, it nevertheless has recorded within its borders approximately 5% of the global diversity of mammals, thus making it one of the worlds megadiverse countries. Over the past ten years, 22 species have been added to the countrys inventory, bringing the total number known as here documented to 271; Chiroptera account for ten of these, having grown to 124 from 114; rodents have increased by eight species, from 47 to 55, with the caveat that we include three invasive species of Muridae that have gone feral. In contrast, the number of orders has decreased by one, by Artiodactyla incorporating the former Cetacea. Notes are provided for all taxonomic novelties since the last update. Since the first taxonomic compendium of the mammals of Costa Rica in 1869, the number of known species has grown by approximately 1.22 species year-1 (R2 = 0.96). Since 1983 however, this growth rate has been 1.64 species year-1 (R2 = 0.98). Despite this strong growth, an asymptote in the number of known species has not been reached. Conservation remains a primary need: over 60% of the countrys mammal species show population trends that are decreasing (13%), unknown (37%), or not assessed (11%), based on IUCN criteria. These analyses suggest that much remains to be known regarding the number of mammal species living in Costa Rica, but also that much more remains to be done to safeguard Costa Ricas exceptional biodiversity heritage.

Multivariate analyses of skull morphology inform the taxonomy and evolution of geomyoid rodents

Morphological analyses are critical to quantify phenotypic variation, identify taxa, inform phylogenetic relationships, and shed light on evolutionary patterns. This work is particularly important in groups that display great morphological disparity. Such is the case in geomyoid rodents, a group that includes 2 of the most species-rich families of rodents in North America: the Geomyidae (pocket gophers) and the Heteromyidae (kangaroo rats, pocket mice, and their relatives). We assessed variation in skull morphology (including both shape and size) among geomyoids to test the hypothesis that there are statistically significant differences in skull measurements at the family, genus, and species levels. Our sample includes 886 specimens representing all geomyoid genera and 39 species. We used the geometric mean to compare size across taxa. We used 14 measurements of the cranium and lower jaw normalized for size to compare shape among and within taxa. Our results show that skull measurements enable the distinction of geomyoids at the family, genus, and species levels. There is a larger amount of size variation within Geomyidae than within Heteromyidae. Our phylomorphospace analysis shows that the skull shape of the common ancestor of all geomyoids was more similar to the common ancestor of heteromyids than that of geomyids. Geomyid skulls display negative allometry whereas heteromyid skulls display positive allometry. Within heteromyids, dipodomyines, and non-dipodomyines show significantly different allometric patterns. Future analyses including fossils will be necessary to test our evolutionary hypotheses.

RM. The earliest dipodomyine heteromyid in North America and the phylogenetic relationships of geomorph rodents

Dipodomyine heteromyids (kangaroo rats and mice) are a diverse group of arid-adapted ricochetal rodents of North America. Here, a new genus and species of a large dipodomyine is reported from early Miocene-aged deposits of the John Day Formation in Oregon that represents the earliest record of the subfamily. The taxon is known from a single specimen consisting of a nearly complete skull, dentary, partial pes, and caudal vertebra. The specimen is characterized by a mosaic of ancestral and highly derived cranial features of heteromyids. Specifically, the dental morphology and some cranial characteristics are similar to early heteromyids, but other aspects of morphology, including the exceptionally inflated auditory bullae, are more similar to known dipodomyines. This specimen was included in a phylogenetic analysis comprising 96 characters and the broadest sampling of living and extinct geomorph rodents of any morphological phylogenetic analysis to date. Results support the monophyly of crown-group Heteromyidae exclusive of Geomyidae and place the new taxon within Dipodomyinae. The new heteromyid is the largest known member of the family. Analyses suggest that large body size evolved several times within Heteromyidae. Overall, the morphology of the new heteromyid supports a mosaic evolution of the open-habitat adaptations that characterize kangaroo rats and mice, with the inflation of the auditory bulla appearing early in the group, and bipedality/ricochetal locomotion appearing later. We hypothesize that cooling and drying conditions in the late Oligocene and early Miocene favored adaptations for life in more open habitats, resulting in increased locomotor specialization in this lineage over time from a terrestrial ancestor.

A chromosome-length reference genome for the endangered Pacific pocket mouse reveals recent inbreeding in a historically large population

High-quality reference genomes are fundamental tools for understanding population history, and can provide estimates of genetic and demographic parameters relevant to the conservation of biodiversity. The federally endangered Pacific pocket mouse (PPM), which persists in three small, isolated populations in southern California, is a promising model for studying how demographic history shapes genetic diversity, and how diversity in turn may influence extinction risk. To facilitate these studies in PPM, we combined PacBio HiFi long reads with Omni-C and Hi-C data to generate a de novo genome assembly, and annotated the genome using RNAseq. The assembly comprised 28 chromosome-length scaffolds (N50 = 72.6 MB) and the complete mitochondrial genome, and included a long heterochromatic region on chromosome 18 not represented in the previously available short-read assembly. Heterozygosity was highly variable across the genome of the reference individual, with 18% of windows falling in runs of homozygosity (ROH) >1 MB, and nearly 9% in tracts spanning >5 MB. Yet outside of ROH, heterozygosity was relatively high (0.0027), and historical N_e estimates were large. These patterns of genetic variation suggest recent inbreeding in a formerly large population. Currently the most contiguous assembly for a heteromyid rodent, this reference genome provides insight into the past and recent demographic history of the population, and will be a critical tool for management and future studies of outbreeding depression, inbreeding depression, and genetic load.

Biomechanical adaptations for burrowing in the incisor enamel microstructure of Geomyidae and Heteromyidae (Rodentia: Geomyoidea)

The enamel microstructure of fossil and extant Geomyoidea (Geomyidae, Heteromyidae) lower incisors incorporates three- or two-layered schmelzmusters with uniserial, transverse Hunter-Schreger bands having parallel and perpendicular or exclusively perpendicular oriented interprismatic matrix. Phylogenetically, these schmelzmusters are regarded as moderately (enamel type 2) to highly derived (enamel type 3). Our analysis detected a zone of modified radial enamel close to the enamel-dentine junction. Modified radial enamel shows a strong phylogenetic signal within the clade Geomorpha as it is restricted to fossil and extant Geomyoidea and absent in Heliscomyidae, Florentiamyidae, and Eomyidae. This character dates back to at least the early Oligocene (early Arikareean, 29 Ma), where it occurs in entoptychine gophers. We contend that this specialized incisor enamel architecture developed as a biomechanical adaptation to regular burrowing activities including chisel-tooth digging and a fiber-rich diet and was probably present in the common ancestor of the clade. We regard the occurrence of modified radial enamel in lower incisors of scratch-digging Geomyidae and Heteromyidae as the retention of a plesiomorphic character that is selectively neutral. The shared occurrence of modified radial enamel is a strong, genetically anchored argument for the close phylogenetic relationship of Geomyidae and Heteromyidae on the dental microstructure level.

A NEW SPECIES OF EIMERIA (APICOMPLEXA: EIMERIIDAE) FROM THE OLIVE-BACKED POCKET MOUSE, PEROGNATHUS FASCIATUS (RODENTIA: HETEROMYIDAE: PEROGNATHINAE), FROM WYOMING

Forty-nine olive-backed pocket mice, Perognathus fasciatus were collected during 2011 and 2012 from 4 sites in Wyoming and examined for coccidian parasites. Fifteen (31%) were found to be passing oocysts of a new species of Eimeria Schneider, 1875. Sporulated oocysts of Eimeria fasciata n. sp. are ellipsoidal to ovoidal, 23.3 × 20.7 (19-27 × 17-25) μm, with a shape index of 1.1; they typically contain a single, smooth, bubble-like oocyst residuum. Oocysts possess 1-2 polar granules, lack a micropyle, and are bilayered with a thickness of 1.3 μm. Sporocysts are ovoidal, 10.0 × 8.2 (8-12 × 7-10) μm, with a shape index of 1.2; they contain a sporocyst residuum that appears similar to a cluster of 1-8 grapes. The Stieda body is small, appearing flattened to knobby, and there are no subStieda or paraStieda bodies. This new eimerian represents the only coccidian, to date, reported from P. fasciatus, as well as the only species from any heteromyid rodent in Wyoming.

New species of Trichuris (Nematoda: Trichuridae) parasitizing Heteromys salvini (Rodentia: Heteromyidae) from Costa Rica, with a key to Trichuris species described from Heteromyidae

Trichuris guanacastei n sp., a parasite isolated from the Salvin' spiny mouse Heteromys salvini, collected from the Guanacaste Conservation Area, Costa Rica, during February 1996 is described. The new species was compared morphologically with the 29 known species that parasitize rodents distributed in 12 families in North and South America; T. guanacastei n. sp. it is characterized by the following set of traits: presence of a spicular tube (measuring 0.72-0.99); thick proximal cloacal tube and a short distal cloacal tube with a total length of 0.72-1.36; eggs 0.03-0.05 long and a semi-protrusible vulva. The new species represents the first one described in the genus in Costa Rica and the fifth one described in the Americas that parasite Heteromyidae.

Geographic variation and evolutionary history of Dipodomys nitratoides (Rodentia: Heteromyidae), a species in severe decline

We examined geographic patterns of diversification in the highly impacted San Joaquin kangaroo rat, Dipodomys nitratoides, throughout its range in the San Joaquin Valley and adjacent basins in central California. The currently recognized subspecies were distinct by the original set of mensural and color variables used in their formal diagnoses, although the Fresno kangaroo rat (D. n. exilis) is the most strongly differentiated with sharp steps in character clines relative to the adjacent Tipton (D. n. nitratoides) and short-nosed (D. n. brevinasus) races. The latter two grade more smoothly into one another but still exhibit independent, and different, character clines within themselves. At the molecular level, as delineated by mtDNA cytochrome b sequences, most population samples retain high levels of diversity despite significant retraction in the species range and severe fragmentation of local populations in recent decades due primarily to landscape conversion for agriculture and secondarily to increased urbanization. Haplotype apportionment bears no relationship to morphologically defined subspecies boundaries. Rather, a haplotype network is shallow, most haplotypes are single-step variants, and the time to coalescence is substantially more recent than the time of species split between D. nitratoides and its sister taxon, D. merriami. The biogeographic history of the species within the San Joaquin Valley appears tied to mid-late Pleistocene expansion following significant drying of the valley resulting from the rain shadow produced by uplift of the Central Coastal Ranges.

Rediscovery of Heteromys nelsoni in its type locality after over a century

We rediscovered a population of Nelson’s spiny pocket mouse (Heteromys nelsoni; Merriam, 1902) in the type locality of Pinabeto in the Mexican state of Chiapas, 121 years after it was last collected. We describe five topotype specimens according to their morphology and external measurements, and we confirm its identity at the species level in the Basic Local Alignment Search Tool (BLAST) of GenBank. As the population of H. nelsoni in Pinabeto is isolated, it is likely to be susceptible to extinction. There is a need to carry out additional scientific studies of this microendemic species in order to obtain more information regarding its biology, ecology and evolutionary history, and to be able to influence environmental policy to protect and conserve this species, as well as the region’s cloud forests.

Biodiversity and Topographic Complexity: Modern and Geohistorical Perspectives

Topographically complex regions on land and in the oceans feature hotspots of biodiversity that reflect geological influences on ecological and evolutionary processes. Over geologic time, topographic diversity gradients wax and wane over millions of years, tracking tectonic or climatic history. Topographic diversity gradients from the present day and the past can result from the generation of species by vicariance or from the accumulation of species from dispersal into a region with strong environmental gradients. Biological and geological approaches must be integrated to test alternative models of diversification along topographic gradients. Reciprocal illumination among phylogenetic, phylogeographic, ecological, paleontological, tectonic, and climatic perspectives is an emerging frontier of biogeographic research.

Nematofauna of Rodents of the Families Heteromyidae and Cricetidae from the Mexican Plateau

As a part of an ongoing project to inventory the helminth parasites of rodents in Mexico, 85 specimens of 2 families of rodents were collected from the Mexican Plateau: Cricetidae ( Neotoma sp., Neotoma leucodon , Onychomys arenicola , Peromyscus sp., Peromyscus eremicus , and Reithrodontomys sp.) and Heteromyidae ( Chaetodipus sp., Chaetodipus eremicus , Chaetodipus hispidus , Dipodomys merriami , Dipodomys ordii , Dipodomys ornatus, Dipodomys spectabilis , Liomys irroratus , Perognathus sp., and Perognathus flavus ). A total of 13 taxa of helminths were found: Heteromyoxyuris longejector, Heteromyoxyuris otomii, Heteromyoxyuris sp., Onchocercidae gen. sp. 1 and sp. 2, Physalopteridae gen. sp., Protospirura dipodomis, Pterygodermatites dipodomis, Subulura sp., Syphacia sp., Trichuris dipodomis, Vexillata liomyos, and Vexillata armande. The highest species richness was recorded in D. merriami (7 taxa). This study is the first report of nematodes from O. arenicola (Physalopteridae gen. sp.) and C. eremicus (H. longejector) and for V. liomyos from D. merriami . All reports of these species of nematodes represent new collection localities in Mexico.

Intercontinental-wide consequences of compromise-breaking adaptations: the case of desert rodents

Desert rodent assemblages from around the world provide convergent, but independent crucibles for testing theory and deducing general ecological principles. The heteromyid rodents of North America and the gerbils of the Middle East and their predators provide such an example. Both sets of rodents face predation from owls and vipers, but the North American species possess unique traits that may represent macroevolutionary breakthroughs: rattlesnakes have infra-red sensitive sensory pits, and heteromyids have cheek pouches. To test their significance, we brought together two gerbils (Middle East), two heteromyid rodents (a kangaroo rat and a pocket mouse; North America) in a common setting (a vivarium in the Negev Desert), and quantified the “opinions” of the rodents towards the North American sidewinder rattlesnake and the Middle Eastern Saharan horned viper and the foraging behavior of each in the face of these snake predators plus owl predators. Gerbils are fairly evenly matched in their anti-predator abilities, while the heteromyids differ widely, and these seem to match well with and may determine the types of mechanisms of species coexistence that operate in the communities of each continent. Evolutionary history, macroevolutionary traits, and risk management therefore combine to determine the characteristics of the organisms and the organization of their communities.

Transcriptomic characterization of the immunogenetic repertoires of heteromyid rodents

Background

When populations evolve under disparate environmental conditions, they experience different selective pressures that shape patterns of sequence evolution and gene expression. These may be manifested in genetic and phenotypic differences such as a diverse immunogenetic repertoire in species from tropical latitudes that have greater and/or different parasite burdens than more temperate species. To test this idea, we compared the transcriptomes of one tropical species (Heteromys desmarestianus) and two species from temperate latitudes (Dipodomys spectabilis and Chaetodipus baileyi) from the Heteromyidae. We did so in a search for positive selection on sequences and/or differential expression, while controlling for phylogenetic history in our choice of species.

Results

We identified 127,812 contigs and annotated 34,878 of these, identifying immune genes associated with interleukins, cytokines, and the production of mast cells. We identified 632 genes that were upregulated in H. desmarestianus (8.7% of genes tested) and 492 (6.7%) that were downregulated. Gene ontology terms including “immune response” were associated with 31 (4.9%) of the 632 upregulated genes. We found preliminary evidence for positive selection on three genes (Palmitoyltransferase ZDHHC5 Ubiquitin-conjugating enzyme E2 N, Krueppel-like factor 10, and Spindle and kinetochore-associated protein 1) along the H. desmarestianus lineage.

Conclusions

Overall our findings pinpoint genes in species from disparate environments that are on different evolutionary trajectories in terms of expression levels and/or nucleotide sequence. Our data indicate there are significant differences in the expression of genes among the spleen transcriptomes of these species and that a number of these differentially expressed genes do not show the same pattern of differential expression in another tissue type. This points to the possibility of expression differences between these species specific to the spleen transcriptome.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-929) contains supplementary material, which is available to authorized users.

Natural selection and the genetic basis of osmoregulation in heteromyid rodents as revealed by RNA-seq

One adaptation of ecological and evolutionary interest is the extraordinary ability of desert rodents to retain water during waste production. Much is known regarding the unique kidney physiology of kangaroo rats (Dipodomys spp.) and their ability to retain water during waste production, yet the genetic basis of these physiological adaptations is relatively unknown. Herein, we utilized RNA-seq data to conduct a comparative study to identify osmoregulatory genes expressed in heteromyid rodents. We sequenced kidney tissue from two temperate desert species (Dipodomys spectabilis and Chaetodipus baileyi) from two separate subfamilies of the Heteromyidae and compared these transcriptomes to a tropical mesic species (Heteromys desmarestianus) from a third subfamily. The evolutionary history of these subfamilies provided a robust phylogenetic control that allowed us to separate shared evolutionary history from convergence. Using two methods to detect differential expression (DE), we identified 1890 genes that showed consistent patterns of DE between the arid and mesic species. A three-species reciprocal BLAST analysis revealed 3511 sets of putative orthologues that, upon comparison to known Mus musculus sequences, revealed 323 annotated and full-length genic regions. Selection tests displayed evidence of positive selection (dn/ds > 1) on six genes in the two desert species and remained significant for one of these genes after correction for multiple testing. Thus, our data suggest that both the coding sequence and expression of genes have been shaped by natural selection to provide the genetic architecture for efficient osmoregulation in desert-adapted heteromyid rodents.

Phylogenetic analysis of mammalian maximal oxygen consumption during exercise

We compiled published values of mammalian maximum oxygen consumption during exercise (VO2max) and supplemented these data with new measurements of VO2max for the largest rodent (capybara), 20 species of smaller-bodied rodents, two species of weasels, and one small marsupial. Many of the new data were obtained with running-wheel respirometers instead of the treadmill systems used in most previous measurements of mammalian VO2max. We used both conventional and phylogenetically informed allometric regression models to analyze VO2max of 77 ‘species’ (including subspecies or separate populations within species) in relation to body size, phylogeny, diet, and measurement method. Both body mass and allometrically mass-corrected VO2max showed highly significant phylogenetic signal (i.e., related species tended to resemble each other). The Akaike Information Criterion corrected for sample size was used to compare 27 candidate models predicting VO2max (all of which included body mass). In addition to mass, the two best-fitting models (cumulative Akaike weight = 0.93) included dummy variables coding for three species previously shown to have high VO2max (pronghorn, horse, and a bat), and incorporated a transformation of the phylogenetic branch lengths under an Ornstein-Uhlenbeck model of residual variation (thus indicating phylogenetic signal in the residuals). We found no statistical difference between wheel- and treadmill-elicited values, and diet had no predictive ability for VO2max. Averaged across all models, the allometric scaling exponent was 0.839, with 95% confidence limits of 0.795 and 0.883, which does not provide support for a scaling exponent of 0.67, 0.75 or unity.

Co-occurrence of small mammals in a tropical dry deciduous forest: comparisons of communities and individual species in Colima, Mexico

Species co-occurrence is an important ecological research area. Mark-and-recapture studies of mammals allow identification of coexisting species, a necessary step in determining mechanisms enabling habitat sharing. Using data from five 1-ha grids in January 2004 in tropical dry deciduous forest of coastal Colima, Mexico, we detected significantly more interspecific overlap than expected between seven species pairs. Oryzomys couesi shared more stations than expected with Sigmodon mascotensis, Baiomys musculus and Peromyscus perfulvus. Baiomys musculus was associated positively with S. mascotensis and Reithrodontomys fulvescens. Heteromys pictus shared fewer stations than expected with O. couesi and S. mascotensis. For species collectively, there was non-random community structuring, with two grids displaying more species aggregation than expected. While two grids had non-random co-occurrence patterns, three grids did not differ from random, which differs from that reported for mammalian taxa on average. Other small-mammal studies have documented species segregation, while this study detected more positive than negative associations. Similarities in preference and habitat use (or diet) are likely explanations for interspecific overlap patterns at stations and co-occurrence patterns among grids. Simultaneously evaluating associations of species pairs and all species on a grid collectively is novel methodology as applied to mammals, adding to understanding of species co-occurrence.

Genetic diversity, recombination, and divergence in animal associated Penicillium dipodomyis

Penicillium dipodomyis is thought to be an exclusively asexual fungus associated with Kangaroo Rats, Dipodomys species, and is unique among Penicillium species in growing at 37°C but producing no known toxins. Lack of recombination within P. dipodomyis would result in limited adaptive flexibility but possibly enhance local adaptation and host selection via maintenance of favourable genotypes. Here, analysis of DNA sequence data from five protein-coding genes shows that recombination occurs within P. dipodomyis on a small spatial scale. Furthermore, detection of mating-type alleles supports outcrossing and a sexual cycle in P. dipodomyis. P. dipodomyis was a weaker competitor in in vitro assays with other Penicillium species found in association with Kanagaroo rats. Bayesian species level analysis suggests that the P. dipodomyis lineage diverged from closely related species also found in cheek pouches of Kangaroo Rats and their stored seeds about 11 million years ago, a similar divergence time as Dipodomys from its sister rodent taxa.

Phylogeography of the dark kangaroo mouse, Microdipodops megacephalus: cryptic lineages and dispersal routes in North America's Great Basin

AIM: The rodent genus Microdipodops (kangaroo mice) includes two sand-obligate endemics of the Great Basin Desert: M. megacephalus and M. pallidus. The dark kangaroo mouse, M. megacephalus, is distributed throughout the Great Basin and our principal aims were to formulate phylogenetic hypotheses for this taxon and make phylogeographical comparisons with its congener. LOCATION: The Great Basin Desert of western North America. METHODS: DNA sequence data from three mitochondrial genes were examined from 186 individuals of M. megacephalus, representing 47 general localities. Phylogenetic inference was used to analyse the sequence data. Directional analysis of phylogeographical patterns was used to examine haplotype sharing patterns and recover routes of gene exchange. Haplotype-area curves were constructed to evaluate the relationship between genetic variation and distributional island size for M. megacephalus and M. pallidus. RESULTS: Microdipodops megacephalus is a rare desert rodent (trapping success was 2.67%). Temporal comparison of trapping data shows that kangaroo mice are becoming less abundant in the study area. The distribution has changed slightly since the 1930s but many northern populations now appear to be small, fragmented, or locally extinct. Four principal phylogroups (the Idaho isolate and the western, central and eastern clades) are evident; mean sequence divergence between phylogroups for cytochrome b is c. 8%. Data from haplotype sharing show two trends: a north-south trend and a web-shaped trend. Analyses of haplotype-area curves reveal significant positive relationships. MAIN CONCLUSIONS: The four phylogroups of M. megacephalus appear to represent morphologically cryptic species; in comparison, a companion study revealed two cryptic lineages in M. pallidus. Estimated divergence times of the principal clades of M. megacephalus (c. 2-4 Ma) indicate that these kangaroo mice were Pleistocene invaders into the Great Basin coincident with the formation of sandy habitats. The north-south and web patterns from directional analyses reveal past routes of gene flow and provide evidence for source-sink population regulation. The web pattern was not seen in the companion study of M. pallidus. Significant haplotype-area curves indicate that the distributional islands are now in approximate genetic equilibrium. The patterns described here are potentially useful to conservation biologists and wildlife managers and may serve as a model for other sand-obligate organisms of the Great Basin.