Functional genomics of adaptation to hypoxic cold-stress in high-altitude deer mice: transcriptomic plasticity and thermogenic performance

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 6,000 m. The final elevation of 6,000 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 genetically based local adaptation, including evolved changes in plasticity. Transcriptomic and physiological measurements identify evolved changes in the acclimation response to hypoxia that contribute to the enhanced thermogenic capacity of highland natives.

Metabolic recovery from submaximal exercise in hypoxia acclimated high altitude deer mice (Peromyscus maniculatus)

Animals living at high-altitude are faced with unremitting low oxygen availability. This can make it difficult to perform daily tasks that require increases in aerobic metabolism. An activity important for survival is aerobic locomotion, and the rapid recovery of muscle metabolism post exercise. Past work shows that hypoxia acclimated high-altitude mice (Peromyscus maniculatus) have a greater reliance on carbohydrates to power exercise than low altitude mice. However, it is unclear how quickly after aerobic exercise these mice can recovery and replenish muscle glycogen stores. The gastrocnemius muscle of high-altitude deer mice has a more aerobic phenotype and a greater capacity to oxidize lipids than low altitude deer mice. This suggests that high altitude mice may recover more rapidly from exercise than their lowland counterparts due to a greater capacity to support glycogen replenishment using intramuscular triglycerides (IMTG). To explore this possibility, we used low- and high-altitude native deer mice born and raised in common lab conditions and acclimated to chronic hypoxia. We determined changes in oxygen consumption following 15 min of aerobic exercise in 12% O2 and sampled skeletal muscles and liver at various time points during recovery to examine changes in key metabolites, including glycogen and IMTG. We found depletion in glycogen stores during exercise only in lowlanders, which returned to resting levels following 90 min of recovery. In contrast, IMTG did not change significantly with exercise or during recovery in either population. These data suggest that exercise recovery is influenced by altitude ancestry in deer mice.

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.

Genomic and Phenotypic Adaptations of Rattus tanezumi to Cold Limit Its Further Northward Expansion and Range Overlap with R. norvegicus

Global climate change has led to shifts in the distribution ranges of many terrestrial species, promoting their migration from lower altitudes or latitudes to higher ones. Meanwhile, successful invaders have developed genetic adaptations enabling the colonization of new environments. Over the past 40 years, Rattus tanezumi (RT) has expanded into northern China (Northwest and North China) from its southern origins. We studied the cold adaptation of RT and its potential for northward expansion by comparing it with sympatric Rattus norvegicus (RN), which is well adapted to cold regions. Through population genomic analysis, we revealed that the invading RT rats have split into three distinct populations: the North, Northwest, and Tibetan populations. The first two populations exhibited high genetic diversity, while the latter population showed remarkably low genetic diversity. These rats have developed various genetic adaptations to cold, arid, hypoxic, and high-UV conditions. Cold acclimation tests revealed divergent thermoregulation between RT and RN. Specifically, RT exhibited higher brown adipose tissue activity and metabolic rates than did RN. Transcriptome analysis highlighted changes in genes regulating triglyceride catabolic processes in RT, including Apoa1 and Apoa4, which were upregulated, under selection and associated with local adaptation. In contrast, RN showed changes in carbohydrate metabolism genes. Despite the cold adaptation of RT, we observed genotypic and phenotypic constraints that may limit its ability to cope with severe low temperatures farther north. Consequently, it is less likely that RT rats will invade and overlap with RN rats in farther northern regions.

Impacts of dietary fat on multi tissue gene expression in the desert-adapted cactus mouse

Understanding the relationship between dietary fat and physiological responses is crucial in species adapted to arid environments where water scarcity is common. In this study, we present a comprehensive exploration of gene expression across five tissues (kidney, liver, lung, gastrointestinal tract, and hypothalamus) and 19 phenotypic measurements, investigating the effects of dietary fat in the desert-adapted cactus mouse ( Peromyscus eremicus ). We show impacts on immune function, circadian gene regulation, and mitochondrial function for mice fed a lower-fat diet compared to mice fed a higher-fat diet. In arid environments with severe water scarcity, even subtle changes in organismal health and water balance can affect physical performance, potentially impacting survival and reproductive success. The study sheds light on the complex interplay between diet, physiological processes, and environmental adaptation, providing valuable insights into the multifaceted impacts of dietary choices on organismal well-being and adaptation strategies in arid habitats.

Evolved changes in phenotype across skeletal muscles in deer mice native to high altitude

The cold and hypoxic conditions at high altitude necessitate high metabolic O2 demands to support thermogenesis while hypoxia reduces O2 availability. Skeletal muscles play key roles in thermogenesis, but our appreciation of muscle plasticity and adaptation at high altitude has been hindered by past emphasis on only a small number of muscles. We examined this issue in deer mice (Peromyscus maniculatus). Mice derived from both high-altitude and low-altitude populations were born and raised in captivity and then acclimated as adults to normoxia or hypobaric hypoxia (12 kPa O2 for 6-8 wk). Maximal activities of citrate synthase (CS), cytochrome c oxidase (COX), β-hydroxyacyl-CoA dehydrogenase (HOAD), hexokinase (HK), pyruvate kinase (PK), and lactate dehydrogenase (LDH) were measured in 20 muscles involved in shivering, locomotion, body posture, ventilation, and mastication. Principal components analysis revealed an overall difference in muscle phenotype between populations but no effect of hypoxia acclimation. High-altitude mice had greater activities of mitochondrial enzymes and/or lower activities of PK or LDH across many (but not all) respiratory, limb, core and mastication muscles compared with low-altitude mice. In contrast, chronic hypoxia had very few effects across muscles. Further examination of CS in the gastrocnemius showed that population differences in enzyme activity stemmed from differences in protein abundance and mRNA expression but not from population differences in CS amino acid sequence. Overall, our results suggest that evolved increases in oxidative capacity across many skeletal muscles, at least partially driven by differences in transcriptional regulation, may contribute to high-altitude adaptation in deer mice.NEW & NOTEWORTHY Most previous studies of muscle plasticity and adaptation in high-altitude environments have focused on a very limited number of skeletal muscles. Comparing high-altitude versus low-altitude populations of deer mice, we show that a large number of muscles involved in shivering, locomotion, body posture, ventilation, and mastication exhibit greater mitochondrial enzyme activities in the high-altitude population. Therefore, evolved increases in mitochondrial oxidative capacity across skeletal muscles contribute to high-altitude adaptation.

Coping with extremes: Alternations in diet, gut microbiota, and hepatic metabolic functions in a highland passerine

