Influence of High Temperatures and Heat Wave on Thermal Biology, Locomotor Performance, and Antioxidant System of High-Altitude Frog Nanorana pleskei Endemic to Qinghai-Tibet Plateau

Investigating how highland amphibians respond to changes in ambient temperature may be of great significance for their fate prediction and effective conservation in the background of global warming. Here, using field individuals as the control group, we investigated the influence of high temperatures (20.5 and 25.5°C) and heat wave (15–26.6°C) on the thermal preference, critical thermal limits, locomotor performance, oxidative stress, and antioxidant enzyme activities in high-altitude frog Nanorana pleskei (3,490 m) endemic to the Qinghai-Tibet Plateau (QTP). After 2 weeks of acclimation to high temperatures and heat wave, the thermal preference (Tpref), critical thermal maximum (CTmax), and range of tolerable temperature significantly increased, while the critical thermal minimum (CTmin) was significantly decreased. The total time of jump to exhaustion significantly decreased, and burst swimming speed significantly increased in frogs acclimated in the high temperature and heat wave groups compared with the field group. In the high temperature group, the level of H2O2 and lipid peroxide (malondialdehyde, MDA), as well as the activities of glutathione peroxidase (GPX), glutathione reductase (GR), catalase (CAT), superoxide dismutase (SOD), and total antioxidant capacity (T-AOC) significantly increased in the liver or muscle. However, in the heat wave group, the MDA content significantly decreased in the liver, and antioxidants activities decreased in the liver and muscle except for CAT activities that were significantly increased in the liver. These results indicated that N. pleskei could respond to the oxidative stress caused by high temperatures by enhancing the activity of antioxidant enzymes. The heat wave did not appear to cause oxidative damage in N. pleskei, which may be attributed to the fact that they have successfully adapted to the dramatic temperature fluctuations on the QTP.

Investigating how highland amphibians respond to changes in ambient temperature may be of great significance for their fate prediction and effective conservation in the background of global warming. Here, using field individuals as the control group, we investigated the influence of high temperatures (20.5 and 25.5 • C) and heat wave (15-26.6 • C) on the thermal preference, critical thermal limits, locomotor performance, oxidative stress, and antioxidant enzyme activities in high-altitude frog Nanorana pleskei (3,490 m) endemic to the Qinghai-Tibet Plateau (QTP). After 2 weeks of acclimation to high temperatures and heat wave, the thermal preference (T pref ), critical thermal maximum (CT max ), and range of tolerable temperature significantly increased, while the critical thermal minimum (CT min ) was significantly decreased. The total time of jump to exhaustion significantly decreased, and burst swimming speed significantly increased in frogs acclimated in the high temperature and heat wave groups compared with the field group. In the high temperature group, the level of H 2 O 2 and lipid peroxide (malondialdehyde, MDA), as well as the activities of glutathione peroxidase (GPX), glutathione reductase (GR), catalase (CAT), superoxide dismutase (SOD), and total antioxidant capacity (T-AOC) significantly increased in the liver or muscle. However, in the heat wave group, the MDA content significantly decreased in the liver, and antioxidants activities decreased in the liver and muscle except for CAT activities that were significantly increased in the liver. These results indicated that N. pleskei could respond to the oxidative stress caused by high temperatures by enhancing the activity of antioxidant enzymes. The heat wave did not appear to cause oxidative damage in N. pleskei, which may be attributed to the fact that they have successfully adapted to the dramatic temperature fluctuations on the QTP.

