Extreme High-Elevation Mammal Surveys Reveal Unexpectedly High Upper Range Limits of Andean Mice.

AbstractIn the world's highest mountain ranges, uncertainty about the upper elevational range limits of alpine animals represents a critical knowledge gap regarding the environmental limits of life and presents a problem for detecting range shifts in response to climate change. Here we report results of mountaineering mammal surveys in the Central Andes, which led to the discovery of multiple species of mice living at extreme elevations that far surpass previously assumed range limits for mammals. We livetrapped small mammals from ecologically diverse sites spanning >6,700 m of vertical relief, from the desert coast of northern Chile to the summits of the highest volcanoes in the Andes. We used molecular sequence data and whole-genome sequence data to confirm the identities of species that represent new elevational records and to test hypotheses regarding species limits. These discoveries contribute to a new appreciation of the environmental limits of vertebrate life.

Diversification, Introgression, and Rampant Cytonuclear Discordance in Rocky Mountains Chipmunks (Sciuridae: Tamias).

Evidence from natural systems suggests that hybridization between animal species is more common than traditionally thought, but the overall contribution of introgression to standing genetic variation within species remains unclear for most animal systems. Here, we use targeted exon capture to sequence thousands of nuclear loci and complete mitochondrial genomes from closely related chipmunk species in the Tamias quadrivittatus group that are distributed across the Great Basin and the central and southern Rocky Mountains of North America. This recent radiation includes six overlapping, ecologically distinct species (Tamias canipes, Tamias cinereicollis, Tamias dorsalis, T. quadrivittatus, Tamias rufus, and Tamias umbrinus) that show evidence for widespread introgression across species boundaries. Such evidence has historically been derived from a handful of markers, typically focused on mitochondrial loci, to describe patterns of introgression; consequently, the extent of introgression of nuclear genes is less well characterized. We conducted a series of phylogenomic and species-tree analyses to resolve the phylogeny of six species in this group. In addition, we performed several population-genomic analyses to characterize nuclear genomes and infer coancestry among individuals. Furthermore, we used emerging quartets-based approaches to simultaneously infer the species tree (SVDquartets) and identify introgression (HyDe). We found that, in spite of rampant introgression of mitochondrial genomes between some species pairs (and sometimes involving up to three species), there appears to be little to no evidence for nuclear introgression. These findings mirror other genomic results where complete mitochondrial capture has occurred between chipmunk species in the absence of appreciable nuclear gene flow. The underlying causes of recurrent massive cytonuclear discordance remain unresolved in this group but mitochondrial DNA appears highly misleading of population histories as a whole. Collectively, it appears that chipmunk species boundaries are largely impermeable to nuclear gene flow and that hybridization, while pervasive with respect to mtDNA, has likely played a relatively minor role in the evolutionary history of this group. [Cytonuclear discordance; hyridization; introgression, phylogenomics; SVDquartets; Tamias.].

Molecular phylogenetics and historical biogeography amid shifting continents in the cockles and giant clams (Bivalvia: Cardiidae).

Reconstructing historical biogeography of the marine realm is complicated by indistinct barriers and, over deeper time scales, a dynamic landscape shaped by plate tectonics. Here we present the most extensive examination of model-based historical biogeography among marine invertebrates to date. We conducted the largest phylogenetic and molecular clock analyses to date for the bivalve family Cardiidae (cockles and giant clams) with three unlinked loci for 110 species representing 37 of the 50 genera. Ancestral ranges were reconstructed using the dispersal-extinction-cladogenesis (DEC) method with a time-stratified paleogeographic model wherein dispersal rates varied with shifting tectonics. Results were compared to previous classifications and the extensive paleontological record. Six of the eight prior subfamily groupings were found to be para- or polyphyletic. Cardiidae originated and subsequently diversified in the tropical Indo-Pacific starting in the Late Triassic. Eastern Atlantic species were mainly derived from the tropical Indo-Mediterranean region via the Tethys Sea. In contrast, the western Atlantic fauna was derived from Indo-Pacific clades. Our phylogenetic results demonstrated greater concordance with geography than did previous phylogenies based on morphology. Time-stratifying the DEC reconstruction improved the fit and was highly consistent with paleo-ocean currents and paleogeography. Lastly, combining molecular phylogenetics with a rich and well-documented fossil record allowed us to test the accuracy and precision of biogeographic range reconstructions.

Local adaptation, plasticity, and evolved resistance to hypoxic cold stress in high-altitude deer mice.

