Molecular Convergent Evolution of the MYBPC2 Gene Among Three High-Elevation Amphibian Species

Potential Effects of High Temperature and Heat Wave on Nanorana pleskei Based on Transcriptomic Analysis

In the context of climate change, understanding how indigenous amphibians of the Qinghai-Tibet plateau react to stresses and their coping mechanisms could be crucial for predicting their fate and successful conservation. A liver transcriptome for Nanorana pleskei was constructed using high-throughput RNA sequencing, and its gene expression was compared with frogs acclimated under either room temperature or high temperature and also heat wave exposed ones. A total of 126,465 unigenes were produced, with 66,924 (52.92%) of them being annotated. Up to 694 genes were found to be differently regulated as a result of abnormal temperature acclimatization. Notably, genes belonging to the heat shock protein (HSP) family were down-regulated in both treated groups. Long-term exposure to high-temperature stress may impair the metabolic rate of the frog and trigger the body to maintain a hypometabolic state in an effort to survive challenging times. During heat waves, unlike the high-temperature group, mitochondrial function was not impaired, and the energy supply was largely normal to support the highly energy-consuming metabolic processes. Genes were more transcriptionally suppressed when treated with high temperatures than heat waves, and the body stayed in low-energy states for combating these long-term adverse environments to survive. It might be strategic to preserve initiation to executive protein activity under heat wave stress. Under both stress conditions, compromising the protection of HSP and sluggish steroid activity occurred in frogs. Frogs were more affected by high temperatures than by heat waves.

Molecular convergent and parallel evolution among four high-elevation anuran species from the Tibetan region

BACKGROUND

To date, evidence for the relative prevalence or rarity of molecular convergent and parallel evolution is conflicting, and understanding of how these processes contribute to adaptation is limited. We compared four high-elevation anuran species (Bufo tibetanus, Nanorana parkeri, Rana kukunoris and Scutiger boulengeri) from the Tibetan region, and examined convergent and parallel amino acid substitutions between them and how they may have contributed to high-elevation adaptation.

RESULTS

Genomic data of the four high-elevation species and eight of their low-elevation close relatives were gathered. A total of 1098 orthologs shared by all species were identified. We first conducted pairwise comparisons using Zhang and Kumar's test. Then, the R_conv index was calculated and convergence/divergence correlation plotting was conducted. Furthermore, genes under positive selection and with elevated evolutionary rate were examined. We detected a large number of amino acid sites with convergent or parallel substitutions. Several pairs of high-elevation species, in particular, R. kukunoris vs N. parkeri and B. tibetanus vs S. boulengeri, had excessive amounts of convergent substitutions compared to neutral expectation. Nevertheless, these sites were mostly concentrated in a small number of genes (3-32), and no genome-wide convergence was detected. Furthermore, the majority of these convergent genes were neither under detectable positive selection nor had elevated evolutionary rates, although functional prediction analysis suggested some of the convergent genes could potentially contribute to high-elevation adaptation.

CONCLUSIONS

There is a substantial amount of convergent evolution at the amino-acid level among high-elevation amphibians, although these sites are concentrated in a few genes, not widespread across the genomes. This may attribute to the fact that all the target species are from the same environment. The relative prevalence of convergent substitutions among high-elevation amphibians provides an excellent opportunity for further study of molecular convergent evolution.