Genomic structural variation contributes to evolved changes in gene expression in high-altitude Tibetan sheep.

Tibetan sheep were introduced to the Qinghai Tibet plateau roughly 3,000 B.P., making this species a good model for investigating genetic mechanisms of high-altitude adaptation over a relatively short timescale. Here, we characterize genomic structural variants (SVs) that distinguish Tibetan sheep from closely related, low-altitude Hu sheep, and we examine associated changes in tissue-specific gene expression. We document differentiation between the two sheep breeds in frequencies of SVs associated with genes involved in cardiac function and circulation. In Tibetan sheep, we identified high-frequency SVs in a total of 462 genes, including EPAS1, PAPSS2, and PTPRD. Single-cell RNA-Seq data and luciferase reporter assays revealed that the SVs had cis-acting effects on the expression levels of these three genes in specific tissues and cell types. In Tibetan sheep, we identified a high-frequency chromosomal inversion that exhibited modified chromatin architectures relative to the noninverted allele that predominates in Hu sheep. The inversion harbors several genes with altered expression patterns related to heart protection, brown adipocyte proliferation, angiogenesis, and DNA repair. These findings indicate that SVs represent an important source of genetic variation in gene expression and may have contributed to high-altitude adaptation in Tibetan sheep.

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.

Genomic structural variation is associated with hypoxia adaptation in high-altitude zokors.

Zokors, an Asiatic group of subterranean rodents, originated in lowlands and colonized high-elevational zones following the uplift of the Qinghai-Tibet plateau about 3.6 million years ago. Zokors live at high elevation in subterranean burrows and experience hypobaric hypoxia, including both hypoxia (low oxygen concentration) and hypercapnia (elevated partial pressure of CO2). Here we report a genomic analysis of six zokor species (genus Eospalax) with different elevational ranges to identify structural variants (deletions and inversions) that may have contributed to high-elevation adaptation. Based on an assembly of a chromosome-level genome of the high-elevation species, Eospalax baileyi, we identified 18 large inversions that distinguished this species from congeners native to lower elevations. Small-scale structural variants in the introns of EGLN1, HIF1A, HSF1 and SFTPD of E. baileyi were associated with the upregulated expression of those genes. A rearrangement on chromosome 1 was associated with altered chromatin accessibility, leading to modified gene expression profiles of key genes involved in the physiological response to hypoxia. Multigene families that underwent copy-number expansions in E. baileyi were enriched for autophagy, HIF1 signalling and immune response. E. baileyi show a significantly larger lung mass than those of other Eospalax species. These findings highlight the key role of structural variants underlying hypoxia adaptation of high-elevation species in Eospalax.

To what extent do physiological tolerances determine elevational range limits of mammals?

A key question in biology concerns the extent to which distributional range limits of species are determined by intrinsic limits of physiological tolerance. Here, we use common-garden data for wild rodents to assess whether species with higher elevational range limits typically have higher thermogenic capacities in comparison to closely related lowland species. Among South American leaf-eared mice (genus Phyllotis), mean thermogenic performance is higher in species with higher elevational range limits, but there is little among-species variation in the magnitude of plasticity in this trait. In the North American rodent genus Peromyscus, highland deer mice (Peromyscus maniculatus) have greater thermogenic maximal oxygen uptake ( V ̇ O 2 max ${\dot V_{{{\mathrm{O}}_{\mathrm{2}}}{\mathrm{max}}}}$ ) than lowland white-footed mice (Peromyscus leucopus) at a level of hypoxia that matches the upper elevational range limit of the former species. In highland deer mice, the enhanced thermogenic V ̇ O 2 max ${\dot V_{{{\mathrm{O}}_{\mathrm{2}}}{\mathrm{max}}}}$ in hypoxia is attributable to a combination of evolved and plastic changes in physiological pathways that govern the transport and utilization of O2 and metabolic substrates. Experiments with Peromyscus mice also demonstrate that exposure to hypoxia during different stages of development elicits plastic changes in cardiorespiratory traits that improve thermogenic V ̇ O 2 max ${\dot V_{{{\mathrm{O}}_{\mathrm{2}}}{\mathrm{max}}}}$ via distinct physiological mechanisms. Evolved differences in thermogenic capacity provide clues about why some species are able to persist in higher-elevation habitats that lie slightly beyond the tolerable limits of other species. Such differences in environmental tolerance also suggest why some species might be more vulnerable to climate change than others.

