Nanorana Parkeri Research Articles

Effects of hibernation on two important contractile tissues in tibetan frogs, Nanorana parkeri: a perspective from transcriptomics and metabolomics approaches

BackgroundIn response to seasonal cold and food shortage, the Xizang plateau frogs, Nanorana parkeri (Anura: Dicroglossidae), enter a reversible hypometabolic state where heart rate and oxygen consumption in skeletal muscle are strongly suppressed. However, the effect of winter hibernation on gene expression and metabolic profiling in these two tissues remains unknown. In the present study, we conducted transcriptomic and metabolomic analyses of heart and skeletal muscle from summer- and winter-collected N. parkeri to explore mechanisms involved in seasonal hibernation.ResultsWe identified 2407 differentially expressed genes (DEGs) in heart and 2938 DEGs in skeletal muscle. Enrichment analysis showed that shared DEGs in both tissues were enriched mainly in translation and metabolic processes. Of these, the expression of genes functionally categorized as “response to stress”, “defense mechanisms”, or “muscle contraction” were particularly associated with hibernation. Metabolomic analysis identified 24 and 22 differentially expressed metabolites (DEMs) in myocardium and skeletal muscle, respectively. In particular, pathway analysis showed that DEMs in myocardium were involved in the pentose phosphate pathway, glycerolipid metabolism, pyruvate metabolism, citrate cycle (TCA cycle), and glycolysis/gluconeogenesis. By contrast, DEMs in skeletal muscle were mainly involved in amino acid metabolism.ConclusionsIn summary, natural adaptations of myocardium and skeletal muscle in hibernating N. parkeri involved transcriptional alterations in translation, stress response, protective mechanisms, and muscle contraction processes as well as metabolic remodeling. This study provides new insights into the transcriptional and metabolic adjustments that aid winter survival of high-altitude frogs N. parkeri.

Evolutionary and functional conservation of IRF7 in the Tibetan frog Nanorana parkeri.

The production of interferons (IFNs) is essential for the control of viral infections, and interferon regulatory factor 7 (IRF7) is considered as a vital regulator for the transcription of type I IFNs. Amphibians appear to possess a highly expanded type I IFN repertoire, consisting of intron-containing genes as observed in fish, and intronless genes as in other higher vertebrates. However, the knowledge on transcriptional regulatory mechanism of these two types of type I IFN genes is rather scarce in amphibians. A IRF7 gene named as Np-IRF7 was identified in Tibetan frog (Nanorana parkeri), and bioinformatic analysis revealed that the predicted protein of Np-IRF7 contains several important structural features known in IRF7. Expression analysis showed that Np-IRF7 gene was widely expressed and rapidly induced by poly(I:C) in different organs/tissues. Interestingly, luciferase reporter assay revealed that intronless IFN promoters were more effectively activated than intron-containing IFN promoter in Np-IRF7-transfected cells. Moreover, the overexpression of Np-IRF7 could induce the expression of ISGs and suppress the replication of FV3 in A6 cells. Np-IRF7 is indeed the ortholog of known IRF7, and IRF7 is structurally conserved in different lineages of vertebrates. Np-IRF7 played distinct roles in the activation of intron-containing and intronless type I IFN promoters, thus inducing the expression of interferon-stimulated antiviral effectors and providing a protection against ranavirus infection. The present research thus contributes to a better understanding of regulatory function of IRF7 in the IFN-mediated antiviral response of anuran amphibians.

Surviving winter on the Qinghai-Xizang Plateau: Extensive reversible protein phosphorylation plays a dominant role in regulating hypometabolism in hibernating Nanorana parkeri.

