Adaptive tail-length evolution in deer mice is associated with differential Hoxd13 expression in early development.

Variation in the size and number of axial segments underlies much of the diversity in animal body plans. Here we investigate the evolutionary, genetic and developmental mechanisms driving tail-length differences between forest and prairie ecotypes of deer mice (Peromyscus maniculatus). We first show that long-tailed forest mice perform better in an arboreal locomotion assay, consistent with tails being important for balance during climbing. We then identify six genomic regions that contribute to differences in tail length, three of which associate with caudal vertebra length and the other three with vertebra number. For all six loci, the forest allele increases tail length, indicative of the cumulative effect of natural selection. Two of the genomic regions associated with variation in vertebra number contain Hox gene clusters. Of those, we find an allele-specific decrease in Hoxd13 expression in the embryonic tail bud of long-tailed forest mice, consistent with its role in axial elongation. Additionally, we find that forest embryos have more presomitic mesoderm than prairie embryos and that this correlates with an increase in the number of neuromesodermal progenitors, which are modulated by Hox13 paralogues. Together, these results suggest a role for Hoxd13 in the development of natural variation in adaptive morphology on a microevolutionary timescale.

Pair bonding and disruption impact lung transcriptome in monogamous Peromyscus californicus.

Social interactions affect physiological and pathological processes, yet their direct impact in peripheral tissues remains elusive. Recently we showed that disruption of pair bonds in monogamous Peromyscus californicus promotes lung tumorigenesis, pointing to a direct effect of bonding status in the periphery (Naderi et al., 2021). Here we show that lung transcriptomes of tumor-free Peromyscus are altered in a manner that depends on pair bonding and superseding the impact of genetic relevance between siblings. Pathways affected involve response to hypoxia and heart development. These effects are consistent with the profile of the serum proteome of bonded and bond-disrupted Peromyscus and were extended to lung cancer cells cultured in vitro, with sera from animals that differ in bonding experiences. In this setting, the species' origin of serum (deer mouse vs FBS) is the most potent discriminator of RNA expression profiles, followed by bonding status. By analyzing the transcriptomes of lung cancer cells exposed to deer mouse sera, an expression signature was developed that discriminates cells according to the history of social interactions and possesses prognostic significance when applied to primary human lung cancers. The results suggest that present and past social experiences modulate the expression profile of peripheral tissues such as the lungs, in a manner that impacts physiological processes and may affect disease outcomes. Furthermore, they show that besides the direct effects of the hormones that regulate bonding behavior, physiological changes influencing oxygen metabolism may contribute to the adverse effects of bond disruption.

Gene regulatory changes underlie developmental plasticity in respiration and aerobic performance in highland deer mice.

Phenotypic plasticity can play an important role in the ability of animals to tolerate environmental stress, but the nature and magnitude of plastic responses are often specific to the developmental timing of exposure. Here, we examine changes in gene expression in the diaphragm of highland deer mice (Peromyscus maniculatus) in response to hypoxia exposure at different stages of development. In highland deer mice, developmental plasticity in diaphragm function may mediate changes in several respiratory traits that influence aerobic metabolism and performance under hypoxia. We generated RNAseq data from diaphragm tissue of adult deer mice exposed to (1) life-long hypoxia (before conception to adulthood), (2) post-natal hypoxia (birth to adulthood), (3) adult hypoxia (6-8 weeks only during adulthood) or (4) normoxia. We found five suites of co-regulated genes that are differentially expressed in response to hypoxia, but the patterns of differential expression depend on the developmental timing of exposure. We also identified four transcriptional modules that are associated with important respiratory traits. Many of the genes in these transcriptional modules bear signatures of altitude-related selection, providing an indirect line of evidence that observed changes in gene expression may be adaptive in hypoxic environments. Our results demonstrate the importance of developmental stage in determining the phenotypic response to environmental stressors.

The discovery of multiple active mys-related LTR-retroelements within the Neotominae subfamily of cricetid rodents.

