Metabolic characteristics and response to high altitude in Phrynocephalus erythrurus (Lacertilia: Agamidae), a lizard dwell at altitudes higher than any other living lizards in the world

Unraveling molecular mechanistic disparities in pathogenic visceral Leishmania resistance between reptiles and mammals through comparative transcriptomic analyses

Leishmaniasis is one of the most important neglected tropical parasitic diseases, manifesting various clinical forms depending on the parasite species and the genetic background of the host. In order to elucidate the underlying mechanisms of reptilian defense against pathogenic Leishmania species and to delineate the global gene expression profile alterations during host-pathogen interaction, we established experimental animal and cell models using both heterothermic lizards (Phrynocephalus przewalskii) and homothermic mammals (BALB/c mice) infected with pathogenic Leishmania infantum (high virulence HCZ strain) and Leishmania donovani (low virulence 801 strain). Overall, the lizards didn't show any obvious clinical symptoms or immune responses in vivo. Using RNA-seq methodology, differentially expressed genes identified in the HCZ and 801-comparison groups of P. przewalskii were primarily associated with arginine biosynthesis, the MAPK signaling pathway and the PI3K-Akt signaling pathway. In contrast, higher parasite loads, exacerbated hepatic inflammatory lesions and enhanced immune responses were observed in BALB/c mice, with DEGs predominantly associated with immunological diseases, innate and adaptive immune responses. By integrating transcriptional data from reptile and mammalian hosts, we elucidated the pivotal role of amino acid metabolism and lipid metabolism in parasite control. In contrast to findings from animal experiments, Leishmania parasites effectively infected peritoneal macrophages of lizards in vitro, demonstrating a high infection rate. Furthermore, we used RT-qPCR to detect changes in cytokine expression in macrophages and found that Th1-type cytokines were significantly upregulated in lizards, facilitating the clearance of the HCZ strain 24 hours post-infection. Conversely, cytokine expression was generally suppressed in BALB/c mice, allowing immune evasion by the parasites.

Selection and Gene Duplication Associated With High-Elevation Diversification in Pristimantis, the Largest Terrestrial Vertebrate Genus

The genus Pristimantis diversified in the tropical Andes mountains and is the most speciose genus of terrestrial vertebrates. Pristimantis are notable among frogs in that they thrive at high elevations (>2,000 m) and are direct developers without a tadpole stage. Despite their ecological significance, little is known about the genetic and physiological traits enabling their success. We conducted transcriptomic analysis on seven Pristimantis species sampled across elevations in the Ecuadorean Andes to explore three hypotheses for their success: (i) unique genes are under selection relative to all other frogs, (ii) common selection occurs across all direct developers, or (iii) common selection occurs across all high-elevation frog clades. Comparative analysis with 34 frog species revealed unique positive selection in Pristimantis genes related to aerobic respiration, hemostasis, signaling, cellular transportation of proteins and ions, and immunity. Additionally, we detected positive selection across all direct developers for genes associated with oxygenase activity and metal ion binding. While many genes under selection in Pristimantis were not positively selected in other high-elevation frog species, we identified some shared genes and pathways linked to lipid metabolism, innate immunity, and cellular redox processes. We observed more positive selection in duplicated- versus single-copy genes, while relaxed purifying selection was prevalent in single-copy genes. Notably, copy number of an innate immunity complement gene was positively correlated with Pristimantis species elevation. Our findings contribute novel insights into the genetic basis of adaptation in Pristimantis and provide a foundation for future studies on the evolutionary mechanisms leading to direct development and coping with high elevations.

Effects of warming at embryonic and larval stages on tadpole fitness in high-altitude Rana kukunoris