In wild animals, diet and gut microbiota interactions are critical moderators of metabolic functions and are highly contingent on habitat conditions. Challenged by the extreme conditions of high-altitude environments, the strategies implemented by highland animals to adjust their diet and gut microbial composition and modulate their metabolic substrates remain largely unexplored. By employing a typical human commensal species, the Eurasian tree sparrow (Passer montanus, ETS), as a model species, we studied the differences in diet, digestive tract morphology and enzyme activity, gut microbiota, and metabolic energy profiling between highland (the Qinghai-Tibet Plateau, QTP; 3230 m) and lowland (Shijiazhuang, Hebei; 80 m) populations. Our results showed that highland ETSs had enlarged digestive organs and longer small intestinal villi, while no differences in key digestive enzyme activities were observed between the two populations. The 18S rRNA sequencing results revealed that the dietary composition of highland ETSs were more animal-based and less plant-based than those of the lowland ones. Furthermore, 16S rRNA sequencing results suggested that the intestinal microbial communities were structurally segregated between populations. PICRUSt metagenome predictions further indicated that the expression patterns of microbial genes involved in material and energy metabolism, immune system and infection, and xenobiotic biodegradation were strikingly different between the two populations. Analysis of liver metabolomics revealed significant metabolic differences between highland and lowland ETSs in terms of substrate utilization, as well as distinct sex-specific alterations in glycerophospholipids. Furthermore, the interplay between diet, liver metabolism, and gut microbiota suggests a dietary shift resulting in corresponding changes in gut microbiota and metabolic functions. Our findings indicate that highland ETSs have evolved to optimize digestion and absorption, rely on more protein-rich foods, and possess gut microbiota tailored to their dietary composition, likely adaptive physiological and ecological strategies adopted to cope with extreme highland environments.

Whole-genome sequencing reveals adaptations of hairy-footed jerboas (Dipus, Dipodidae) to diverse desert environments

BACKGROUND

Environmental conditions vary among deserts across the world, spanning from hyper-arid to high-elevation deserts. However, prior genomic studies on desert adaptation have focused on desert and non-desert comparisons overlooking the complexity of conditions within deserts. Focusing on the adaptation mechanisms to diverse desert environments will advance our understanding of how species adapt to extreme desert environments. The hairy-footed jerboas are well adapted to diverse desert environments, inhabiting high-altitude arid regions, hyper-arid deserts, and semi-deserts, but the genetic basis of their adaptation to different deserts remains unknown.

RESULTS

Here, we sequenced the whole genome of 83 hairy-footed jerboas from distinct desert zones in China to assess how they responded under contrasting conditions. Population genomics analyses reveal the existence of three species in hairy-footed jerboas distributed in China: Dipus deasyi, Dipus sagitta, and Dipus sowerbyi. Analyses of selection between high-altitude desert (elevation ≥ 3000m) and low-altitude desert (< 500m) populations identified two strongly selected genes, ATR and HIF1AN, associated with intense UV radiation and hypoxia in high-altitude environments. A number of candidate genes involved in energy and water homeostasis were detected in the comparative genomic analyses of hyper-arid desert (average annual precipitation < 70mm) and arid desert (< 200mm) populations versus semi-desert (> 360mm) populations. Hyper-arid desert animals also exhibited stronger adaptive selection in energy homeostasis, suggesting water and resource scarcity may be the main drivers of desert adaptation in hairy-footed jerboas.

CONCLUSIONS

Our study challenges the view of deserts as homogeneous environments and shows that distinct genomic adaptations can be found among desert animals depending on their habitats.

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.

Markhor-derived Introgression of a Genomic Region Encompassing PAPSS2 Confers High-altitude Adaptability in Tibetan Goats

Understanding the genetic mechanism of how animals adapt to extreme conditions is fundamental to determine the relationship between molecular evolution and changing environments. Goat is one of the first domesticated species and has evolved rapidly to adapt to diverse environments, including harsh high-altitude conditions with low temperature and poor oxygen supply but strong ultraviolet radiation. Here, we analyzed 331 genomes of domestic goats and wild caprid species living at varying altitudes (high > 3000 m above sea level and low < 1200 m), along with a reference-guided chromosome-scale assembly (contig-N50: 90.4 Mb) of a female Tibetan goat genome based on PacBio HiFi long reads, to dissect the genetic determinants underlying their adaptation to harsh conditions on the Qinghai-Tibetan Plateau (QTP). Population genomic analyses combined with genome-wide association studies (GWAS) revealed a genomic region harboring the 3'-phosphoadenosine 5'-phosphosulfate synthase 2 (PAPSS2) gene showing strong association with high-altitude adaptability (PGWAS = 3.62 × 10-25) in Tibetan goats. Transcriptomic data from 13 tissues revealed that PAPSS2 was implicated in hypoxia-related pathways in Tibetan goats. We further verified potential functional role of PAPSS2 in response to hypoxia in PAPSS2-deficient cells. Introgression analyses suggested that the PAPSS2 haplotype conferring the high-altitude adaptability in Tibetan goats originated from a recent hybridization between goats and a wild caprid species, the markhor (Capra falconeri). In conclusion, our results uncover a hitherto unknown contribution of PAPSS2 to high-altitude adaptability in Tibetan goats on QTP, following interspecific introgression and natural selection.

Modified Metabolism and Response to UV Radiation: Gene Expression Variations Along an Elevational Gradient in the Asiatic Toad (Bufo gargarizans)

High-elevation adaptation provides an excellent system for examining adaptive evolution, and adaptive variations may manifest at gene expression or any other phenotypic levels. We examined gene expression profiles of Asiatic toads (Bufo gargarizans) along an elevational gradient from both wild and common-garden acclimated populations. Asiatic toads originated from high altitudes have distinctive gene expression patterns. We identified 18 fixed differentially expressed genes (DEGs), which are different in both wild and acclimated samples, and 1217 plastic DEGs, which are different among wild samples. The expression levels of most genes were linearly correlated with altitude gradient and down-regulated in high-altitude populations. Expression variations of several genes associated with metabolic process are fixed, and we also identified a co-expression module that is significantly different between acclimated populations and has functions related to DNA repair. The differential expression of the vast majority genes, however, are due to phenotypic plasticity, revealing the highly plastic nature of gene expression variations. Expression modification of some specific genes related to metabolism and response to UV radiation play crucial role in adaptation to high altitude for Asiatic toads. Common-garden experiments are essential for evaluating adaptive evolution of natural populations.

Genetic variation in HIF-2α attenuates ventilatory sensitivity and carotid body growth in chronic hypoxia 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, 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.

ABSTRACT

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 are still poorly understood. We investigated the influence of the high-altitude variant of Epas1 in North American deer mice (Peromyscus maniculatus) on 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 physiological effects of Epas1 alleles (Epas1^H and Epas1^L , respectively) could be examined on an admixed genomic background. In general, chronic hypoxia (4 weeks at 12 kPa O₂ ) enhanced ventilatory chemosensitivity (assessed as the acute ventilatory response to hypoxia), increased total ventilation and arterial O₂ saturation during progressive poikilocapnic hypoxia, and increased haematocrit and blood haemoglobin content across genotypes. However, effects of chronic hypoxia on ventilatory chemosensitivity were attenuated in mice that were homozygous for the high-altitude Epas1 allele (Epas1^H/H ). Carotid body growth and glomus cell hyperplasia, which was strongly induced in Epas1^L/L mice in chronic hypoxia, was not observed in Epas1^H/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. Abstract figure legend ventilation and carotid body growth in hypoxia, so we investigated the role genetic variants in Epas1 in highaltitude deer mice on the control of breathing. In the lab, hybrids between high- and lowaltitude populations of deer mice were created to disrupt linkages between genetic loci so physiological effects of Epas1 alleles (Epas1H and Epas1L, respectively) could be examined on an admixed genomic background. The high-altitude variant was associated with reduced ventilatory chemosensitivity and carotid body growth after 4 weeks of chronic hypoxia, compared to mice homozygous for the low-altitude allele (Epas1LL). These results help us better understand the genetic basis for the unique physiological phenotype of high-altitude natives. This article is protected by copyright. All rights reserved.