INTRODUCTION
In the background of global warming, the annual average air temperature in Qinghai Tibet Plateau (QTP) is increasing, with increase in frequency and intensity of extreme temperature events (Stocker et al., 2014;Zhao et al., 2014;Yi et al., 2018). Poikilothermic vertebrates, such as amphibians and nonavian reptiles, although have behavioral and weak physiological thermoregulatory mechanisms, their growth, reproduction, foraging, immunity, and competitiveness are easily affected by environmental temperatures (Panov and Mcqueen, 1998;Cherkasov et al., 2006;Pörtner and Farrell, 2008;Peng et al., 2020). Since 1980, 122 species have gone extinct, and more than 32% of amphibian species are vulnerable, endangered, or critically endangered because of climate change, habitat loss and so on (Collins and Storfer, 2003;Stuart et al., 2004;Mendelson et al., 2006;Wake and Vredenburg, 2008). Han et al. (2016) have demonstrated that temperature significantly affects the daily food intake, sprint speed, and length of continuous locomotion in Qinghai Plateau toad-headed lizard Phrynocephalus vlangalii. Peng et al. (2020) found that high temperature accelerates the development of tadpoles and reduces their body size during metamorphosis in Tibetan brown frog Rana kukunoris. The body mass of the corn snake Pantherophis guttatus is decreased significantly in a simulated heat wave (Stahlschmidt et al., 2017). Therefore, understanding how highland amphibians respond to climate warming and heat wave events may have great significance for their effective conservation.
The thermal biology of ectotherms comprise preferred body temperature (T pref ) and thermal tolerance (critical thermal maximum CT max , and critical thermal minimum, CT min ) (Angilletta et al., 2010b). The T pref , CT max , and CT min of ectotherms are affected by many factors, such as thermal acclimation trials (Avalos et al., 2020), light cycle (Lapwong et al., 2020), and ontogeny (Wollmuth et al., 1987). For example, T pref is higher in the warm/wet season than in the cold/dry season in the terrestrial toad Rhinella icterica (César et al., 2018). Furthermore, the common coqui frog Eleutherodactylus coqui, acclimated to a cold temperature (16 • C, 30 days), has CT min lower than that of those acclimated to a high temperature (25 • C) (Haggerty, 2016). The plasticity of thermal biology of amphibians can help them to adapt to the warming ecological environment and meet the needs of different physiology (Hernandez-Sandoval et al., 2018). The locomotion performance of individuals largely determine their adaptability, such as foraging, reproduction, and avoiding predators (Clusella-Trullas et al., 2010), while the locomotion performance of ectotherms is closely related to ambient temperature (Chown and Nicolson, 2004;Stillwell and Fox, 2005). The individual locomotion ability gradually decreases once the ambient temperature goes beyond optimal temperature range (Hertz et al., 1983;Bennett, 1990). However, some locomotion properties, i.e., stamina and burst performance, show plasticity after acclimation to different temperatures (Clusella-Trullas et al., 2010). For instance, the tropical clawed frog Xenopus tropicalis, acclimated to lower temperature (24 • C, 2 months), has stronger stamina than individuals acclimated to high temperature (29 • C, 2 months) (Padilla et al., 2019). The soft-shelled turtle Pelodiscus sinensis, acclimated to higher temperature (30 • C, 4 weeks), has faster swimming speed than those acclimated to a low temperature (10 and 20 • C, 4 weeks) (Wu et al., 2013). Nevertheless, there are limited studies that reported whether high-altitude amphibians also exhibit plasticity in thermal biology and locomotion performance after acclimation to high temperatures and heat wave.
The oxidative stress and antioxidant defense ability of organism are closely related to the ambient temperature . A high-temperature environment can increase metabolism and accelerate oxygen consumption, which will increase the production of reactive oxygen species. For example, hydrogen peroxide (H 2 O 2 ) content significantly increases in the liver of the giant spiny frog (Quasipaa spinosa) under heat stresses (Liu et al., 2018). The antioxidant defense system includes enzymatic antioxidants, such as catalase (CAT), glutathione peroxidase (GPX), glutathione reductase (GR), superoxide dismutase (SOD), and some non-enzymatic antioxidants, such as glutathione (GSH; Cheng et al., 2015;Paital et al., 2016). Liu et al. (2018) found that the antioxidant enzyme system in the liver of giant spiny frogs (Q. spinosa) was enhanced to cope with oxidative stress when they were exposed at 30 • C for 48 h. Furthermore, the antioxidant system does not change in corn snake (P. guttatus) when acclimated to a simulated heat wave, although its oxidative damage level significantly decreases (Stahlschmidt et al., 2017). The telomere length of desert lizard (P. przewalskii) shortens while SOD activity did not change in a stimulated heat wave (Zhang et al., 2018).These research studies suggest that high temperature and heat wave may have a different impact on individuals.
Nanorana pleskei is endemic to the Qinghai Tibet Plateau and belongs to the Dicroglossinae family. It is distributed at 2,800 to 5,100 m above sea level and is considered to be one of the amphibians with highest altitude distribution in the world (Zhang et al., 2012). Niu et al. (2018) investigated oxidative stress and antioxidant defense in hibernating N. parkeri, and the results showed that the level of oxidative stress and oxidative damage was significantly increased, but that the total antioxidant capacity was significantly decreased compared with the summer frogs. Zhang et al. (2021) found that long-term cold acclimation would inhibit the antioxidant defense system of N. pleskei. However, there is little know about how N. pleskei responds to the warming environment at the organism's performance and physiological and molecular levels.
Herein, to investigate the influence of high temperatures and heat wave on the high-altitude N. pleskei in the background of global and plateau warming, we measured the thermal biology, locomotor performance, oxidative stress, and antioxidant defense system in liver and muscle of frogs acclimated to high temperatures and heat waves, and compared the results with that of the field individuals. Our results provide a reference for studies on the physiological and biochemical mechanisms of high-altitude amphibians adapted to high temperatures, and for survival prediction and protection of high-altitude amphibians.