Retrotransposon families in the rodent family Cricetidae have been understudied in contrast to Muridae, both taxa classified within the superfamily Muroidea. Therefore, we carried out a study to advance our knowledge of the unique mys LTR-retroelement identified in Peromyscus leucopus, by incorporating intra-ORF PCR, quantitative dot blots, DNA and protein library screens, the generation of molecular phylogenies, and analyses of orthologous LTR-retroelement loci. These analyses led to the discovery of three additional related families of LTR-retroelements, which include a 2900 bp full-length element of mys-related sequences (mysRS), an 8000 bp element containing the mys ORF1 sequence (mORF1) with ERV-related sequences downstream in the reverse orientation, as well as an 1800 bp element primarily consisting of mys ORF2 (mORF2) related sequences flanked by LTRs. Our data revealed only a few full-length mys elements among genera of the Neotominae subfamily of cricetid rodents, most existing as partial copies. The mysRS and mORF1 elements are also limited to the genomes of the Neotominae subfamily, whereas mORF2 appears to be restricted to the Peromyscus genus. Molecular phylogenies demonstrating concerted evolution along with an assessment of orthologous loci in Peromyscus for the presence or absence of elements are consistent with activity of these novel LTR-retroelement families within this genus. Together with known activity of various families of non-LTR retroelements in Peromyscus species, we propose that retrotransposons have been continually contributing to the dynamics of Peromyscus genomes promoting genomic diversity and may be correlated with the evolution of more than 50 identified Peromyscus species.

Transposable Element Interactions Shape the Ecology of the Deer Mouse Genome.

The genomic landscape of transposable elements (TEs) varies dramatically across species, with some TEs demonstrating greater success in colonizing particular lineages than others. In mammals, long interspersed nuclear element (LINE) retrotransposons are typically more common than any other TE. Here, we report an unusual genomic landscape of TEs in the deer mouse, Peromyscus maniculatus. In contrast to other previously examined mammals, long terminal repeat elements occupy more of the deer mouse genome than LINEs (11% and 10%, respectively). This pattern reflects a combination of relatively low LINE activity and a massive invasion of lineage-specific endogenous retroviruses (ERVs). Deer mouse ERVs exhibit diverse origins spanning the retroviral phylogeny suggesting they have been host to a wide range of exogenous retroviruses. Notably, we trace the origin of one ERV lineage, which arose ∼5-18 million years ago, to a close relative of feline leukemia virus, revealing inter-ordinal horizontal transmission. Several lineage-specific ERV subfamilies have very high copy numbers, with the top five most abundant accounting for ∼2% of the genome. We also observe a massive amplification of Kruppel-associated box domain-containing zinc finger genes, which likely control ERV activity and whose expansion may have been facilitated by ectopic recombination between ERVs. Finally, we find evidence that ERVs directly impacted the evolutionary trajectory of LINEs by outcompeting them for genomic sites and frequently disrupting autonomous LINE copies. Together, our results illuminate the genomic ecology that shaped the unique deer mouse TE landscape, shedding light on the evolutionary processes that give rise to variation in mammalian genome structure.

Phylogeography within the Peromyscus maniculatus species group: Understanding past distribution of genetic diversity and areas of refugia in western North America.

The effects of anthropogenic climate change on biodiversity have been recognized on every continent, ocean, and across different taxonomic groups. Here, we study the range dynamics and demography of a cosmopolitan species: the deer mouse, Peromyscus maniculatus. We generated a multilocus SNP dataset using the ddRADseq protocol for 218 individuals across the geographic range within three western North American lineages of this species group. We evaluated population structure using several methods and explored the correlation between geographic and genetic distances. We modeled the demographic history using a site frequency spectrum approach and used a machine learning algorithm to infer current and past (Last Glacial Maximum; LGM) environmental suitability. Lastly, we explored the origin of population expansion for the identified lineages. The genome-wide SNP dataset was able to identify-three regionally distinct groups- 1) P. m. gambelii (southern California); 2) P. keeni (Pacific Northwest); 3) P. m. sonoriensis (a broad population spanning the Pacific Northwest through central California and across the Rocky Mountains into the Great Plains). Demographic analysis indicated the splits between the three populations occurred within the last 500 thousand years, with one very recent (late Holocene) split. Ecological niche models for each of these lineages predicted suitable environment present throughout their known ranges for current conditions, and a severe reduction of northern habitat in the past. The deer mouse has responded to past climate changes by expanding its range during interglacial periods and contracting its range during glacial periods leading to strong population differentiation. But lower magnitude climate change or other processes within the Holocene interglacial period led to population differentiation as well, which is likely still ongoing today given the substantial anthropogenic climate change and other landscape transformations caused by humans during the Anthropocene. By understanding the historical processes that led to the contemporary geographic distribution of biodiversity, we can determine the relative importance of different factors that shape biodiversity, now and into the future.