Challenging ecogeographical rules: Phenotypic variation in the Mountain Treeshrew (Tupaia montana) along tropical elevational gradients

Bergmann’s and Allen’s rules were defined to describe macroecological patterns across latitudinal gradients. Bergmann observed a positive association between body size and latitude for endothermic species while Allen described shorter appendages as latitude increases. Almost two centuries later, there is still ongoing discussion about these patterns. Temperature, the common variable in these two rules, varies predictably across both latitude and elevation. Although these rules have been assessed extensively in mammals across latitude, particularly in regions with strong seasonality, studies on tropical montane mammals are scarce. We here test for these patterns and assess the variation of several other locomotory, diet-associated, body condition, and thermoregulatory traits across elevation in the Mountain Treeshrew (Tupaia montana) on tropical mountains in Borneo. Based on morphological measurements from both the field and scientific collections, we found a complex pattern: Bergmann’s rule was not supported in our tropical mountain system, since skull length, body size, and weight decreased from the lowest elevations (<1000 m) to middle elevations (2000–2500 m), and then increased from middle elevations to highest elevations. Allen’s rule was supported for relative tail length, which decreased with elevation, but not for ear and hindfoot length, with the former remaining constant and the latter increasing with elevation. This evidence together with changes in presumed diet-related traits (rostrum length, zygomatic breadth and upper tooth row length) along elevation suggest that selective pressures other than temperature, are playing a more important role shaping the morphological variation across the distribution of the Mountain Treeshrew. Diet, food acquisition, predation pressure, and/or intra- and inter-specific competition, are some of the potential factors driving the phenotypic variation of this study system. The lack of variation in body condition might suggest local adaptation of this species across its elevational range, perhaps due to generalist foraging strategies. Finally, a highly significant temporal effect was detected in several traits but not in others, representing the first phenotypic variation temporal trends described on treeshrews.

Phenotypic plasticity in the energy metabolism of a small Andean rodent: Effect of short-term thermal acclimation and developmental conditions

The study of phenotypic variation within species in response to different environments is a central issue in evolutionary and ecological physiology. Particularly, ambient temperature is one of the most important factors modulating interactions between animals and their environment. Phyllotis xanthopygus, a small Andean rodent, exhibits intraspecific differences along an altitudinal gradient in traits relevant to energy balance that persist after acclimation to common experimental temperatures. Therefore, we aim to explore geographic variations in energetic traits of P. xanthopygus and to assess the contribution of phenotypic plasticity to population differences. We compared metabolic rate and thermal conductance in response to different acclimation temperatures in animals collected at distinct altitudes (F0 generation) and in their offspring, born and raised under common-garden conditions (F1 generation). We found intraspecific differences in resting metabolic rate (RMR) of animals collected at different altitudes that were no longer evident in the F1 generation. Furthermore, although both generations showed the same pattern of RMR flexibility in response to acclimation temperature, its magnitude was lower for the F1 individuals. This suggests that developmental conditions affect the short-term acclimation capacity of this trait during adulthood. On the other hand, thermal conductance (C) showed irreversible plasticity, as animals raised in the laboratory at stable warm conditions had a relatively higher C than the animals from the field, showing no adjustments to thermal acclimation during adulthood in either group. In sum, our results support the hypothesis that the developmental environment shapes energetic traits, emphasizing the relevance of incorporating ontogeny in physiological studies.

Altitude acclimatization, hemoglobin-oxygen affinity, and circulatory oxygen transport in hypoxia

In mammals and other air-breathing vertebrates that live at high altitude, adjustments in convective O₂ transport via changes in blood hemoglobin (Hb) content and/or Hb-O₂ affinity can potentially mitigate the effects of arterial hypoxemia. However, there are conflicting views about the optimal values of such traits in hypoxia, partly due to the intriguing observation that hypoxia-induced acclimatization responses in humans and other predominantly lowland mammals are frequently not aligned in the same direction as evolved phenotypic changes in high-altitude natives. Here we review relevant theoretical and empirical results and we highlight experimental studies of rodents and humans that provide insights into the combination of hematological changes that help attenuate the decline in aerobic performance in hypoxia. For a given severity of hypoxia, experimental results suggest that optimal values for hematological traits are conditional on the states of other interrelated phenotypes that govern different steps in the O₂-transport pathway.

A Multilevel Assessment of Plasticity in Response to High-Altitude Environment for Agama Lizards

Upslope range shifting has been documented in diverse species in response to global warming. Plasticity, which refers to the ability of organisms to alter their phenotypes in changing environments, is crucial for the survival of those that newly migrated to a high-altitude environment. The scope and mechanisms of plasticity across biological levels, however, have rarely been examined. We used two agama lizards (genus Phrynocephalus) as model systems and a transplant experiment to comprehensively assess their plasticity on multiple organization levels. Two low-altitude (934 m) agama species, Phrynocephalus axillaris (oviparous) and P. forsythii (viviparous), were transplanted to a high-altitude site (3,400 m). After acclimation for 6 weeks in seminatural enclosures, plasticity was measured from bite force, tail display behavior, gene expression, and metabolome. Both lizards were capable of acclimating to the high-altitude environment without sacrificing their performance in bite force, but they also showed high plasticity in tail display behavior by either decreasing the intensity of a specific display component (P. forsythii) or by the trade-off between display components (P. axillaris). Genes and metabolites associated with lipids, especially fatty acid metabolism, exhibited significant differentiation in expression, compared to individuals from their native habitats. Improved fatty acid storage and metabolism appeared to be a common response among animals at high altitudes. Despite distinct reproductive modes that may differ in response to physiological pressure, the two lizards demonstrated high concordance in plasticity when they faced a novel environment at high altitudes. Taken together, lizards likely acclimate to high-altitude environments by reducing behavioral activity and increasing energy efficiency after range shifting. Our results provide new insights into our understanding of phenotypic plasticity and its importance in today’s changing climate.

The genomic basis of high-elevation adaptation in wild house mice (Mus musculus domesticus) from South America

Understanding the genetic basis of environmental adaptation in natural populations is a central goal in evolutionary biology. The conditions at high elevation, particularly the low oxygen available in the ambient air, impose a significant and chronic environmental challenge to metabolically active animals with lowland ancestry. To understand the process of adaptation to these novel conditions and to assess the repeatability of evolution over short timescales, we examined the signature of selection from complete exome sequences of house mice (Mus musculus domesticus) sampled across two elevational transects in the Andes of South America. Using phylogenetic analysis, we show that house mice colonized high elevations independently in Ecuador and Bolivia. Overall, we found distinct responses to selection in each transect and largely nonoverlapping sets of candidate genes, consistent with the complex nature of traits that underlie adaptation to low oxygen availability (hypoxia) in other species. Nonetheless, we also identified a small subset of the genome that appears to be under parallel selection at the gene and SNP levels. In particular, three genes (Col22a1, Fgf14, and srGAP1) bore strong signatures of selection in both transects. Finally, we observed several patterns that were common to both transects, including an excess of derived alleles at high elevation, and a number of hypoxia-associated genes exhibiting a threshold effect, with a large allele frequency change only at the highest elevations. This threshold effect suggests that selection pressures may increase disproportionately at high elevations in mammals, consistent with observations of some high-elevation diseases in humans.

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.