Ethics Statement
This study was carried out under the approval of the Ethics Committee of Animal Experiments at Lanzhou University and on the basis of principles from the China Council on Animal Care. Every effort was made to minimize the numbers used and any suffering experienced by the animals in the experiment.

Animal Acclimation
Adult N. pleskei frogs (n = 84) were captured in Maqu County (33.77 • N, 101.74 • E, 3,490 m asl), Gansu province, China, at the end of July, 2020. They were transported to the research station of alpine meadow and wetland ecosystems of Lanzhou University in Maqu, Azi Branch Station (33.67 • N, 101.87 • E, 3,512 m asl). The frogs used in the experiment were captured on postbreeding period, and the sex ratio was 1:1 in the measurement of biochemical indexes. Acclimation to high temperatures and heat wave, and the measurement of thermal biology and locomotor performance were all performed in the research station 15 km from the sampling site. The frogs were randomly separated into four groups (n = 21 for each group). The thermal trials included three treatments consisting of two high-temperature treatments and one heat wave treatment. The frogs captured directly from the field were assigned to the control group. We set two high temperature groups in this study, ETem20 (20.5 ± 0.5 • C) and ETem25 (25 ± 0.5 • C), according to the preferred body temperature of field frogs (21 ± 0.5 • C) and the monthly maximum (22.5 ± 0.5 • C) in August in Maqu County from 2005 to 2017, respectively ( Figure 1A). To further simulate the temperature fluctuation experienced by frogs during the active period, the heat wave group (HWG) was set based on the monthly maximum (12.1 ± 0.5 • C) and minimum (1.2 ± 0.5 • C) average temperatures in August in Maqu County from 2005 to 2017 ( Figure 1A). The temperatures of the ETem20, ETem25, and heat wave groups were controlled 2 2 weeks by an automatic temperature control system, and three temperature data loggers (iButton DS1922G; CA, United States) were fixed in the bottom of each aquarium to record temperature changes in real-time ( Figure 1B). According to our records of temperature data loggers during the acclimations, the average temperature of the heat wave group was 19.3 • C, and the range of diurnal temperature fluctuation was 15-26.6 • C ( Figure 1B).
The frogs (n = 21 for each group) were placed in glass aquaria (48.7 cm × 35.5 cm × 28.5 cm) containing water and grass from the wild. The animals were fed with insects that were captured near the habitat of frogs, such as mosquitoes, screwworms, and locust larvae. The water and grass were replaced every 3 days. The photoperiod was 12 L: 12 D. Before the experiment, the cloacal temperature, snout-to-vent length (SVL), and body mass of all the frogs were measured. Body mass index (BMI, body mass × SVL −2 × 100%) was calculated before and after the acclimation treatments. Then, 13 frogs from each group were used to measure thermal biology and locomotor performance parameters, and 8 frogs from each group were sacrificed to collect the whole liver and hind leg muscle tissues. The hepatosomatic Changes in temperature in ETem20, ETem25, and HWG, and changes in indoor temperature and water temperature of Nanorana pleskei for 2 consecutive days. ETem20, ETem25, HWG, and water temperature represent the group that elevated temperature to 20 • C, the group that elevated temperature to 25 • C, water temperature at sampling site heat, and wave group, respectively.
index (HSI, liver mass × body mass −1 × 100%) and muscle index (MI, muscle mass × body mass −1 × 100%) were also calculated. All tissues were quick-frozen in liquid nitrogen immediately and stored at -80 • C until the measurement of biochemical parameters.