Genetic Analysis of the Stereotypic Phenotype in Peromyscus maniculatus (deer mice).

Peromyscus maniculatus, including the laboratory stock BW, have been used as a model organism for autism spectrum disorder and obsessive-compulsive disorder because of the high occurrence of stereotypy. Several studies have identified neurological and environmental components of the phenotype; however, the heritability of the phenotype has not been examined. This study characterizes the incidence and heritability of vertical jumping stereotypy (VS) and backflipping (BF) behavior in the BW stock of the Peromyscus Genetic Stock Center, which are indicative of autism spectrum disorders. In addition, interspecies crosses between P. maniculatus and P. polionotus were also performed to further dissect genetically stereotypic behavior. The inheritance pattern of VS suggests that multiple genes result in a quantitative trait with low VS being dominant over high VS. The inheritance pattern of BF suggests that fewer genes are involved, with one allele causing BF in a dominant fashion. An association analysis in BW could reveal the underlying genetic loci associated with stereotypy in P. maniculatus, especially for the BF behavior.

Sexual selection drives the coevolution of male and female reproductive traits in Peromyscus mice.

When females mate with multiple partners within a single reproductive cycle, sperm from rival males may compete for fertilization of a limited number of ova, and females may bias the fertilization of their ova by particular sperm. Over evolutionary timescales, these two forms of selection shape both male and female reproductive physiology when females mate multiply, yet in monogamous systems, post-copulatory sexual selection is weak or absent. Here, we examine how divergent mating strategies within a genus of closely related mice, Peromyscus, have shaped the evolution of reproductive traits. We show that in promiscuous species, males exhibit traits associated with increased sperm production and sperm swimming performance, and females exhibit traits that are predicted to limit sperm access to their ova including increased oviduct length and a larger cumulus cell mass surrounding the ova, compared to monogamous species. Importantly, we found that across species, oviduct length and cumulus cell density are significantly correlated with sperm velocity, but not sperm count or relative testes size, suggesting that these female traits may have coevolved with increased sperm quality rather than quantity. Taken together, our results highlight how male and female traits evolve in concert and respond to changes in the level of post-copulatory sexual selection.

Chromosomal inversion polymorphisms shape the genomic landscape of deer mice.

Chromosomal inversions are an important form of structural variation that can affect recombination, chromosome structure and fitness. However, because inversions can be challenging to detect, the prevalence and hence the significance of inversions segregating within species remains largely unknown, especially in natural populations of mammals. Here, by combining population-genomic and long-read sequencing analyses in a single, widespread species of deer mouse (Peromyscus maniculatus), we identified 21 polymorphic inversions that are large (1.5-43.8 Mb) and cause near-complete suppression of recombination when heterozygous (0-0.03 cM Mb-1). We found that inversion breakpoints frequently occur in centromeric and telomeric regions and are often flanked by long inverted repeats (0.5-50 kb), suggesting that they probably arose via ectopic recombination. By genotyping inversions in populations across the species' range, we found that the inversions are often widespread and do not harbour deleterious mutational loads, and many are likely to be maintained as polymorphisms by divergent selection. Comparisons of forest and prairie ecotypes of deer mice revealed 13 inversions that contribute to differentiation between populations, of which five exhibit significant associations with traits implicated in local adaptation. Taken together, these results show that inversion polymorphisms have a significant impact on recombination, genome structure and genetic diversity in deer mice and likely facilitate local adaptation across the widespread range of this species.

An enhancer of Agouti contributes to parallel evolution of cryptically colored beach mice.