Genomic insights into the mystery of mouse mummies on the summits of Atacama volcanoes.

Our understanding of the limits of animal life is continually revised by scientific exploration of extreme environments. Here we report the discovery of mummified cadavers of leaf-eared mice, Phyllotis vaccarum, from the summits of three different Andean volcanoes at elevations 6,029-6,233 m above sea level in the Puna de Atacama in Chile and Argentina. Such extreme elevations were previously assumed to be completely uninhabitable by mammals. In combination with a live-captured specimen of the same species from the nearby summit of Volcán Llullaillaco (6,739 m)1, the summit mummies represent the highest altitude physical records of mammals in the world. We also report a chromosome-level genome assembly for P. vaccarum that, in combination with a whole-genome re-sequencing analysis and radiocarbon dating analysis, provides insights into the provenance and antiquity of the summit mice. Radiocarbon data indicate that the most ancient of the mummies are, at most, a few centuries old. Genomic polymorphism data revealed a high degree of continuity between the summit mice and conspecifics from lower elevations in the surrounding Altiplano. Genomic data also revealed equal numbers of males and females among the summit mice and evidence of close kinship between some individuals from the same summits. These findings bolster evidence for resident populations of Phyllotis at elevations >6,000 m and challenge assumptions about the environmental limits of vertebrate life and the physiological tolerances of small mammals.

Elevational range extension of the Puna Mouse, Punomys (Cricetidae), with the first record of the genus from Chile.

We report an elevational record for the Andean sigmodontine Puna Mouse Punomys, which is also the first record of the genus in Chile. The record is based on a mummified specimen that we discovered at an elevation of 5,461 m (17,917 feet) in the caldera of Volcán Acamarachi, Región de Antofagasta, Chile. Results of a morphological assessment suggest that the specimen can be provisionally referred to the species P. lemminus. This new record also extends the known geographic distribution of the genus by 700 km to the south and brings the known Chilean mammal richness to a total of 170 living species and 88 genera. This finding highlights the need for increased survey efforts in more remote, high-elevation regions and demonstrates that there is still much to be learned about the mammal fauna of the Andean Altiplano.

Se reporta un registro altitudinal para el roedor sigmodontino Punomys, el cual corresponde a su vez al primer hallazgo del género para Chile. Este se basa en un espécimen momificado encontrado a una elevación de 5,461 m en la caldera del Volcán Acamarachi, Región de Antofagasta, Chile. Los caracteres morfológicos sugieren que el espécimen puede ser referido provisionalmente a la especie P. lemminus. Este nuevo registro amplía la distribución geográfica conocida del género en 700 km al sur, y eleva la riqueza de mamíferos vivientes chilenos a un nuevo total de 170 especies y 88 géneros. Este hallazgo resalta la necesidad de aumentar los esfuerzos de prospección en las regiones más remotas y de mayor altitud y demuestra que aún queda mucho por aprender sobre el ensamble de los mamíferos del Altiplano andino.

Extreme high-elevation mammal surveys reveal unexpectedly high upper range limits of Andean mice.

In 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 live-trapped small mammals from ecologically diverse sites spanning >6700 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.

Blood Variation Implicates Respiratory Limits on Elevational Ranges of Andean Birds.