A fundamental question in evolutionary biology concerns the relative contributions of phenotypic plasticity vs. local adaptation (genotypic specialization) in enabling wide-ranging species to inhabit diverse environmental conditions. Here we conduct a long-term hypoxia acclimation experiment to assess the relative roles of local adaptation and plasticity in enabling highland and lowland deer mice ( Peromyscus maniculatus ) to sustain aerobic thermogenesis at progressively increasing elevations. We assessed the relative physiological performance capacities of highland and lowland natives as they were exposed to progressive, stepwise increases in hypoxia, simulating the gradual ascent from sea level to an elevation of 6000 m. The final elevation of 6000 m far exceeds the highest attainable elevations within the species' range, and therefore tests the animals' ability to tolerate levels of hypoxia that surpass the prevailing conditions within their current distributional limits. Our results demonstrate that highland natives exhibit superior thermogenic capacities at the most severe levels of hypoxia, suggesting that the species' broad fundamental niche and its ability to inhabit such a broad range of elevational zones is attributable to a combination of genetically based local adaptation and plasticity. Transcriptomic and physiological measurements identify evolved changes in the acclimation response to hypoxia that contribute to the enhanced thermogenic capacity of highland natives. SIGNIFICANCE STATEMENT: In species that are distributed across steep environmental gradients, the ability to inhabit a broad range of conditions may be attributable to local adaptation and/or a generalized acclimatization capacity (phenotypic plasticity). By experimentally acclimating highland and lowland deer mice ( Peromyscus maniculatus ) to progressively increasing levels of hypoxia during a simulated ascent to 6000 m, we assessed the relative roles of evolved and plastic changes in thermogenic capacity. At especially severe levels of hypoxia, the superior thermogenic performance of highland natives relative to lowland conspecifics suggests that the broad fundamental niche of deer mice is largely attributable to local adaptation to different elevational zones, including evolved plasticity in gene expression and respiratory traits.

Genomic resolution of cryptic species diversity in chipmunks.

Discovery of cryptic species is essential to understand the process of speciation and assessing the impacts of anthropogenic stressors. Here, we used genomic data to test for cryptic species diversity within an ecologically well-known radiation of North American rodents, western chipmunks (Tamias). We assembled a de novo reference genome for a single species (Tamias minimus) combined with new and published targeted sequence-capture data for 21,551 autosomal and 493 X-linked loci sampled from 121 individuals spanning 22 species. We identified at least two cryptic lineages corresponding with an isolated subspecies of least chipmunk (T. minimus grisescens) and with a restricted subspecies of the yellow-pine chipmunk (Tamias amoenus cratericus) known only from around the extensive Craters of the Moon lava flow. Additional population-level sequence data revealed that the so-called Crater chipmunk is a distinct species that is abundant throughout the coniferous forests of southern Idaho. This cryptic lineage does not appear to be most closely related to the ecologically and phenotypically similar yellow-pine chipmunk but does show evidence for recurrent hybridization with this and other species.

Expression Differentiation Is Constrained to Low-Expression Proteins over Ecological Timescales.

Protein expression level is one of the strongest predictors of protein sequence evolutionary rate, with high-expression protein sequences evolving at slower rates than low-expression protein sequences largely because of constraints on protein folding and function. Expression evolutionary rates also have been shown to be negatively correlated with expression level across human and mouse orthologs over relatively long divergence times (i.e., ∼100 million years). Long-term evolutionary patterns, however, often cannot be extrapolated to microevolutionary processes (and vice versa), and whether this relationship holds for traits evolving under directional selection within a single species over ecological timescales (i.e., <5000 years) is unknown and not necessarily expected. Expression is a metabolically costly process, and the expression level of a particular protein is predicted to be a tradeoff between the benefit of its function and the costs of its expression. Selection should drive the expression level of all proteins close to values that maximize fitness, particularly for high-expression proteins because of the increased energetic cost of production. Therefore, stabilizing selection may reduce the amount of standing expression variation for high-expression proteins, and in combination with physiological constraints that may place an upper bound on the range of beneficial expression variation, these constraints could severely limit the availability of beneficial expression variants. To determine whether rapid-expression evolution was restricted to low-expression proteins owing to these constraints on highly expressed proteins over ecological timescales, we compared venom protein expression levels across mainland and island populations for three species of pit vipers. We detected significant differentiation in protein expression levels in two of the three species and found that rapid-expression differentiation was restricted to low-expression proteins. Our results suggest that various constraints on high-expression proteins reduce the availability of beneficial expression variants relative to low-expression proteins, enabling low-expression proteins to evolve and potentially lead to more rapid adaptation.