Thermal Biology
Preferred body temperature (T pref ) and ambient temperature (T amb ) were measured in a homemade cuboid-shaped glass tank (120 cm × 15 cm × 30 cm), which was covered up with a shading cloth in all the trials to create a completely dark environment (César et al., 2018). A layer of moist gauze was spread on the bottom of the tank to ensure an adequately humid environment in the trials. One side of the tank was warmed using a ceramic lamp placed inside, and the other side was cooled using ice packs placed outside. Three temperature data loggers were placed on both ends and the middle to record temperature changes in real time. Our results showed that the range of temperature gradient was 14.6 ± 0.5 to 39.1 ± 0.5 • C, and that the temperature in the middle was 20.7 ± 0.5 • C. Frogs (n = 13 for each group) were placed in the temperature gradient for about 1 h to adapt before each trial. After 1 h, the cloacal temperatures of the frogs were measured every 30 min using a handheld electronic thermometer (TM6801b, China). Each frog was tested three times, and the average was defined as T pref . The ambient temperature (T amb ) in the location of each frog was also measured three times using the same thermometer, and the average was defined as T amb .
The frogs used for measuring T pref and T amb were used further to measure thermal tolerances after resting for 24 h. The thermal tolerances of frogs comprise the non-lethal critical thermal maximum (CT max ) and minimum (CT min ), which were the cloacal temperatures at which frogs lose their ability of righting reflex at high and low temperatures (Lowe and Vance, 1955). The time interval between the measure of CT max and CT min was 48 h. Low-temperature thermostat bath (HMDC-2006, China) was used to control temperatures in the measurement of CT max and CT min . In the process of measuring CT max and CT min , water temperature was gradually increased or decreased at a rate of.33 • C min −1 . In the measurement of thermal tolerances, the frogs were checked every 5 min at the beginning, and then checked every 1 min when their movement slows down until individuals lose their ability of righting reflex, and then the cloacal temperature was recorded as the CT max and CT min . The measurement time of each group was not more than 30 min to avoid short-term adaptation.

Stamina and Burst Swimming Speed
After 48 h of measuring CT max and CT min , the stamina of frogs was measured in a 2.7-m homemade runway with a water temperature of about 16 • C. An individual was chased by a glass stick until it was exhausted and lost the ability of righting reflex, and then the total distance and time of jumping to exhaustion were recorded to characterize the stamina of frogs. The burst swimming speed of frogs was further measured after 48 h of rest. The frogs were placed in a custom-made cuboid-shaped glass tank (120 cm × 15 cm × 30 cm) containing 10-cm deep water collected from the field, and a ruler with an accuracy of.1 cm was fixed at the bottom of the tank to record the distance of swimming. A digital camera (Canon A610; Canon, Japan) was fixed above the tank to videotape the swimming movements of the frogs when they were being stimulated with a glass stick on the tail. Then, the videos were analyzed to measure the swimming distance during the 1 s using the Adobe Premiere Pro 2020 frame-by-frame. Each video was analyzed 10 times, and the maximum distance during the 1 s was regarded as burst swimming speed (cm/s).

Measurement of Oxidative Damage and Antioxidant Activity Indices
The liver and muscle tissues from the eight frogs in each group were accurately weighed, and then we immediately added a quantitative pre-cooled.8% saline solution in a 1:4 ratio in 1.5-ml centrifugal tubes to thoroughly grind, and then the supernatant fluid was collected by centrifuging the tissue fluid at 3,000 rpm and 4 • C for 10 min. Afterward, the supernatant fluid was diluted in proportion to 10, 1, 0.3, and.1% tissue solution with the pre-cooled.8% saline solution. Chemical assay kits were used to determine the total protein content, MDA, H 2 O 2 , GPX, SOD, CAT, GR, T-AOC, GSH, and GSSG (Nanjing Jiancheng Bioengineering Institute, Nanjing, China) according to the instructions. The GSH-eq was calculated as the sum of GSH and 2 × GSSG.