AbstractThe extent to which species ranges reflect intrinsic physiological tolerances is a major question in evolutionary ecology. To date, consensus has been hindered by the limited tractability of experimental approaches across most of the tree of life. Here, we apply a macrophysiological approach to understand how hematological traits related to oxygen transport shape elevational ranges in a tropical biodiversity hot spot. Along Andean elevational gradients, we measured traits that affect blood oxygen-carrying capacity-total and cellular hemoglobin concentration and hematocrit, the volume percentage of red blood cells-for 2,355 individuals of 136 bird species. We used these data to evaluate the influence of hematological traits on elevational ranges. First, we asked whether the sensitivity of hematological traits to changes in elevation is predictive of elevational range breadth. Second, we asked whether variance in hematological traits changed as a function of distance to the nearest elevational range limit. We found that birds showing greater hematological sensitivity had broader elevational ranges, consistent with the idea that a greater acclimatization capacity facilitates elevational range expansion. We further found reduced variation in hematological traits in birds sampled near their elevational range limits and at high absolute elevations, patterns consistent with intensified natural selection, reduced effective population size, or compensatory changes in other cardiorespiratory traits. Our findings suggest that constraints on hematological sensitivity and local genetic adaptation to oxygen availability promote the evolution of the narrow elevational ranges that underpin tropical montane biodiversity.

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.

Evolution and molecular basis of a novel allosteric property of crocodilian hemoglobin.

The extraordinary breath-hold diving capacity of crocodilians has been ascribed to a unique mode of allosterically regulating hemoglobin (Hb)-oxygenation in circulating red blood cells. We investigated the origin and mechanistic basis of this novel biochemical phenomenon by performing directed mutagenesis experiments on resurrected ancestral Hbs. Comparisons of Hb function between the common ancestor of archosaurs (the group that includes crocodilians and birds) and the last common ancestor of modern crocodilians revealed that regulation of Hb-O2 affinity via allosteric binding of bicarbonate ions represents a croc-specific innovation that evolved in combination with the loss of allosteric regulation by ATP binding. Mutagenesis experiments revealed that evolution of the novel allosteric function in crocodilians and the concomitant loss of ancestral function were not mechanistically coupled and were caused by different sets of substitutions. The gain of bicarbonate sensitivity in crocodilian Hb involved the direct effect of few amino acid substitutions at key sites in combination with indirect effects of numerous other substitutions at structurally disparate sites. Such indirect interaction effects suggest that evolution of the novel protein function was conditional on neutral mutations that produced no adaptive benefit when they first arose but that contributed to a permissive background for subsequent function-altering mutations at other sites. Due to the context dependence of causative substitutions, the unique allosteric properties of crocodilian Hb cannot be easily transplanted into divergent homologs of other species.

Epistatic Effects Between Amino Acid Insertions and Substitutions Mediate Toxin resistance of Vertebrate Na+,K+-ATPases.

The recurrent evolution of resistance to cardiotonic steroids (CTS) across diverse animals most frequently involves convergent amino acid substitutions in the H1-H2 extracellular loop of Na+,K+-ATPase (NKA). Previous work revealed that hystricognath rodents (e.g., chinchilla) and pterocliform birds (sandgrouse) have convergently evolved amino acid insertions in the H1-H2 loop, but their functional significance was not known. Using protein engineering, we show that these insertions have distinct effects on CTS resistance in homologs of each of the two species that strongly depend on intramolecular interactions with other residues. Removing the insertion in the chinchilla NKA unexpectedly increases CTS resistance and decreases NKA activity. In the sandgrouse NKA, the amino acid insertion and substitution Q111R both contribute to an augmented CTS resistance without compromising ATPase activity levels. Molecular docking simulations provide additional insight into the biophysical mechanisms responsible for the context-specific mutational effects on CTS insensitivity of the enzyme. Our results highlight the diversity of genetic substrates that underlie CTS insensitivity in vertebrate NKA and reveal how amino acid insertions can alter the phenotypic effects of point mutations at key sites in the same protein domain.

Evidence of a population of leaf-eared mice Phyllotis vaccarum above 6,000 m in the Andes and a survey of high-elevation mammals.

Biologists have long pondered the extreme limits of life on Earth, including the maximum elevation at which species can live and reproduce. Here we review evidence of a self-sustaining population of mice at an elevation that exceeds that of all previously reported for mammals. Five expeditions over 10 years to Volcán Llullaillaco on the Argentina/Chile border observed and collected mice at elevations ranging from 5,070 m at the mountain's base to the summit at 6,739 m (22,110 feet). Previously unreported evidence includes observations and photographs of live animals and mummified remains, environmental DNA, and a soil microbial community reflecting animal activity that are evaluated in combination with previously reported video recordings and capture of live mice. All of the evidence identifies the mouse as the leaf-eared mouse Phyllotis vaccarum, and it robustly places the population within a haplotype group containing individuals from the Chilean Atacama Desert and nearby regions of Argentina. A critical review of the literature affirms that this population is not only an elevational record for mammals but for all terrestrial vertebrates to date, and we further find that many extreme elevations previously reported for mammals are based on scant or dubious evidence.