Global warming may affect the early developmental stages of high-altitude amphibians, thereby influencing their later fitness. Yet, this has been largely unexplored. To investigate whether and how the temperatures experienced by embryonic and larval stages affect their fitness at later developmental stages, we designed two experiments in which the embryos and larvae were treated with three temperatures (24, 18 and 12 °C), respectively. Then, the life history traits of the tadpoles during the metamorphotic climax in all treatments were evaluated, including growth rate, survival rate, morphology, thermal physiology, swimming performance, standard metabolic rate (SMR), oxidative and antioxidative system, and metabolic enzyme activities. The results revealed that elevated temperature accelerated metamorphosis but decreased body size at metamorphosis. Additionally, warming during the embryonic and larval stages decreased the thermal tolerance range and induced increased oxidative stress. Furthermore, high embryonic temperature significantly decreased the hatching success, but had no significant effect on swimming performance and SMR. Warming during larval periods was harmful to the survival and swimming performance of tadpoles. The effect size analysis revealed that the negative impacts of embryonic temperature on certain physiological traits, such as growth and development, survival and swimming performance, were more pronounced than those of larval temperature. Our results highlight the necessity for particular attention to be paid to the early stages of amphibians, notably the embryonic stages when evaluating the impact of global warming on their survival.

Effects of microhabitat features on the intraspecific variability of the distribution and functional traits in a highest elevational distributed lizard

Exploring the microhabitat determinants of organisms distribution and functional traits differences can help us better understand the importance of intraspecific variations in ecological niches. Investigations on animals functional niche primarily focused on differences among species and tended to neglect the potential variability within species, despite the fact that the ecological and evolutionary importance of intraspecific variations was widely recognized. In this study, we examined the influence of microhabitat features on the intraspecific variability of the distribution and functional traits of a highest elevational distributed lizard species Phrynocephalus erythrurus. To do so, field work was conducted between July and August, 2020 and August and September, 2021 in Namtso watershed in central Xizang, China. Specifically, 11 transects were sampled for P. erythrurus individuals, which were measured for a set of 10 morphological traits. Moreover, 11 microhabitat variables that potentially affect the distribution of lizards were also measured for each transect. Our results indicated that juveniles, males, and females exhibited different functional traits, allowing them to occupy distinct functional space. The distribution of juveniles, males, and females was determined by different microhabitat variables such as illuminance and air temperature. More importantly, these variables also determined the intraspecific functional traits variability in this lizard species. All of these results supported previous claims that intraspecific traits variation should be incorporated into functional ecological studies, and diverse microhabitat features should be conserved to maintain high intraspecific diversity. Future studies can focus on the food analysis to explore the linkage between functional traits and resources utilization within animal populations.

Comparative metabolomics analysis reveals high-altitude adaptations in a toad-headed viviparous lizard, Phrynocephalus vlangalii

Extreme environmental conditions at high altitude, such as hypobaric hypoxia, low temperature, and strong UV radiation, pose a great challenge to the survival of animals. Although the mechanisms of adaptation to high-altitude environments have attracted much attention for native plateau species, the underlying metabolic regulation remains unclear. Here, we used a multi-platform metabolomic analysis to compare metabolic profiles of liver between high- and low-altitude populations of toad-headed lizards, Phrynocephalus vlangalii, from the Qinghai-Tibet Plateau. A total of 191 differential metabolites were identified, consisting of 108 up-regulated and 83 down-regulated metabolites in high-altitude lizards as compared with values for low-altitude lizards. Pathway analysis revealed that the significantly different metabolites were associated with carbohydrate metabolism, amino acid metabolism, purine metabolism, and glycerolipid metabolism. Most intermediary metabolites of glycolysis and the tricarboxylic acid cycle were not significantly altered between the two altitudes, but most free fatty acids as well as β-hydroxybutyric acid were significantly lower in the high-altitude population. This may suggest that high-altitude lizards rely more on carbohydrates as their main energy fuel rather than lipids. Higher levels of phospholipids occurred in the liver of high-altitude populations, suggesting that membrane lipids may undergo adaptive remodeling in response to low-temperature stress at high altitude. In summary, this study demonstrates that metabolic profiles differ substantially between high- and low-altitude lizard populations, and that these differential metabolites and metabolic pathways can provide new insights to reveal mechanisms of adaptation to extreme environments at high altitude.

Oxygen environment and metabolic oxygen demand predictably interact to affect thermal behavior in a lizard, Sceloporus occidentalis