Statistical Analysis
All the data were represented as means ± standard error and analyzed on the SPSS 20.0 (IBM SPSS statistics 20.0). The data were checked by Kolmogorov-Smirnov test for normality. The statistical differences (P < 0.05) of each measured parameter among the high temperatures, heat wave, and field groups were analyzed by one-way ANOVA. The results of homogeneity variance were showed in the one-way ANOVA analysis, and LSD post hoc test for data in homogeneity and Tamhane's T2 test for data in heterogeneity were used in the post hoc multiple comparisons to detect the differences between every two groups (P < 0.05). Furthermore, a paired t-test was performed to compare the significant difference in measured parameters before and after acclimation of each group. All graphs were generated with Origin 2018 64Bit (OriginLab, American).

DISCUSSION
This study investigated the effect of high temperatures and heat wave on thermal biology, locomotion performance, oxidative stress, and antioxidant defense systems in the high-altitude frog N. pleskei. The results indicated that thermal biology and locomotion performance showed plasticity for N. pleskei acclimated for 2 weeks at high temperatures and heat wave. The level of oxidative stress, oxidative damage, and antioxidant activity were remarkably enhanced in frogs acclimated at high temperatures. However, heat wave nearly has no significant influence on the oxidative stress of N. pleskei and even decreased the oxidative damage and activities of mostly antioxidant enzymes.

Preferred Body Temperature and Thermal Tolerances
The increasing ambient temperature may reduce the activity time of poikilotherms, and further affect their fitness-related activities, such as foraging and mating, which will directly lead to the decline in population size and even lead to species extinction (Gvozdík, 2012). Although behavioral thermoregulation enables them to maintain their body temperature within an optimal range, it also comes at the cost of increased energy consumption and predation risk (Sartorius et al., 2002). T pref is closely related to the optimal temperature range of biochemical and physiological activities of organisms (Angilletta et al., 2010a). In this study, the T pref of N. pleskei was significantly increased after acclimation to high temperatures and heat wave, compared with field frogs. The significant effect of thermal acclimation on T pref has been confirmed in other ectotherms, such as terrestrial toad (R. icterica) (César et al., 2018), fresh-water prawn (Macrobrachium tenellum) (Hernandez-Sandoval et al., 2018), river shrimp (Cryphiops caementarius) (Avalos et al., 2020), and soft-shelled turtle (P. sinensis) (Wu et al., 2013). The great ability of thermal acclimation makes them remain active at high temperature, and expand the time window for activity, so as to partially or completely compensate for the impact of climate change (Gvozdík, 2012). CT max and CT min are usually defined as the upper limit and lower limit of individual survival, and animals will experience movement disorder, spasm, and loss of righting reflex or balance when the ambient temperature reaches or exceeds the hightemperature and low-temperature limits (Lowe and Vance, 1955). We found that both the CT max (31.62 • C) and CT min (3.17 • C) of unacclimated frog N. pleskei were higher than the annual maximum (24.7 • C) and minimum(−24.4 • C) temperatures in their distribution range. Ectothermic animals have thermal tolerance after a short period of sublethal temperature, which is an acute response, enabling them to cope with short-term fluctuations in ambient temperature (Hong et al., 2009). That was aligned with the results, 2 weeks of acclimation to high temperatures significantly increased the temperature tolerance range of N. pleskei (CT max was significantly increased, and CT min was significantly decreased). The plasticity of thermal tolerance has also been demonstrated in other exothermic animals, such as Galaxias zebratus (Olsen et al., 2021), Argyrosomus regius (Kır et al., 2017), and Eleutherodactylus coqui (Haggerty, 2016).
The T pref , CT max , and CT min of frogs acclimated in heat wave also showed plasticity. The T amb and CT max increased and CT min decreased significantly, although there was no significant difference in the T pref compared with field frogs. Brusch et al. (2016) investigated the effects of microhabitat type, species, and SVL on thermal tolerances of 10 lizards from Costa Rica and revealed that microhabitat type is an important factor to predict the maximum tolerated temperature regardless of size or evolutionary relationship. The climate conditions in Tibetan plateau are complex and changeable, such as windy and dry in spring and winter, rainy in summer and autumn, and great circadian temperature difference (Zhang et al., 2006). In addition, Zhao et al. (2014) found that the summer days and daily maximum and minimum air temperature have all risen, while the freezing and frost days have all declined in the QTP from 1963 to 2012. The warming trend of lowest temperature was more striking than that of highest temperature. Chen et al. (2020) predicted that the number of heat wave days in the Qinghai Tibet Plateau will significantly increase in the future. Therefore, we suggested that the plasticity of thermal biology in frogs acclimated to high temperatures and heat wave may be closely related to the climate of microhabitat of N. pleskei on the plateau, and that it would show great adaptability to the variability of habitat climate. The thermal tolerance of highaltitude N. pleskei may be due to their adaptation mechanism to the changing thermal environment, enabling them to survive in the short-term exposure at high temperatures. A similar phenomenon has also been found in the anuran species that have transformed from aquatic to terrestrial or semi-terrestrial during ontogenesis, and their increased thermal tolerance is also important for their survival, as they may be exposed to a wide range of diurnal temperature fluctuations (Sherman, 1980;Rome et al., 1992).