Durante mucho tiempo los biólogos han reflexionado sobre los límites extremos de altura a la que las especies pueden vivir y reproducirse. Aquí presentamos nueva evidencia sobre la existencia de una población de ratones establecida a una elevación que supera todos los reports previos para mamíferos. Durante 10 años fueron realizadas 5 expediciones al Volcán Llullaillaco, ubicado en la frontera entre Argentina y Chile; observando y colectando ratones en elevaciones que van desde los 5,070 m hasta la cima de 6,739 m (22,110 feet). La nueva evidencia incluye fotografías de restos momificados, ADN ambiental y la actividad microbiana del suelo que confirman la presencia del animal, la cual fue analizada junto a videos reportados anteriormente y la captura de ejemplares vivos. Toda esta información indica que dicha población corresponde al ratón orejudo amarillento Phyllotis vaccarum y lo posicionan dentro de un grupo de haplotipos compuesto por individuos del Desierto de Atacama y regiones cercanas en Argentina. La revisión crítica de la literatura demostró que esta población no solo es un récord de elevación para los mamíferos, sino para todos los vertebrados terrestres; igualmente, que los reportes de elevaciones extremas reportados para mamíferos se derivan de evidencias escasas y dudosas.

Constraints on the evolution of toxin-resistant Na,K-ATPases have limited dependence on sequence divergence.

A growing body of theoretical and experimental evidence suggests that intramolecular epistasis is a major determinant of rates and patterns of protein evolution and imposes a substantial constraint on the evolution of novel protein functions. Here, we examine the role of intramolecular epistasis in the recurrent evolution of resistance to cardiotonic steroids (CTS) across tetrapods, which occurs via specific amino acid substitutions to the α-subunit family of Na,K-ATPases (ATP1A). After identifying a series of recurrent substitutions at two key sites of ATP1A that are predicted to confer CTS resistance in diverse tetrapods, we then performed protein engineering experiments to test the functional consequences of introducing these substitutions onto divergent species backgrounds. In line with previous results, we find that substitutions at these sites can have substantial background-dependent effects on CTS resistance. Globally, however, these substitutions also have pleiotropic effects that are consistent with additive rather than background-dependent effects. Moreover, the magnitude of a substitution's effect on activity does not depend on the overall extent of ATP1A sequence divergence between species. Our results suggest that epistatic constraints on the evolution of CTS-resistant forms of Na,K-ATPase likely depend on a small number of sites, with little dependence on overall levels of protein divergence. We propose that dependence on a limited number sites may account for the observation of convergent CTS resistance substitutions observed among taxa with highly divergent Na,K-ATPases (See S1 Text for Spanish translation).

Time Domains of Hypoxia Responses and -Omics Insights.

The ability to respond rapidly to changes in oxygen tension is critical for many forms of life. Challenges to oxygen homeostasis, specifically in the contexts of evolutionary biology and biomedicine, provide important insights into mechanisms of hypoxia adaptation and tolerance. Here we synthesize findings across varying time domains of hypoxia in terms of oxygen delivery, ranging from early animal to modern human evolution and examine the potential impacts of environmental and clinical challenges through emerging multi-omics approaches. We discuss how diverse animal species have adapted to hypoxic environments, how humans vary in their responses to hypoxia (i.e., in the context of high-altitude exposure, cardiopulmonary disease, and sleep apnea), and how findings from each of these fields inform the other and lead to promising new directions in basic and clinical hypoxia research.

Genomic signatures of high-altitude adaptation and chromosomal polymorphism in geladas.