The climate crisis necessitates predicting how organisms respond to changing environments, but this requires understanding the mechanisms underlying thermal tolerance. The Hierarchical Mechanisms of Thermal Limitation (HMTL) hypothesis proposes that respiratory capacity and marginal stability of proteins and membranes interact hierarchically to determine thermal performance and limits. An untested prediction of the HMTL hypothesis is that behavioral anapyrexia (i.e., reduced body temperature in hypoxia) is exacerbated when metabolic demand is high. We tested this prediction by manipulating western fence lizards' (Sceloporus occidentalis) metabolic demand and oxygen environment, then measuring selected body temperatures. Lizards with elevated metabolic demand selected lower body temperatures at higher oxygen concentrations than resting lizards, but this occurred only at oxygen concentrations <12% O₂ , suggesting thermal limits are unaffected by naturally-occurring oxygen variation. Given our results and the ubiquity of behavioral anapyrexia, the HMTL hypothesis may generally explain how oxygen and temperature interactively affect reptile performance.

Physiological and Biochemical Adaptations to High Altitude in Tibetan Frogs, Nanorana parkeri

The Xizang plateau frog, N. parkeri (Anura: Dicroglossidae), is endemic to the Tibetan Plateau, ranging from 2,850 to 5,100 m above sea level. The present study explores physiological and biochemical adaptations to high altitude in this species with a particular emphasis on parameters of hematology, oxidative stress, and antioxidant defense in adult and juvenile N. parkeri collected from high (4,600 m a.s.l) and low (3,400 m a.s.l) altitudes. Hematological results showed that hemoglobin concentration ([Hb]), hematocrit (Hct), and red blood cell (RBC) counts were significantly higher in high-altitude N. parkeri. High-altitude juveniles had lower RBC sizes than low-altitude juveniles. Higher levels of GSH and GSSG were indicated only in juveniles from high altitude, not in adults. High-altitude individuals also showed lower oxidative damage, assessed as malondialdehyde (MDA) and carbonyl groups (CG) in the liver. High-altitude adults also showed higher activities of superoxide dismutase (SOD), catalase (CAT), glutathione peroxidase (GPX), and glutathione-S-transferase (GST) as well as total antioxidant capacity (T-AOC) in the liver as compared to low-altitude adults. Moreover, higher GPX activity and T-AOC were observed in the heart and brain of high-altitude adults. Liver CAT, GPX, and T-AOC showed significant increases in high-altitude juveniles. Vitamin C content was also higher in the heart of high-altitude frogs compared to low-altitude individuals. In summary, the high-altitude population of N. parkeri showed more robust hematological parameters, less oxidative damage, and stronger antioxidant defenses than the low-altitude population, all contributing to increased protection for survival in high-altitude environments.

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.

A Re-Assessment of Positive Selection on Mitochondrial Genomes of High-Elevation Phrynocephalus Lizards

Due to their integral roles in oxidative phosphorylation, mitochondrially encoded proteins represent common targets of selection in response to altitudinal hypoxia across high-altitude taxa. While previous studies revealed evidence of positive selection on mitochondrial genomes of high-altitude Phrynocephalus lizards, their conclusions were restricted by out-of-date phylogenies and limited taxonomic sampling. Using topologies derived from both nuclear and mitochondrial DNA phylogenies, we re-assessed the evidence of positive selection on the mitochondrial genomes of high-altitude Phrynocephalus. We sampled representative species from all four main lineages and sequenced the mitochondrial genome of P. maculatus, a putative sister taxon to the high-altitude group. Positive selection was assessed through two widely used branch-site tests: the branch-site model in PAML and BUSTED in HyPhy. No evidence of positive selection on mitochondrial genes was detected on branches leading to two most recent common ancestors of high-altitude species; however, we recovered evidence of positive selection on COX1 on the P. forsythii branch, which represents a reversal from high- to low-elevation environments. A positively selected site therein marked a threonine to valine substitution at position 419. We suggest this bout of selection occurred as the ancestors of P. forsythii re-colonized lower altitude environments north of the Tibetan Plateau. Despite their role in oxidative phosphorylation, we posit that mitochondrial genes are unlikely to have represented historical targets of selection for high-altitude adaptation in Phrynocephalus. Consequently, future studies should address the roles of nuclear genes and differential gene expression.