Stamina and Burst Swimming Speed
As found in most frogs and toads (Knowles and Weigl, 1990;Wilson and Franklin, 1999;Wilson et al., 2000;Padilla et al., 2019), the locomotor performance of highland N. pleskei was also affected by the acclimation temperature. Compared with the field group, the total time of jumping to exhaustion significantly decreased while the burst swimming speed of the frogs was significantly increased after the acclimation to high temperatures and heat wave. The locomotor performance of the animals, especially the burst performance, is important for avoiding predators and catching prey. However, stamina may be more important in finding and copulating with mates. Probably there is a trade-off between these two locomotor properties, as they depend on different underlying physiological mechanisms (Padilla et al., 2019). Our results demonstrated that N. pleskei was able to jump faster until exhaustion and swim faster after the acclimation to high temperature and heat wave, and that this may be more likely to contribute to their foraging and avoiding predators successfully to maintain normal physical condition (Preest and Pough, 2003). Accordingly, we speculated that the strengthening of burst ability and the attenuation of endurance may be beneficial to the survival of frogs in a short time, but with the continuous rising of global temperatures, further research is needed to assess the impact of extreme high temperature on the survival of plateau amphibians.

Oxidative Stress and Oxidative Damage
The oxygen consumption and metabolism of animals increase with rise in temperature within a certain range, resulting in increased generation of reactive oxygen species (ROS; Paital et al., 2016). The overproduction of ROS will damage DNA, proteins, lipids, and other biomolecules (Cheng et al., 2015). The thermal acclimation significantly increased the level of H 2 O 2 (Figure 4B), and significantly decreased GSH content and GSH/GSSG ratio in the liver and muscle in N. pleskei (Table 3). This suggests that high-temperature-acclimated N. pleskei is in a state of oxidative stress compared to the field frogs. The body mass index (BMI) and muscle index (MI) values significantly decreased in the frogs after acclimation in high temperatures compared with those before acclimation, and the hepatosomatic index (HSI) was lowest in ETem25 among all the groups (Table 1). These may be the result of the energy budget among antioxidation, growth, and development for N. pleskei in response to high temperatures (Ferral et al., 2020). Niu et al. (2018) and Zhang et al. (2021) have demonstrated that the content of GSH-eq and GSSG in N. parkeri was markedly decreased in winter and long-term cold acclimation. On the contrary, the content of GSH-eq and GSSG obviously increased in the liver and muscle in N. pleskei acclimated at high temperatures and heat wave. This inconsistency may be related to food intake. The frogs were unfed during winter and cold acclimation, but food was plenty during thermal acclimation in our experiment. Therefore, it is possible that the de novo synthesis of glutathione increased in the frogs during our thermal acclimation, so the GSH-eq content in the liver and muscle significantly increased (Zhang et al., 2012;Niu et al., 2018). GSH is one of the important nonenzymatic antioxidants in cells to inhibit the production of ROS (Meister and Anderson, 1983). In our study, GSH is greatly oxidized to GSSG by the catalysis of GPX to inhibit the further production of ROS, because the activity of GPX was significantly enhanced in N. pleskei acclimated at high temperatures and heat wave ( Figure 5C). This maintains the redox status in hepatocytes and muscle cells of frogs under thermal acclimation, which may be the main reason for the increase in GSSG content and decrease in GSH/GSSG ratio in the liver and muscle.
The acclimation to heat wave did not significantly influence the level of H 2 O 2 in the liver and muscle, but it remarkably reduced the lipid oxidative damage in the liver (Figure 4). This may be due to the increase in CAT activity (Figure 5A), because it decomposes H 2 O 2 into H 2 O and oxygen (Lledías et al., 1998). The GSH content in the liver and muscle, and GSH/GSSG ratio in the muscle significantly decreased after the acclimation to heat wave. These results correspond with the study on corn snake (P. guttatus), in which the acclimation to heat wave reduced its oxidative stress (Stahlschmidt et al., 2017). However, the telomere length and overwinter survival of desert toad-headed agama (Phrynocephalus przewalskii) were shortened after exposure to a stimulated heat wave (Zhang et al., 2018). In Asian yellow pond turtle (Mauremys mutica), heat wave does not significantly affect the level of MDA in the liver and muscle (Li et al., 2021). Therefore, heat wave has different physiological effects on different species, positive or negative. Our study indicated that simulated heat wave does not induce oxidative stress, and even significantly decreased oxidative damage in frogs. In addition, the BMI did not significantly change in the frogs acclimated in heat wave compared with that before the acclimation to heat wave. Therefore, we suggested that heat wave may have positive effects on N. pleskei. The possible reason for the non-significant or positive effects of heat wave is that the thermal regimes of the heat wave group may be too mild for frogs and they can recover from heat stress in the lowtemperature period.