Primates have adapted to numerous environments and lifestyles but very few species are native to high elevations. Here we investigated high-altitude adaptations in the gelada (Theropithecus gelada), a monkey endemic to the Ethiopian Plateau. We examined genome-wide variation in conjunction with measurements of haematological and morphological traits. Our new gelada reference genome is highly intact and assembled at chromosome-length levels. Unexpectedly, we identified a chromosomal polymorphism in geladas that could potentially contribute to reproductive barriers between populations. Compared with baboons at low altitude, we found that high-altitude geladas exhibit significantly expanded chest circumferences, potentially allowing for greater lung surface area for increased oxygen diffusion. We identified gelada-specific amino acid substitutions in the alpha-chain subunit of adult haemoglobin but found that gelada haemoglobin does not exhibit markedly altered oxygenation properties compared with lowland primates. We also found that geladas at high altitude do not exhibit elevated blood haemoglobin concentrations, in contrast to the normal acclimatization response to hypoxia in lowland primates. The absence of altitude-related polycythaemia suggests that geladas are able to sustain adequate tissue-oxygen delivery despite environmental hypoxia. Finally, we identified numerous genes and genomic regions exhibiting accelerated rates of evolution, as well as gene families exhibiting expansions in the gelada lineage, potentially reflecting altitude-related selection. Our findings lend insight into putative mechanisms of high-altitude adaptation while suggesting promising avenues for functional hypoxia research.

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 O2 transport via changes in blood hemoglobin (Hb) content and/or Hb-O2 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 O2-transport pathway.

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.

Changes in hemoglobin function and isoform expression during embryonic development in the American alligator, Alligator mississippiensis.

In the developing embryos of egg-laying vertebrates, O2 flux takes place across a fixed surface area of the eggshell and the chorioallantoic membrane. In the case of crocodilians, the developing embryo may experience a decrease in O2 flux when the nest becomes hypoxic, which may cause compensatory adjustments in blood O2 transport. However, whether the switch from embryonic to adult hemoglobin isoforms (isoHbs) plays some role in these adjustments is unknown. Here, we provide a detailed characterization of the developmental switch of isoHb synthesis in the American alligator, Alligator mississippiensis. We examined the in vitro functional properties and subunit composition of purified alligator isoHbs expressed during embryonic developmental stages in normoxia and hypoxia (10% O2). We found distinct patterns of isoHb expression in alligator embryos at different stages of development, but these patterns were not affected by hypoxia. Specifically, alligator embryos expressed two main isoHbs: HbI, prevalent at early developmental stages, with a high O2 affinity and high ATP sensitivity, and HbII, prevalent at later stages and identical to the adult protein, with a low O2 affinity and high CO2 sensitivity. These results indicate that whole blood O2 affinity is mainly regulated by ATP in the early embryo and by CO2 and bicarbonate from the late embryo until adult life, but the developmental regulation of isoHb expression is not affected by hypoxia exposure.

Whole-Genome Duplications and the Diversification of the Globin-X Genes of Vertebrates.

Globin-X (GbX) is an enigmatic member of the vertebrate globin gene family with a wide phyletic distribution that spans protostomes and deuterostomes. Unlike canonical globins such as hemoglobins and myoglobins, functional data suggest that GbX does not have a primary respiratory function. Instead, evidence suggests that the monomeric, membrane-bound GbX may play a role in cellular signaling or protection against the oxidation of membrane lipids. Recently released genomes from key vertebrates provide an excellent opportunity to address questions about the early stages of the evolution of GbX in vertebrates. We integrate bioinformatics, synteny, and phylogenetic analyses to characterize the diversity of GbX genes in nonteleost ray-finned fishes, resolve relationships between the GbX genes of cartilaginous fish and bony vertebrates, and demonstrate that the GbX genes of cyclostomes and gnathostomes derive from independent duplications. Our study highlights the role that whole-genome duplications (WGDs) have played in expanding the repertoire of genes in vertebrate genomes. Our results indicate that GbX paralogs have a remarkably high rate of retention following WGDs relative to other globin genes and provide an evolutionary framework for interpreting results of experiments that examine functional properties of GbX and patterns of tissue-specific expression. By identifying GbX paralogs that are products of different WGDs, our results can guide the design of experimental work to explore whether gene duplicates that originate via WGDs have evolved novel functional properties or expression profiles relative to singleton or tandemly duplicated copies of GbX.