Changes in protein abundance and reversible protein phosphorylation (RPP) play important roles in regulating hypometabolism but have never been documented in overwintering frogs at high altitudes. To test the hypothesis that protein abundance and phosphorylation change in response to winter hibernation, we conducted a comprehensive and quantitative proteomic and phosphoproteomic analysis of the liver of the Xizang plateau frog, Nanorana parkeri, living on the Qinghai-Xizang (Tibet) Plateau (QTP). In total, 5 170 proteins and 5 695 phosphorylation sites in 1 938 proteins were quantified. Based on proteomic analysis, 674 differentially expressed proteins (438 up-regulated, 236 down-regulated) were screened in hibernating N. parkeri versus summer individuals. Functional enrichment analysis revealed that higher expressed proteins in winter were significantly enriched in immune-related signaling pathways, whereas lower expressed proteins were mainly involved in metabolic processes. A total of 4 251 modified sites (4 147 up-regulated, 104 down-regulated) belonging to 1 638 phosphoproteins (1 555 up-regulated, 83 down-regulated) were significantly changed in the liver. During hibernation, RPP regulated a diverse array of proteins involved in multiple functions, including metabolic enzymatic activity, ion transport, protein turnover, signal transduction, and alternative splicing. These changes contribute to enhancing protection, suppressing energy-consuming processes, and inducing metabolic depression. Moreover, the activities of phosphofructokinase, glutamate dehydrogenase, and ATPase were all significantly lower in winter compared to summer. In conclusion, our results support the hypothesis and demonstrate the importance of RPP as a regulatory mechanism when animals transition into a hypometabolic state.

Evolutionary conservation and antiviral effect of interferon regulatory factor 1 (IRF1) in Tibetan frog Nanorana parkeri

Interferon regulatory factor 1 (IRF1) is considered to be a vital regulator of antiviral immune responses, and the knowledge on functions of IRF1 in antiviral defense is rather limited in amphibians, which have undergone worldwide declines due to the prevalence of ranavirus infections. In the present study, an IRF1 gene named as Np-IRF1 was identified and characterized in Tibetan frog (Nanorana parkeri). Bioinformatic analysis revealed that Np-IRF1 and IRF1 from different species of vertebrates have a conserved genomic organization of eight exons and nine introns, and locate a conserved genomic region. Expression analysis showed that Np-IRF1 gene was widely expressed and induced by poly(I:C) in different organs/tissues. Moreover, the overexpression of Np-IRF1 could induce the expression of interferon-stimulated genes, and reduce the CPE and viral loads in ranavirus-infected A6 cells. The present research thus contributes to a better understanding of regulatory function of IRF1 in the antiviral immune responses of anuran amphibians.

Hepatic transcriptome and gut microbiome provide insights into freeze tolerance in the high-altitude frog, Nanorana parkeri

Among amphibians, freeze tolerance is a low-temperature survival strategy that has been well studied in several species. One influence on animal health and survival under adverse conditions is the gut microbiome. Gut microbes can be greatly affected by temperature fluctuations but, to date, this has not been addressed in high-altitude species. Nanorana parkeri (Anura: Dicroglossidae) lives at high altitudes on the Tibetan plateau and shows a good freeze tolerance. In the present study, we addressed two goals: (1) analysis of the effects of whole body freezing on the liver transcriptome, and (2) assess modifications of the gut microbiome as a consequence of freezing. We found that up-regulated genes in liver were significantly enriched in lipid and fatty acid metabolism that could contribute to accumulating the cryoprotectant glycerol and raising levels of unsaturated fatty acids. The results suggest the crucial importance of membrane adaptations and fuel reserves for freezing survival of these frogs. Down-regulated genes were significantly enriched in the immune response and inflammatory response, suggesting that energy-consuming processes are inhibited to maintain metabolic depression during freezing. Moreover, freezing had a significant effect on intestinal microbiota. The abundance of bacteria in the family Lachnospiraceae was significantly increased after freezing exposure, which likely supports freezing survival of N. parkeri. The lower abundance of bacteria in the family Peptostreptococcaceae in frozen frogs may be associated with the hypometabolic state and decreased immune response. In summary, these findings provide insights into the regulatory mechanisms of freeze tolerance in a high-altitude amphibian at the level of gene expression and gut microbiome, and contribute to enhancing our understanding of the adaptations that support frog survival in high-altitude extreme environments.