Coping with extremes: High-altitude sparrows enhance metabolic and thermogenic capacities in the pectoralis muscle and suppress in the liver relative to their lowland counterparts

Animals living at high altitudes are challenged by the extreme environmental conditions of cold temperature and hypobaric hypoxia. It is not well understood how high-altitude birds enhance the capacity of metabolic thermogenesis and allocate metabolic capacity in different organs to maximize survival in extreme conditions of a cold winter. The Qinghai-Tibet Plateau (QTP) is the largest and highest plateau globally, offering a natural laboratory for investigating coping mechanisms of organisms inhabiting extreme environments. To understand the adaptive strategies in the morphology and physiology of small songbirds on the QTP, we compared plasma triiodothyronine (T₃), pectoralis muscle mitochondrial cytochrome c oxidase (COX) and state IV capacities, the expression of peroxisome proliferator-activated receptor γ coactivator α (PGC-1α), adenine nucleotide translocase (ANT), uncoupling protein (UCP), and adenosine monophosphate-dependent kinase (AMPK) α1 mRNA in the pectoralis and liver of Eurasian tree sparrows (Passer montanus) from high-altitude (3,230 m), medium-altitude (1400 m), and low-altitude (80 m) regions. Our results showed that high-altitude sparrows had greater body masses, longer wings and tarsometatarsi, but comparable bill lengths relative to medium- and low-altitude individuals. High-altitude sparrows had higher plasma T₃ levels and pectoralis muscle mitochondrial COX capacities than their lowland counterparts. They also upregulated the pectoralis muscle mRNA expression of UCP, PGC-1α, and ANT proteins relative to low-altitude sparrows. Unlike pectoralis, high-altitude sparrows significantly down-regulated hepatic AMPKα1 and ANT protein expression as compared with their lowland counterparts. Our results contribute to understanding the morphological, biochemical, and molecular adaptations in free-living birds to cope with the cold seasons in the extreme environment of the QTP.

Increased Reliance on Carbohydrates for Aerobic Exercise in Highland Andean Leaf-Eared Mice, but Not in Highland Lima Leaf-Eared Mice

Exercise is an important performance trait in mammals and variation in aerobic capacity and/or substrate allocation during submaximal exercise may be important for survival at high altitude. Comparisons between lowland and highland populations is a fruitful approach to understanding the mechanisms for altitude differences in exercise performance. However, it has only been applied in very few highland species. The leaf-eared mice (LEM, genus Phyllotis) of South America are a promising taxon to uncover the pervasiveness of hypoxia tolerance mechanisms. Here we use lowland and highland populations of Andean and Lima LEM (P. andium and P. limatus), acclimated to common laboratory conditions, to determine exercise-induced maximal oxygen consumption (V˙O₂max), and submaximal exercise metabolism. Lowland and highland populations of both species showed no difference in V˙O₂max running in either normoxia or hypoxia. When run at 75% of V˙O₂max, highland Andean LEM had a greater reliance on carbohydrate oxidation to power exercise. In contrast, highland Lima LEM showed no difference in exercise fuel use compared to their lowland counterparts. The higher carbohydrate oxidation seen in highland Andean LEM was not explained by maximal activities of glycolytic enzymes in the gastrocnemius muscle, which were equivalent to lowlanders. This result is consistent with data on highland deer mouse populations and suggests changes in metabolic regulation may explain altitude differences in exercise performance.

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.

Genomic analysis unveils mechanisms of northward invasion and signatures of plateau adaptation in the Asian house rat

The Asian house rat (AHR), Rattus tanezumi, has recently invaded the northern half of China. The AHR is a highly adaptive rat species that has also successfully conquered the Qinghai-Tibet Plateau (QTP) and replaced the brown rat (BR), R. norvegicus, at the edge of the QTP. Here, we assembled a draft genome of the AHR and explored the mechanisms of its northward invasion and the genetic basis underlying plateau adaptation in this species. Population genomic analyses revealed that the northwardly invasive AHRs consisted of two independent and genetically distinct populations which might result from multiple independent primary invasion events. One invasive population exhibited reduced genetic diversity and distinct population structure compared with its source population, while the other displayed preserved genetic polymorphisms and little genetic differentiation from its source population. Genes involved in G-protein coupled receptors and carbohydrate metabolism may contribute to the local adaptation of northern AHRs. In particular, RTN4 was identified as a key gene for AHRs in the QTP that favours adaptation to high-altitude hypoxia. Coincidently, the physiological performance and transcriptome profiles of hypoxia-exposed rats both showed better hypoxia adaptation in AHRs than in BRs that failed to colonize the heart of the QTP, which may have facilitated the replacement of the BR population by the invading AHRs at the edge of the QTP. This study provides profound insights into the multiple origins of the northwardly invasive AHR and the great tolerance to hypoxia in this species.

Disentangling environmental drivers of circadian metabolism in desert-adapted mice

Metabolism is a complex phenotype shaped by natural environmental rhythms, as well as behavioral, morphological and physiological adaptations. Metabolism has been historically studied under constant environmental conditions, but new methods of continuous metabolic phenotyping now offer a window into organismal responses to dynamic environments, and enable identification of abiotic controls and the timing of physiological responses relative to environmental change. We used indirect calorimetry to characterize metabolic phenotypes of the desert-adapted cactus mouse (Peromyscus eremicus) in response to variable environmental conditions that mimic their native environment versus those recorded under constant warm and constant cool conditions, with a constant photoperiod and full access to resources. We found significant sexual dimorphism, with males being more prone to dehydration than females. Under circadian environmental variation, most metabolic shifts occurred prior to physical environmental change and the timing was disrupted under both constant treatments. The ratio of CO2 produced to O2 consumed (the respiratory quotient) reached greater than 1.0 only during the light phase under diurnally variable conditions, a pattern that strongly suggests that lipogenesis contributes to the production of energy and endogenous water. Our results are consistent with historical descriptions of circadian torpor in this species (torpid by day, active by night), but reject the hypothesis that torpor is initiated by food restriction or negative water balance.

Distinct Mechanisms Underlie Developmental Plasticity and Adult Acclimation of Thermogenic Capacity in High-Altitude Deer Mice

Developmental plasticity can elicit phenotypic adjustments that help organisms cope with environmental change, but the relationship between developmental plasticity and plasticity in adult life (e.g., acclimation) remains unresolved. We sought to examine developmental plasticity and adult acclimation in response to hypoxia of aerobic capacity (V̇O_2max) for thermogenesis in deer mice (Peromyscus maniculatus) native to high altitude. Deer mice were bred in captivity and exposed to normoxia or one of four hypoxia treatments (12 kPa O₂) across life stages: adult hypoxia (6–8 weeks), post-natal hypoxia (birth to adulthood), life-long hypoxia (before conception to adulthood), and parental hypoxia (mice conceived and raised in normoxia, but parents previously exposed to hypoxia). Hypoxia during perinatal development increased V̇O_2max by a much greater magnitude than adult hypoxia. The amplified effect of developmental hypoxia resulted from physiological plasticity that did not occur with adult hypoxia – namely, increases in lung ventilation and volume. Evolved characteristics of deer mice enabled developmental plasticity, because white-footed mice (P. leucopus; a congener restricted to low altitudes) could not raise pups in hypoxia. Parental hypoxia had no persistent effects on V̇O_2max. Therefore, developmental plasticity can have much stronger phenotypic effects and can manifest from distinct physiological mechanisms from adult acclimation.