Identifying the genetic basis of repeatedly evolved traits provides a way to reconstruct their evolutionary history and ultimately investigate the predictability of evolution. Here, we focus on the oldfield mouse (Peromyscus polionotus), which occurs in the southeastern United States, where it exhibits considerable color variation. Dorsal coats range from dark brown in mainland mice to near white in mice inhabiting sandy beaches; this light pelage has evolved independently on Florida's Gulf and Atlantic coasts as camouflage from predators. To facilitate genomic analyses, we first generated a chromosome-level genome assembly of Peromyscus polionotus subgriseus. Next, in a uniquely variable mainland population (Peromyscus polionotus albifrons), we scored 23 pigment traits and performed targeted resequencing in 168 mice. We find that pigment variation is strongly associated with an ∼2-kb region ∼5 kb upstream of the Agouti signaling protein coding region. Using a reporter-gene assay, we demonstrate that this regulatory region contains an enhancer that drives expression in the dermis of mouse embryos during the establishment of pigment prepatterns. Moreover, extended tracts of homozygosity in this Agouti region indicate that the light allele experienced recent and strong positive selection. Notably, this same light allele appears fixed in both Gulf and Atlantic coast beach mice, despite these populations being separated by >1,000 km. Together, our results suggest that this identified Agouti enhancer allele has been maintained in mainland populations as standing genetic variation and from there, has spread to and been selected in two independent beach mouse lineages, thereby facilitating their rapid and parallel evolution.

Genetic variation in HIF-2α attenuates ventilatory sensitivity and carotid body growth in chronic hypoxia in high-altitude deer mice.

The gene encoding HIF-2α, Epas1, has experienced a history of natural selection in many high-altitude taxa, but the functional role of mutations in this gene is still poorly understood. We investigated the influence of the high-altitude variant of Epas1 in North American deer mice (Peromyscus maniculatus) on the control of breathing and carotid body growth during chronic hypoxia. We created hybrids between high- and low-altitude populations of deer mice to disrupt linkages between genetic loci so that the physiological effects of Epas1 alleles (Epas1H and Epas1L , respectively) could be examined on an admixed genomic background. In general, chronic hypoxia (4 weeks at 12 kPa O2 ) enhanced ventilatory chemosensitivity (assessed as the acute ventilatory response to hypoxia), increased total ventilation and arterial O2 saturation during progressive poikilocapnic hypoxia, and increased haematocrit and blood haemoglobin content across genotypes. However, the effects of chronic hypoxia on ventilatory chemosensitivity were attenuated in mice that were homozygous for the high-altitude Epas1 allele (Epas1H/H ). Carotid body growth and glomus cell hyperplasia, which was strongly induced in Epas1L/L mice in chronic hypoxia, was not observed in Epas1H/H mice. Epas1 genotype also modulated the effects of chronic hypoxia on metabolism and body temperature depression in hypoxia, but had no effects on haematological traits. These findings confirm the important role of HIF-2α in modulating ventilatory sensitivity and carotid body growth in chronic hypoxia, and show that genetic variation in Epas1 is responsible for evolved changes in the control of breathing and metabolism in high-altitude deer mice. KEY POINTS: High-altitude natives of many species have experienced natural selection on the gene encoding HIF-2α, Epas1, including high-altitude populations of deer mice. HIF-2α regulates ventilation and carotid body growth in hypoxia, and so the genetic variants in Epas1 in high-altitude natives may underlie evolved changes in control of breathing. Deer mice from controlled crosses between high- and low-altitude populations were used to examine the effects of Epas1 genotype on an admixed genomic background. The high-altitude variant was associated with reduced ventilatory chemosensitivity and carotid body growth in chronic hypoxia, but had no effects on haematology. The results help us better understand the genetic basis for the unique physiological phenotype of high-altitude natives.

A chromosomal inversion contributes to divergence in multiple traits between deer mouse ecotypes.

How locally adapted ecotypes are established and maintained within a species is a long-standing question in evolutionary biology. Using forest and prairie ecotypes of deer mice (Peromyscus maniculatus), we characterized the genetic basis of variation in two defining traits-tail length and coat color-and discovered a 41-megabase chromosomal inversion linked to both. The inversion frequency is 90% in the dark, long-tailed forest ecotype; decreases across a habitat transition; and is absent from the light, short-tailed prairie ecotype. We implicate divergent selection in maintaining the inversion at frequencies observed in the wild, despite high levels of gene flow, and explore fitness benefits that arise from suppressed recombination within the inversion. We uncover a key role for a large, previously uncharacterized inversion in the evolution and maintenance of classic mammalian ecotypes.

Burrowing deeper to uncover the genetic architecture of behavioral evolution.

Allele-specific expression analysis of hybrid mice provides new insights into the genetic substrates of behavioral evolution. As a complement to QTL mapping, this approach, described by Hu et al. in this issue of Cell Reports, holds promise for identifying causative regulatory loci that influence species-specific behavior.

cis-Regulatory changes in locomotor genes are associated with the evolution of burrowing behavior.