Analysis of evolutionary and functional features of the bullfrog SULT1 family

We identified the bullfrog Rana catesbeiana sulfotransferase 1 (SULT1) family from the BLAST search tool of the public databases based on the SULT1 families of Nanorana parkeri, Xenopus laevis, and Xenopus tropicalis as queries, revealing the characteristics of the anuran SULT1 family. The results showed that the anuran SULT1 family comprises six subfamilies, four of which were related to the mammalian SULT1 subfamily. Additionally, the bullfrog has two SULT1Cc subfamily members that are consistent with the characteristics of the expanded Xenopus SULT1C subfamily. Several members of the bullfrog SULT1 family were suggested to play important roles in sulfation during metamorphosis. Among these, cDNAs encoding SULT1Cc1 and SULT1Y1 were cloned, and the sulfation activity was analyzed using recombinant proteins. The affinity for 2-naphthol and 3′-phosphoadenosine 5′-phosphosulfate (PAPS) and the enzymatic reaction rate were higher in SULT1Cc1 than in SULT1Y1. Both the enzymes showed inhibitory effect of many thyroid hormones (THs) analogs on the sulfation of 2-naphthol. The potency of sulfation activities of SULT1Cc1 and SULT1Y1 against T4 indicated their possible role in the intracellular T4 clearance during metamorphosis.

Phenotypic trait variations in the frog Nanorana parkeri: differing adaptive strategies to altitude between sexes

In many animals, changes in altitude drive adaptive variation in body size. However, how other phenotypic traits change when faced with different environments has been little studied in ectotherms. In this study, we selected the high Himalaya frog Nanorana parkeri as a model species for investigating the adaptive evolution of phenotypic traits that respond to altitude in both sexes. First, we found that body mass in populations at higher altitudes was lower than at lower altitudes in females, with no difference observed in males. Second, we found significant differences in fresh liver mass, fresh heart mass, and the ratio of liver mass to body mass with increasing altitude, while hindlimb length decreased with altitude in both sexes. Third, snout-urostyle length, hindlimb length, fresh heart mass and the ratio of heart mass to body mass showed significant negative correlations with increasing altitude in both sexes. In contrast, body mass showed a significant correlation with altitude in females but not males. On the other hand, the ratio of liver mass to body mass showed a significant correlation with altitude in males but not in females. Thus, the species displayed sex-specific organ-size variation along elevation gradients, which may trade-off in life history strategies among populations. We speculate that selection favours a larger heart and liver mass to maintain a higher respiratory rate and energy consumption as an adaptation to high-altitude environments.

Integrated analysis of transcriptome and metabolome data reveals insights for molecular mechanisms in overwintering Tibetan frogs, Nanorana parkeri.

Nanorana parkeri (Anura, Dicroglossidae) is a unique frog living at high altitude on the Tibetan plateau where they must endure a long winter dormancy at low temperatures without feeding. Here, we presented a comprehensive transcriptomic and metabolomic analysis of liver tissue from summer-active versus overwintering N. parkeri, providing the first broad analysis of altered energy metabolism and gene expression in this frog species. We discovered that significantly up-regulated genes (2,397) in overwintering frogs mainly participated in signal transduction and immune responses, phagosome, endocytosis, lysosome, and autophagy, whereas 2,169 down-regulated genes were mainly involved in metabolic processes, such as oxidation-reduction process, amino acid metabolic process, fatty acid metabolic process, and TCA cycle. Moreover, 35 metabolites were shown to be differentially expressed, including 22 down-regulated and 13 up-regulated in winter. These included particularly notable reductions in the concentrations of most amino acids. These differentially expressed metabolites were mainly involved in amino acid biosynthesis and metabolism. To sum up, these findings suggest that gene expression and metabolic processes show adaptive regulation in overwintering N. parkeri, that contributes to maintaining homeostasis and enhancing protection in the hypometabolic state. This study has greatly expanded our understanding of the winter survival mechanisms in amphibians.