Elevation as a selective force on mitochondrial respiratory chain complexes of the Phrynocephalus lizards in the Tibetan plateau

Animals living in extremely high elevations have to adapt to low temperatures and low oxygen availability (hypoxia), but the underlying genetic mechanisms associated with these adaptations are still unclear. The mitochondrial respiratory chain can provide >95% of the ATP in animal cells, and its efficiency is influenced by temperature and oxygen availability. Therefore, the respiratory chain complexes (RCCs) could be important molecular targets for positive selection associated with respiratory adaptation in high-altitude environments. Here, we investigated positive selection in 5 RCCs and their assembly factors by analyzing sequences of 106 genes obtained through RNA-seq of all 15 Chinese Phrynocephalus lizard species, which are distributed from lowlands to the Tibetan plateau (average elevation >4,500 m). Our results indicate that evidence of positive selection on RCC genes is not significantly different from assembly factors, and we found no difference in selective pressures among the 5 complexes. We specifically looked for positive selection in lineages where changes in habitat elevation happened. The group of lineages evolving from low to high altitude show stronger signals of positive selection than lineages evolving from high to low elevations. Lineages evolving from low to high elevation also have more shared codons under positive selection, though the changes are not equivalent at the amino acid level. This study advances our understanding of the genetic basis of animal respiratory metabolism evolution in extreme high environments and provides candidate genes for further confirmation with functional analyses.

Changes in intestinal microbiota across an altitudinal gradient in the lizard Phrynocephalus vlangalii

High altitude is an important driving force in animal evolution. However, the effect of altitude on gut microbial communities in reptiles has not been examined in detail. Here, we investigated the intestinal microbiota of three populations of the lizard Phrynocephalus vlangalii living at different altitudes using 16S rRNA gene sequencing. Bacteroidetes, Firmicutes, and Proteobacteria were the most abundant phyla. Bacteroides, Odoribacter, and Parabacteroides were the most abundant genera. Significant differences in the intestinal microbiota composition were found among the three populations from different altitudes. The proportions of Verrucomicrobia and Akkermansia decreased, whereas Bacteroides increased significantly with altitude. Greater abundance of Bacteroides at higher altitude led to the fractional increase in the phylum Bacteroides relative to other phyla. Hypoxia may be the main factor that caused intestinal microbiota variation in P. vlangalii along the altitude gradient. Overall, our study suggested that the community composition and structure of intestinal microbiota of the lizard P. vlangalii varied along altitudes, and such differences likely play a certain role in highland adaptation. Our findings warrant a further study that would determine whether ambient and body temperatures play a key role in the modulation of intestinal microbiota in reptiles.

Evolutionary analysis of mitochondrially encoded proteins of toad-headed lizards, Phrynocephalus, along an altitudinal gradient

Background

Animals living at high altitude must adapt to environments with hypoxia and low temperatures, but relatively little is known about underlying genetic changes. Toad-headed lizards of the genus Phrynocephalus cover a broad altitudinal gradient of over 4000 m and are useful models for studies of such adaptive responses. In one of the first studies to have considered selection on mitochondrial protein-coding regions in an ectothermic group distributed over such a wide range of environments, we analysed nineteen complete mitochondrial genomes from all Chinese Phrynocephalus (including eight genomes sequenced for the first time). Initial analyses used site and branch-site model (program: PAML) approaches to examine nonsynonymous: synonymous substitution rates across the mtDNA tree.

Results

Ten positively selected sites were discovered, nine of which corresponded to subunits ND2, ND3, ND4, ND5, and ND6 within the respiratory chain enzyme mitochondrial Complex I (NADH Coenzyme Q oxidoreductase). Four of these sites showed evidence of general long-term selection across the group while the remainder showed evidence of episodic selection across different branches of the tree. Some of these branches corresponded to increases in altitude and/or latitude. Analyses of physicochemical changes in protein structures revealed that residue changes at sites that were under selection corresponded to major functional differences. Analyses of coevolution point to coevolution of selected sites within the ND4 subunit, with key sites associated with proton translocation across the mitochondrial membrane.

Conclusions

