# Rapid Adjustments in Thermal Tolerance and the Metabolome to Daily Environmental Changes – A Field Study on the Arctic Seed Bug Nysius groenlandicus

Natasja Krog Noer, Mathias Hamann Sørensen, H Colinet, D Renault, Simon Bahrndorff, Torsten Nygaard Kristensen

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Natasja Krog Noer, Mathias Hamann Sørensen, H Colinet, D Renault, Simon Bahrndorff, et al.. Rapid Adjustments in Thermal Tolerance and the Metabolome to Daily Environmental Changes – A Field Study on the Arctic Seed Bug Nysius groenlandicus. Frontiers in Physiology, 2022, 13, pp.818485. 10.3389/fphys.2022.818485. [hal-03585727](https://hal-03585727v1)

HAL Id: hal-03585727

https://univ-rennes.hal.science/hal-03585727v1

Submitted on 6 May 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_3.png" /></figure>

Check for
updates

# Rapid Adjustments in Thermal Tolerance and the Metabolome to Daily Environmental Changes – A Field Study on the Arctic Seed Bug *Nysius groenlandicus*

Natasja Krog Noer¹*, Mathias Hamann Sørensen¹, Hervé Colinet², David Renault²,³, Simon Bahrndorff¹† and Torsten Nygaard Kristensen¹†

¹ Department of Chemistry and Bioscience, Aalborg University, Aalborg, Denmark, ² UMR 6553, CNRS, Ecosystèmes, Biodiversité, Évolution, University of Rennes 1, Rennes, France, ³ Institut Universitaire de France, Paris, France

Laboratory investigations on terrestrial model-species, typically of temperate origin, have demonstrated that terrestrial ectotherms can cope with daily temperature variations through rapid hardening responses. However, few studies have investigated this ability and its physiological basis in the field. Especially in polar regions, where the temporal and spatial temperature variations can be extreme, are hardening responses expected to be important. Here, we examined diurnal adjustments in heat and cold tolerance in the Greenlandic seed bug *Nysius groenlandicus* by collecting individuals for thermal assessment at different time points within and across days. We found a significant correlation between observed heat or cold tolerance and the ambient microhabitat temperatures at the time of capture, indicating that *N. groenlandicus* continuously and within short time-windows respond physiologically to thermal changes and/or other environmental variables in their microhabitats. Secondly, we assessed underlying metabolomic fingerprints using GC-MS metabolomics in a subset of individuals collected during days with either low or high temperature variation. Concentrations of metabolites, including sugars, polyols, and free amino acids varied significantly with time of collection. For instance, we detected elevated sugar levels in animals caught at the lowest daily field temperatures. Polyol concentrations were lower in individuals collected in the morning and evening and higher at midday and afternoon, possibly reflecting changes in temperature. Additionally, changes in concentrations of metabolites associated with energetic metabolism were observed across collection times. Our findings suggest that in these extreme polar environments hardening responses are marked and likely play a crucial role for coping with microhabitat temperature variation on a daily scale, and that metabolite levels are actively altered on a daily basis.

**Keywords:** arctic, climate change, diurnal environmental variation, GC-MS metabolomics, insects, phenotypic plasticity, temperature variation, thermal tolerance

**OPEN ACCESS**

**Edited by:**

Frank Seebacher,
The University of Sydney, Australia

**Reviewed by:**

Vladimír Koštál,
Centre for Biology, Academy
of Sciences of the Czech Republic
(ASCR), Czechia
Folco Giomi,
Independent Researcher, Padua, Italy

*Correspondence:

Natasja Krog Noer
nkn@bio.aau.dk

† These authors share last authorship

**Specialty section:**

This article was submitted to
Invertebrate Physiology,
a section of the journal
Frontiers in Physiology

**Received:** 22 November 2021

**Accepted:** 03 January 2022

**Published:** 16 February 2022

**Citation:**

Noer NK, Sørensen MH, Colinet H, Renault D, Bahrndorff S and Kristensen TN (2022) *Rapid Adjustments in Thermal Tolerance and the Metabolome to Daily Environmental Changes – A Field Study on the Arctic Seed Bug Nysius groenlandicus*. *Front. Physiol.* 13:818485.
doi: 10.3389/fphys.2022.818485
# INTRODUCTION

Terrestrial ectotherms are subject to large spatial and temporal variability in their thermal environment (Kearney and Porter, 2009). In terrestrial ecosystems, daily changes in temperature can be substantial, and vary greatly with microhabitat characteristics such as topography and orientation, vegetation cover, shading and more (Sears et al., 2019; Kearney et al., 2020; Lembrechts and Lenoir, 2020). Some of the most extreme environments are found in the polar regions where the winters are long and cold, and the summers short and periodically hot (Bahrndorff et al., 2021b). During the arctic summer, daily temperatures can vary by >30°C and reach subzero temperatures at night (Convey et al., 2018; Davey et al., 2021). Organisms, including insects, living in these environments must therefore be able to survive and reproduce over a wide range of temperatures (Deere et al., 2006; Bahrndorff et al., 2021a). This can be achieved by evolutionary adaptation to the local thermal conditions across generations, or by fast adjustments of the physiology within the lifetime of an organism *via* phenotypic plasticity (Scheiner, 1993; Fusco and Minelli, 2010; Kristensen et al., 2020). Evolutionary adaptation to changing and periodically stressful temperatures can be slow, and are sometimes constrained by genetic trade-offs or lack of adaptive genetic variation (Araújo et al., 2013; Hoffmann et al., 2013). Conversely, rapid plastic adaptive changes can rescue individuals exposed to biotic and abiotic challenges at a shorter timescale, including daily environmental fluctuations (Colinet and Hoffmann, 2012; Noer et al., 2022). Plastic changes might therefore be particularly relevant for arctic species exposed to unpredictable and rapid changes in the environment.

Organisms can respond plastically to short-term exposure to sub-optimal temperatures through hardening or by acclimation at longer term exposures (Colinet and Hoffmann, 2012; Schou et al., 2017). Hardening responses to extreme or acute temperatures are thought to counter rapid thermal stress, such as daily temperature extremes and stochastic events (Koveos, 2001; Hoffmann et al., 2003; Kelty, 2007; Overgaard and Sørensen, 2008; Jensen et al., 2019). The other form of more gradual acclimation includes seasonal acquisition of cold or heat tolerance that is induced by changes in temperature and photoperiod interacting with other abiotic factors (reviewed by Chown and Terblanche, 2006; Teets and Denlinger, 2013). There are several published examples of cold acclimation and rapid cold hardening in arctic arthropods (e.g., Bahrndorff et al., 2007; Everatt et al., 2013), but very few studies have investigated physiological acclimation of polar terrestrial arthropods to high temperatures (Sørensen et al., 2019; Bahrndorff et al., 2021b). Traditionally, thermal plasticity of insects has been investigated using model-organisms kept and hardened/acclimated to constant controlled temperatures in the laboratory (Angilletta, 2009; Colinet et al., 2015; Javal et al., 2016; Ketola and Kristensen, 2017). However, recent work on the impacts of temperature variability on thermal tolerance have emphasized that thermal performance based on constant temperatures do not always accurately predict performance under variable conditions in the laboratory (reviewed by Colinet et al., 2015; Vázquez et al., 2017), nor in the field (see e.g.,

Kingsolver and Nagle, 2007; Loeschcke and Hoffmann, 2007; Kristensen et al., 2008; Ketola and Kristensen, 2017). This potential mismatch in the conclusions arising from investigations based on constant *versus* fluctuating temperatures partly results from the non-linear impact of temperatures on thermal performance (Jensen's inequality) (Ruel and Ayres, 1999; Colinet et al., 2015), time-by-temperature interactions (Foray et al., 2013; Kingsolver et al., 2015), and methodology (Chown et al., 2009; Mitchell and Hoffmann, 2010; Bahrndorff et al., 2016). Based on such results, the potential for transferring the knowledge obtained from the laboratory to field conditions, and thus forecast reliable predictions of the effects of climate change on the responses and geographic distribution of insects, is being increasingly questioned (Fischer et al., 2011; Kingsolver et al., 2015; Kinzner et al., 2019; Taylor et al., 2021).

The physiological and molecular mechanisms enabling arthropods to tolerate temperature stress has previously focused on controlled laboratory studies (but see Tomanek and Somero, 1999, 2002; Buckley et al., 2001; Gracey et al., 2008; Kristensen et al., 2012; Vasquez et al., 2019). Studies on temperate and polar species suggest (causation is typically lacking in such studies) that metabolites such as sugars, free amino acids and polyols can be associated with changes in cold tolerance measures (Zachariassen, 1985; Fields et al., 1998; Somme, 1999; Holmstrup et al., 2002; Michaud and Denlinger, 2007, 2010; Overgaard et al., 2007, 2014). Changes in polyols on the other hand have been associated with changes in heat tolerance (Hendrix and Salvucci, 1998; Wolfe et al., 1998; Salvucci et al., 2000; Benoit et al., 2009). However, few have attempted to describe how these metabolites are affected by dynamic and fluctuating temperatures as encountered in nature (but see Kristensen et al., 2012; Noer et al., 2020; Sheldon et al., 2020).

In this study, we examined the effects of daily variation in the microhabitat temperatures on plastic adjustments of heat and cold tolerance of the Greenlandic seed bug *Nysius groenlandicus* (Zetterstedt) during summer in Southern Greenland. *Nysius groenlandicus* is a univoltine species, and widespread and abundant in Arctic and sub-Arctic regions. Previous work on the species have revealed that it can rapidly increase heat tolerance when exposed to high and stressful temperatures under laboratory conditions (Sørensen et al., 2019), and thus *N. groenlandicus* represents a valuable polar insect model for field-based description of daily changes in individuals' thermal tolerance. Further, we examined the metabolic fingerprints within days with high and low temperature variation using a quantitative targeted gas chromatography-mass spectrometry (GC-MS) approach. We hypothesized that the ability to tolerate low and high temperatures is constantly fine-tuned to respond to temporally fluctuating temperatures in field-collected individuals of *N. groenlandicus* as an adaptation to the highly variable environmental conditions within and across days. Thus, we expected specimens collected in early morning and late evening to have the highest cold tolerance, while those collected at midday exhibiting the highest heat tolerance. Finally, we expected to see daily changes in metabolites known to improve cold tolerance (sugars and amino acids) in individuals collected in early morning and late evening, while metabolites enhancing
heat tolerance (polyols) would show higher concentrations in individuals sampled during the warm periods of the day.

# MATERIALS AND METHODS

## Field Collection and Experimental Design

The field work was conducted in July–August 2018 in Narsarsuaq (Southern Greenland, 61.160°N, 45.424°W). This region is characterized by cool temperatures, long winters and short and thermally variable summers (Böcher and Nachman, 2001; Bahrndorff et al., 2021b). The study site was a heath-like, grass-covered area, where adult *N. groenlandicus* were collected from the grasses using a sweep net (Figure 1). All individuals used in the experiment were caught in a 50 m × 50 m area less than 100 m from the laboratory facilities.

Adult individuals of *N. groenlandicus* were collected at four time points (08:00 am, 12:00 pm, 04:00 pm, and 08:00 pm) during each of 5 days (depending on weather conditions and the abundance of *N. groenlandicus*; see Supplementary Table 1 for

exact collection times and dates), thereafter referred to as day 1 to day 5. At the time of field-collection, each individual was placed in a 4 mL screw-cap glass vial (45 mm × 14.7 mm) and placed in the shade on the ground to prevent abrupt changes in the thermal conditions. The sex of individuals was then assessed by eye, and the vials transferred to the laboratory within 30–45 min of collection. Immediately after returning, the heat and cold tolerances were scored using 20 females and 20 males for each assay (see next section).

Additional individuals were collected for subsequent metabolomics analysis at the same four sampling times at two dates (22/08/18 and 27/08/18, representing day 4 and 5, respectively, see Supplementary Table 1), representing days with either high or low observed temperature variation. At each collection time, eight samples of five females (only females were used for the metabolomics studies) were collected, transferred directly into ice cold RNAlater, and stored at –20°C for approximately 1 week. Then the samples were transferred to our laboratory in Denmark where they were stored at –80°C and later used for metabolomic fingerprinting. Collection of

<figure><img src="image_4.png" /><figcaption>**FIGURE 1** | Microhabitat temperature (°C) averaged in 1-h bins across the five experimental days. The temperatures were measured at 15 cm above the soil surface in the shade every 5 min. Samples for GC-MS analysis were collected on day 4–5 (black lines). Shaded columns indicate the timespans where specimen are collected.</figcaption></figure>
samples in RNAlater is amenable for downstream metabolomics
analysis if stored at subzero temperatures (van Eijsden et al.,
2013; Harris, 2018).

The air temperature at the collection site was continuously recorded with 5-min intervals in the shade using Easylog USB data loggers (LASCAR Electronics, EL-USB-2+). The temperature was measured in the shade to avoid warm temperature-spikes in the measurements caused by direct solar radiation (Maclean et al., 2021) and the loggers were placed 15 cm above the soil surface to reflect the thermal environment at the top of the grasses where *N. groenlandicus* was most abundant, and was caught with the sweeping net. Based on these recordings, the mean temperature was calculated for the 1-h timespan prior to testing thermal tolerances. The mean temperature immediately prior to thermal assessment has been shown to be highly correlated with the heat tolerance in a range of insect species collected in the field (Noer et al., 2022).

# Heat and Cold Tolerance Assays

## Heat Knockdown Time

Heat tolerance was measured as heat knockdown time (HKDT), i.e., a measure of the time before individuals go into a heat-induced coma/die at a high stressful temperature (Terblanche et al., 2007; Bak et al., 2020). For *N. groenlandicus* HKDT constitutes a direct measure of heat tolerance as we register the time individuals can withstand before they die. The vials containing field-collected individuals were mounted to a rack and submerged into a temperature-controlled water bath (PolyScience MX Immersion Circulator: MX-CA12E) maintained at 48°C. The individuals were then observed and stimulated with flashes of light and gentle taps on the vial caps with a metal rod. The time until movement ceased was noted for each individual. The chosen HKDT temperature was based on experiences from previous work on the species (Sørensen et al., 2019), and on unpublished preliminary results showing that *N. groenlandicus* individuals went into coma within 20-40 min at 48°C. Earlier results suggest a fast heat hardening response for this species and since we wanted a measure of their acute heat tolerance, and aimed for reducing hardening responses induced while testing, this HKDT was found relevant.

## Chill Coma Recovery

Cold tolerance was measured as the temperature at which the bugs regained the ability to move after being knocked down by exposure to a low temperature ($T_{recovery}$) following a modified procedure of the method described in Overgaard et al. (2011). Thus, we used a proxy rather than direct measures of cold tolerance. The glass vials containing the individuals were mounted to a rack and submerged into a glycol-water solution that was kept at -3°C using a thermostat (LAUDA Proline Edition X RP 1845-C, LAUDA DR. R. WOBSER GMBH & CO., KG, Germany). This temperature was based on the lower critical temperature (CTmin) which induces a cold coma for the species (Bahrndorff et al., 2021b). Immediately after submersion, the temperature of the solution was increased at a rate of 0.2°C/min. Pilot studies showed that this temperature (-3°C) was sufficient to induce chill coma within a few minutes

with full survival upon returning to room temperature (data not presented). Following the HKDT procedure, the individuals were observed and provoked using light flashes and gentle taps, and the temperature at which individuals first moved any body part was noted as their chill coma recovery temperature ($T_{recovery}$). A low $T_{recovery}$ is interpreted as high cold tolerance.

# Metabolomic Fingerprinting

We adapted the methods detailed in Thiébaut et al. (2020) for detecting and quantifying the metabolite content from whole-body extracts of female *N. groenlandicus*. For each field sampling time, eight replicates were used, each consisting of five pooled females to obtain sufficient biomass (~4 mg dry mass per sample). Each sample was vacuum dried (Speed Vac Concentrator, miVac; Genevac Ltd., Ipswich, England) and weighed (Mettler Toledo UMX2, accurate to 0.001 mg) before extractions so that metabolite concentrations could be reported according to dry mass. The samples were first homogenized for 90 s with two tungsten beads in 450 µL of a solution of ice cold methanol-chloroform (ratio 2:1, v:v) using a bead beater (Qiagen MM301; Retsch GmbH, Haan, Germany) set at 25 Hz. To separate the homogenate in two distinct phases (lipid-rich containing phase, and aqueous phase containing metabolites), 300 µL of cold ultrapure water was added to each tube before centrifugation (Sigma 2-16K, Sigma GmbH, Harz, Germany) for 10 min at 4000 g at 4°C. The supernatants containing metabolites were transferred to new tubes and stored at -80°C until analysis by GC-MS. Prior to the analysis, 120 µL of the metabolite extract was transferred to glass vials and vacuum dried.

The derivatization of the samples (dry residues) was automatized with a CTC CombiPAL autosampler (CTC Analytics AG, Zwingen, Switzerland). Dried samples were re-suspended in 30 µL of 25 mg mL⁻¹ methoxyamine-hydrochloride (CAS Number: 593-56-6, SIGMA-ALDRICH, St. Louis, MO, United States) in pyridine prior to incubation under orbital shaking at 40°C for 60 min. Following incubation, 30 µL of N-methyl-bis(trifluoroacetamine) (BSTFA, CAS Number: 685-27-8) was added, and derivatization was conducted at 40°C for 60 min under agitation. An Agilent 7890B gas chromatograph coupled to a 5977B mass spectrometer was used for the separation and detection of the metabolites. For each sample, 1 µL was injected (Injector temperature: 250°C; split mode with a split ratio of 25:1); the temperature of the oven increased from 70 to 170°C at a rate of 5°C/min, from 170 to 280°C at 7°C/min, and from 280 to 320°C at 15°C/min, and then remained at 320°C for 4 min. We used a 30-m fused silica column (HP5 MS 30 m, I.D. 0.25 mm, thickness 0.25 µm, 5% Diphenyl/95% Dimethylpolysiloxane, Agilent Technologies), and the gas carrier (Helium) had a flow of 1 mL per min. The temperatures of the transfer line and ion source were 280 and 230°C, respectively. Metabolite fragmentation and ionization were carried out by electronic impact (electron energy: 70 eV) and detected with the full scan mode. The detected peaks were identified and annotated with MassHunter (Agilent). Most detected metabolites were identified, and calibration curves were used with pure compounds for calculating the concentration of each metabolite.
# Statistical Analysis

Differences in thermal tolerance with sampling time were examined using two-way ANOVAs for male and female *N. groenlandicus* separately with “time to knockdown” or “chill coma recovery temperature” as dependent variables, and “time of day” and “day” as independent variables. The tests were run on individuals. To ensure normal distribution, the data were transformed using rank inverse transformation and the models were run on both transformed and non-transformed data for validation. In addition, we used Pearson's correlations to examine if microhabitat temperature affected thermal tolerances rather than sampling time alone. The Pearson's correlations were calculated using the mean field temperatures observed in the 1 h preceding each collection round and the mean HKDT and T<sub>recovery</sub> for each assay and sex (n = 20).

We used one-way ANOVAs to examine if individual metabolite concentrations within each day differed according to collection time. The concentrations of all quantified metabolites were then scaled and mean centered. For each sampling day, between-class PCA (R-package “ade4”), were run to explore the daily temporal structure of the metabolomic profiles. Monte Carlo tests (1000 permutations) were used to examine the significance of differences in metabolite profiles among classes of individuals from the four sampling times. Further, the metabolites that contributed the most to separation of groups were identified and ranked by their correlation to the principal components that described most of the inertia in the data. All analyses were carried out using the software “R” (R Core Team, 2020). Raw temperature and GC-MS files used for the analyses are available in the **Supplementary Material**.

# RESULTS

## Microhabitat Temperatures

Observed microhabitat temperatures of the different sampling periods are shown in **Figure 1**. The largest daily amplitude recorded was 28.8°C on day 1, while the lowest of 11.6°C was on day 4. The warmest average temperature in the 1-h timespan prior to thermal tolerance tests was 26.3°C (see **Supplementary Table 1** for exact temperatures) and was recorded at midday on day 2; the coldest average temperature prior to tests was 1.0°C and was recorded in the morning of day 5.

The microhabitat temperatures varied markedly between the two experimental days where samples were collected for metabolomic profiling (day 4 and 5; **Figure 1**). Within day 4, the temperatures were relatively constant with only 4.9°C difference in the average temperatures at the four daily collection times. Day 4 was also characterized by the lowest recorded daily amplitude. Day 5 was characterized by a high temperature variation with an amplitude of 25°C; the average temperatures at the four collection times differed by 21°C (**Figure 1** and **Supplementary Table 1**).

## Thermal Tolerance

Thermal tolerances varied significantly within and between days for both male and female individuals (**Supplementary Figure 2** and ANOVA **Supplementary Table 2**). The average HKDT at

each collection time and day was correlated with the average field temperature observed 1 h prior to heat tolerance assessment. The correlation was significant for females but not for males (**Figure 2A**). The regression slope, representing the change in HKDT per °C change in field temperature, was 0.36 min/°C for females and 0.18 min/°C for males. In addition, HKDT was overall higher for females than for males at any given field temperature with a difference of 5 min at the intercept. The maximal difference in mean HKDT measured across all temperatures was 23 min for females and 12 min for males. The relation between mean T<sub>recovery</sub> and the average microhabitat temperature preceding cold assays revealed a significant positive relationship for both females and males (**Figure 2B**). The increase in recovery temperature per °C change in field temperature were 0.09 and 0.07°C for females and males, respectively. The largest difference in T<sub>recovery</sub> found across all days was 3.3°C for females and 3.4°C for males.

For the days where females were collected for metabolomic fingerprinting, the HKDT and T<sub>recovery</sub> varied markedly across the thermally variable day (day 5), and less so on the thermally stable day (day 4) (**Supplementary Figure 2**). Thus, the correlation between the field temperature and heat tolerance of the females collected only at the time points used for metabolomic profiling was strong and highly significant (R = 0.94, *p* = <0.001). The relationship between field temperature and T<sub>recovery</sub> of females collected on day 4 and 5 was also positive and directional, however, not significant (R = 0.52, *p* = 0.19).

## Metabolomic Profiling

### Description of the Quantified Metabolites in the Samples

The full scan monitoring allowed us to identify and quantify 33 metabolites with the quadrupole GC-MS platform (**Supplementary Table 3**; and see raw metabolite concentrations in **Supplementary Table 4**). We quantified 13 free amino acids, five sugars, six polyols, six metabolic intermediates, and three other metabolites. Xylitol and ethanolamine concentrations were below the quantification limit (i.e., below the lowest concentration of the calibration curve) in many samples and were therefore excluded from the quantitative profiling. The most abundant metabolites across all sampling times and days were phosphoric acid, proline, glutamate, and tyrosine.

### Effects of Collection Time on Metabolic Fingerprints

Individual metabolite concentrations for each sampling time are displayed in **Figure 3**. We observed that a large number of metabolites varied in concentrations across collection times on day 5. There was a significant effect of “time of day” on the levels of eight metabolites (**Supplementary Table 5**). Conversely, the metabolite concentrations varied less on collection day 4, though there was a significant effect of “time of day” on concentration for seven of these (**Supplementary Table 5**). To examine whether the concentration levels varied similarly across collection times for the two collection days, we ran 2-way ANOVAs for each metabolite. Six metabolites varied distinctively across the four collection times for the 2 days, and these are depicted as
<figure><img src="image_5.png" /></figure>

**FIGURE 2** | Scatter plot showing the mean (A) heat knockdown time in minutes, and (B) chill coma recovery temperature ($T_{recovery}$) of females (circles) and males (triangles), as a function of the average microhabitat temperature 1 h prior to assay start. At eight sampling times, females were collected for metabolomic fingerprints simultaneously with individuals used for the thermal assays (blue circles). Black points are data points with no metabolomic fingerprints associated. Solid (female) and dotted (male) lines represent regressions between field temperature and thermal tolerances. Bars are standard errors of the mean.

the interaction between collection day and time on **Figure 3**. Finally, the sugars glucose, fructose, and galactose and the polyol glycerol occurred in larger concentrations in individuals collected on day 5.

Between-class PCAs were run separately on the metabolite concentrations from individuals sampled on day 4 (**Figure 4A**) and day 5 (**Figure 4B**). PC1 and PC2 cumulated 80 and 83% inertia on day 4 and 5, respectively; thus, on both days, the first two PCs explained most of the between-class variation. All classes showed a clear-cut separation meaning that metabolomic fingerprints differed among the different sampling times. The separation appeared stronger on day 5 (**Figure 4B**) than on day 4 (**Figure 4A**). Monte-Carlo randomization tests confirmed differences in metabolomic profiles among classes on day 5 ($p < 0.001$) and day 4 ($p = 0.026$) and the significance levels underpins that metabolomic profiles differed more markedly between sampling times on the thermally variable day 5 compared to on the less thermal variable day 4, as evidenced by much lower ellipses overlap.

## Effects of Low Daily Thermal Fluctuations on Metabolic Fingerprints

On day 4, which was characterized by low temperature variation, the metabolomic profiles of individuals collected in the morning and evening differed from the profiles from midday and afternoon along PC1 (**Figure 4A**). Further, the metabolic fingerprints of individuals collected during midday and afternoon separated along PC2.

The metabolites most positively correlated to PC1 were the sugars galactose and glucose (**Figure 4A**). These two sugars

were more abundant in individuals sampled in the morning and evening compared to those sampled at midday and afternoon. Lactic acid and citrulline were the metabolites that were most negatively correlated to PC1, and they were more abundant in individuals sampled at midday and afternoon. The polyol glycerol-3-phosphate and succinic acid had the highest positive associations with PC2 (**Supplementary Figure 3**) and were characterized by higher concentrations in individuals sampled in the morning and afternoon. Only one metabolite, glucose-6-phosphate, were negatively associated with PC2 and thus more abundant in the morning.

## Effects of High Daily Thermal Fluctuations on Metabolic Fingerprints

On day 5, metabolomic profiles of individuals collected in the morning separated strongly from individuals sampled in the afternoon (**Figure 4B**). These two collection times were also the two thermal extremes of the day (i.e., the lowest and highest temperatures of the day, see **Figure 1** and **Supplementary Table 1**); this observation suggests that PC1 explained metabolic changes correlated to high daily thermal variation. Further, the metabolic fingerprints of the individuals collected at these time points separated from the fingerprints of individuals collected at midday and in the evening along PC2.

The metabolites that were positively associated with PC1 on day 5 were phenylalanine, glutamic acid, citric acid, and to a lesser extent other amino acids, including valine, leucine, serine, and alanine (**Figure 4B**). These were more abundant in individuals sampled in the morning compared to individuals from the other collection times. The metabolites that were most negatively
<figure><img src="image_6.png" /><figcaption>**FIGURE 3** | Individual metabolite concentrations (nmol mg⁻¹ dry mass) measured in whole-body extracts of female *N. groenlandicus* collected from the field at four consecutive time points (morning, midday, afternoon, and evening) during day 4 (blue) and day 5 (red). Each point represent mean concentration (n = 8) and bars are standard errors of the mean. Differences between concentrations (log-transformed) and collection days (d), collection times (t), and the interaction between day and time (d:t) was investigated using 2-way ANOVAs and significance *p*-values are shown on each plot if significant. The metabolites are classified from A to E in the header according to functional group; A, amino acids; B, sugars; C, polyols; D, metabolic intermediates; and E, other metabolites.</figcaption></figure>

correlated to PC1 were succinic acid, fumaric acid, glyceric acid and ornithine. Further, negative correlations included several polyols (adonitol, glycerol-3-phosphate, arabitol) whose concentrations were higher in individuals collected in the afternoon than in the morning. The sugars fructose, glucose and galactose, and the amino acids isoleucine and lysine, as well as succinic acid and glycerol were positively associated with PC2 (Supplementary Figure 3), and were thus more abundant in individuals sampled in the evening. Negative associations were few, but the most negatively associated metabolite was the sugar glucose-6-phosphate which was more abundant in the morning and afternoon.

# DISCUSSION

## Diel Variations in the Thermal Tolerance of Field-Sampled Insects

In our study, we showed a linear relationship between ambient microhabitat temperature and measures of both cold and heat tolerance of field-collected specimens of *N. groenlandicus*. Given the temperatures observed in the field during summer in

Narsarsuaq includes subzero night temperatures and peak day temperatures above 40°C (Sørensen et al., 2019; Bahrndorff et al., 2021b) we advocate that the ability to withstand and remain active at high temperatures and recover fast from low temperature coma is ecologically important, especially for a univoltine species such as *N. groenlandicus*. This allows the species to, e.g., forage and mate in a transient environment and short summer season. Thus, our results point to plasticity in thermal tolerances being of strong importance for survival of insects in this region. Recently, the heat hardening capacity of *N. groenlandicus* has been examined in the laboratory, showing that heat tolerance can be increased within 45 min when the insects are exposed to high temperatures (Sørensen et al., 2019). It was also found that the heat hardening effect was reverted within 2 h when the insects were returned to cooler temperatures. Our results are consistent with these former observations, and the rapid adjustments that we observed in both heat and cold tolerance from field-sampled *N. groenlandicus* additionally support the assumption that hardening responses are important for coping with rapid changes in ambient temperature in the field. This is in contrast to the slower hardening responses observed in some temperate
<figure><img src="image_7.png" /><figcaption>**FIGURE 4** | Between-class Principal Component Analyses (top) and metabolite correlations to PC1 (bottom) based on GC-MS quantification of metabolites from whole-body extracts of female *N. groenlandicus* sampled in the field at four consecutive sampling time points (8:00 am, 12:00 pm, 4:00 pm, 8:00 pm) on (A) the thermally stable sampling day (day 4) and (B) the thermally variable sampling day (day 5). Lines represent individual sample positions respective to centroids (n = 8).</figcaption></figure>

insect species and indicates that high thermal variability of the environment can be a selective agent for rapid plastic responses (Dahlgaard et al., 1998; Bahrndorff et al., 2009; Alemu et al., 2017).

Interestingly, the adjustments occurred not only at extreme temperatures, but also at temperatures that are not considered as stressful or sub-lethal to the species. The body temperatures experienced by *N. groenlandicus*, like other insect species, might differ from the temperatures measured in the shade due to behavioral thermoregulation (e.g., Stevenson, 1985; Kearney et al., 2009), such as seeking microhabitat temperatures that deviate from air temperatures (Stevenson, 1985; Böcher and Nachman, 2001; Danks, 2004; Kearney et al., 2009). However, as sampled individuals are collected in the same microhabitat as where air temperatures are measured, we do argue that the difference between air and body temperatures is likely to be minor. The adult life-stage of *N. groenlandicus* seemingly has a high preferred body temperature of approximately 30°C (Böcher and Nachman, 2001) and an extreme thermal tolerance breadth spanning from critical lower limits (CTmin) of -3.4°C to critical upper thermal limits (CTmax) of 48.5-52°C (Böcher and Nachman, 2001; Sørensen et al., 2019; Bahrndorff et al., 2021b). Thus, only the lowest temperatures recorded in the field approximated sub-lethal conditions for the species. Typically, hardening responses are described as being induced

by stressful conditions (Angilletta, 2009). For example it has been described from laboratory and field studies performed on *Drosophila* spp. that hardening temperatures ca. 10-15°C above optimal rearing temperatures are needed to induce adaptive increases in heat tolerance, and cause upregulation of heat shock proteins (Sørensen et al., 2003; King and MacRae, 2015). Here, despite exposure to temperatures well within their thermal comfort zone, we show that heat and cold tolerance changed daily in *N. groenlandicus*. This finding might represent an evolutionary adaptation to the extreme climatic variations of arctic environments, but may also suggest that temperature variation act in concert with changes in air humidity and/or other climate variables to affect thermal tolerances as found for several other polar and sub-polar species (Block et al., 1994; Hodkinson et al., 1996; Benoit et al., 2009; Everatt et al., 2015).

Another important discovery was that the patterns of plastic changes in cold tolerance were similar for males and females, while distinct patterns were seen for heat tolerance for the two sexes (Figure 2). Females had a higher HKDT (+5 min HKDT) compared to males and tended to exhibit a stronger plasticity for that trait when field temperatures varied. Often, studies on thermal plasticity in insects find that the variation in upper thermal limits is constrained, and less plastic, compared to lower thermal limits (reviewed by Hoffmann et al., 2013; see also Chown, 2001; Overgaard et al., 2011; Alford et al., 2012). This
is likely resulting from the difficulty of insects to seek shelter from cold temperatures, thus resulting in a stronger selection pressure for plasticity of cold tolerance (Hoffmann et al., 2013). Our findings suggest that the selection pressure for cold tolerance may have been similar in male and female *N. groenlandicus* because no differences were observed in T<sub>recovery</sub>. Conversely, heat stress is often countered by behavioral thermoregulation in ectotherms, for instance by seeking shadow or migrating below-ground (Huey and Tewksbury, 2009; Kearney et al., 2009). The higher heat tolerance and plasticity for this trait in females could be explained by the univoltine life history and the short arctic summers requiring females to seek out warm temperatures to rapidly complete their life cycle (Bahrndorff et al., 2021a). Further, our results indicate a trade-off between heat and cold tolerance. Thus, individuals sampled at middays and afternoons are overall more heat tolerant and less cold tolerant compared to individuals sampled during mornings and evenings. Similar results have been found for thermal tolerance of the fruit fly *Drosophila melanogaster* kept under natural and semi-natural conditions (Overgaard and Sørensen, 2008; Schou et al., 2015). A consequence might be maladaptive plastic responses to environmental cues, as climatic conditions are prospected to become more unpredictably variable in the future (Kingsolver and Huey, 1998; Huey et al., 1999; Manenti et al., 2014).

## Daily Thermal Variations and Metabolic Fingerprints

The separation of the metabolic fingerprints from field-sampled *N. groenlandicus* was much stronger when individuals were collected during the thermally variable day (day 5), as compared with the less temperature variable day (day 4). This finding supports our hypothesis that microhabitat environmental conditions have a strong impact on diurnal changes of the physiology of adult *N. groenlandicus*. It also suggests that the observed diurnal variation in metabolic fingerprints cannot be explained by circadian clock regulations alone because similar patterns would be expected on the two collection days if that was the case.

On the less variable day, the average temperature prior to testing the insects was rather similar for the four collection periods (maximal temperature difference of 4.9°C; **Figure 1**). The measured changes in metabolite concentrations on this day were thus mostly independent of temperature, and rather reflected adjustments in energetic metabolism over the day or circadian regulated responses independent of temperature. This assumption is supported by the grouping pattern of the metabolomic profiles of individuals collected at the four different time points of day 4. Metabolomic profiles in the morning and evening were more similar and separated from those of individuals collected on the midday and afternoon. This pattern may reflect that the activity of the individuals was higher during midday and afternoon and in turn increased the energetic needs and metabolism in general. Consistently, sugars (glucose, fructose, galactose), some metabolic intermediates (citric and fumaric acid), and a range of free amino acids (phenylalanine, valine, serine, glutamine, and tyrosine), all being important

substrates for glycolysis and Krebs cycle, varied in rhythmic patterns on both day 4 and 5 (**Figure 3**). These patterns might constitute circadian clock mechanisms that are regulated independently of temperature, humidity and other variable abiotic factors as seen, e.g., in *D. melanogaster* (Rhoades et al., 2018).

On the temperature variable day, the pattern of separation was markedly different than the one reported for day 4. The groups separating the strongest and explaining most of the variation in the data belonged to the individuals collected in the morning and afternoon, representing the time points with the lowest and the highest temperatures of day 5. Throughout this day, the sugars fructose, glucose and galactose, occurred in higher quantities compared to day 4 and especially fructose and glucose accumulated 2-3 fold in the evening in the individuals sampled on day 5, despite the temperature not being different from the temperature on day 4 at this time point (**Figures 1, 3**). Typically, sugar accumulation is associated with cold shock responses (Jagdale et al., 2005; Lalouette et al., 2007; Michaud and Denlinger, 2007; Overgaard et al., 2007; Holmstrup et al., 2010; Teets et al., 2011) and seasonal preparation for diapause (Koštál et al., 2001; Watanabe, 2002; Vasquez et al., 2019). Sugars have osmoprotective properties that may play important protective roles in cold tolerance possibly through stabilization of cell membranes and macromolecular structures even at low concentration (Gekko and Timasheff, 1981; Yancey, 2005; Koštál et al., 2016) or by maintaining haemolymph osmolality despite low [Na<sup>+</sup>] and [K<sup>+</sup>] due to cold exposure (MacMillan et al., 2015). Thus, sugars might be accountable for the higher cold tolerance observed in individuals on days characterized by large temperature amplitude.

Polyols accumulated during the warmest periods of the days and especially inositol, fluctuated on day 5. Accumulation of the polyols mannitol and sorbitol in whiteflies and aphids with daily warm peaks has been associated with increased heat tolerance under natural and semi-natural conditions (Hendrix and Salvucci, 1998; Wolfe et al., 1998; Salvucci et al., 2000). This could indicate that polyols contribute to regulation of heat tolerance in *N. groenlandicus*.

It is possible that oscillations of sugars and polyols were affected by temperature-dependent activity patterns such as feeding, mating and general metabolism, and this might confound the effects of temperature and humidity alone on thermal tolerance. Foraging or feeding rates are partly governed by upper and lower activity thresholds of organisms (Everatt et al., 2013). Our own unpublished data on the activity of *N. groenlandicus* show that the species is virtually inactive at temperatures below 15°C and activity peaks at 35-40°C. This might suggest that feeding is constricted to the warmest periods of the day (typically between 20 and 30°C at the given study site). *Nysius* species feed on phloem sap and plant seeds (Böcher, 1972; Broadle et al., 1986; Böcher et al., 2015; Tiwari and Wratten, 2019; Maharjan et al., 2020), and therefore ingest large quantities of sucrose, which is produced by photosynthesis in plants. In other hemipteran phloem-feeders, sucrose is hydrolyzed to glucose and fructose when ingested and rapidly converted to trehalose or polyols, which are less toxic compounds to store
in the hemolymph at high concentrations (Becker et al., 1996; Hendrix and Salvucci, 1998). Thus, there might be a direct link between feeding behavior and sugar and polyol levels. For instance, trehalose concentrations increased more in the afternoon/evening on the warm day when feeding rates are expected to be higher (Figure 3).

Additionally, selective feeding on protein- and lipid rich diets impact on thermal tolerance in several insect species. For instance, the dung beetle *Thorectes lusitanicus* has been found to selectively supplement its diet with acorn, and experiments showed that beetles that were fed on acorn had a hemolymph supercooling point that was 5°C lower than individuals fed on cow-dung (Verdú et al., 2010). The shift was associated with alterations in hemolymph cryoprotectant content. Likewise, Rho and Lee (2017) showed that the beetle *Tenebrio molitor* selectively chose a carbohydrate rich diet at cold and warm temperatures, opposed to a more balanced protein-carbohydrate diet at intermediate temperatures. Switches in feeding preference with cold stress have also been found for *Drosophila* species (Brankatschk et al., 2018; Strassburger and Teleman, 2018). These aspects should be examined further in future studies.

Whether the oscillations of sugars and polyols constitute protective responses rather than consequences of altered energetic metabolism or temperature-dependent feeding is not evident from our results. However, the contrasting patterns of sugar and polyol oscillation may indicate that polyols are converted to sugars during the coldest times of the thermally variable day, maybe as a protective response. It might also explain the negative trade-off observed between heat and cold tolerance in this study. However, this is speculative and should be examined further in future studies along with common garden experiments revealing the importance of circadian rhythm regulation of *N. groenlandicus* behavior and physiology.

## CONCLUSION

Here, we showed that thermal tolerance was correlated with ambient microhabitat temperature in the Greenlandic seed bug, *N. groenlandicus*. Thus, we show that hardening responses observed under constant laboratory conditions in a previous study (Sørensen et al., 2019) occur also in field settings for this species. Interestingly, we found that the plastic adjustments of thermal tolerance occurred at relative benign temperatures contrasting experimental evidence from laboratory studies on primarily model species, suggesting that these responses occur at sub- or supra lethal temperatures. Thus, plasticity of thermal tolerance is likely affected by multiple factors under natural conditions and constitute an example of evolutionary adaption to the extreme and variable Arctic and sub-Arctic habitats of *N. groenlandicus*. Further, we showed that field heat hardening causes increased heat tolerance and reduced cold tolerance and vice versa with cold acclimation, suggesting a trade-off. GC-MS investigation of field collected individuals revealed candidate metabolites that are regulated according to thermal and other abiotic and biotic conditions that vary on a diurnal basis. Distinct metabolic fingerprints associated with temperature on

the thermally variable day (day 5) were not observed at the thermal stable day (day 4). This suggests that our results cannot be explained by circadian clock regulated mechanisms alone, but that these metabolites are partly regulated by temperature variability (or variation in correlated environmental or physiological variables) and constitute important physiological mechanisms controlling diurnal variability in thermal tolerances in the species. The strong plastic responses observed in heat and cold tolerance and in the metabolomic profiles in our study suggests ongoing strong selection for plasticity in this highly fluctuating polar environment. The genetic architecture of these traits are mainly investigated in model organisms where studies on *D. melanogaster* suggest significant heritable variation for plasticity of thermal stress tolerance traits and metabolite profiles (Gerken et al., 2015; Hangartner and Hoffmann, 2015; Ørsted et al., 2017; Rohde et al., 2021).

## DATA AVAILABILITY STATEMENT

All data, including raw field temperature files and GC-MS data, are presented in this article/Supplementary Material.

## AUTHOR CONTRIBUTIONS

SB, TNK, and NKN conceived the ideas and designed the methodology. MHS, SB, TNK, and NKN collected the data. DR, HC, and NKN processed samples for metabolomics analysis and analyzed the metabolomics data. TNK and NKN led the writing of the manuscript. All authors contributed to the drafts and gave final approval for publication.

## FUNDING

This work was supported by a grant from the Danish Council for Independent Research (DFF8021-00014B), the European co-funded Partnership BiodivClim-191 ASICS (0156-00024B), and Carlsberg (Carlsbergfondet, CF17-0415). Travel grants was provided by North2North and SAFT (Selskab for Arktisk Forskning og Teknologi). Metabolomics analyses were supported by IRP CNRS "PRICES" and IPEV 136 "Subanteco."

## ACKNOWLEDGMENTS

We thank Finn Thomsen and Storch Lund for access to field sites and laboratory facilities in Narsarsuaq, Greenland. Further, we are grateful to EcoChimie Platform (EcoChim) from UMS OSUR 3343 for target quantitative metabolomics analyses.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2022.818485/full#supplementary-material
REFERENCES

Alemu, T., Alemneh, T., Pertoldi, C., Ambelu, A., and Bahrdorff, S. (2017). Costs and benefits of heat and cold hardening in a soil arthropod. *Biol. J. Linn. Soc.* 122, 765–773. doi: 10.1093/biolinnean/blx092

Alford, L., Blackburn, T. M., and Bale, J. S. (2012). Effect of latitude and acclimation on the lethal temperatures of the peach-potato aphid *Myzus persicae*. *Agric. For. Entomol.* 14, 69–79. doi: 10.1111/j.1461-9563.2011.00553.x

Angilletta, M. J. (2009). "Thermal acclimation," in *Thermal Adaptation: A Theoretical and Empirical Synthesis*, ed. M. J. Angilletta (New York, NY: Oxford University Press), 126–156.

Araújo, M. B., Ferri-Yáñez, F., Bozinovic, F., Marquet, P. A., Valladares, F., and Chown, S. L. (2013). Heat freezes niche evolution. *Ecol. Lett.* 16, 1206–1219. doi: 10.1111/ele.12155

Bahrdorff, S., Gertsen, S., Pertoldi, C., and Kristensen, T. N. (2016). Investigating thermal acclimation effects before and after a cold shock in *Drosophila melanogaster* using behavioural assays. *Biol. J. Linn. Soc.* 117, 241–251. doi: 10.1111/bij.12659

Bahrdorff, S., Lauritzen, J. M. S., Sørensen, M. H., Noer, N. K., and Kristensen, T. N. (2021b). Responses of terrestrial polar arthropods to high and increasing temperatures. *J. Exp. Biol.* 224:jeb230797. doi: 10.1242/jeb.230797

Bahrdorff, S., Alemu, T., Kristensen, T. N., Sørensen, M. H., Höye, T., and Holmstrup, M. (2021a). Thermal adaptations of adults and eggs in the arctic insect, *Nysius groenlandicus*. *Polar Biol.* 44, 491–498. doi: 10.1007/s00300-021-02807-6

Bahrdorff, S., Mariën, J., Loeschcke, V., and Ellers, J. (2009). Dynamics of heat-induced thermal stress resistance and hsp70 expression in the springtail, *Orchesella cincta*. *Funct. Ecol.* 23, 233–239. doi: 10.1111/j.1365-2435.2009.01541.x

Bahrdorff, S., Petersen, S. O., Loeschcke, V., Overgaard, J., and Holmstrup, M. (2007). Differences in cold and drought tolerance of high arctic and sub-arctic populations of *Megaphorura arctica* Tullberg 1876 (Onychiuridae: Collembola). *Cryobiology* 55, 315–323. doi: 10.1016/j.cryobiol.2007.09.001

Bak, C. W., Bahrdorff, S., Noer, N. K., Jørgensen, L. B., Overgaard, J., and Kristensen, T. N. (2020). Comparison of static and dynamic assays when quantifying thermal plasticity of drosophilids. *Insects* 11:537. doi: 10.3390/insects11080537

Becker, A., Schlöder, P., Steele, J. E., and Wegener, G. (1996). The regulation of trehalose metabolism in insects. *Experientia* 52, 433–439.

Benoit, J. B., Lopez-martinez, G., Elnitsky, M. A., Lee, R. E., and Denlinger, D. L. (2009). Dehydration-induced cross tolerance of *Belgica antarctica* larvae to cold and heat is facilitated by trehalose accumulation. *Comp. Biochem. Physiol. A Mol. Integr. Physiol.* 152, 518–523. doi: 10.1016/j.cbpa.2008.12.009

Block, W., Webb, N. R., Coulson, S., Hodkinson, I. D., and Worland, M. R. (1994). Thermal adaptation in the Arctic collembolan *Onychiurus arcticus* (Tullberg). *J. Insect Physiol.* 40, 715–722. doi: 10.1016/0022-1910(94)90099-X

Böcher, J. (1972). Feeding biology of *Nysius groenlandicus* (Heteroptera: Miridae) in Greenland. With a note on oviposition in relation to food-source and dispersal of the species. *Medd. Grønland* 191, 1–41.

Böcher, J., Kristensen, N. P., Pape, T., and Vilhelmsen, L. (2015). *The Greenland Entomofauna: An Identification Manual of Insects, Spiders and Their Allies* (Fauna Entomologica Scandinavica). Leiden: Brill. doi: 10.1163/ej.9789004157705.i-265.29

Böcher, J., and Nachman, G. (2001). Temperature and humidity responses of the arctic-alpine seed bug *Nysius groenlandicus*. *Entomol. Exp. Appl.* 99, 319–330. doi: 10.1046/j.1570-7458.2001.00831.x

Brankatschk, M., Gutmann, T., Knittelfelder, O., Palladini, A., Prince, E., Grzybek, M., et al. (2018). A temperature-dependent switch in feeding preference improves *Drosophila* development and survival in the cold. *Dev. Cell* 46, 781–793.e4. doi: 10.1016/j.DEVCEL.2018.05.028

Broadle, R., Simpson, B., and Beavis, C. (1986). Damage by "Nysius" spp. (Hemiptera: Lygaeidae) in non-stressed sunflower ('Helianthus annus' L.) crops. *Gen. Appl. Entomol.* 18, 17–24.

Buckley, B. A., Owen, M. E., and Hofmann, G. E. (2001). Adjusting the thermostat: the threshold induction temperature for the heat-shock response in intertidal mussels (genus *Mytilus*) changes as a function of thermal history. *J. Exp. Biol.* 204, 3571–3579.

Chown, S. L. (2001). Physiological variation in insects: hierarchical levels and implications. *J. Insect Physiol.* 47, 649–660. doi: 10.1016/S0022-1910(00)00016-3

Chown, S. L., Jumbam, K. R., Sørensen, J. G., and Terblanche, J. S. (2009). Phenotypic variance, plasticity and heritability estimates of critical thermal limits depend on methodological context. *Funct. Ecol.* 23, 133–140. doi: 10.1111/j.1365-2435.2008.01481.x

Chown, S. L., and Terblanche, J. S. (2006). Physiological diversity in insects: ecological and evolutionary contexts. *Adv. Insect Phys.* 33, 50–152. doi: 10.1016/S0065-2806(06)33002-0

Colinet, H., and Hoffmann, A. A. (2012). Comparing phenotypic effects and molecular correlates of developmental, gradual and rapid cold acclimation responses in *Drosophila melanogaster*. *Funct. Ecol.* 26, 84–93. doi: 10.1111/j.1365-2435.2011.01898.x

Colinet, H., Sinclair, B. J., Vernon, P., and Renault, D. (2015). Insects in fluctuating thermal environments. *Ann. Rev. Entomol.* 60, 123–140. doi: 10.1146/annurev-ento-010814-021017

Convey, P., Coulson, S. J., Worland, M. R., and Sjöblom, A. (2018). The importance of understanding annual and shorter-term temperature patterns and variation in the surface levels of polar soils for terrestrial biota. *Polar Biol.* 41, 1587–1605.

Dahlgaard, J., Loeschcke, V., Michalak, P., and Justesen, J. (1998). Induced thermotolerance and associated expression of the heat-shock protein Hsp70 in adult *Drosophila melanogaster*. *Funct. Ecol.* 12, 786–793. doi: 10.1046/j.1365-2435.1998.00246.x

Danks, H. V. (2004). Seasonal adaptations in Arctic insects. *Integr. Comp. Biol.* 44, 85–94. doi: 10.1093/ICB/44.2.85

Davey, M. C., Pickup, J., and Block, W. (2021). Temperature variation and its biological significance in fellfield habitats on a maritime Antarctic island. *Antarct. Sci.* 4, 383–388. doi: 10.1017/S0954102092000567

Deere, J. A., Sinclair, B. J., Marshall, D. J., and Chown, S. L. (2006). Phenotypic plasticity of thermal tolerances in five oribatid mite species from sub-Antarctic Marion Island. *J. Insect Physiol.* 52, 693–700. doi: 10.1016/j.jinsphys.2006.03.009

Everatt, M. J., Bale, J. S., Convey, P., Worland, M. R., and Hayward, S. A. L. (2013). The effect of acclimation temperature on thermal activity thresholds in polar terrestrial invertebrates. *J. Insect Physiol.* 59, 1057–1064. doi: 10.1016/j.jinsphys.2013.08.003

Everatt, M. J., Convey, P., Bale, J. S., Worland, M. R., and Hayward, S. A. (2015). Responses of invertebrates to temperature and water stress: a polar perspective. *J. Therm. Biol.* 54, 118–132. doi: 10.1016/j.jtherbio.2014.05.004

Fields, P. G., Fleurat-Lessard, F., Lavenseau, L., Febvay, G., Peypelut, L., and Bonnot, G. (1998). The effect of cold acclimation and deacclimation on cold tolerance, trehalose and free amino acid levels in *Sitophilus granarius* and *Cryptolestes ferrugineus* (Coleoptera). *J. Insect Physiol.* 44, 955–965. doi: 10.1016/S0022-1910(98)00055-9

Fischer, K., Kölzow, N., Höltje, H., and Karl, I. (2011). Assay conditions in laboratory experiments: is the use of constant rather than fluctuating temperatures justified when investigating temperature-induced plasticity? *Oecologia* 166, 23–33. doi: 10.1007/s00442-011-1917-0

Foray, V., Desouhant, E., Voituron, Y., Larvor, V., Renault, D., Colinet, H., et al. (2013). Does cold tolerance plasticity correlate with the thermal environment and metabolic profiles of a parasitoid wasp? *Comp. Biochem. Physiol. A Mol. Integr. Physiol.* 164, 77–83. doi: 10.1016/j.cbpa.2012.10.018

Fusco, G., and Minelli, A. (2010). Phenotypic plasticity in development and evolution: facts and concepts. *Philos. Trans. R. Soc. B* 365, 547–556. doi: 10.1098/rstb.2009.0267

Gekko, K., and Timasheff, S. N. (1981). Mechanism of protein stabilization by glycerol: preferential hydration in glycerol-water mixtures. *Biochemistry* 20, 4667–4676. doi: 10.1021/bi00519a023

Gerken, A. R., Eller, O. C., Hahn, D. A., and Morgan, T. J. (2015). Constraints, independence, and evolution of thermal plasticity: probing genetic architecture of long-and short-term thermal acclimation. *Proc. Natl. Acad. Sci. U.S.A.* 112, 4399–4404.

Gracey, A. Y., Chaney, M. L., Boomhower, J. P., Tyburczy, W. R., Connor, K., and Somero, G. N. (2008). Rhythms of gene expression in a fluctuating intertidal environment. *Curr. Biol.* 14, 1501–1507. doi: 10.1016/j.cub.2008.08.049

Hangartner, S., and Hoffmann, A. A. (2015). Evolutionary potential of multiple measures of upper thermal tolerance in *Drosophila melanogaster*. *Funct. Ecol.* 30, 442–452. doi: 10.1111/1365-2435.12499

Harris, D. T. (2018). Biobanking and omics. *Front. Biol.* 13:287–292. doi: 10.1007/s11515-018-1505-3

Hendrix, D. L., and Salvucci, M. E. (1998). Polyol metabolism in homopterans at high temperatures: accumulation of mannitol in aphids (Aphididae:
Homoptera) and sorbitol in whiteflies (Aleyrodidae: Homoptera). *Comp. Biochem. Physiol. A Mol. Integr. Physiol.* **120**, 487–494. doi: **10.1016/S1095-6433(98)10058-2**

Hodkinson, I. D., Coulson, S. J., Webb, N. R., and Block, W. (1996). Can high Arctic soil microarthropods survive elevated summer temperatures? *Funct. Ecol.* **10**, 314–321. doi: **10.2307/2390278**

Hoffmann, A. A., Chown, S. L., and Clusella-Trullas, S. (2013). Upper thermal limits in terrestrial ectotherms: how constrained are they? *Funct. Ecol.* **27**, 934–949. doi: **10.1111/j.1365-2435.2012.02036.x**

Hoffmann, A. A., Sørensen, J. G., and Loeschcke, V. (2003). Adaptation of *Drosophila* to temperature extremes: bringing together quantitative and molecular approaches. *J. Therm. Biol.* **28**, 175–216. doi: **10.1016/S0306-4565(02)00057-8**

Holmstrup, M., Bayley, M., Pedersen, S. A., and Zachariassen, K. E. (2010). "Interactions between cold, desiccation and environmental toxin," in *Low Temperature Biology of Insects*, eds D. L. Denlinger and R. E. Lee (Cambridge: Cambridge University Press), 166–190.

Holmstrup, M., Hedlund, K., and Boriss, H. (2002). Drought acclimation and lipid composition in *Folsomia candida*: implications for cold shock, heat shock and acute desiccation stress. *J. Insect Physiol.* **48**, 961–970. doi: **10.1016/S0022-1910(02)00175-0**

Huey, R. B., Berrigan, D., Gilchrist, G. W., and Herron, J. C. (1999). Testing the adaptive significance of acclimation: a strong inference approach. *Am. Zool.* **39**, 323–336. doi: **10.1093/icb/39.2.323**

Huey, R. B., and Tewksbury, J. J. (2009). Can behavior douse the fire of climate warming? *Proc. Natl. Acad. Sci. U.S.A.* **106**, 3647–3648. doi: **10.1073/pnas.0900934106**

Jagdale, G. B., Grewal, P. S., and Salminen, S. O. (2005). Both heat-shock and cold-shock influence trehalose metabolism in an entomopathogenic nematode. *J. Parasitol.* **91**, 988–994. doi: **10.1645/GE-504R.1**

Javal, M., Renault, D., and Colinet, H. (2016). Impact of fluctuating thermal regimes on *Drosophila melanogaster* survival to cold stress. *Anim. Biol.* **66**, 427–444. doi: **10.1163/15707563-00002510**

Jensen, A., Alemu, T., Alemneh, T., Pertoldi, C., and Bahrndorff, S. (2019). Thermal acclimation and adaptation across populations in a broadly distributed soil arthropod. *Funct. Ecol.* **33**, 833–845. doi: **10.1111/1365-2435.13291**

Kearney, M., Gillingham, P. K., Bramer, I., Duffy, J. P., and MacLean, I. M. D. (2020). A method for computing hourly, historical, terrain-corrected microclimate anywhere on earth. *Methods Ecol. Evol.* **11**, 38–43. doi: **10.1111/2041-210X.13330**

Kearney, M., and Porter, W. (2009). Mechanistic niche modelling: combining physiological and spatial data to predict species' ranges. *Ecol. Lett.* **12**, 334–350. doi: **10.1111/j.1461-0248.2008.01277.x**

Kearney, M., Shine, R., and Porter, W. P. (2009). The potential for behavioral thermoregulation to buffer “cold-blooded” animals against climate warming. *Proc. Natl. Acad. Sci. U.S.A.* **106**, 3835–3840. doi: **10.1073/pnas.0808910310**

Kelty, J. (2007). Rapid cold-hardening of *Drosophila melanogaster* in a field setting. *Physiol. Entomol.* **32**, 343–350. doi: **10.1111/j.1365-3032.2007.00584.x**

Ketola, T., and Kristensen, T. N. (2017). Experimental approaches for testing if tolerance curves are useful for predicting fitness in fluctuating environments. *Front. Ecol. Evol.* **5**:129. doi: **10.3389/fevo.2017.00129**

King, A. M., and MacRae, T. H. (2015). Insect heat shock proteins during stress and diapause. *Ann. Rev. Entomol.* **60**, 59–75. doi: **10.1146/annurev-ento-011613-162107**

Kingsolver, J. G., Higgins, J. K., and Augustine, K. E. (2015). Fluctuating temperatures and ectotherm growth: distinguishing non-linear and time-dependent effects. *J. Exp. Biol.* **218**, 2218–2225. doi: **10.1242/jeb.120733**

Kingsolver, J. G., and Huey, R. B. (1998). Evolutionary analyses of morphological and physiological plasticity in thermally variable environments. *Am. Zool.* **38**, 545–560. doi: **10.1093/icb/38.3.545**

Kingsolver, J. G., and Nagle, A. (2007). Evolutionary divergence in thermal sensitivity and diapause of field and laboratory populations of *Manduca sexta*. *Physiol. Biochem. Zool.* **80**, 473–479. doi: **10.1086/519962**

Kinzner, M. T., Kinzner, M. C., Kaufmann, R., Hoffmann, A. A., Arthofer, W., Schlick-Steiner, B. C., et al. (2019). Is temperature preference in the laboratory ecologically relevant for the field? The case of *Drosophila nigrosparsa*. *Glob. Ecol. Conserv.* **18**:e00638. doi: **10.1016/j.gecco.2019.e00638**

Koštál, V., Korbelová, J., Poupardin, R., Moos, M., and Šimek, P. (2016). Arginine and proline applied as food additives stimulate high freeze tolerance in larvae of *Drosophila melanogaster*. *J. Exp. Biol.* **219**, 2358–2367. doi: **10.1242/jeb.142158**

Koštál, V., Šlachta, M., and Šimek, P. (2001). Cryoprotective role of polyols independent of the increase in supercooling capacity in diapausing adults of *Pyrrhocoris apterus* (Heteroptera: Insecta). *Comp. Biochem. Physiol. B Biochem. Mol. Biol.* **130**, 365–374. doi: **10.1016/S1096-4959(01)00441-9**

Koveos, D. S. (2001). Rapid cold hardening in the olive fruit fly *Bactrocera oleae* under laboratory and field conditions. *Entomol. Exp. Appl.* **101**, 257–263. doi: **10.1046/j.1570-7458.2001.00910.x**

Kristensen, T. N., Hoffmann, A., Overgaard, J., Sørensen, J. G., Hallas, R., and Loeschcke, V. (2008). Costs and benefits of cold acclimation in field-released *Drosophila*. *Proc. Natl. Acad. Sci. U.S.A.* **105**, 216–221. doi: **10.1073/pnas.0708074105**

Kristensen, T. N., Ketola, T., and Kronholm, I. (2020). Adaptation to environmental stress at different timescales. *Ann. N. Y. Acad. Sci.* **1476**, 5–12. doi: **10.1111/nyas.13974**

Kristensen, T. N., Overgaard, J., Hoffmann, A. A., Nielsen, N. C., and Malmendal, A. (2012). Inconsistent effects of developmental temperature acclimation on low-temperature performance and metabolism in *Drosophila melanogaster*. *Evol. Ecol. Res.* **14**, 821–837.

Lalouette, L., Kostal, V., Colinet, H., Gagneul, D., and Renault, D. (2007). Cold exposure and associated metabolic changes in adult tropical beetles exposed to fluctuating thermal regimes. *FEBS J.* **274**, 1759–1767. doi: **10.1111/j.1742-4658.2007.05723.x**

Lembrechts, J. J., and Lenoir, J. (2020). Microclimatic conditions anywhere at any time! *Glob. Change Biol.* **26**, 337–339. doi: **10.1111/gcb.14942**

Loeschke, V., and Hoffmann, A. A. (2007). Consequences of heat hardening on a field fitness component in *Drosophila* depend on environmental temperature. *Am. Nat.* **169**, 175–183. doi: **10.1086/510632**

Maclean, I. M. D., Duffy, J. P., Haesen, S., Govaert, S., De Frenne, P., Vanneste, T., et al. (2021). On the measurement of microclimate. *Methods Ecol. Evol.* **12**, 1397–1410. doi: **10.1111/2041-210x.13627**

MacMillan, H. A., Ferguson, L. V., Nicolai, A., Donini, A., Staples, J. F., and Sinclair, B. J. (2015). Parallel ionoregulatory adjustments underlie phenotypic plasticity and evolution of *Drosophila* cold tolerance. *J. Exp. Biol.* **218**, 423–432. doi: **10.1242/jeb.115790**

Maharjan, R., Yoon, Y., Jang, Y., Jeong, M., Jung, T. W., Ha, T. J., et al. (2020). Oviposition and development response of perilla seed bugs (*Nysius* sp.) (Heteroptera: Lygaeidae) to five crop seeds. *J. Appl. Entomol.* **144**, 806–816. doi: **10.1111/JEN.12814**

Manenti, T., Sørensen, J. G., Moghadam, N. N., and Loeschcke, V. (2014). Predictability rather than amplitude of temperature fluctuations determines stress resistance in a natural population of *Drosophila simulans*. *J. Evol. Biol.* **27**, 2113–2122. doi: **10.1111/jeb.12463**

Michaud, M. R., and Denlinger, D. L. (2007). Shifts in the carbohydrate, polyol, and amino acid pools during rapid cold-hardening and diapause-associated cold-hardening in flesh flies (*Sarcophaga crassipalpis*): a metabolomic comparison. *J. Comp. Physiol. B Biochem. Syst. Environ. Physiol.* **177**, 753–763. doi: **10.1007/s00360-007-0172-5**

Michaud, R., and Denlinger, D. L. (2010). "Genomics, proteomics and metabolomics: finding the other players in insect cold-tolerance," in *Low Temperature Biology of Insects*, eds D. L. Denlinger and R. E. Lee (New York, NY: Cambridge University Press), 91–115.

Mitchell, K. A., and Hoffmann, A. A. (2010). Thermal ramping rate influences evolutionary potential and species differences for upper thermal limits in *Drosophila*. *Funct. Ecol.* **24**, 694–700. doi: **10.1111/j.1365-2435.2009.01666.x**

Noer, N. K., Ørsted, M., Schiffer, M., Hoffmann, A. A., Bahrndorff, S., and Kristensen, T. N. (2022). Into the wild – a field study on the evolutionary and ecological importance of thermal plasticity in ectotherms across temperate and tropical regions. *Philos. Trans. R. Soc. B* **376**:20210004. doi: **10.1098/rstb.2021.0004**

Noer, N. K., Pagter, M., Bahrndorff, S., Malmendal, A., and Kristensen, T. N. (2020). Impacts of thermal fluctuations on heat tolerance and its metabolomic basis in *Arabidopsis thaliana*, *Drosophila melanogaster*, and *Orchesella cincta*. *PLoS One* **15**:e0237201. doi: **10.1371/journal.pone.0237201**
Ørsted, M., Rohde, P. D., Hoffmann, A. A., Sørensen, P., and Kristensen, T. N. (2017). Environmental variation partitioned into separate heritable components. *Evolution* 72, 136–152. doi: 10.1111/evo.13391

Overgaard, J., Kristensen, T. N., Mitchell, K. A., and Hoffmann, A. A. (2011). Thermal tolerance in widespread and tropical *Drosophila* species: does phenotypic plasticity increase with latitude? *Am. Nat.* 178, S80–S96. doi: 10.1086/661780

Overgaard, J., Malmendal, A., Sørensen, J. G., Bundy, J. G., Loeschcke, V., Nielsen, N. C., et al. (2007). Metabolomic profiling of rapid cold hardening and cold shock in *Drosophila melanogaster*. *J. Insect Physiol.* 53, 1218–1232. doi: 10.1016/j.jinsphys.2007.06.012

Overgaard, J., and Sørensen, J. G. (2008). Rapid thermal adaptation during field temperature variations in *Drosophila melanogaster*. *Cryobiology* 56, 159–162. doi: 10.1016/j.cryobiol.2008.01.001

Overgaard, J., Sørensen, J. G., Com, E., and Colinet, H. (2014). The rapid cold hardening response of *Drosophila melanogaster*: complex regulation across different levels of biological organization. *J. Insect Physiol.* 62, 46–53. doi: 10.1016/j.jinsphys.2014.01.009

R Core Team (2020). *R: A Language and Environment for Statistical Computing*. Vienna: R Foundation for Statistical Computing.

Rho, M. S., and Lee, K. P. (2017). Temperature-driven plasticity in nutrient use and preference in an ectotherm. *Oecologia* 185, 401–413. doi: 10.1007/S00442-017-3959-4/FIGURES/4

Rhoades, S. D., Nayak, K., Zhang, S. L., Sehgal, A., and Weljie, A. M. (2018). Circadian- and light-driven metabolic rhythms in *Drosophila melanogaster*. *J. Biol. Rhythms* 33, 126–136. doi: 10.1177/0748730417753003

Rohde, P. D., Kristensen, T. N., Sarup, P., Muñoz, J., and Malmendal, A. (2021). Prediction of complex phenotypes using the *Drosophila melanogaster* metabolome. *Heredity* 126, 717–732. doi: 10.1038/s41437-021-00404-1

Ruel, J. J., and Ayres, M. P. (1999). Jensen's inequality predicts effects of environmental variation. *Trends Ecol. Evol.* 14, 361–366. doi: 10.1016/S0169-5347(99)01664-X

Salvucci, M. E., Stecher, D. S., and Henneberry, T. J. (2000). Heat shock proteins in whiteflies, an insect that accumulates sorbitol in response to heat stress. *J. Therm. Biol.* 25, 363–371. doi: 10.1016/S0306-4565(99)00108-4

Scheiner, S. M. (1993). Genetics and evolution of phenotypic plasticity. *Ann. Rev. Ecol. Evol. Syst.* 24, 35–68. doi: 10.1146/annurev.es.24.110193.000343

Schou, M. F., Kristensen, T. N., Pedersen, A., Karlsson, G., Loeschcke, V., and Malmendal, A. (2017). Metabolic and functional characterization of effects of developmental temperature in *Drosophila melanogaster*. *Am. J. Physiol. Regul. Integr. Comp. Physiol.* 312, R211–R222. doi: 10.1152/ajpregu.00268.2016

Schou, M. F., Loeschcke, V., and Kristensen, T. N. (2015). Strong costs and benefits of winter acclimatization in *Drosophila melanogaster*. *PLoS One* 10:e0130307. doi: 10.1371/journal.pone.0130307

Sears, M. W., Riddell, E. A., Rusch, T. W., and Angilletta, M. J. (2019). The world still is not flat: lessons learned from organismal interactions with environmental heterogeneity in terrestrial environments. *Integr. Comp. Biol.* 59, 1049–1058. doi: 10.1093/icb/icz130

Sheldon, K. S., Padash, M., Carter, A. W., and Marshall, K. E. (2020). Different amplitudes of temperature fluctuation induce distinct transcriptomic and metabolomic responses in the dung beetle *Phanaeus vindex*. *J. Exp. Biol.* 223(Pt 23):jeb233239. doi: 10.1242/jeb.233239

Sømme, L. (1999). The physiology of cold hardiness in terrestrial arthropods. *Eur. J. Entomol.* 96, 1–10.

Sørensen, J. G., Kristensen, T. N., and Loeschcke, V. (2003). The evolutionary and ecological role of heat shock proteins. *Ecol. Lett.* 6, 1025–1037. doi: 10.1046/j.1461-0248.2003.00528.x

Sørensen, M. H., Kristensen, T. N., Lauritzen, J. M. S., Noer, N. K., Høye, T. T., and Bahrdorff, S. (2019). Rapid induction of the heat hardening response in an Arctic insect. *Biol. Lett.* 15:20190613. doi: 10.1098/rsbl.2019.0613

Stevenson, R. D. (1985). The relative importance of behavioral and physiological adjustments controlling body temperature in terrestrial ectotherms. *Am. Nat.* 126, 362–386. doi: 10.1086/284423

Strassburger, K., and Teleman, A. A. (2018). Flies eat their veggies to survive the cold. *Dev. Cell* 46, 671–672. doi: 10.1016/j.devcel.2018.05.030

Taylor, E. N., Diele-Viegas, L. M., Gangloff, E. J., Hall, J. M., Halpern, B., Massey, M. D., et al. (2021). The thermal ecology and physiology of reptiles and amphibians: a user's guide. *J. Exp. Zool. A Ecol. Integr. Physiol.* 335, 13–44. doi: 10.1002/jez.2396

Teets, N. M., and Denlinger, D. L. (2013). Physiological mechanisms of seasonal and rapid cold-hardening in insects. *Physiol. Entomol.* 38, 105–116. doi: 10.1111/phen.12019

Teets, N. M., Kawarasaki, Y., Lee, R. E., and Denlinger, D. L. (2011). Survival and energetic costs of repeated cold exposure in the Antarctic midge, *Belgica antarctica*: a comparison between frozen and supercooled larvae. *J. Exp. Biol.* 214, 806–814. doi: 10.1242/jeb.051912

Terblanche, J. S., Deere, J. A., Clusella-Trullas, S., Janion, C., and Chown, S. L. (2007). Critical thermal limits depend on methodological context. *Proc. R. Soc. B Biol. Sci.* 274, 2935–2942. doi: 10.1098/rspb.2007.0985

Thièbaut, G., Tarayre, M., Jambon, O., Le Bris, N., Colinet, H., and Renault, D. (2020). Variation of thermal plasticity for functional traits between populations of an invasive aquatic plant from two climatic regions. *Hydrobiologia* 848, 2077–2091. doi: 10.1007/s10750-020-04452-2

Tiwari, S., and Wratten, S. D. (2019). Profile biology and management of the New Zealand endemic wheat bug, *Nysius huttoni* (Hemiptera: Lygaeidae). *J. Integr. Pest Manag.* 10, 34–35. doi: 10.1093/jipm/pmz032

Tomanek, L., and Somero, G. N. (1999). Evolutionary and acclimation-induced variation in the heat-shock responses of congeneric marine snails (genus *Tegula*) from different thermal habitats: implications for limits of thermotolerance and biogeography. *J. Exp. Biol.* 202, 2925–2936.

Tomanek, L., and Somero, G. N. (2002). Interspecific- and acclimation-induced variation in levels of heat-shock proteins 70 (hsp70) and 90 (hsp90) and heat-shock transcription factor-1 (HSF1) in congeneric marine snails (genus *Tegula*): implications for regulation of hsp gene expression. *J. Exp. Biol.* 205, 677–685. doi: 10.1242/jeb.205.5.677

van Eijsden, R. G. E., Stassen, C., Daenen, L., Van Mulders, S. E., Bapat, P. M., Siewers, V., et al. (2013). A universal fixation method based on quaternary ammonium salts (RNALater) for omics-technologies: *Saccharomyces cerevisiae* as a case study. *Biotechnol. Lett.* 35, 891–900. doi: 10.1007/s10529-013-1163-0

Vasquez, M. C., Lippert, M. R., White, C., Walter, R. K., and Tomanek, L. (2019). Proteomic changes across a natural temperature gradient in a marine gastropod. *Mar. Environ. Res.* 149, 137–147. doi: 10.1016/j.marenvres.2019.06.002

Vázquez, D. P., Gianoli, E., Morris, W. F., and Bozinovic, F. (2017). Ecological and evolutionary impacts of changing climatic variability. *Biol. Rev.* 92, 22–42. doi: 10.1111/brv.12216

Verdu, J. R., Casas, J. L., Lobo, J. M., and Numa, C. (2010). Dung beetles eat acorns to increase their ovarian development and thermal tolerance. *PLoS One* 5:e10114. doi: 10.1371/journal.pone.0010114

Watanabe, M. (2002). Cold tolerance and myo-inositol accumulation in overwintering adults of a lady beetle, *Harmonia axyridis* (Coleoptera: Coccinellidae). *Eur. J. Entomol.* 99, 5–9. doi: 10.14411/eje.2002.002

Wolfe, G. R., Hendrix, D. L., and Salvucci, M. E. (1998). A thermoprotective role for sorbitol in the silverleaf whitefly, *Bemisia argentifolii*. *J. Insect Physiol.* 44, 597–603. doi: 10.1016/s0022-1910(98)00035-3

Yancey, P. H. (2005). Organic osmolytes as compatible, metabolic and counteracting cytoprotectants in high osmolarity and other stresses. *J. Exp. Biol.* 208, 2819–2830. doi: 10.1242/jeb.01730

Zachariassen, K. E. (1985). Physiology of cold tolerance in insects. *Physiol. Rev.* 65, 799–832. doi: 10.1152/physrev.1985.65.4.799

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

**Publisher's Note:** All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Copyright © 2022 Noer, Sørensen, Colinet, Renault, Bahrdorff and Kristensen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.