The highest-elevation frog provides insights into mechanisms and evolution of defenses against high UV radiation

Defense against ultraviolet (UV) radiation exposure is essential for survival, especially in high-elevation species. Although some specific genes involved in UV response have been reported, the full view of UV defense mechanisms remains largely unexplored. Herein, we used integrated approaches to analyze UV responses in the highest-elevation frog, Nanorana parkeri. We show less damage and more efficient antioxidant activity in skin of this frog than those of its lower-elevation relatives after UV exposure. We also reveal genes related to UV defense and a corresponding temporal expression pattern in N. parkeri. Genomic and metabolomic analysis along with large-scale transcriptomic profiling revealed a time-dependent coordinated defense mechanism in N. parkeri. We also identified several microRNAs that play important regulatory roles, especially in decreasing the expression levels of cell cycle genes. Moreover, multiple defense genes (i.e., TYR for melanogenesis) exhibit positive selection with function-enhancing substitutions. Thus, both expression shifts and gene mutations contribute to UV adaptation in N. parkeri. Our work demonstrates a genetic framework for evolution of UV defense in a natural environment.

Effects of acute heat exposure on oxidative stress and antioxidant defenses in overwintering frogs, Nanorana parkeri

Climate warming is intensifying on the Tibetan Plateau and poses a serious threat to amphibians that live there. Although hibernation physiology and ecology of Nanorana parkeri, a frog species native to the Tibetan plateau, has been well studied, little information is available about the physiological and biochemical responses to acute rising temperature. Here, we conducted an acute warming experiment comparing hibernating N. parkeri (acclimated at 4 °C) and frogs acutely warmed to 10 °C for 12 h, comparing indicators of oxidative stress and antioxidant defense between the two groups. Acute warming led to a significant increase in the content of oxidized glutathione and the ratio of oxidized:reduced glutathione in liver and muscle, indicating that rapid warming causing oxidative stress could be a negative factor for frogs inhabiting the Tibetan plateau. Malondialdehyde content increased by 57% only in muscle but decreased significantly in three tissues. Protein carbonyl groups rose by 15% in brain, 28% in liver, and 21% in muscle after heat exposure but vitamin C content in heart decreased by 44%. Acute heat exposure also induced tissue-specific upregulation of antioxidant enzyme activities. Catalase activity increased by 27% in heat-exposed frogs, as compared to controls, and glutathione peroxidase activity rose by 12% in brain, 30% in liver, and 12% in muscle. Glutathione-S-transferase activity was also enhanced in heart and muscle of heat-exposed frogs. Acute warming also resulted in a rise in total antioxidant capacity in all tissues examined except kidney, relative to controls. In summary, our findings show that acute heat exposure to hibernating N. parkeri causes a tissue-specific increase in oxidative damage and antioxidant defenses, with skeletal muscle being the most affected tissue. These results reveal the physiological responses to acute temperature change in overwintering N. parkeri.

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.

Physiological Ecology of Winter Hibernation by the High-Altitude Frog Nanorana parkeri.

AbstractThe Xizang plateau frog, Nanorana parkeri (Anura: Dicroglossidae), enters a dormant state in the winter in response to seasonal cold and lack of food. To investigate the physiological and ecological characteristics of overwintering in this species, we measured habitat conditions (hibernacula temperatures, body temperature, and water quality variables), morphology, metabolite concentrations, total antioxidant capacity (T-AOC), and bacteria-killing ability (BKA) of plasma during summer and winter. We found that N. parkeri hibernates underwater at the bottom of ponds (10-20-cm depth). Dissolved oxygen content in the water decreases significantly (by 12%) in the winter compared with summer, suggesting that overwintering N. parkeri may experience hypoxia. Body mass, body mass index, hepatosomatic index, and hepatic glycogen concentration all increased significantly in winter-collected frogs as compared to summer-collected individuals, indicating that overwintering N. parkeri accumulates high fuel/energy reserves to support prolonged periods of hibernation. A significant reduction in glucose, urea, and lactate concentrations in most organs may be closely related to metabolic depression in overwintering N. parkeri. Liver lactate concentration rose significantly in winter-collected frogs, suggesting that anaerobic metabolism dominates when this species overwinters. The T-AOC of plasma showed a significant reduction in winter, suggesting a reduced need for antioxidant defenses. Oppositely, the BKA of plasma increased significantly in winter versus summer, indicating that innate immunity was enhanced during overwintering. In summary, these behavioral (migrating to caves), physiological, and biochemical adjustments may be key for the successful overwintering of this high-altitude frog.