Our results identify mitochondrial Complex I as a target for environment-mediated selection in this group of lizards, a complex that frequently appears to be under selection in other organisms. This makes these lizards good candidates for more detailed future studies of molecular evolution.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4569-1) contains supplementary material, which is available to authorized users.

Gene expression variations in high-altitude adaptation: a case study of the Asiatic toad (Bufo gargarizans)

Background

Genome-wide investigation of molecular mechanisms for high-altitude adaptation has attracted great attention in the last few years. In order to understand the contribution of gene expression level variations to high-altitude adaptation in Asiatic toads (Bufo gargarizans), we implemented a reciprocal transplant experiment between low- and high-altitude sites and sequenced 12 transcriptomes from brain, heart, and liver tissues.

Results

A large number of genes with expression differences (DEGs) between high- and low-altitude individuals (193 fixed and 844 plastic) were identified, and the majority of them were tissue specific. Heart displayed the largest number of DEGs, both plastic and fixed. Fixed DEGs were particularly concentrated in functions associated with muscle contraction, and the majority of them were down-regulated in high-altitude individuals. Plastic DEGs were highly concentrated in several energy metabolism related functional categories, and the majority of them were also down-regulated at high-altitude environments. In liver samples, genes associated with nutrient metabolism experienced a broad-scale expression down-regulation in high-altitude toads.

Conclusions

These broadly suppressed expression patterns at high altitudes are in strong contrast to those of endothermic homeotherms, suggesting poikilothermic vertebrates may have adopted different strategies at high altitudes. Our results strongly support that both genotypic specialization and phenotypic plasticity play crucial role in adaptation to high altitude for Asiatic toads. Poikilothermic vertebrates are among the most hypoxia-tolerant animals known, and many molecular mechanisms remain elusive. We hope that our results will provide useful directions for future research.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0529-z) contains supplementary material, which is available to authorized users.

Non-Mammalian Vertebrates: Distinct Models to Assess the Role of Ion Gradients in Energy Expenditure

Animals store metabolic energy as electrochemical gradients. At least 50% of mammalian energy is expended to maintain electrochemical gradients across the inner mitochondrial membrane (H⁺), the sarcoplasmic reticulum (Ca⁺⁺), and the plasma membrane (Na⁺/K⁺). The potential energy of these gradients can be used to perform work (e.g., transport molecules, stimulate contraction, and release hormones) or can be released as heat. Because ectothermic species adapt their body temperature to the environment, they are not constrained by energetic demands that are required to maintain a constant body temperature. In fact, ectothermic species expend seven to eight times less energy than similarly sized homeotherms. Accordingly, ectotherms adopt low metabolic rates to survive cold, hypoxia, and extreme bouts of fasting that would result in energy wasting, lactic acidosis and apoptosis, or starvation in homeotherms, respectively. Ectotherms have also evolved unique applications of ion gradients to allow for localized endothermy. Endothermic avian species, which lack brown adipose tissue, have been integral in assessing the role of H⁺ and Ca⁺⁺ cycling in skeletal muscle thermogenesis. Accordingly, the diversity of non-mammalian vertebrate species allows them to serve as unique models to better understand the role of ion gradients in heat production, metabolic flux, and adaptation to stressors, including obesity, starvation, cold, and hypoxia.

Genetic signals of high-altitude adaptation in amphibians: a comparative transcriptome analysis

BACKGROUND

High-altitude adaptation provides an excellent system for studying how organisms cope with multiple environmental stressors and interacting genetic modifications. To explore the genetic basis of high-altitude adaptation in poikilothermic animals, we acquired transcriptome sequences from a high-altitude population and a low-altitude population of the Asiatic toad (Bufo gargarizans). Transcriptome data from another high-altitude amphibian, Rana kukunoris and its low-altitude relative R. chensiensis, which are from a previous study, were also incorporated into our comparative analysis.

RESULTS

