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

Does brain size of Asiatic toads (Bufo gargarizans) trade-off with other energetically expensive organs along altitudinal gradients?

Brain size variation is often attributed to energetic trade-offs with other metabolically expensive tissues and organs, which is a prediction of the expensive brain hypothesis (EBH). Here we examine Asiatic toads (Bufo gargarizans) along altitudinal gradients and test size trade-offs between the brain and four visceral organs (heart, liver, alimentary tract, and kidney) with altitude. Body size and scaled mass index (a proxy for total energy intake) decline with altitude, implying stronger energetic constraints at high altitudes. Relative brain size decreases along altitudinal gradients, while visceral organs mostly increase in relative sizes. Using structural equation modeling, a significant negative relationship between brain size and a latent variable "budget," which represents the energy allocation to the four visceral organs, is detected among high-altitudinal toads. Heart appears to have the largest and most consistent response to changes in energy allocation. No such relationships are observed among toads at middle- and low-altitudes, where high energy intake may allow individuals to forego energetic trade-offs. When applying EBH to poikilotherms, a great emphasis should be placed on total energy intake in addition to energy allocation. Future research on EBH will benefit from more intra-specific comparisons and the evaluation of fitness consequences beyond energy limitation.

Genomic Insights Into Genetic Basis of Evolutionary Conservatism and Innovation in Frogs

Examining closely related species evolving in similar environments offers valuable insights into the mechanisms driving phylogenetic conservatism and evolutionary lability. This can elucidate the intricate relationship between inheritance and environmental factors. Nonetheless, the precise genomic dynamics and molecular underpinnings of this process remain enigmatic. This study explores the evolutionary conservatism and adaptation exhibited by two closely related high-altitude frog species: Nanorana parkeri and N. pleskei. We assembled a high-quality genome for Tibetan N. pleskei and compared it to the genomes of N. parkeri and their lowland relatives. Our findings reveal that these two Tibetan frog species diverged approximately 16.6 million years ago, pointing to a possible ancestral colonization of high-elevation habitats. Following this colonization, significant adaptive evolution occurred in both coding and non-coding regions of the ancestral lineage. This evolution led to notable phenotypic alterations, as evidenced by the reduced body size. Also, due to purifying selection, most ancestral adaptive features persisted in descendant species, indicating a strong element of evolutionary conservatism. However, descendant species evolved novel adaptations to exacerbated environmental challenges in the Tibet Plateau, mainly related to hypoxia response. Furthermore, our analysis underscores the critical role of regulatory variations in descendant adaptive evolution. Notably, hub genes in networks, such as EGLN3, accumulated more variations in regulatory regions as they were transmitted from ancestors to descendants. In sum, our study sheds light on the profound and lasting impact of genetic heritage on species' adaptive evolution.

Ecological adaptations of amphibians to environmental changes along an altitudinal gradient (Case Study: Bufo gargarizans) from phenotypic and genetic perspectives

BACKGROUND

Organisms have evolved a range of phenotypic and genetic adaptations to live in different environments along an altitudinal gradient. Herein, we studied the widely distributed Chinese toad, Bufo gargarizans, as a model and used an integrated phenotype-genotype approach to assess adaptations to different altitudinal environments.

RESULTS

Comparison of populations from four altitudes (50 m, 1200 m, 2300 m, and 3400 m) showed more effective defenses among high-altitude toads. These included thickened epidermis, more epidermal capillaries and granular glands, greater gland size in skin, and higher antioxidant enzyme activities in plasma. High-altitude toads also showed increased erythrocytes and hematocrit and elevated hemoglobin concentration, potentially improving oxygen delivery. Elevated altitude led to a metabolic shift from aerobic to anaerobic metabolism, and high-altitude populations favored carbohydrates over fatty acids to fuel for energy metabolism. Differentially expressed genes were associated with adaptive phenotypic changes. For instance, expression of genes associated with fatty acid metabolism showed greater suppression at high altitude (3400 m), consistent with decreased flux of β-hydroxybutyric acid and lower free fatty acids levels. Moreover, down-regulation of genes involved in carbon metabolism processes at high altitude (3400 m) were coincident with reduced TCA cycle flux. These results suggest that high-altitude toads adopt a metabolic suppression strategy for survival under harsh environmental conditions. Moreover, the hypoxia-inducible factor signaling cascade was activated at high altitude.

CONCLUSIONS

Collectively, these results advance our comprehension of adaptation to high-altitude environments by revealing physiological and genetic mechanisms at work in Chinese toads living along altitudinal gradients.

Context-dependent effects of thermal acclimation on physiological correlates of animal personality in Asiatic toads