A chromosomal level genome sequence for Quasipaa spinosa (Dicroglossidae) reveals chromosomal evolution and population diversity.

Quasipaa spinosa is an Asian commercial Dicroglossidae species noted for its spiny chest found in adult males. Here, we report the first chromosomal level Q. spinosa genome employing PacBio long read sequencing and high-resolution chromosome conformation capture (Hi-C) technology. The total length of the final assembled genome was 2,839,292,578bp, with contig N50 of 3.79Mb and scaffold N50 of 327.44Mb. Approximately 99.30% of the length of the assembled genome sequences were anchored to 13 chromosomes with the assistance of Hi-C reads. A total of 26,173 protein-coding genes were predicted, and 95.98% of the genes were functionally annotated. The annotated genes covered a total of 92.10% of the complete vertebrate core gene set according to the BUSCO pipeline evaluation. Approximately 41million years ago, Q. spinosa began to diverge from its dicroglossid sister taxon Nanorana parkeri. The Q. spinosa genome revealed obvious chromosomal fissions compared with Xenopus tropicalis, which probably represented a specific chromosome evolutionary history within frogs. Population analysis showed that Chinese Q. spinosa could be divided into eastern and western genetic clusters, with the western population showing higher diversity than the eastern population. The effective population size of Q. spinosa showed a continuously decreasing trend from one million years ago to 10,000years ago. In summary, this study sheds light on Q. spinosa evolution and population differentiation, providing a valuable genomic resource for further biological and genetic studies on this species, and other closely related frog taxa.

Antioxidant and non-specific immune defenses in partially freeze-tolerant Xizang plateau frogs, Nanorana parkeri.

The Xizang plateau frog Nanorana parkeri can tolerate brief and partial freezing of their body. To determine the significant role of antioxidant defense and non-specific immune defense in freezing survival of this species, we assayed parameters of oxidative damage, antioxidant defense and non-specific immune enzymes during freezing exposure (-2°C for 12h) in five organs (heart, brain, liver, kidney, and skeletal muscle). The results showed that freezing led to a significant rise in the content of malondialdehyde (MDA) and carbonyl groups (CG) in brain, liver and kidney tissues. The activities of superoxide dismutase (SOD), catalase (CAT) and glutathione peroxidase (GPX) increased significantly in brain and liver tissues with an augmentation of total antioxidant capacity (T-AOC). Apparent increments in muscle SOD activity and liver GST activity were also observed during freezing exposure. Vitamin C content significantly decreased in liver and kidney but a significant increase occurred in brain. Activities of non-specific immune enzymes, acid phosphatase (ACP) and alkaline phosphatase (AKP), were also assessed. ACP activity was significantly reduced in all five tissues tested whereas AKP activity decreased significantly in four tissues but rose in brain. In summary, freezing is accompanied by oxidative stress in the high-altitude frog, N. parkeri, as documented by increases in the content of MDA and CG in tissues. Freezing exposure also induced tissue-specific changes in the antioxidant defenses showing that activation of antioxidant systems is a part of the survival strategy of this in a high-altitude frog during freezing. Such up-regulation of antioxidant enzymes suggests a particularly important role for them in the liver and brain, serving as an anticipatory mechanism to deal with the ROS challenge during freeze/thaw episodes. Our findings contribute to extending the current understanding of the mechanisms of freeze tolerance in high-altitude frogs.