How evolution modifies complex, innate behaviors is largely unknown. Divergence in many morphological traits, and some behaviors, is linked to cis-regulatory changes in gene expression. Given this, we compare brain gene expression of two interfertile sister species of Peromyscus mice that show large and heritable differences in burrowing behavior. Species-level differential expression and allele-specific expression in F1 hybrids indicate a preponderance of cis-regulatory divergence, including many genes whose cis-regulation is affected by burrowing behavior. Genes related to locomotor coordination show the strongest signals of lineage-specific selection on burrowing-induced cis-regulatory changes. Furthermore, genetic markers closest to these candidate genes associate with variation in burrow shape in a genetic cross, suggesting an enrichment for loci affecting burrowing behavior near these candidate locomotor genes. Our results provide insight into how cis-regulated gene expression can depend on behavioral context and how this dynamic regulatory divergence between species may contribute to behavioral evolution.

The utility of a closed breeding colony of Peromyscus leucopus for dissecting complex traits.

Deermice of the genus Peromyscus are well suited for addressing several questions of biologist interest, including the genetic bases of longevity, behavior, physiology, adaptation, and their ability to serve as disease vectors. Here, we explore a diversity outbred approach for dissecting complex traits in Peromyscus leucopus, a nontraditional genetic model system. We take advantage of a closed colony of deer-mice founded from 38 individuals and subsequently maintained for ∼40-60 generations. From 405 low-pass short-read sequenced deermice we accurate impute genotypes at 16 million single nucleotide polymorphisms. Conditional on observed genotypes simulations were conducted in which three different sized quantitative trait loci contribute to a complex trait under three different genetic models. Using a stringent significance threshold power was modest, largely a function of the percent variation attributable to the simulated quantitative trait loci, with the underlying genetic model having only a subtle impact. We additionally simulated 2,000 pseudo-individuals, whose genotypes were consistent with those observed in the genotyped cohort and carried out additional power simulations. In experiments employing more than 1,000 mice power is high to detect quantitative trait loci contributing greater than 2.5% to a complex trait, with a localization ability of ∼100 kb. We finally carried out a Genome-Wide Association Study on two demonstration traits, bleeding time and body weight, and uncovered one significant region. Our work suggests that complex traits can be dissected in founders-unknown P. leucopus colony mice and similar colonies in other systems using easily obtained genotypes from low-pass sequencing.

Methylation studies in Peromyscus: aging, altitude adaptation, and monogamy.

DNA methylation-based biomarkers of aging have been developed for humans and many other mammals and could be used to assess how stress factors impact aging. Deer mice (Peromyscus) are long-living rodents that have emerged as an informative model to study aging, adaptation to extreme environments, and monogamous behavior. In the present study, we have undertaken an exhaustive, genome-wide analysis of DNA methylation in Peromyscus, spanning different species, stocks, sexes, tissues, and age cohorts. We describe DNA methylation-based estimators of age for different species of deer mice based on novel DNA methylation data generated on highly conserved mammalian CpGs measured with a custom array. The multi-tissue epigenetic clock for deer mice was trained on 3 tissues (tail, liver, and brain). Two human-Peromyscus clocks accurately measure age and relative age, respectively. We present CpGs and enriched pathways that relate to different conditions such as chronological age, high altitude, and monogamous behavior. Overall, this study provides a first step towards studying the epigenetic correlates of monogamous behavior and adaptation to high altitude in Peromyscus. The human-Peromyscus epigenetic clocks are expected to provide a significant boost to the attractiveness of Peromyscus as a biological model.

Genetic variation in haemoglobin is associated with evolved changes in breathing in high-altitude deer mice.

Physiological systems often have emergent properties but the effects of genetic variation on physiology are often unknown, which presents a major challenge to understanding the mechanisms of phenotypic evolution. We investigated whether genetic variants in haemoglobin (Hb) that contribute to high-altitude adaptation in deer mice (Peromyscus maniculatus) are associated with evolved changes in the control of breathing. We created F2 inter-population hybrids of highland and lowland deer mice to test for phenotypic associations of α- and ß-globin variants on a mixed genetic background. Hb genotype had expected effects on Hb-O2 affinity that were associated with differences in arterial O2 saturation in hypoxia. However, high-altitude genotypes were also associated with breathing phenotypes that should contribute to enhancing O2 uptake in hypoxia. Mice with highland α-globin exhibited a more effective breathing pattern, with highland homozygotes breathing deeper but less frequently across a range of inspired O2, and this difference was comparable to the evolved changes in breathing pattern in deer mouse populations native to high altitude. The ventilatory response to hypoxia was augmented in mice that were homozygous for highland ß-globin. The association of globin variants with variation in breathing phenotypes could not be recapitulated by acute manipulation of Hb-O2 affinity, because treatment with efaproxiral (a synthetic drug that acutely reduces Hb-O2 affinity) had no effect on breathing in normoxia or hypoxia. Therefore, adaptive variation in Hb may have unexpected effects on physiology in addition to the canonical function of this protein in circulatory O2 transport.