A test of altitude-related variation in aerobic metabolism of Andean birds

Endotherms at high altitude face the combined challenges of cold and hypoxia. Cold increases thermoregulatory costs, and hypoxia may limit both thermogenesis and aerobic exercise capacity. Consequently, in comparisons between closely related highland and lowland taxa, we might expect to observe consistent differences in basal metabolic rate (BMR), maximal metabolic rate (MMR) and aerobic scope. Broad-scale comparative studies of birds reveal no association between BMR and native elevation, and altitude effects on MMR have not been investigated. We tested for altitude-related variation in aerobic metabolism in 10 Andean passerines representing five pairs of closely related species with contrasting elevational ranges. Mass-corrected BMR and MMR were significantly higher in most highland species relative to their lowland counterparts, but there was no uniform elevational trend across all pairs of species. Our results suggest that there is no simple explanation regarding the ecological and physiological causes of elevational variation in aerobic metabolism.

Limited Evidence for Parallel Evolution Among Desert-Adapted Peromyscus Deer Mice

Warming climate and increasing desertification urge the identification of genes involved in heat and dehydration tolerance to better inform and target biodiversity conservation efforts. Comparisons among extant desert-adapted species can highlight parallel or convergent patterns of genome evolution through the identification of shared signatures of selection. We generate a chromosome-level genome assembly for the canyon mouse (Peromyscus crinitus) and test for a signature of parallel evolution by comparing signatures of selective sweeps across population-level genomic resequencing data from another congeneric desert specialist (Peromyscus eremicus) and a widely distributed habitat generalist (Peromyscus maniculatus), that may be locally adapted to arid conditions. We identify few shared candidate loci involved in desert adaptation and do not find support for a shared pattern of parallel evolution. Instead, we hypothesize divergent molecular mechanisms of desert adaptation among deer mice, potentially tied to species-specific historical demography, which may limit or enhance adaptation. We identify a number of candidate loci experiencing selective sweeps in the P. crinitus genome that are implicated in osmoregulation (Trypsin, Prostasin) and metabolic tuning (Kallikrein, eIF2-alpha kinase GCN2, APPL1/2), which may be important for accommodating hot and dry environmental conditions.

Plasticity of non-shivering thermogenesis and brown adipose tissue in high-altitude deer mice

High altitude environments challenge small mammals with persistent low ambient temperatures that require high rates of aerobic heat production in face of low O2 availability. An important component of thermogenic capacity in rodents is non-shivering thermogenesis (NST) mediated by uncoupled mitochondrial respiration in brown adipose tissue (BAT). NST is plastic, and capacity for heat production increases with cold acclimation. However, in lowland native rodents, hypoxia inhibits NST in BAT. We hypothesize that highland deer mice (Peromyscus maniculatus) overcome the hypoxic inhibition of NST through changes in BAT mitochondrial function. We tested this hypothesis using lab born and raised highland and lowland deer mice, and a lowland congeneric (Peromyscus leucopus), acclimated to either warm normoxia (25°C, 760 mmHg) or cold hypoxia (5°C, 430 mmHg). We determined the effects of acclimation and ancestry on whole-animal rates of NST, the mass of interscapular BAT (iBAT), and uncoupling protein (UCP)-1 protein expression. To identify changes in mitochondrial function, we conducted high-resolution respirometry on isolated iBAT mitochondria using substrates and inhibitors targeted to UCP-1. We found that rates of NST increased with cold hypoxia acclimation but only in highland deer mice. There was no effect of cold hypoxia acclimation on iBAT mass in any group, but highland deer mice showed increases in UCP-1 expression and UCP-1-stimulated mitochondrial respiration in response to these stressors. Our results suggest that highland deer mice have evolved to increase the capacity for NST in response to chronic cold hypoxia, driven in part by changes in iBAT mitochondrial function.

Lipid oxidation during thermogenesis in high-altitude deer mice (Peromyscus maniculatus)

When at their maximum thermogenic capacity (cold-induced V̇o_2max), small endotherms reach levels of aerobic metabolism as high, or even higher, than running V̇o_2max. How these high rates of thermogenesis are supported by substrate oxidation is currently unclear. The appropriate utilization of metabolic fuels that could sustain thermogenesis over extended periods may be important for survival in cold environments, like high altitude. Previous studies show that high capacities for lipid use in high-altitude deer mice may have evolved in concert with greater thermogenic capacities. The purpose of this study was to determine how lipid utilization at both moderate and maximal thermogenic intensities may differ in high- and low-altitude deer mice, and strictly low-altitude white-footed mice. We also examined the phenotypic plasticity of lipid use after acclimation to cold hypoxia (CH), conditions simulating high altitude. We found that lipids were the primary fuel supporting both moderate and maximal rates of thermogenesis in both species of mice. Lipid oxidation increased threefold in mice from 30°C to 0°C, consistent with increases in oxidation of [¹³C]palmitic acid. CH acclimation led to an increase in [¹³C]palmitic acid oxidation at 30°C but did not affect total lipid oxidation. Lipid oxidation rates at cold-induced V̇o_2max were two- to fourfold those at 0°C and increased further after CH acclimation, especially in high-altitude deer mice. These are the highest mass-specific lipid oxidation rates observed in any land mammal. Uncovering the mechanisms that allow for these high rates of oxidation will aid our understanding of the regulation of lipid metabolism.

Mitochondrial Short-Term Plastic Responses and Long-Term Evolutionary Dynamics in Animal Species

How do species respond or adapt to environmental changes? The answer to this depends partly on mitochondrial epigenetics and genetics, new players in promoting adaptation to both short- and long-term environmental changes. In this review, we explore how mitochondrial epigenetics and genetics mechanisms, such as mtDNA methylation, mtDNA-derived noncoding RNAs, micropeptides, mtDNA mutations, and adaptations, can contribute to animal plasticity and adaptation. We also briefly discuss the challenges in assessing mtDNA adaptive evolution. In sum, this review covers new advances in the field of mitochondrial genomics, many of which are still controversial, and discusses processes still somewhat obscure, and some of which are still quite speculative and require further robust experimentation.

Adaptive Shifts in Gene Regulation Underlie a Developmental Delay in Thermogenesis in High-Altitude Deer Mice

Aerobic performance is tied to fitness as it influences an animal's ability to find food, escape predators, or survive extreme conditions. At high altitude, where low O2 availability and persistent cold prevail, maximum metabolic heat production (thermogenesis) is an aerobic performance trait that is closely linked to survival. Understanding how thermogenesis evolves to enhance survival at high altitude will yield insight into the links between physiology, performance, and fitness. Recent work in deer mice (Peromyscus maniculatus) has shown that adult mice native to high altitude have higher thermogenic capacities under hypoxia compared with lowland conspecifics, but that developing high-altitude pups delay the onset of thermogenesis. This finding suggests that natural selection on thermogenic capacity varies across life stages. To determine the mechanistic cause of this ontogenetic delay, we analyzed the transcriptomes of thermoeffector organs-brown adipose tissue and skeletal muscle-in developing deer mice native to low and high altitude. We demonstrate that the developmental delay in thermogenesis is associated with adaptive shifts in the expression of genes involved in nervous system development, fuel/O2 supply, and oxidative metabolism pathways. Our results demonstrate that selection has modified the developmental trajectory of the thermoregulatory system at high altitude and has done so by acting on the regulatory systems that control the maturation of thermoeffector tissues. We suggest that the cold and hypoxic conditions of high altitude force a resource allocation tradeoff, whereby limited energy is allocated to developmental processes such as growth, versus active thermogenesis, during early development.