Metabolic responses of plasma to extreme environments in overwintering Tibetan frogs Nanorana parkeri: a metabolome integrated analysis

Many animals lower their metabolic rate in response to low temperatures and scarcity of food in the winter in phenomena called hibernation or overwintering. Living at high altitude on the Tibetan Plateau where winters are very cold, the frog Nanorana parkeri, survives in one of the most hostile environments on Earth but, to date, relatively little is known about the biochemical and physiological adjustments for overwintering by this species. The present study profiled changes in plasma metabolites of N. parkeri between winter and summer using UHPLC-QE-MS non-target metabolomics in order to explore metabolic adaptations that support winter survival. The analysis showed that, in total, 11 metabolites accumulated and 95 were reduced in overwintering frogs compared with summer-active animals. Metabolites that increased included some that may have antioxidant functions (canthaxanthin, galactinol), act as a metabolic inhibitor (mono-ethylhexylphthalate), or accumulate as a product of anaerobic metabolism (lactate). Most other metabolites in plasma showed reduced levels in winter and were generally involved in energy metabolism including 11 amino acids (proline, isoleucine, leucine, valine, phenylalanine, tyrosine, arginine, tryptophan, methionine, threonine and histidine) and 4 carbohydrates (glucose, citrate, succinate, and malate). Pathway analysis indicated that aminoacyl-tRNA biosynthesis, phenylalanine, tyrosine and tryptophan biosynthesis, and nitrogen metabolism were potentially the most prominently altered pathways in overwintering frogs. Changes to these pathways are likely due to fasting and global metabolic depression in overwintering frogs. Concentrations of glucose and urea, commonly used as cryoprotectants by amphibians that winter on land, were significantly reduced during underwater hibernation in N. parkeri. In conclusion, winter survival of the high-altitude frog, N. parkeri was accompanied by substantial changes in metabolomic profiles and this study provides valuable information towards understanding the special adaptive mechanisms of N. parkeri to winter stresses.

Evolution of hes gene family in vertebrates: the hes5 cluster genes have specifically increased in frogs

Backgroundhes genes are chordate homologs of Drosophila genes, hairy and enhancer of split, which encode a basic helix-loop-helix (bHLH) transcriptional repressor with a WRPW motif. Various developmental functions of hes genes, including early embryogenesis and neurogenesis, have been elucidated in vertebrates. However, their orthologous relationships remain unclear partly because of less conservation of relatively short amino acid sequences, the fact that the genome was not analyzed as it is today, and species-specific genome duplication. This results in complicated gene names in vertebrates, which are not consistent in orthologs. We previously revealed that Xenopus frogs have two clusters of hes5, named “the hes5.1 cluster” and “the hes5.3 cluster”, but the origin and the conservation have not yet been revealed.ResultsHere, we elucidated the orthologous and paralogous relationships of all hes genes of human, mouse, chicken, gecko, zebrafish, medaka, coelacanth, spotted gar, elephant shark and three species of frogs, Xenopus tropicalis (X. tropicalis), X. laevis, Nanorana parkeri, by phylogenetic and synteny analyses. Any duplicated hes5 were not found in mammals, whereas hes5 clusters in teleost were conserved although not as many genes as the three frog species. In addition, hes5 cluster-like structure was found in the elephant shark genome, but not found in cyclostomata.ConclusionThese data suggest that the hes5 cluster existed in the gnathostome ancestor but became a single gene in mammals. The number of hes5 cluster genes were specifically large in frogs.

Analysis of convergent and parallel amino acid substitutions in the HSP90AA1 gene among high-elevation anurans