Linking genetic, morphological, and behavioural divergence between inland island and mainland deer mice.

The island syndrome hypothesis (ISH) stipulates that, as a result of local selection pressures and restricted gene flow, individuals from island populations should differ from individuals within mainland populations. Specifically, island populations are predicted to contain individuals that are larger, less aggressive, more sociable, and that invest more in their offspring. To date, tests of the ISH have mainly compared oceanic islands to continental sites, and rarely smaller spatial scales such as inland watersheds. Here, using a novel set of genome-wide SNP markers in wild deer mice (Peromyscus maniculatus) we conducted a genomic assessment of predictions underlying the ISH in an inland riverine island system: analysing island-mainland population structure, and quantifying heritability of phenotypes thought to underlie the ISH. We found clear genomic differentiation between the island and mainland populations and moderate to high marker-based heritability estimates for overall variation in traits previously found to differ in line with the ISH between mainland and island locations. FST outlier analyses highlighted 12 loci associated with differentiation between mainland and island populations. Together these results suggest that the island populations examined are on independent evolutionary trajectories, the traits considered have a genetic basis (rather than phenotypic variation being solely due to phenotypic plasticity). Coupled with the previous results showing significant phenotypic differentiation between the island and mainland groups in this system, this study suggests that the ISH can hold even on a small spatial scale.

Genomic variation in captive deer mouse (Peromyscus maniculatus) populations.

BACKGROUND: Deer mice (genus Peromyscus) are the most common rodents in North America. Despite the availability of reference genomes for some species, a comprehensive database of polymorphisms, especially in those maintained as living stocks and distributed to academic investigators, is missing. In the present study we surveyed two populations of P. maniculatus that are maintained at the Peromyscus Genetic Stock Center (PGSC) for polymorphisms across their 2.5 × 109 bp genome. RESULTS: High density of variation was identified, corresponding to one SNP every 55 bp for the high altitude stock (SM2) or 207 bp for the low altitude stock (BW) using snpEff (v4.3). Indels were detected every 1157 bp for BW or 311 bp for SM2. The average Watterson estimator for the BW and SM2 populations is 248813.70388 and 869071.7671 respectively. Some differences in the distribution of missense, nonsense and silent mutations were identified between the stocks, as well as polymorphisms in genes associated with inflammation (NFATC2), hypoxia (HIF1a) and cholesterol metabolism (INSIG1) and may possess value in modeling pathology. CONCLUSIONS: This genomic resource, in combination with the availability of P. maniculatus from the PGSC, is expected to promote genetic and genomic studies with this animal model.

Broad Concordance in the Spatial Distribution of Adaptive and Neutral Genetic Variation across an Elevational Gradient in Deer Mice.

When species are continuously distributed across environmental gradients, the relative strength of selection and gene flow shape spatial patterns of genetic variation, potentially leading to variable levels of differentiation across loci. Determining whether adaptive genetic variation tends to be structured differently than neutral variation along environmental gradients is an open and important question in evolutionary genetics. We performed exome-wide population genomic analysis on deer mice sampled along an elevational gradient of nearly 4,000 m of vertical relief. Using a combination of selection scans, genotype-environment associations, and geographic cline analyses, we found that a large proportion of the exome has experienced a history of altitude-related selection. Elevational clines for nearly 30% of these putatively adaptive loci were shifted significantly up- or downslope of clines for loci that did not bear similar signatures of selection. Many of these selection targets can be plausibly linked to known phenotypic differences between highland and lowland deer mice, although the vast majority of these candidates have not been reported in other studies of highland taxa. Together, these results suggest new hypotheses about the genetic basis of physiological adaptation to high altitude, and the spatial distribution of adaptive genetic variation along environmental gradients.