High-Altitude Adaptation: Mechanistic Insights from Integrated Genomics and Physiology

Population genomic analyses of high-altitude humans and other vertebrates have identified numerous candidate genes for hypoxia adaptation, and the physiological pathways implicated by such analyses suggest testable hypotheses about underlying mechanisms. Studies of highland natives that integrate genomic data with experimental measures of physiological performance capacities and subordinate traits are revealing associations between genotypes (e.g., hypoxia-inducible factor gene variants) and hypoxia-responsive phenotypes. The subsequent search for causal mechanisms is complicated by the fact that observed genotypic associations with hypoxia-induced phenotypes may reflect second-order consequences of selection-mediated changes in other (unmeasured) traits that are coupled with the focal trait via feedback regulation. Manipulative experiments to decipher circuits of feedback control and patterns of phenotypic integration can help identify causal relationships that underlie observed genotype–phenotype associations. Such experiments are critical for correct inferences about phenotypic targets of selection and mechanisms of adaptation.

Phenotypic plasticity, genetic assimilation, and genetic compensation in hypoxia adaptation of high-altitude vertebrates

Important questions about mechanisms of physiological adaptation concern the role of phenotypic plasticity and the extent to which acclimatization responses align with genetic responses to selection. Such questions can be addressed in experimental studies of high-altitude vertebrates by investigating how mechanisms of acclimatization to hypoxia in lowland natives may influence genetic adaptation to hypoxia in highland natives. Evidence from high-altitude mammals suggest that evolved changes in some physiological traits involved canalization of the ancestral acclimatization response to hypoxia (genetic assimilation), a mechanism that results in an evolved reduction in plasticity. In addition to cases where adaptive plasticity may have facilitated genetic adaptation, evidence also suggests that some physiological changes in high-altitude natives are the result of selection to mitigate maladaptive plastic responses to hypoxia (genetic compensation). Examples of genetic compensation involve the attenuation of hypoxic pulmonary hypertension in Tibetan humans and other mammals native to high altitude. Here we discuss examples of adaptive physiological phenotypes in high-altitude natives that may have evolved by means of genetic assimilation or genetic compensation.

Physiological Genomics of Adaptation to High-Altitude Hypoxia

Population genomic studies of humans and other animals at high altitude have generated many hypotheses about the genes and pathways that may have contributed to hypoxia adaptation. Future advances require experimental tests of such hypotheses to identify causal mechanisms. Studies to date illustrate the challenge of moving from lists of candidate genes to the identification of phenotypic targets of selection, as it can be difficult to determine whether observed genotype-phenotype associations reflect causal effects or secondary consequences of changes in other traits that are linked via homeostatic regulation. Recent work on high-altitude models such as deer mice has revealed both plastic and evolved changes in respiratory, cardiovascular, and metabolic traits that contribute to aerobic performance capacity in hypoxia, and analyses of tissue-specific transcriptomes have identified changes in regulatory networks that mediate adaptive changes in physiological phenotype. Here we synthesize recent results and discuss lessons learned from studies of high-altitude adaptation that lie at the intersection of genomics and physiology.

Experimental study of hypoxia-induced changes in gene expression in an Asian pika, Ochotona dauurica

Acclimation to environmental changes driven by alterations in gene expression will serve as an important response for some species facing rapid Anthropogenic climate change. Pikas, genus Ochotona, are particularly vulnerable to climate change and current trends suggest that only the highest, coldest elevations within their ranges may remain suitable habitat for these species. In this study we aimed to assess the role of changes in gene expression in potentially facilitating elevational movements in pikas by measuring gene expression in the only known captive pika population, Ochotona dauurica, in response to hypoxic conditions. Using a controlled experiment, we exposed four male pikas to oxygen concentrations characteristic of sea-level, 2,000 m, and 4,000 m for 5 days each. Using blood samples collected after each treatment, we used RNAseq to determine if candidate pathways were undergoing significant changes in gene expression at different levels of oxygen (~100%, ~77%, and ~61% of sea-level oxygen concentrations). Gene set enrichment analyses showed that gene sets associated with the oxidative phosphorylation pathway and electron transport chain were significantly enriched for up-regulated genes in the 4,000 m samples compared to samples from the same individuals at lower-elevation conditions. Up-regulation of these pathways is consistent with known mechanisms of oxygen compensation. Our results suggest that these pikas have the acclimation capacity to tolerate oxygen concentrations characteristic of any elevation within their species range and that gene expression can be changed in a matter of days to accommodate drastically different oxygen concentrations. Thus, rapid and radical elevational movements that may be required of some pika species to avoid warmer temperatures in the Anthropocene will likely not be limited by hypoxic stress.

Naked mole-rats suppress energy metabolism and modulate membrane cholesterol in chronic hypoxia

Naked mole-rats (NMRs) are mammalian champions of hypoxia tolerance that enter metabolic suppression to survive in low oxygen environments. Common physiological mechanisms used by animals to suppress metabolic rate include downregulating energy metabolism (ATP supply) as well as ion pumps (primary cellular ATP consumers). A recent goldfish study demonstrated that remodeling of membrane lipids may mediate these responses, but it is unknown if NMR employs the same strategies; therefore, we aimed to test the hypotheses that these fossorial mammals 1) downregulate the activity of key enzymes of glycolysis, tricarboxylic acid (TCA) cycle, and β-oxidation, 2) inhibit sodium-potassium-ATPase, and 3) alter membrane lipids in response to chronic hypoxia. We found that NMRs exposed to 11% oxygen for 4 wk had a lower metabolic rate by 34%. This suppression occurs concurrently with tissue-specific 25-99% decreases in metabolic enzymes activities, a 77% decrease in brain sodium/potassium-ATPase activity, and widespread changes in membrane cholesterol abundance. By reducing glycolytic and β-oxidation fluxes, NMRs decrease the supply of acetyl-CoA to the TCA cycle. By contrast, there is a 94% upregulation of citrate synthase in the heart, possibly to support circulation and thus oxygen supply to other organs. Taken together, these responses may reflect a coordinated physiological response to hypoxia, but a clear functional link between changes in membrane composition and enzyme activities could not be established. Nevertheless, this is the first demonstration that hypometabolic NMRs alter the lipid composition of their membranes in response to chronic in vivo exposure to hypoxia.

Phenotypic plasticity as a long-term memory easing readaptations to ancestral environments

Phenotypic plasticity refers to environment-induced phenotypic changes without mutation and is present in all organisms. The role of phenotypic plasticity in organismal adaptations to novel environments has attracted much attention, but its role in readaptations to ancestral environments is understudied. To address this question, we use the reciprocal transplant approach to investigate the multitissue transcriptomes of chickens adapted to the Tibetan Plateau and adjacent lowland. While many genetic transcriptomic changes had occurred in the forward adaptation to the highland, plastic changes largely transform the transcriptomes to the preferred state when Tibetan chickens are brought back to the lowland. The same trend holds for egg hatchability, a key component of the chicken fitness. These findings, along with highly similar patterns in comparable experiments of guppies and Escherichia coli, demonstrate that organisms generally "remember" their ancestral environments via phenotypic plasticity and reveal a mechanism by which past experience affects future evolution.