A significant amount of convergent and parallel amino acid substitutions in the HSP90AA1 gene has been detected among four species of high-elevation anurans: Bufo tibetanus, Scutiger boulengeri, Rana kukunoris, and Nanorana parkeri. As HSP proteins are involved in response to environmental stress, it is possible these mutations play a role in high-elevation adaptation. In this study, I investigated the functional consequences of these substitutions and inferred their potential links to adaptation. I examined HSP90AA1 sequences of 13 anuran species previously studied. Using PROVEAN, I isolated three deleterious mutations: P65S, K195A, and _199I, each shared between two of the high-elevation species. I further analyzed the protein structure, stability change, and structural damage using model predictions. Based on its buried location and cavity expansion, P65S was predicted to most likely alter protein function. Furthermore, I examined HSP90AA1 sequences of over 100 other animal species available from public databases and found that serine at site 65 is ubiquitously present in cold-water fish, suggesting the substitution is related to cold adaptation. Alanine at site 195 and isoleucine at site 199 were not found in any other species, but these substitutions also might impact protein function as they are predicted to be destabilizing and their ancestral residues have reported post-translational modifications in orthologs. Tests of protein function and an investigation of more sequences from high-elevation species would help to further link these substitutions to adaptation, particularly P65S. Identifying mutations that contribute to high-elevation adaptation would aid in uncovering the molecular mechanisms of adaptation.

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

BackgroundTo 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.ResultsGenomic 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 Rconv 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.ConclusionsThere 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.

Freeze tolerance and the underlying metabolite responses in the Xizang plateau frog, Nanorana parkeri.

The frog Nanorana parkeri (Dicroglossidae) is endemic to the Tibetan Plateau, and overwinters shallow pond within damp caves for up to 6 months. Herein, we investigate the freeze tolerance of this species and profile changes in liver and skeletal muscle metabolite levels using an untargeted LC-MS-based metabolomic approach to investigate molecular mechanisms that may contribute to freezing survival. We found that three of seven specimens of N. parkeri could survive after being frozen for 12h at -2.0°C with 39.91% ± 5.4% (n = 7) of total body water converted to ice. Freezing exposure induced partial dehydration of the muscle, which contributed to decreasing the amount of freezable water within the muscle and could be protective for the myocytes themselves. A comparative metabolomic analysis showed that freezing elicited significant responses, and a total of 33 and 36 differentially expressed metabolites were identified in the liver and muscle, respectively. These metabolites mainly participate in alanine, aspartic acid and glutamic acid metabolism, arginine and proline metabolism, and D-glutamine and D-glutamate metabolism. After freezing exposure, the contents of ornithine, melezitose, and maltotriose rose significantly; these may act as cryoprotectants. Additionally, the content of 8-hydroxy-2-deoxyguanine, 7-Ketocholesterol and hypoxanthine showed a marked increase, suggesting that freezing induced oxidative stress in the frogs. In summary, N. parkeri can tolerate a brief and partial freezing of their body, which was accompanied by substantial changes in metabolomic profiles after freezing exposure.

Genomewide analysis of microsatellite markers based on sequenced database in two anuran species

Eukaryotic and prokaryotic cell genomes exhibit multiple microsatellites. In this study, we characterized microsatellites in genomes and genes of Nanorana parkeri and Xenopus laevis. This characterization was used for gene ontology (GO) analysis of coding sequences (CDS). Compared to the genome of N. parkeri, the genome of X. laevis is larger and contains more number of microsatellites, but the diversity of both species are similar. Trinucleotide repeats in the genome of N. parkeri and dinucleotide and tetranucleotide repeats inthe genome of X. laevis were the most diverse. In both the species, diversity of microsatellites was highest in intergenic regions, followed by intron and exon regions, and lowest in coding regions. Microsatellites in CDS are thus subject to higher selective pressure. Many microsatellites are concentrated upstream and downstream of genes in both species, suggesting suppression of repeats in the middle of protein-CDS. Repeats are enriched in regions near gene termini purely due to the biophysical constraints of protein structure. In GO analysis, two and five unique GO terms, only found in N. parkeri and X. laevis, respectively, indicate advantageous mutations during species evolution. Biological process, cellular component and molecular function ontology reflected in the GO analysis predicted that the microsatellites located in CDS can alter protein function and may provide a molecular basis for species adaptation to new and changing environments.