More than 40,000 transcripts were obtained from each transcriptome, and 5107 one-to-one orthologs were identified among the four taxa for comparative analysis. A total of 29 (Bufo) and 33 (Rana) putative positively selected genes were identified for the two high-altitude species, which were mainly concentrated in nutrient metabolism related functions. Using SNP-tagging and F_ST outlier analysis, we further tested 89 other nutrient metabolism related genes for signatures of natural selection, and found that two genes, CAPN2 and ITPR1, were likely under balancing selection. We did not detect any positively selected genes associated with response to hypoxia.

CONCLUSIONS

Amphibians clearly employ different genetic mechanisms for high-altitude adaptation compared to endotherms. Modifications of genes associated with nutrient metabolism feature prominently while genes related to hypoxia tolerance appear to be insignificant. Poikilotherms represent the majority of animal diversity, and we hope that our results will provide useful directions for future studies of amphibians as well as other poikilotherms.

Comparative transcriptomic analysis revealed adaptation mechanism of Phrynocephalus erythrurus, the highest altitude Lizard living in the Qinghai-Tibet Plateau

BACKGROUND

Organisms living at high altitudes must overcome three major environmental challenges: hypoxia, cold, and intense UV radiation. The molecular mechanisms that enable these challenges to be overcome have mainly been studied in endothermic organisms; relatively little attention has been paid to poikilothermic species. Here, we present deep transcriptome sequencing in two closely related lizards, the high altitude-dwelling Phrynocephalus erythrurus and the lowland-dwelling P. putjatia, to identify candidate genes under positive selection and to explore the convergent evolutionary adaptation of poikilothermic animals to high altitude life.

RESULTS

More than 70 million sequence reads were generated for each species via Illumina sequencing. De novo assembly produced 56,845 and 63,140 transcripts for P. erythrurus and P. putjatia, respectively. P. erythrurus had higher Ka/Ks ratios than P. putjatia, implying an accelerated evolutionary rate in the high altitude lizard lineage. 206 gene ontology (GO) categories with accelerated evolutionary rates and 43 candidate positively selected genes were detected along the P. erythrurus lineage. Some of these GO categories have functions associated with responses to hypoxia, energy metabolism and responses to UV damage. We also found that the high-altitude ranid frog R. kukunoris had higher Ka/Ks ratios than the closely related low-altitude frog R. chensinensis, and that the functional categories with accelerated evolutionary rates in R. kukunoris overlapped extensively with those detected along the P. erythrurus lineage.

CONCLUSIONS

The mechanisms of high altitude adaptation in P. erythrurus were tentatively inferred. By comparing two pairs of low- and high-altitude poikilothermic species, we found that similar functional categories had undergone positive selection in high altitude-dwelling Phrynocephalus and Rana lineages, indicating that similar mechanisms of adaptation to high altitude might have evolved in both genera. Our findings provide important guidance for future functional studies on high altitude adaptation in poikilothermic animals.

Differences in Hematological Traits between High- and Low-Altitude Lizards (Genus Phrynocephalus)

Phrynocephalus erythrurus (Lacertilia: Agamidae) is considered to be the highest living reptile in the world (about 4500-5000 m above sea level), whereas Phrynocephalus przewalskii inhabits low altitudes (about 1000-1500 m above sea level). Here, we report the differences in hematological traits between these two different Phrynocephalus species. Compared with P. przewalskii, the results indicated that P. erythrurus own higher oxygen carrying capacity by increasing red blood cell count (RBC), hemoglobin concentration ([Hb]) and hematocrit (Hct) and these elevations could promote oxygen carrying capacity without disadvantage of high viscosity. The lower partial pressure of oxygen in arterial blood (PaO2) of P. erythrurus did not cause the secondary alkalosis, which may be attributed to an efficient pulmonary system for oxygen (O2) loading. The elevated blood-O2 affinity in P. erythrurus may be achieved by increasing intrinsic O2 affinity of isoHbs and balancing the independent effects of potential heterotropic ligands. We detected one α-globin gene and three β-globin genes with 1 and 33 amino acid substitutions between these two species, respectively. Molecular dynamics simulation results showed that amino acids substitutions in β-globin chains could lead to the elimination of hydrogen bonds in T-state Hb models of P. erythrurus. Based on the present data, we suggest that P. erythrurus have evolved an efficient oxygen transport system under the unremitting hypobaric hypoxia.