Antioxidant Defense System
Under normal conditions, the production and elimination of ROS maintain a dynamic balance in cell metabolism.
Thermal stress could impair the balance between oxidant and antioxidant defense mechanisms, causing oxidative stress and oxidative damage in tissues of organisms (Maud et al., 2006;Tu et al., 2012). As we have mentioned above, higher temperature acclimation induced oxidative stress and oxidative damage in N. pleskei (Figure 4). To offset the adverse effects of peroxides, the organism usually activates the antioxidant enzymes and some non-enzymatic antioxidants to resist oxidative stress and prevent or repair oxidative damage (Tu et al., 2012;Duan et al., 2016). In this study, compared to the field group, the activities of GPX, GR, CAT, SOD, and T-AOC in the liver were significantly increased in ETem20 (Figure 5). In addition, the activities of GPX and SOD and T-AOC in the liver, and GPX and GR activity in the muscle were significantly increased in ETem25 (Figure 5). These results indicated that the overall antioxidant capacity of N. pleskei is enhanced after high-temperature acclimation. Similar inference was found in giant spiny frog (Quasipaa spinosa) (Liu et al., 2018), puffer fish (Takifugu obscurus) (Cheng et al., 2015), and ark shell (Scapharca broughtonii) (An and Choi, 2010).
In terms of the effect of heat wave on antioxidant defense, our results showed that the activities of GR, GPX, and SOD decreased in the liver and muscle, but that the activity of CAT in the liver increased in high-altitude frog, N. pleskei (Figure 5). This may be the reason for the significant decrease in lipid peroxide level. Similar results were found in P. guttatus in the acclimation to simulated heat wave, in which the antioxidative capacity did not change but the oxidative damage level obviously decreased (Stahlschmidt et al., 2017). On the contrary, although the activity of T-SOD and CAT, and T-AOC are not affected by heat wave, the GPX activity in the muscle in M. mutica significantly increased (Li et al., 2021). The activities of antioxidant enzymes were attenuated for N. pleskei in the acclimation to simulated heat wave, but its oxidative damage level dramatically decreased. In the process of long-term adaptation to the plateau environment, N. pleskei may have evolved a unique mechanism to deal with temperature fluctuations. Therefore, their oxidative damage and antioxidant enzyme activity did not significantly change and even attenuated with simulated heat wave treatment.

CONCLUSION
In conclusion, thermal preference, thermal tolerance, stamina, and burst performance showed plasticity under thermal acclimation and simulated heat wave. In addition, thermal acclimation disturbs the redox state, which is manifested as oxidative stress and oxidative damage in the liver and muscle of the high-altitude frog N. pleskei. The antioxidant capacity was enhanced significantly after high-temperatures acclimation. Unexpectedly, the heat wave did not cause oxidative stress and oxidative damage, and even decreased the activity of antioxidant enzymes. In summary, the impact of heat wave on amphibians is not as great as we expected, while the short-term extreme high temperatures have a greater impact on them. These results can provide a reference for studies on the adaptation of other high-altitude amphibians and protection of the QTP ecosystem.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.

ETHICS STATEMENT
The animal study was reviewed and approved by the Ethics Committee of Animal Experiments at Lanzhou University.