Ontogenesis of evolved changes in respiratory physiology in deer mice native to high altitude

High-altitude environments are cold and hypoxic, and many high-altitude natives have evolved changes in respiratory physiology that improve O2 uptake in hypoxia as adults. Altricial mammals undergo a dramatic metabolic transition from ectothermy to endothermy in early post-natal life, which may influence the ontogenetic development of respiratory traits at high altitude. We examined the developmental changes in respiratory and haematological traits in deer mice (Peromyscus maniculatus) native to high altitude, comparing the respiratory responses to progressive hypoxia between highland and lowland deer mice. Among adults, highlanders exhibited higher total ventilation and a more effective breathing pattern (relatively deeper tidal volumes), for mice that were caught and tested at their native altitudes and those lab-raised in normoxia. Lab-raised progeny of each population were also tested at post-natal day (P)7, 14, 21 and 30. Highlanders developed an enhanced hypoxic ventilatory response by P21, concurrent with the full maturation of the carotid bodies, and their more effective breathing pattern arose by P14; these ages correspond to critical benchmarks in the full development of homeothermy in highlanders. However, highlanders exhibited developmental delays in ventilatory sensitivity to hypoxia, hyperplasia of type I cells in the carotid body and increases in blood haemoglobin content compared with lowland mice. Nevertheless, highlanders maintained consistently higher arterial O2 saturation in hypoxia across development, in association with increases in blood-O2 affinity that were apparent from birth. We conclude that evolved changes in respiratory physiology in high-altitude deer mice become expressed in association with the post-natal development of endothermy.

Faraway, so close. The comparative method and the potential of non-model animals in mitochondrial research

Inference from model organisms has been the engine for many discoveries in life science, but indiscriminate generalization leads to oversimplifications and misconceptions. Model organisms and inductive reasoning are irreplaceable: there is no other way to tackle the complexity of living systems. At the same time, it is not advisable to infer general patterns from a restricted number of species, which are very far from being representative of the diversity of life. Not all models are equal. Some organisms are suitable to find similarities across species, other highly specialized organisms can be used to focus on differences. In this opinion piece, we discuss the dominance of the mechanistic/reductionist approach in life sciences and make a case for an enhanced application of the comparative approach to study processes in all their various forms across different organisms. We also enlist some rising animal models in mitochondrial research, to exemplify how non-model organisms can be chosen in a comparative framework. These taxa often do not possess implemented tools and dedicated methods/resources. However, because of specific features, they have the potential to address still unanswered biological questions. Finally, we discuss future perspectives and caveats of the comparative method in the age of 'big data'. This article is part of the theme issue 'Linking the mitochondrial genotype to phenotype: a complex endeavour'.

Strong phenotypic divergence in spite of low genetic structure in the endemic Mangrove Warbler subspecies (Setophaga petechia xanthotera) of Costa Rica

Despite the enormous advances in genetics, links between phenotypes and genotypes have been made for only a few nonmodel organisms. However, such links can be essential to understand mechanisms of ecological speciation. The Costa Rican endemic Mangrove Warbler subspecies provides an excellent subject to study differentiation with gene flow, as it is distributed along a strong precipitation gradient on the Pacific coast with no strong geographic barriers to isolate populations. Mangrove Warbler populations could be subject to divergent selection driven by precipitation, which influences soil salinity levels, which in turn influences forest structure and food resources. We used single nucleotide polymorphisms (SNPs) and morphological traits to examine the balance between neutral genetic and phenotypic divergence to determine whether selection has acted on traits and genes with functions related to specific environmental variables. We present evidence showing: (a) associations between environmental variables and SNPs, identifying candidate genes related to bill morphology (BMP) and osmoregulation, (b) absence of population genetic structure in neutrally evolving markers, (c) divergence in bill size across the precipitation gradient, and (d) strong phenotypic differentiation (P ST) which largely exceeds neutral genetic differentiation (F ST) in bill size. Our results indicate an important role for salinity, forest structure, and resource availability in maintaining phenotypic divergence of Mangrove Warblers through natural selection. Our findings add to the growing body of literature identifying the processes involved in phenotypic differentiation along environmental gradients in the face of gene flow.

Physiological and genomic evidence that selection on the transcription factor Epas1 has altered cardiovascular function in high-altitude deer mice

Evolutionary adaptation to extreme environments often requires coordinated changes in multiple intersecting physiological pathways, but how such multi-trait adaptation occurs remains unresolved. Transcription factors, which regulate the expression of many genes and can simultaneously alter multiple phenotypes, may be common targets of selection if the benefits of induced changes outweigh the costs of negative pleiotropic effects. We combined complimentary population genetic analyses and physiological experiments in North American deer mice (Peromyscus maniculatus) to examine links between genetic variation in transcription factors that coordinate physiological responses to hypoxia (hypoxia-inducible factors, HIFs) and multiple physiological traits that potentially contribute to high-altitude adaptation. First, we sequenced the exomes of 100 mice sampled from different elevations and discovered that several SNPs in the gene Epas1, which encodes the oxygen sensitive subunit of HIF-2α, exhibited extreme allele frequency differences between highland and lowland populations. Broader geographic sampling confirmed that Epas1 genotype varied predictably with altitude throughout the western US. We then discovered that Epas1 genotype influences heart rate in hypoxia, and the transcriptomic responses to hypoxia (including HIF targets and genes involved in catecholamine signaling) in the heart and adrenal gland. Finally, we used a demographically-informed selection scan to show that Epas1 variants have experienced a history of spatially varying selection, suggesting that differences in cardiovascular function and gene regulation contribute to high-altitude adaptation. Our results suggest a mechanism by which Epas1 may aid long-term survival of high-altitude deer mice and provide general insights into the role that highly pleiotropic transcription factors may play in the process of environmental adaptation.

Developmental delay in shivering limits thermogenic capacity in juvenile high-altitude deer mice (Peromyscus maniculatus)

Many endotherms native to cold and hypoxic high-altitude (HA) environments have evolved a highly vascularized and aerobic skeletal muscle. This specialized muscle phenotype contributes via shivering to an enhanced capacity for aerobic thermogenesis (cold-induced V̇ _O2,max). However, it is unclear how selection at HA for shivering thermogenesis acts early in the development of small altricial mammals, which are born with immature skeletal muscles and without the capacity for homeothermic endothermy. We have previously shown that postnatal maturation of brown adipose tissue and non-shivering thermogenesis is delayed in HA native deer mouse pups (Peromyscus maniculatus). To assess whether HA adaptation has also altered the developmental program of skeletal muscle and shivering thermogenesis, we used laboratory-reared descendants of deer mice native to low altitude (LA, 430 m a.s.l.) and HA (4350 m a.s.l.) and a LA congeneric outgroup (P. leucopus). We found that LA juveniles were able to shiver robustly at 2 weeks after birth. However, HA juveniles were unlikely able to shiver at this point, resulting in a 30% lower capacity for thermoregulation compared with lowlanders. It was only at 27 days after birth that HA juveniles had established the aerobic muscle phenotype characteristic of HA adults and a superior cold-induced V̇ _O2,max compared with LA mice of the same age. The capacity for shivering may be delayed in HA mice to allow energy to be allocated to other important processes such as growth.