Persistent individual variation in behaviour, or 'personality', is a widespread phenomenon in animals, and understanding the evolution of animal personality is a key task of current biology. Natural selection has been proposed to promote the integration of personality with animal 'intrinsic states', such as metabolic or endocrine traits, and this integration varies with ecological conditions. However, these external ecological modulatory effects have rarely been examined. Here, we investigate the effects of thermal acclimation on between-individual covariations between physiology and behaviour in Asiatic toads (Bufo gargarizans) along an altitudinal gradient. Our results reveal that the thermal modulatory effects on the covariations depend on the altitudinal population. Specifically, at low altitudes, between-individual covariations are highly plastic, with risk-taking behaviour covarying with baseline glucocorticoids (GCs) under warm acclimation, but risk-taking and exploration behaviour covarying with resting metabolic rate (RMR) under cold acclimation. In contrast, between-individual covariations are relatively fixed at high altitudes, with risk-taking behaviour consistently covarying with baseline GCs. Furthermore, at low altitudes, changes in covariations between RMR and personality are associated with adjustment of energy management models. Evidently, animal physiological states that determine or covary with personality can adapt according to the seasonal thermal environment and the thermal evolutionary background of populations. Our findings highlight the importance of a multi-system physiological approach to understand the evolution of animal personality.

Evolutionary Insights into the Relationship of Frogs, Salamanders, and Caecilians and Their Adaptive Traits, with an Emphasis on Salamander Regeneration and Longevity

The extant amphibians have developed uncanny abilities to adapt to their environment. I compared the genes of amphibians to those of other vertebrates to investigate the genetic changes underlying their unique traits, especially salamanders' regeneration and longevity. Using the well-supported Batrachia tree, I found that salamander genomes have undergone accelerated adaptive evolution, especially for development-related genes. The group-based comparison showed that several genes are under positive selection, rapid evolution, and unexpected parallel evolution with traits shared by distantly related species, such as the tail-regenerative lizard and the longer-lived naked mole rat. The genes, such as EEF1E1, PAFAH1B1, and OGFR, may be involved in salamander regeneration, as they are involved in the apoptotic process, blastema formation, and cell proliferation, respectively. The genes PCNA and SIRT1 may be involved in extending lifespan, as they are involved in DNA repair and histone modification, respectively. Some genes, such as PCNA and OGFR, have dual roles in regeneration and aging, which suggests that these two processes are interconnected. My experiment validated the time course differential expression pattern of SERPINI1 and OGFR, two genes that have evolved in parallel in salamanders and lizards during the regeneration process of salamander limbs. In addition, I found several candidate genes responsible for frogs' frequent vocalization and caecilians' degenerative vision. This study provides much-needed insights into the processes of regeneration and aging, and the discovery of the critical genes paves the way for further functional analysis, which could open up new avenues for exploiting the genetic potential of humans and improving human well-being.

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.

Comparative transcriptomes of three different skin sites for the Asiatic toad (Bufo gargarizans)

Toads release toxic dry secretions from glands in their skin. Toxin possesses a wide range of biological effects, but little is known about its specific gene expression pattern and regulatory mechanisms. The Asiatic toad (Bufo gargarizans) is widely used to produce toxin. Here, we explored the gene expression of 30 tissue samples from three different skin sites (parotoid gland, dorsal skin, and abdomen skin) of B. gargarizans. After de novo assembly, 783,130 unigenes with an average length of 489 bp (N50 = 556 bp) were obtained. A total of 9,248 significant differentially expressed genes (DEGs) were detected. There were 8,819 DEGs between the parotoid gland and abdomen skin and 1,299 DEGs between the dorsal skin and abdomen skin, while only 1,283 DEGs were obtained between the parotoid gland and dorsal skin. Through enrichment analysis, it was found that the detected differential gene expressions corresponded to the different functions of different skin sites. Our key findings were the genetic expression of toxin secretion, the protection function of skin, and the related genes such as HSD3B, Cyp2c, and CAT, LGALS9. In conclusion, we provide useful transcript resources to study the gene expression and gene function of B. gargarizans and other amphibians. The detected DEGs between different sites of the skin provided better insights into the genetic mechanisms of toxin secretion and the protection function of skin for amphibians.

Comparison of hematological traits and oxygenation properties of hemoglobins from highland and lowland Asiatic toad (Bufo gargarizans)

The Asiatic toad (Bufo gargarizans) belonging to the family of Bufonidae (Anura: Amphibia) is successfully residing on the Qinghai-Tibetan Plateau (QTP). To investigate whether the oxygen delivery undergoes adaptive adjustments to high-altitude environments in Asian toads inhabiting the QTP (Zoige County, 3446 m), choosing low-altitude populations (Chengdu City, 500 m) as control, we measured hematological traits, O₂ affinities of whole blood, Hb-O₂ affinities of purified Hbs, their sensitivities to temperature, and allosteric effectors (H⁺, Cl^- and ATP). Our results showed that high-altitude Asiatic toads possessed significantly increased hemoglobin concentration, hematocrit, and red blood cell count, but significantly decreased erythrocyte volume compared with low-altitude toads. The whole blood and purified Hbs of high-altitude Asiatic toads both exhibited significantly higher O₂ affinities compared with low-altitude toads. Substantially increased intrinsic Hb-O₂ affinities of high-altitude Asiatic toads Hbs are likely to be the main reason for its elevated Hb-O₂ affinities given the anionic cofactor sensitivities of high- and low-altitude toads were similar. The Hbs of high-altitude toads were also characterized by distinctly strong Bohr effects at the low temperature and low-temperature sensitivities. The adaptive adjustments of hematological traits could enhance the blood-O₂ carrying capacity of high-altitude Asiatic toads. The increased Hb-O₂ affinities could safeguard the pulmonary O₂ uploading under hypoxia. The strong Bohr effects at the low temperature could help the release of O₂ in metabolic tissues and cold limbs, while low-temperature sensitivity could minimize the effect of temperature fluctuation on the Hb-O₂ affinity.