Thermal cues drive plasticity of desiccation resistance in montane salamanders with implications for climate change

Organisms rely upon external cues to avoid detrimental conditions during environmental change. Rapid water loss, or desiccation, is a universal threat for terrestrial plants and animals, especially under climate change, but the cues that facilitate plastic responses to avoid desiccation are unclear. We integrate acclimation experiments with gene expression analyses to identify the cues that regulate resistance to water loss at the physiological and regulatory level in a montane salamander (Plethodon metcalfi). Here we show that temperature is an important cue for developing a desiccation-resistant phenotype and might act as a reliable cue for organisms across the globe. Gene expression analyses consistently identify regulation of stem cell differentiation and embryonic development of vasculature. The temperature-sensitive blood vessel development suggests that salamanders regulate water loss through the regression and regeneration of capillary beds in the skin, indicating that tissue regeneration may be used for physiological purposes beyond replacing lost limbs.

Life Ascending: Mechanism and Process in Physiological Adaptation to High-Altitude Hypoxia

To cope with the reduced availability of O₂ at high altitude, air-breathing vertebrates have evolved myriad adjustments in the cardiorespiratory system to match tissue O₂ delivery with metabolic O₂ demand. We explain how changes at interacting steps of the O₂ transport pathway contribute to plastic and evolved changes in whole-animal aerobic performance under hypoxia. In vertebrates native to high altitude, enhancements of aerobic performance under hypoxia are attributable to a combination of environmentally induced and evolved changes in multiple steps of the pathway. Additionally, evidence suggests that many high-altitude natives have evolved mechanisms for attenuating maladaptive acclimatization responses to hypoxia, resulting in counter-gradient patterns of altitudinal variation for key physiological phenotypes. For traits that exhibit counteracting environmental and genetic effects, evolved changes in phenotype may be cryptic under field conditions and can only be revealed by rearing representatives of high-and low-altitude populations under standardized environmental conditions to control for plasticity.

Evolution of physiological performance capacities and environmental adaptation: insights from high-elevation deer mice (Peromyscus maniculatus)

Analysis of variation in whole-animal performance can shed light on causal connections between specific traits, integrated physiological capacities, and Darwinian fitness. Here, we review and synthesize information on naturally occurring variation in physiological performance capacities and how it relates to environmental adaptation in deer mice (Peromyscus maniculatus). We discuss how evolved changes in aerobic exercise capacity and thermogenic capacity have contributed to adaptation to high elevations. Comparative work on deer mice at high and low elevations has revealed evolved differences in aerobic performance capacities in hypoxia. Highland deer mice have consistently higher aerobic performance capacities under hypoxia relative to lowland natives, consistent with the idea that it is beneficial to have a higher maximal metabolic rate (as measured by the maximal rate of O₂ consumption, VO_2max) in an environment characterized by lower air temperatures and lower O₂ availability. Observed differences in aerobic performance capacities between highland and lowland deer mice stem from changes in numerous subordinate traits that alter the flux capacity of the O₂-transport system, the oxidative capacity of tissue mitochondria, and the relationship between O₂ consumption and ATP synthesis. Many such changes in physiological phenotype are associated with hypoxia-induced changes in gene expression. Research on natural variation in whole-animal performance forms a nexus between physiological ecology and evolutionary biology that requires insight into the natural history of the study species.

Genome-Wide Microarray Analysis Suggests Transcriptomic Response May Not Play a Major Role in High- to Low-Altitude Acclimation in Harvest Mouse (Micromys minutus)

Simple Summary

Micromys minutus is a small rodent species that has a wide range of vertical distribution in Taiwan. By comparing the gene expression profile of the skeletal muscle tissues taken from individuals native to the high-altitude environment and those transferred to the low-altitude captive site, the Tnfrsf12a gene was demonstrated to have a differential expression pattern. Although this finding may be correlated with the altitude acclimation, the observation of only one gene transcript with significant alteration leads us to suggest that genetic response may not play a major role in altitude acclimation in M. minutus. Future comparative functional genomics studies involving other organ systems (in addition to skeletal muscles) and alarger sample size are warranted for better insight into the altitude acclimation of this small rodent species.

Abstract

The harvest mouse (Micromys minutus) is a small rodent species with a wide range of vertical distribution in Taiwan, extending from the sea level to 3100 m altitude. This species has recently suffered from habitat loss in high-altitude areas due to orchard cultivation, which may have resulted in mouse migration from high to low altitude. To investigate whether there is any physiological mechanism involved in altitude acclimation, rat cDNA microarray was used to compare transcriptomic patterns of the skeletal muscle tissues taken from individuals native to the high-altitude environment and those transferred to the low-altitude captive site. Of the 23,188 genes being analyzed, 47 (33 up-regulated and 14 down-regulated) were found to have differential expression (fold change > 4 or < −4, ANOVA p < 0.05). However, after multiple testing correction with a false discovery rate (FDR), only the result for Tnfrsf12a was found to be statistically significant (fold change = 13, FDR p < 0.05). The result was confirmed by quantitative polymerase chain reaction (q-PCR). The expression of Tnfrsf12a possibly relates to the skeletal muscle biology and thus can be correlated with altitude acclimation. However, finding only one gene transcript with significant alteration suggests that transcriptomic response may not play a major role in high- to low-altitude acclimation in harvest mouse.

The Mitochondrial Basis for Adaptive Variation in Aerobic Performance in High-Altitude Deer Mice

Mitochondria play a central role in aerobic performance. Studies aimed at elucidating how evolved variation in mitochondrial physiology contributes to adaptive variation in aerobic performance can therefore provide a unique and powerful lens to understanding the evolution of complex physiological traits. Here, we review our ongoing work on the importance of changes in mitochondrial quantity and quality to adaptive variation in aerobic performance in high-altitude deer mice. Whole-organism aerobic capacity in hypoxia (VO2max) increases in response to hypoxia acclimation in this species, but high-altitude populations have evolved consistently greater VO2max than populations from low altitude. The evolved increase in VO2max in highlanders is associated with an evolved increase in the respiratory capacity of the gastrocnemius muscle. This appears to result from highlanders having more mitochondria in this tissue, attributed to a higher proportional abundance of oxidative fiber-types and a greater mitochondrial volume density within oxidative fibers. The latter is primarily caused by an over-abundance of subsarcolemmal mitochondria in high-altitude mice, which is likely advantageous for mitochondrial O2 supply because more mitochondria are situated adjacent to the cell membrane and close to capillaries. Evolved changes in gastrocnemius phenotype appear to be underpinned by population differences in the expression of genes involved in energy metabolism, muscle development, and vascular development. Hypoxia acclimation has relatively little effect on respiratory capacity of the gastrocnemius, but it increases respiratory capacity of the diaphragm. However, the mechanisms responsible for this increase differ between populations: lowlanders appear to adjust mitochondrial quantity and quality (i.e., increases in citrate synthase [CS] activity, and mitochondrial respiration relative to CS activity) and they exhibit higher rates of mitochondrial release of reactive oxygen species, whereas highlanders only increase mitochondrial quantity in response to hypoxia acclimation. In contrast to the variation in skeletal muscles, the respiratory capacity of cardiac muscle does not appear to be affected by hypoxia acclimation and varies little between populations. Therefore, evolved changes in mitochondrial quantity and quality make important tissue-specific contributions to adaptive variation in aerobic performance in high-altitude deer mice.