Metabolic cold adaptation in the Asiatic toad: intraspecific comparison along an altitudinal gradient

The metabolic cold adaptation (MCA) hypothesis predicts an increase in metabolic rate and thermal sensitivity of poikilotherms from cold environments as compared to those from warm environments, when measured under standardized conditions. This compensatory response is also expected to evolve in life history and behavioral traits if the reductions in these phenotypic traits at low temperature involves in a reduction in fitness. We investigated the extent to which the level of energy intake (measured as feeding rate), energy turnover (measured as standard metabolic rate, SMR) and the energy budget (energy allocation to growth and physical activity) are influenced by climatic conditions in three populations of the Asiatic toad (Bufo gargarizans) distributed across an altitudinal gradient of 1350 m in the Qionglai Mountains of Western China. We found a similar thermal reaction norm of SMR at both population and individual levels; therefore, the data did not support the MCA hypothesis. However, there was a co-gradient variation (CoGV) for mass change rate in which the high and medium altitudinal populations displayed slower mass change rates than their counterparts from low altitudes. Moreover, this CoGV pattern was accompanied by a low feeding rate and high physical activity for the high- and medium-altitude populations. Our results highlight that adjustments in energy intake and energy allocation to behaviors, but not energy allocation to metabolism of maintenance, could act as an energetic strategy to accommodate the varied growth efficiency in Asiatic toads along an altitudinal gradient.

Transcriptome analysis reveals plasticity in gene regulation due to environmental cues in Primula sikkimensis, a high altitude plant species

BACKGROUND

Studying plasticity in gene expression in natural systems is crucial, for predicting and managing the effects of climate change on plant species. To understand the contribution of gene expression level variations to abiotic stress compensation in a Himalaya plant (Primula sikkimensis), we carried out a transplant experiment within (Ambient), and beyond (Below Ambient and Above Ambient) the altitudinal range limit of species. We sequenced nine transcriptomes (three each from each altitudinal range condition) using Illumina sequencing technology. We compared the fitness variation of transplants among three transplant conditions.

RESULTS

A large number of significantly differentially expressed genes (DEGs) between below ambient versus ambient (109) and above ambient versus ambient (85) were identified. Transcripts involved in plant growth and development were mostly up-regulated in below ambient conditions. Transcripts involved in signalling, defence, and membrane transport were mostly up-regulated in above ambient condition. Pathway analysis revealed that most of the genes involved in metabolic processes, secondary metabolism, and flavonoid biosynthesis were differentially expressed in below ambient conditions, whereas most of the genes involved in photosynthesis and plant hormone signalling were differentially expressed in above ambient conditions. In addition, we observed higher reproductive fitness in transplant individuals at below ambient condition compared to above ambient conditions; contrary to what we expect from the cold adaptive P. sikkimensis plants.

CONCLUSIONS

We reveal P. sikkimensis's capacity for rapid adaptation to climate change through transcriptome variation, which may facilitate the phenotypic plasticity observed in morphological and life history traits. The genes and pathways identified provide a genetic resource for understanding the temperature stress (both the hot and cold stress) tolerance mechanism of P. sikkimensis in their natural environment.

Whole-genome de novo sequencing reveals unique genes that contributed to the adaptive evolution of the Mikado pheasant

Background

The Mikado pheasant (Syrmaticus mikado) is a nearly endangered species indigenous to high-altitude regions of Taiwan. This pheasant provides an opportunity to investigate evolutionary processes following geographic isolation. Currently, the genetic background and adaptive evolution of the Mikado pheasant remain unclear.

Results

We present the draft genome of the Mikado pheasant, which consists of 1.04 Gb of DNA and 15,972 annotated protein-coding genes. The Mikado pheasant displays expansion and positive selection of genes related to features that contribute to its adaptive evolution, such as energy metabolism, oxygen transport, hemoglobin binding, radiation response, immune response, and DNA repair. To investigate the molecular evolution of the major histocompatibility complex (MHC) across several avian species, 39 putative genes spanning 227 kb on a contiguous region were annotated and manually curated. The MHC loci of the pheasant revealed a high level of synteny, several rapidly evolving genes, and inverse regions compared to the same loci in the chicken. The complete mitochondrial genome was also sequenced, assembled, and compared against four long-tailed pheasants. The results from molecular clock analysis suggest that ancestors of the Mikado pheasant migrated from the north to Taiwan about 3.47 million years ago.

Conclusions

This study provides a valuable genomic resource for the Mikado pheasant, insights into its adaptation to high altitude, and the evolutionary history of the genus Syrmaticus, which could potentially be useful for future studies that investigate molecular evolution, genomics, ecology, and immunogenetics.