Exposure to stressful low temperatures during development can result in the accumulation of deleterious physiological effects called chill injury. Metabolic imbalances, disruptions in ion homeostasis and oxidative stress contribute to the increased mortality of chill-injured insects. Interestingly, survival can be significantly increased when chill-susceptible insects are exposed to a daily warm-temperature pulse during chilling. We hypothesize that warm pulses allow for the repair of damage associated with chill injury. Here, we describe transcriptional responses during exposure to a fluctuating thermal regime, relative to constant chilled temperatures, during pupal development in the alfalfa leafcutting bee, Megachile rotundata, using a combination of RNA-seq and qPCR. Pupae were exposed to either a constant, chilled temperature of 6°C, or 6°C with a daily pulse of 20°C for 7 days. RNA-seq after experimental treatment revealed differential expression of transcripts involved in construction of cell membranes, oxidation–reduction and various metabolic processes. These mechanisms provide support for shared physiological responses to chill injury across taxa. The large number of differentially expressed transcripts observed after 7 days of treatment suggests that the initial divergence in expression profiles between the two treatments occurred upstream of the time point sampled. Additionally, the differential expression profiles observed in this study show little overlap with those differentially expressed during temperature stress in the diapause state of M. rotundata. While the mechanisms governing the physiological response to low-temperature stress are shared, the specific transcripts associated with the response differ between life stages.

In ectothermic animals such as insects, environmental temperatures play a crucial role in the animal's physiology. Exposure to stressful, low temperatures during development can have significant, detrimental effects on adult insects (Bale and Hayward, 2010; Bennett et al., 2013; O'Neill et al., 2011; Whitfield and Richards, 1992; Yocum et al., 1994, 2006). The degree to which insects respond to these temperature stressors is highly dependent upon the timing, duration and amplitude of exposure (Colinet et al., 2015). With global climate change threatening to increase temperature variance (Easterling et al., 2000; Hansen et al., 2012), insects will likely be faced with increased exposure to low, seasonally unpredictable temperatures (Vasseur et al., 2014). In chill-susceptible insects, prolonged exposure to low temperatures can lead to the accumulation of deleterious physiological effects called chill injury (Lee, 2010).

Chill injury causes significant increases in mortality (Colinet et al., 2011; Kostál, 2006; Lee, 2010; Renault et al., 2004). The downstream physiological consequences of chill injury are complex but appear to be caused in part by membrane phase transitions (Lee, 2010). Damage to cell membranes results in metabolic imbalances and perturbations in ion homeostasis (Kostál et al., 2004, 2006) that have been observed in chill-injured insects. Furthermore, temperature stress is associated with oxidative stress (Lalouette et al., 2011; Rojas and Leopold, 1996) and impairment of neuromuscular system function (Bennett et al., 2015; Yocum et al., 1994), both of which are implicated as downstream consequences of cell membrane damage (Lee, 2010).

The decreased survival and cellular damage associated with chill injury can be mitigated by exposure to transient, daily increases in temperature during chilling [termed fluctuating thermal regime (FTR)]. Exposure to these brief, daily temperature increases during chilling results in increases in survival in chill-injured insects (Chen and Denlinger, 1992; Colinet et al., 2006, 2007; Jing et al., 2005; Kostál et al., 2007, 2016; Leopold et al., 1998; Marshall and Sinclair, 2012; Renault et al., 2004, 2011; Yocum et al., 2012). Exposure to FTR during chilling allows for either the protection against or repair of damage caused by chill injury (Colinet et al., 2015; Lee, 2010). While it is possible that brief, daily exposure to warmer temperatures during chilling allows for a physiological preparation (i.e. protection) for chilling, the hypothesis that exposure to warmer temperatures allows reparative mechanisms to act has more empirical support (Colinet et al., 2006; Kostál et al., 2007; Nedvěd et al., 1998; Torson et al., 2015). FTR repairs chill injury damage in many taxa (Boardman et al., 2013; Jing et al., 2005; Kostál et al., 2007; Leopold et al., 1998; Nedvěd et al., 1998; Renault et al., 2004), but it is unclear whether the mechanisms governing the beneficial effects associated with FTR exposure are conserved.

Current management practices for the solitary alfalfa leafcutting bee, Megachile rotundata (Fabricius 1787), a valuable agricultural pollinator, can include stressful low temperature exposure at two life stages: prepupal post-diapause quiescence and pupal development. Megachile rotundata females emerge and provision offspring in the early summer. Larvae complete growth in about a week and enter a prepupal diapause for overwintering. Under commercial management, the overwintering period lasts 7–10 months (Pitts-Singer and Cane, 2011), but can be artificially extended to over a year (Rinehart et al., 2013). However, extended overwintering causes chill injury (Rinehart et al., 2013; Torson et al., 2015). In the spring, managers increased storage temperatures to initiate pupation and adult eclosion. If weather delays peak floral bloom, managers slow pupal development using low temperatures. Low temperature exposure during pupation also leads to chill injury (Rinehart et al., 2011). Diapausing prepupae and developing pupae are very different life stages both physiologically and molecularly (MacRae, 2010). Interestingly, exposure to FTR during low temperature exposure increases survival in both contexts (Rinehart et al., 2011, 2013), and reduces the sub-lethal effects of chill injury (Bennett et al., 2015).

Gene expression comparisons could reveal whether the response to FTR is conserved across different developmental stages of M. rotundata. An RNA-seq study conducted in post-diapause quiescent prepupal M. rotundata provided support for the physiological evidence in other chill-injured insects and offered insight into mechanisms driving the reparative effect of FTR versus exposure to a constant 6°C (Torson et al., 2015). Individuals exposed to FTR had increased expression of transcripts functioning in the maintenance of ion homeostasis, response to oxidative stress and an increase in metabolic function; all congruent with the hypothesis that FTR allows for repair of damage during the warm phase. These transcriptional profiles could represent a conserved mechanism. Given that FTR has similar protective effects on survival in both overwintering prepupae and pupae exposed to cold during development, we hypothesized that a conserved suite of genes would underlie the response. The goals of our study were to (1) identify differential gene expression in individuals exposed to FTR protocol during pupal development and (2) determine whether expression profiles observed in pupal FTR support a conserved reparative mechanism of FTR across life stages in M. rotundata. Developing pupae were exposed to either a constant, low-temperature stress (static thermal regime, STR) or the same low temperature with a daily, 1-h warm-temperature pulse for 7 days (FTR). RNA was harvested from bees exposed to either treatment after the warm pulse on day seven. We hypothesized that, as a result of a conserved mechanism driving the reparative effect of FTR across taxa, the exposure to fluctuating temperatures during low-temperature stress during this life stage would elicit transcriptional responses similar to prepupal responses observed previously (Torson et al., 2015).

Temperature protocols

Megachile rotundata were purchased from JWM Leafcutter (Nampa, ID, USA), and were of Canadian origin. Diapausing prepupae arrived as loose cells and were housed in Percival model I-30BLL reach-in incubators at 6±0.5°C under darkness until initiation of the experiment. To end overwinting and initiate development, prepupae were transferred to a Percival model I-30BLL incubator at 29°C. To assess the progression of pupal development prior to treatment, bees not destined for experimental treatment were dissected out of their brood cells and placed in 24-well plates to function as guide plates. The guide plates were stored in plastic containers with breathable lids and a saturated NaCl solution was used to maintain humidity at 75% (as described in Winston and Bates, 1960). Experimental bees were maintained at 29°C until ∼50% of the bees in the guide plates reached a developmental stage characterized by melanization of the eyes, called the ‘red-eye’ stage. At the red-eye stage, bees were grouped into one of two treatments: FTR or STR (Fig. 1). Bees reared under the FTR protocol were exposed to 6°C with a daily warm pulse of 20°C consisting of a 1 h ramp up to 20°C (0.23°C min−1), a 1 h incubation at 20°C, and a 1 h ramp down back to 6°C for each of the 7 days of treatment. Those reared under the STR treatment were maintained at 6°C for the duration of the 7-day treatment.

Fig. 1.

Temperature regimes and RNA collection time points.Megachile rotundata were reared at 29°C until ‘red-eye’ pigmentation stage during pupal development. Individuals were exposed to either a static thermal regime (STR) of 6°C or a fluctuating thermal regime (FTR) of 6°C with a 1-h pulse to 20°C with 1-h ramp-up and ramp-down times. RNA was harvested at T7 for the RNA-seq and T0T7 for the qPCR time series.

Fig. 1.

Temperature regimes and RNA collection time points.Megachile rotundata were reared at 29°C until ‘red-eye’ pigmentation stage during pupal development. Individuals were exposed to either a static thermal regime (STR) of 6°C or a fluctuating thermal regime (FTR) of 6°C with a 1-h pulse to 20°C with 1-h ramp-up and ramp-down times. RNA was harvested at T7 for the RNA-seq and T0T7 for the qPCR time series.

Bees used in the experiments were collected over two field seasons. The RNA-seq and validation was conducted on bees derived from the 2013 summer field season, which arrived in the laboratory in late March 2014. The qPCR time series was conducted on bees derived from the 2014 summer field season, which arrived in the laboratory in late March 2015. Before arriving at the laboratory, prepupal bees were overwintered at a constant temperature (4–6°C) by the supplier (JWM Leafcutter).

RNA-seq

RNA collection

Total RNA was collected from whole M. rotundata pupae exposed to either the FTR or STR treatment (three biological replicates per treatment; one biological replicate is one individual). To ensure the assessment of response to low temperature exposure and not a warm-temperature exposure, individuals from both treatments were harvested for RNA during the cold phase (6°C; Fig. 1) and maintained on ice until the animals were harvested. Animals were harvested during the cold phase (after FTR warm pulse) to ensure a measurement of low temperature response in both treatments. Samples were harvested shortly after (∼1 h) the warm pulse in the event that the effect was transient. Pupae were transferred directly from ice and homogenized in 500 μl of TRIzol using a polypropylene pestle in a nuclease-free 1.5 ml microcentrifuge tube. The Invitrogen TRIzol protocol (Carlsbad, CA, USA) was followed for RNA extraction. Final product was stored as an ethanol precipitate at −80°C until needed.

RNA-seq library preparation and sequencing

Prior to sequencing, RNA pellets were dissolved in DEPC-treated H2O and shipped at a concentration of no less than 20 ng µl−1 (∼10 µg total RNA) on dry ice overnight to the University of Georgia. Upon arrival to the University of Georgia Genomics Facility (Athens, GA, USA), RNA samples from both FTR and STR treatments were checked for quality using an Aligent Bioanalyzer (Santa Clara, CA, USA). A total of 4 μg of total RNA was used for input in the Illumina Truseq Stranded mRNA kit. Samples were pooled and run on a NextSeq500 using a Version 1 Mid-Output 300 cycle kit with PE150 settings. Demultiplexing and adapter trimming were done by default when data were processed in BaseSpace (Illumina, San Diego, CA, USA).

Differential expression analysis

Raw sequence data (BioProject accession: PRJNA352846) were quality checked using FastQC (version 0.10.1 http://www.bioinformatics.babraham.ac.uk/projects/fastqc/; Babraham Bioinformatics). TopHat (version 2.0.9; Trapnell et al., 2012) was used to align the raw reads against the M. rotundata reference genome (BioProject accession: PRJNA66515). Assembly of mapped reads was conducted using Cufflinks (version 2.0.2; Trapnell et al., 2012) and differential expression analysis was carried out using Cuffdiff (version 2.0.2; Trapnell et al., 2012) via the iPlant Collaborative Discovery Environment (Stanzione, 2011). Transcripts were considered significantly differentially expressed when the q-value (Benjamini–Hochberg-corrected P-value) was below the false discovery rate (FDR) of 0.05. The R package cummeRbund (Trapnell et al., 2012) and the Java-based functional annotation program Blast2GO (Conesa et al., 2005) were used for downstream analysis and the generation of differential expression figures.

Life-stage-specific RNA-seq comparison

To assess whether significant overlap between gene expression profiles exist during low-temperature stress between post-diapause quiescent prepupae (Torson et al., 2015) and developing pupae, Venn diagrams were constructed using the R package VennDiagram (Chen and Boutros, 2011). Statistical significance of the overlap of differentially expressed transcripts between data sets was determined via hypergeometric distribution using the R function phyper. The total number of transcripts for comparison was 11,195, which is the total number of transcripts in the reference transcriptome for M. rotundata.

RNA-seq validation

RNA collection

Residual RNA from each RNA-seq sample (2013 field season; described above) was aliquoted and stored at −80°C under ethanol and used for the qPCR validation of the computational methods described above. The 10 most upregulated transcripts in response to FTR were selected for validation.

Sample preparation

RNA samples were diluted to a concentration of 0.333 μg μl−1 and subsequently treated with Invitrogen DNAase I to remove genomic DNA contaminants. A total of 1 μg RNA for each of the three biological replicates per treatment was used as template for first-strand cDNA synthesis using the Invitrogen Super Script III first strand synthesis system for RT-PCR.

qPCR

Real-time quantitative PCR reactions were conducted on a Roche LightCycler 480 II (Indianapolis, ID, USA) using Roche FastStart Universal SYBR Green I Master Mix with ROX. Primers for all targets (10 most upregulated transcripts from RNA-seq experiment; Table S1) and eight reference gene candidates were designed using the Integrated DNA Technologies IDR program (Coralville, IA, USA). Candidate reference genes were selected using this RNA-seq data set and libraries generated in Torson et al. (2015) (18 libraries in total) to identify stably expressed transcripts. Expression of differentially expressed transcripts in qPCR runs were normalized against two most stably expressed reference targets identified using the algorithm geNorm (Vandesompele et al., 2002).

To ensure quality control between and within plates, three controls were used on each plate: a ‘no template control’ consisting of all enzymatic components without cDNA template, a ‘negative RT control’ lacking reverse transcriptase, and a ‘positive RT control’ using a control cDNA template and control primers for each plate. The positive RT control also served as calibrator reactions in the event that comparisons spanned multiple plates.

Biogazelle qbase+ (Ghent, Belgium) and Prism GraphPad (La Jolla, CA, USA) were used for statistical analysis and graphical representations of the data, respectively. Statistical significance for each of the 10 upregulated and eight non-differentially expressed transcripts was determined by t-test of the calibrated normalized relative quantity (CNRQ) of each transcript from qbase+.

qPCR time series

RNA collection

Megachile rotundata from the 2014 field season were reared at 29°C until eye pigmentation stage and then split into either the FTR or STR protocol described above. Total RNA for the qPCR time series was collected from whole M. rotundata pupae exposed to either the FTR or STR treatment (six biological replicates; one individual per replicate). Before exposure to treatment, pupae were harvested to serve as a baseline for gene expression (T0; Fig. 1). After placement into treatment, individuals were harvested during the cold phase (6°C, shortly after the warm pulse in FTR) for each of the 7 days of treatment (T1T7; Fig. 1). Animals were placed on ice during dissection from brood cells until they were homogenized for RNA extraction to prevent an additional warm-temperature exposure. The Invitrogen TRIzol protocol was followed for RNA extraction. Final product was stored as an ethanol precipitate at −80°C until needed. Sample preparation and qPCR protocols were similar to RNA-seq validation (see previous section). Two-way ANOVA with Bonferroni's correction for multiple comparisons determined statistical significance of differential expression, and relationships across time, within treatments were determined using linear regression.

Assembly

Three biological replicates were harvested from each treatment, totaling six individual transcriptomic libraries (BioProject accession: PRJNA352846). Each library averaged a total of 47,279,330 reads (Table 1). Raw reads were mapped to the M. rotundata reference genome (BioProject accession: PRJNA66515) with an average alignment rate of 69.11%. FastQC (version 0.10.1) was used to determine baseline quality of raw reads, and percentage of reads mapped to the reference genome and mate pair concordance (Table 1) were used as a metric of assembly quality.

Table 1.

Each paired-end read RNA-seq library was mapped against the Megachile rotundata genome

Each paired-end read RNA-seq library was mapped against the Megachile rotundata genome
Each paired-end read RNA-seq library was mapped against the Megachile rotundata genome

Differential expression analysis and validation

Using the Tuxedo protocol (Trapnell et al., 2012), expression analysis revealed a total of 434 differentially expressed transcripts between STR and FTR treatments after 7 days of exposure to treatment (307 downregulated and 127 upregulated in FTR). Real-time quantitative PCR confirmed the expression profiles of eight transcripts upregulated in FTR, relative to STR, and six non-differentially expressed transcripts that we used as reference gene candidates (Table S2). Differential expression was confirmed in all cases with at least one primer pair. We were able to confirm expression of eight of the top 10 FTR-upregulated transcripts (Tables S2 and S3) and six of eight non-differentially expressed transcripts (Table S3).

Gene ontology

A Benjamini–Hochberg-corrected Fisher's exact test (one-tailed, FDR=0.05) of upregulated and downregulated transcripts in response to FTR revealed no overrepresented or underrepresented functional classes relative to the gene ontology (GO) distribution in the reference transcriptome. Owing to the lack of enrichment for GO terms, relative to the reference transcriptome, we used direct counts of GO terms of the differentially expressed transcripts in this study for downstream analyses and discussions. Direct counts of the transcripts differentially expressed between FTR and STR treatments on T7 of treatment revealed classes of transcripts functioning in biological processes including regulation of transcription, oxidation–reduction, signal transduction and various metabolic processes (Fig. S1). Also, a direct count of biological process GO terms (Fig. 2) revealed the upregulation in FTR of transcripts functioning in signal transduction and trans-membrane transport. Furthermore, GO terms associated with the construction of cellular components, including plasma membrane and various membrane components, were upregulated in response to FTR (Fig. 2).

Fig. 2.

Gene ontology (GO) of transcripts upregulated in FTR. The most prominent GO terms for biological processes, molecular function and cellular components of transcripts upregulated in M. rotundata exposed to FTR for 7 days.

Fig. 2.

Gene ontology (GO) of transcripts upregulated in FTR. The most prominent GO terms for biological processes, molecular function and cellular components of transcripts upregulated in M. rotundata exposed to FTR for 7 days.

Life-stage-specific comparisons

When the identities of differentially expressed transcripts observed in this study are compared with those presented during the two time points assessed during post-diapause quiescence (Torson et al., 2015), 30 (6.9%) differentially expressed transcripts in the pupal response (434 transcripts) were also differentially expressed during the early time point (215 transcripts), before mortality diverged in prepupal STR/FTR comparison. This overlap is more than we would expect by chance, based on hypergeometric distribution (P=1.27E-10). Seventeen of the 434 differentially expressed transcripts (3.9%) overlapped with the late time point (after mortality diverged, 64 transcripts) in prepupae, also more than expected by chance (P=1.68E–11). Six (1.4%) transcripts differentially expressed in pupae were differentially expressed in both prepupal time points (Fig. 3A).

Fig. 3.

Comparison of differential expression profiles of pupae and post-diapause quiescent prepupae exposed to FTR treatment. The differential expression profiles in this study (interrupted development) were compared with two time points sampled during long-term FTR treatment during the prepupal life stage (Torson et al., 2015). (A) Comparison of differentially expressed transcripts; (B) transcripts upregulated in FTR; (C) transcripts downregulated in FTR; (D) graphical representation of the direction of regulation of the six shared transcripts from A.

Fig. 3.

Comparison of differential expression profiles of pupae and post-diapause quiescent prepupae exposed to FTR treatment. The differential expression profiles in this study (interrupted development) were compared with two time points sampled during long-term FTR treatment during the prepupal life stage (Torson et al., 2015). (A) Comparison of differentially expressed transcripts; (B) transcripts upregulated in FTR; (C) transcripts downregulated in FTR; (D) graphical representation of the direction of regulation of the six shared transcripts from A.

While differentially expressed transcripts overlap more than expected by chance, when direction of regulation (whether the transcript is upregulated or downregulated in FTR-exposed individuals) is considered, all overlap among the three differential expression profiles disappears (Fig. 3B,C). In fact, the only differential expression patterns conserved between the two life stages are regulated in opposite directions (Fig. 3D).

Exposure to daily fluctuations in temperature (FTR) during chilling can mitigate the deleterious physiological effects associated with chill injury (Colinet et al., 2015). FTR appears to provide a protective benefit across all taxa studied so far, but we lack evidence for a conserved mechanism driving this protective effect. In this study, we exposed developing M. rotundata pupae to an STR, known to lead to chill injury, or to an FTR, which diminishes the negative effects of chilling (Bennett et al., 2015). An RNA-seq analysis after 7 days of treatment revealed 434 transcripts that were differentially expressed between the two low-temperature treatments. These differentially expressed transcripts provide evidence for cellular damage, oxidative stress and repair of damage (Table 2; Table S4). Additionally, we observed a lack of conservation in differential gene expression in response to FTR exposure across life stages in M. rotundata, suggesting a robust physiological response driven by variable transcriptional mechanisms.

Table 2.

Candidate genes involved in cell membrane structure and function

Candidate genes involved in cell membrane structure and function
Candidate genes involved in cell membrane structure and function

Membrane imbalance as a driving factor of chill injury

Plants, microbes, and other animals change the physical composition of cell membranes in response to changes in ambient temperature, namely the proportion of unsaturated fatty acids in the phospholipid bilayer (Cossins, 1983). In insects, direct evidence of shifts in membrane phase transitions have yet to be observed (Hayward et al., 2014; Lee, 2010), but strong correlative evidence has been provided (Boardman et al., 2011; Kostál et al., 2006, 2007; Macmillan et al., 2012). Chill injury in insects has been associated with the collapse of ion gradients presumably caused by damage at the cellular level and possibly resulting in decreased activity of membrane-bound enzymes (Hazel, 1995; MacMillan and Sinclair, 2011). Extracellular increases of potassium and diminished levels of sodium and magnesium have been observed in chill-injured tropical cockroaches (Kostál et al., 2004, 2006, 2007). These changes in hemolymph (extracellular) ion concentrations could explain the loss of neuromuscular coordination observed in chill-injured insects (Andersen et al., 2016; MacMillan and Sinclair, 2011) and may be a direct effect of damage to cell membranes. In this study, transcripts such as sodium potassium calcium exchanger 6, a mitochondrial potassium-dependent sodium–calcium antiporter (Table 2), are upregulated in individuals exposed to the STR protocol. Mitochondrial function decreases in Drosophila melanogaster exposed to constant low temperature (Colinet et al., 2017). Proper mitochondrial function would be heavily dependent on membrane composition, so the differential expression of this sodium–potassium–calcium exchanger suggests that STR-exposed individuals may be attempting to restore proper ion gradients in response to cellular damage. Additionally, the upregulation of transcripts with GO terms associated with the construction of various cellular components, including cell membranes and integral membrane components (Fig. 2, Table 2), in response to FTR exposure indicates that one of the FTR's mechanisms of action is to repair the membrane damage that has accumulated over the duration of low-temperature stress.

The upregulation of aquaporins, water-selective transmembrane channels, has been observed in freeze-tolerant insects (Philip et al., 2011), and the blockage of those channels has been shown to decrease freeze tolerance owing to an inability of the cells to move water out of the cell; subsequently causing damage as a result of osmotic pressure (Philip et al., 2008). We observed the upregulation of aquaporins during constant low-temperature stress (STR), suggesting that these individuals may be attempting to return to proper osmotic conditions. Furthermore, STR-exposed M. rotundata showed increased expression of two isoforms coding for the facilitated trehalose transporter TRET1. This transmembrane protein mediates the movement of trehalose synthesized in the fat body bidirectionally across the cell membrane (Kanamori et al., 2010). The presence of two isoforms for Tret-1 is consistent with orthologs in D. melanogaster (Zhao et al., 2000). The enrichment of Tret-1 in STR-reared individuals not only indicates changes in membrane composition, but suggests that these individuals may be transporting trehalose into the hemolymph because of its known role as a cryoprotectant (Bale and Hayward, 2010). Interestingly, transcripts coding for a trehalase precursor, an enzyme needed to convert trehalose into glucose (Shukla et al., 2015) and a facilitated glucose transporter (Table 2) were upregulated in FTR-reared individuals, indicating a shift in hemolymph composition between the two treatments.

In addition to changes in structural components within membranes, such as aquaporins and trehalose transporters, an increase in unsaturated fatty acids has been observed in response to decreased temperatures across taxa (Cossins, 1983). In our dataset, three different transcripts coding for acyl-delta desaturases, enzymes functioning in the transition between saturated and unsaturated fatty acids, were upregulated in STR-reared individuals (Table 2). The upregulation of transcripts coding for desaturase enzymes has been associated with cold hardiness in the onion maggot, Delia antiqua (Kayukawa et al., 2007), indicating this response is conserved. Increased expression of these transcripts suggests that either bees exposed to STR are trying to compensate for a decrease in temperature or FTR-reared individuals may have less of a strain on their cell membranes because of the temporary increases in temperature.

Antioxidant system activation and repair of damage in the FTR response

Oxidative stress has been implicated in the occurrence of chill injury in many insect species (Joanisse and Storey, 1996, 1998; Lalouette et al., 2011; Rojas and Leopold, 1996; Torson et al., 2015) and has been hypothesized to be a downstream consequence of membrane phase transitions (Lee, 2010). Quiescent post-wandering D. melanogaster larvae exposed to extended low-temperature regimes also show enrichment in transcripts that suggest redox imbalance and possibly oxidative stress (Koštál et al., 2016).

Gene expression indicates that STR-exposed individuals are repairing damage, possibly in response to oxidative stress. The upregulation of a gene coding for a myb domain-containing protein in STR-reared individuals (Table S4) implies that M. rotundata is indeed experiencing a stress response. In plant species, this class of transcription factors is known to respond to abiotic stressors, including decreased temperature. The functional unit of the translated proteins is highly conserved, implying a similar function in insects (Ambawat et al., 2013). Furthermore, increased expression of oxidase-peroxidase (Table S4) is indicative of STR-reared individuals experiencing increased levels of peroxide, a consequence of malfunction in mitochondrial respiration (Prasad et al., 1994). If antioxidant systems cannot combat the accumulation of free radicals, such as peroxide, damage to cellular structures would be likely.

The upregulation of transcripts in STR-reared bees such as myofilin isoform b, sestrin-like and DNA damage-binding protein 1 (Table S4) suggests that STR-reared individuals experience DNA damage, potentially caused by increased levels of reactive oxygen species in the bees during chilling. In Drosophila, expression levels of sestrins increase in response to target of rapamycin (TOR) under increased levels of reactive oxygen species (Lee et al., 2010a). The increased expression of these sestrins has also been linked to genotoxic stress including DNA damage (Hay, 2008; Lee et al., 2010b). While sestrins, a class of highly conserved cytoplasmic proteins, have not been shown to repair DNA damage directly, they may be part of a larger signaling cascade. Direct assessments of DNA damage will be required to assess whether oxidative stress may be playing a significant role in the deleterious effects of chill injury or whether high levels of antioxidant expression in FTR-reared individuals is just an artifact of increased metabolic rate owing to transient increases in ambient temperature.

Conservation of mechanisms across life stages

In addition to the benefits of FTR observed during pupal development in M. rotundata, significant increases in survival have been observed when prepupal bees in post-diapause quiescence are exposed to a similar treatment (Rinehart et al., 2013). For the duration of diapause and post-diapause quiescence in both natural and managed populations, M. rotundata are exposed to low temperatures and have mechanisms to cope with these temperatures. Conversely, other life stages, such as developing pupae, which we have assessed in this study, may not be adapted for low-temperature exposure and may have different mechanisms driving low-temperature stress responses, and subsequently the beneficial effects of FTR. Adult bees observed after STR treatment in both life stages exhibit diminished neuromuscular function, behavioral abnormalities and diminished flight capacity (Bennett et al., 2013, 2015). These phenotypes suggest chill injury accumulation in both contexts. Initially, we predicted that the temperature stress during the pupal stage would elicit transcriptional responses similar to the prepupal responses observed in Torson et al. (2015). This hypothesis was not supported. Markedly different differential expression profiles (Fig. 3) between the two life stages indicates that either (1) the mechanism responsible for a conserved FTR benefit vary throughout development or (2) the benefit of FTR exposure changes throughout development.

Even within the pupal life stage, temperature exposure can have variable effects on gene expression depending on context. Survival benefits of FTR, relative to STR, are robust to field season (Bennett et al., 2015; Rinehart et al., 2016), but gene expression changes at a finer scale. Expression profiles of the most differentially expressed transcripts identified in this RNA-seq study (Table S2) were not maintained in M. rotundata pupae from the following field season (2014; Fig. S2). The direction of regulation, relative to STR, was maintained, but expression was highly variable, even with a doubling of the sample size. High variability between individuals suggests that slight differences in developmental timing are introducing variation or season-specific environmental history is influencing gene expression, or both. However, it should be noted that the individuals used for qPCR were exposed to one additional month of overwintering prior to pupal development. This short extension of the overwintering period, compared with the RNA-seq individuals, could have also influenced gene expression during subsequent pupal development. Gene expression in response to FTR exposure now appears to be highly context dependent (i.e. environmental history and developmental stage) and provides further evidence that the conservation of the benefit of FTR benefit may not be at the level of gene expression, but at a higher level of organization. This result provides the basis for our hypothesis that chill injury repair in response to FTR is a robust physiological response that is driven by life-stage-specific mechanisms.

In biological systems, an animal's ability to maintain a certain characteristic response under perturbations in environmental conditions is known as a robust response (Stelling et al., 2004). With similar physiological responses having been observed across taxa, and, thus far, a lack of evidence of a conserved mechanism at the level of gene expression, we hypothesize that the beneficial effects of chill injury repair during FTR exposure represents a robust physiological response driven by variable gene expression. Redundancy is one basic mechanism of a robust system. Redundancy provides alternative molecular strategies to carry out an action (Hartman et al., 2001). For example, glycosylphosphatidylinositol (GPI)-linked glycoproteins in yeast, termed fungal adhesions, have functionally distinct roles separated temporally and spatially (Guo et al., 2000), but members of this protein family can replace one another when inappropriately expressed in the organism. Similar to our results in M. rotundata, work in yeast has shown little overlap in gene expression profiles in response to stressors (Berry et al., 2011), suggesting that variable downstream effectors (i.e. gene expression response) may enact similar physiological responses (i.e. redundancy or genetic buffering).

Between the larval and pupal life stages in M. rotundata, a complete restructuring of the animal takes place, resulting in entirely different physiologies. Redundancy in transcript function, or ‘buffering’, may explain why variable gene expression patterns are observed in what appear to be similar physiological responses to FTR exposure. However, changes in gene or protein expression and metabolite accumulation might be secondary effects rather than direct or adaptive responses to chilling and/or FTR exposure. As a result of this notion, an alternative hypothesis is that the variable gene expression profiles between life stages in M. rotundata is not evidence of different life-stage-specific mechanisms driving the same response at the physiological level, but actually different responses driven by stimuli at a higher level of organization. Additional functional genomic work to identify potentially redundant genes and physiological assays to verify similar physiological responses to FTR exposure will be necessary to provide support for the hypothesis that redundancy allows for conservation of the FTR response across life stages in M. rotundata.

qPCR time series

Expression of a subset of the most significantly upregulated and downregulated transcripts in the RNA-seq study after 7 days of treatment (Table S2) was assessed across the duration of low-temperature stress (Fig. 1) using qPCR. Generally, direction of regulation between treatments at T7 was maintained between this time series (different field season) and the RNA-seq time point, but the qPCR time series had greater variation in expression between biological replicates within treatment and time point, resulting in a lack of statistical significance even after sample sizes were increased from three to six. Expression of the transcript coding for a zinc finger SCAN domain-containing protein had a significant time effect (F7,70=3.752, P=0.0012), showing an increase in expression in individuals exposed to STR (F1,46=17.04, P=0.0002) over time and a significant difference between treatments at T6 (P=0.0183). Expression profiles of Cabut for STR and FTR both decreased over time (STR, F1,46=6.026, P=0.0179; FTR, F1,46=8.779, P=0.0048), with time (F7,70=3.609, P=0.0022) and treatment (F7,70=7.180, P=0.0231) explaining a significant amount of variation in gene expression, but only T2 was significantly different between treatments (P=0.0085; Fig. S2). The serine protease Stubble was upregulated in FTR-reared individuals at T2 (P=0.0196).

The authors thank Marnie Larson of the USDA-ARS, Fargo, ND, USA, for her technical assistance and advice. We also thank Kendra J. Greenlee of North Dakota State University for her comments on the manuscript.

Author contributions

Conceptualization: A.S.T., G.D.Y., J.P.R., J.H.B.; Methodology: A.S.T., S.A.N.; Validation: A.S.T., S.A.N., K.M.K.; Formal analysis: A.S.T.; Investigation: A.S.T.; Resources: A.S.T.; Data curation: A.S.T.; Writing - original draft: A.S.T.; Writing - review & editing: A.S.T., G.D.Y., J.P.R., S.A.N., K.M.K., J.H.B.; Visualization: A.S.T.; Supervision: G.D.Y., J.P.R., J.H.B.; Project administration: G.D.Y., J.P.R.; Funding acquisition: G.D.Y., J.P.R.

Funding

This research was funded by the United States Department of Agriculture.

Data availability

All raw RNA-seq reads were submitted to the NCBI Sequence Read Archives (SRA) under the SRA BioProject accession number PRJNA352846: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA352846/.

Ambawat
,
S.
,
Sharma
,
P.
,
Yadav
,
N. R.
and
Yadav
,
R. C.
(
2013
).
MYB transcription factor genes as regulators for plant responses: an overview
.
Physiol. Mol. Biol. Plants
19
,
307
-
321
.
Andersen
,
M. K.
,
Folkersen
,
R.
,
MacMillan
,
H. A.
and
Overgaard
,
J.
(
2016
).
Cold-acclimation improves chill tolerance in the migratory locust through preservation of ion balance and membrane potential
.
J. Exp. Biol.
220
,
487
-
496
.
Bale
,
J. S.
and
Hayward
,
S. A. L.
(
2010
).
Insect overwintering in a changing climate
.
J. Exp. Biol.
213
,
980
-
994
.
Bennett
,
M. M.
,
Petersen
,
K.
,
Yocum
,
G.
,
Rinehart
,
J.
,
Kemp
,
W.
and
Greenlee
,
K. J.
(
2013
).
Effects of extended prepupal storage duration on adult flight physiology of the alfalfa leafcutting bee (Hymenoptera: Megachilidae)
.
J. Econ. Entomol.
106
,
1089
-
1097
.
Bennett
,
M. M.
,
Cook
,
K. M.
,
Rinehart
,
J. P.
,
Yocum
,
G. D.
,
Kemp
,
W. P.
and
Greenlee
,
K. J.
(
2015
).
Exposure to suboptimal temperatures during metamorphosis reveals a critical developmental window in the solitary bee, Megachile rotundata
.
Physiol. Biochem. Zool.
88
,
508
-
520
.
Berry
,
D. B.
,
Guan
,
Q.
,
Hose
,
J.
,
Haroon
,
S.
,
Gebbia
,
M.
,
Heisler
,
L. E.
,
Nislow
,
C.
,
Giaever
,
G.
and
Gasch
,
A. P.
(
2011
).
Multiple means to the same end: the genetic basis of acquired stress resistance in yeast
.
PLoS Genet.
7
,
e1002353
.
Boardman
,
L.
,
Terblanche
,
J. S.
and
Sinclair
,
B. J.
(
2011
).
Transmembrane ion distribution during recovery from freezing in the woolly bear caterpillar Pyrrharctia isabella (Lepidoptera: Arctiidae)
.
J. Insect Physiol.
57
,
1154
-
1162
.
Boardman
,
L.
,
Sørensen
,
J. G.
and
Terblanche
,
J. S.
(
2013
).
Physiological responses to fluctuating thermal and hydration regimes in the chill susceptible insect, Thaumatotibia leucotreta
.
J. Insect Physiol.
59
,
781
-
794
.
Chen
,
H.
and
Boutros
,
P. C.
(
2011
).
VennDiagram: a package for the generation of highly-customizable Venn and Euler diagrams in R
.
BMC Bioinform.
12
,
35
.
Chen
,
C.-P.
and
Denlinger
,
D. L.
(
1992
).
Reduction of cold injury in flies using an intermittent pulse of high temperature
.
Cryobiology
29
,
138
-
143
.
Colinet
,
H.
,
Renault
,
D.
,
Hance
,
T.
and
Vernon
,
P.
(
2006
).
The impact of fluctuating thermal regimes on the survival of a cold-exposed parasitic wasp, Aphidius colemani
.
Physiol. Entomol.
31
,
234
-
240
.
Colinet
,
H.
,
Hance
,
T.
,
Vernon
,
P.
,
Bouchereau
,
A.
and
Renault
,
D.
(
2007
).
Does fluctuating thermal regime trigger free amino acid production in the parasitic wasp Aphidius colemani (Hymenoptera: Aphidiinae)
.
Comp. Biochem. Physiol. A. Mol. Integr. Physiol.
147
,
484
-
492
.
Colinet
,
H.
,
Lalouette
,
L.
and
Renault
,
D.
(
2011
).
A model for the time–temperature–mortality relationship in the chill-susceptible beetle, Alphitobius diaperinus, exposed to fluctuating thermal regimes
.
J. Therm. Biol.
36
,
403
-
408
.
Colinet
,
H.
,
Sinclair
,
B. J.
,
Vernon
,
P.
and
Renault
,
D.
(
2015
).
Insects in fluctuating thermal environments
.
Annu. Rev. Entomol.
60
,
123
-
140
.
Colinet
,
H.
,
Renault
,
D.
and
Roussel
,
D.
(
2017
).
Cold acclimation allows Drosophila flies to maintain mitochondrial functioning under cold stress
.
Insect Biochem. Mol. Biol.
80
,
52
-
60
.
Conesa
,
A.
,
Götz
,
S.
,
García-Gómez
,
J. M.
,
Terol
,
J.
,
Talón
,
M.
and
Robles
,
M.
(
2005
).
Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research
.
Bioinformatics
21
,
3674
-
3676
.
Cossins
,
A. R.
(
1983
).
Acclimation to environmental change
. In
Cellular Acclimation to Environmental Change
(ed.
A.
Cossins
and
P.
Sheterline
), pp. 3-32.
Cambridge
:
Cambridge University Press
.
Easterling
,
D. R.
,
Meehl
,
G. A.
,
Parmesan
,
C.
,
Changnon
,
S. A.
,
Karl
,
T. R.
and
Mearns
,
L. O.
(
2000
).
Climate extremes: observations, modeling and impacts
.
Science
289
,
2068
-
2074
.
Guo
,
B.
,
Styles
,
C. A.
,
Feng
,
Q.
and
Fink
,
G. R.
(
2000
).
A Saccharomyces gene family involved in invasive growth, cell-cell adhesion, and mating
.
Proc. Natl. Acad. Sci. USA
97
,
12158
-
12163
.
Hansen
,
J.
,
Sato
,
M.
and
Ruedy
,
R.
(
2012
).
Perception of climate change
.
Proc. Natl. Acad. Sci. USA
109
,
E2415
-
E2423
.
Hartman
,
J. L.
, IV
,
Garvik
,
B.
and
Hartwell
,
L.
(
2001
).
Principles for the buffering of genetic variation
.
Science
291
,
1001
-
1004
.
Hay
,
N.
(
2008
).
p53 Strikes mTORC1 by employing sestrins
.
Cell Metab.
8
,
184
-
185
.
Hayward
,
S. A. L.
,
Manso
,
B.
and
Cossins
,
A. R.
(
2014
).
Molecular basis of chill resistance adaptations in poikilothermic animals
.
J. Exp. Biol.
217
,
6
-
15
.
Hazel
,
J. R.
(
1995
).
Thermal adaptation in biological membranes: Is homeoviscous adaptation the explanation?
Annu. Rev. Physiol.
57
,
19
-
42
.
Jing
,
X.-H.
,
Wang
,
X.-H.
and
Kang
,
L.
(
2005
).
Chill injury in the eggs of the migratory locust, Locusta migratoria (Orthoptera: Acrididae): the time-temperature relationship with high-temperature interruption
.
Insect Sci
. 12,
171
-
178
.
Joanisse
,
D.
and
Storey
,
K.
(
1996
).
Oxidative stress and antioxidants in overwintering larvae of cold-hardy goldenrod gall insects
.
J. Exp. Biol.
199
,
1483
-
1491
.
Joanisse
,
D. R.
and
Storey
,
K. B.
(
1998
).
Oxidative stress and antioxidants in stress and recovery of cold-hardy insects
.
Insect Biochem. Mol. Biol.
28
,
23
-
30
.
Kanamori
,
Y.
,
Saito
,
A.
,
Hagiwara-Komoda
,
Y.
,
Tanaka
,
D.
,
Mitsumasu
,
K.
,
Kikuta
,
S.
,
Watanabe
,
M.
,
Cornette
,
R.
,
Kikawada
,
T.
and
Okuda
,
T.
(
2010
).
The trehalose transporter 1 gene sequence is conserved in insects and encodes proteins with different kinetic properties involved in trehalose import into peripheral tissues
.
Insect Biochem. Mol. Biol.
40
,
30
-
37
.
Kayukawa
,
T.
,
Chen
,
B.
,
Hoshizaki
,
S.
and
Ishikawa
,
Y.
(
2007
).
Upregulation of a desaturase is associated with the enhancement of cold hardiness in the onion maggot, Delia antiqua
.
Insect Biochem. Mol. Biol.
37
,
1160
-
1167
.
Kostál
,
V.
(
2006
).
Eco-physiological phases of insect diapause
.
J. Insect Physiol.
52
,
113
-
127
.
Kostál
,
V.
,
Vambera
,
J.
and
Bastl
,
J.
(
2004
).
On the nature of pre-freeze mortality in insects: water balance, ion homeostasis and energy charge in the adults of Pyrrhocoris apterus
.
J. Exp. Biol.
207
,
1509
-
1521
.
Kostál
,
V.
,
Yanagimoto
,
M.
and
Bastl
,
J.
(
2006
).
Chilling-injury and disturbance of ion homeostasis in the coxal muscle of the tropical cockroach (Nauphoeta cinerea)
.
Comp. Biochem. Physiol. B. Biochem. Mol. Biol.
143
,
171
-
179
.
Kostál
,
V.
,
Renault
,
D.
,
Mehrabianová
,
A.
and
Bastl
,
J.
(
2007
).
Insect cold tolerance and repair of chill-injury at fluctuating thermal regimes: role of ion homeostasis
.
Comp. Biochem. Physiol. A. Mol. Integr. Physiol.
147
,
231
-
238
.
Kostál
,
V.
,
Korbelová
,
J.
,
Štětina
,
T.
,
Poupardin
,
R.
,
Colinet
,
H.
,
Zahradníčková
,
H.
,
Opekarová
,
I.
,
Moos
,
M.
and
Šimek
,
P.
(
2016
).
Physiological basis for low-temperature survival and storage of quiescent larvae of the fruit fly Drosophila melanogaster
.
Sci. Rep.
6
,
32346
.
Lalouette
,
L.
,
Williams
,
C. M.
,
Hervant
,
F.
,
Sinclair
,
B. J.
and
Renault
,
D.
(
2011
).
Metabolic rate and oxidative stress in insects exposed to low temperature thermal fluctuations
.
Comp. Biochem. Physiol. A. Mol. Integr. Physiol.
158
,
229
-
234
.
Lee
,
R. E.
(
2010
).
A primer on insect cold-tolerance
. In
Low Temperature Biology of Insects
. (ed.
D. L.
Denlinger
and
R. E.
Lee
), pp. 3-34.
Cambridge, NY
:
Cambridge University Press
.
Lee
,
J. H.
,
Budanov
,
A. V.
,
Park
,
E. J.
,
Birse
,
R.
,
Kim
,
T. E.
,
Perkins
,
G. A.
,
Ocorr
,
K.
,
Ellisman
,
M. H.
,
Bodmer
,
R.
,
Bier
,
E.
, et al. 
(
2010a
).
Sestrin as a feedback inhibitor of TOR that prevents age-related pathologies
.
Science
327
,
1223
-
1228
.
Lee
,
J. H.
,
Bodmer
,
R.
,
Bier
,
E.
,
Karin
,
M.
(
2010b
).
Sestrins at the crossroad between stress and aging
.
Aging
2
,
369
-
374
.
Leopold
,
R. A.
,
Rojas
,
R. R.
and
Atkinson
,
P. W.
(
1998
).
Post pupariation cold storage of three species of flies: increasing chilling tolerance by acclimation and recurrent recovery periods
.
Cryobiology
36
,
213
-
224
.
MacMillan
,
H. A.
and
Sinclair
,
B. J.
(
2011
).
The role of the gut in insect chilling injury: cold-induced disruption of osmoregulation in the fall field cricket, Gryllus pennsylvanicus
.
J. Exp. Biol.
214
,
726
-
734
.
Macmillan
,
H. A.
,
Williams
,
C. M.
,
Staples
,
J. F.
and
Sinclair
,
B. J.
(
2012
).
Reestablishment of ion homeostasis during chill-coma recovery in the cricket Gryllus pennsylvanicus
.
Proc. Natl. Acad. Sci. USA
109
,
20750
-
20755
.
MacRae
,
T. H.
(
2010
).
Gene expression, metabolic regulation and stress tolerance during diapause
.
Cell. Mol. Life Sci.
67
,
2405
-
2424
.
Marshall
,
K. E.
and
Sinclair
,
B. J.
(
2012
).
The impacts of repeated cold exposure on insects
.
J. Exp. Biol.
215
,
1607
-
1613
.
Nedvěd
,
O.
,
Lavy
,
D.
and
Verhoef
,
H. A.
(
1998
).
Modelling the time-temperature relationship in cold injury and effect of high-temperature interruptions on survival in a chill-sensitive collembolan
.
Funct. Ecol.
12
,
816
-
824
.
O'Neill
,
K. M.
,
O'Neill
,
R. P.
,
Kemp
,
W. P.
and
Delphia
,
C. M.
(
2011
).
Effect of temperature on post-wintering development and total lipid content of alfalfa leafcutting bees
.
Environ. Entomol.
40
,
917
-
930
.
Philip
,
B. N.
,
Yi
,
S.-X.
,
Elnitsky
,
M. A.
and
Lee
,
R. E.
(
2008
).
Aquaporins play a role in desiccation and freeze tolerance in larvae of the goldenrod gall fly, Eurosta solidaginis
.
J. Exp. Biol.
211
,
1114
-
1119
.
Philip
,
B. N.
,
Kiss
,
A. J.
and
Lee
,
R. E.
(
2011
).
The protective role of aquaporins in the freeze-tolerant insect Eurosta solidaginis: functional characterization and tissue abundance of EsAQP1
.
J. Exp. Biol.
214
,
848
-
857
.
Pitts-Singer
,
T. L.
,
Cane
,
J. H.
(
2011
).
The alfalfa leafcutting bee, Megachile rotundata: the world's most intensively managed solitary bee
.
Annu. Rev. Entomol.
56
,
221
-
237
. https://doi.org/10.1146/annurev-ento-120709-144836
Prasad
,
T. K.
,
Anderson
,
M. D.
,
Martin
,
B. A.
and
Stewart
,
C. R.
(
1994
).
Evidence for chilling-induced oxidative stress in maize seedlings and a regulatory role for hydrogen peroxide
.
Plant Cell
6
,
65
-
74
.
Renault
,
D.
,
Nedved
,
O.
,
Hervant
,
F.
and
Veron
,
P.
(
2004
).
The importance of fluctuating thermal regimes for repairing chill injuries in the tropical beetle Alphitobius diaperinus (Coleoptera : Tenebrionidae) during exposure to low temperature
.
Physiol. Entomol.
29
,
139
-
145
.
Rinehart
,
J. P.
,
Yocum
,
G. D.
,
West
,
M.
and
Kemp
,
W. P.
(
2011
).
A fluctuating thermal regime improves survival of cold-mediated delayed emergence in developing Megachile rotundata (Hymenoptera: Megachilidae)
.
J. Econ. Entomol.
104
,
1162
-
1166
.
Rinehart
,
J. P.
,
Yocum
,
G. D.
,
Kemp
,
W. P.
and
Greenlee
,
K. J.
(
2013
).
A fluctuating thermal regime improves long-term survival of quiescent prepupal Megachile rotundata (Hymenoptera : Megachilidae)
.
Apic. Soc. Insects
106
,
1081
-
1088
.
Rinehart
,
J. P.
,
Yocum
,
G. D.
,
Kemp
,
W. P.
and
Bowsher
,
J. H.
(
2016
).
Optimizing fluctuating thermal regime storage of developing Megachile rotundata (Hymenoptera: Megachilidae)
.
J. Econ. Entomol.
109
,
993
-
1000
.
Rojas
,
R. R.
and
Leopold
,
R. A.
(
1996
).
Chilling injury in the housefly: evidence for the role of oxidative stress between pupariation and emergence
.
Cryobiology
33
,
447
-
458
.
Shukla
,
E.
,
Thorat
,
L. J.
,
Nath
,
B. B.
and
Gaikwad
,
S. M.
(
2015
).
Insect trehalase: physiological significance and potential applications
.
Glycobiology
25
,
357
-
367
.
Stanzione
,
D.
(
2011
).
The iPlant collaborative: cyberinfrastructure to feed the world
.
Computer (Long. Beach. Calif).
44
,
44
-
52
.
Stelling
,
J.
,
Sauer
,
U.
,
Szallasi
,
Z.
,
Doyle
,
F. J.
and
Doyle
,
J.
(
2004
).
Robustness of cellular functions
.
Cell
118
,
675
-
685
.
Torson
,
A. S.
,
Yocum
,
G. D.
,
Rinehart
,
J. P.
,
Kemp
,
W. P.
and
Bowsher
,
J. H.
(
2015
).
Transcriptional responses to fluctuating thermal regimes underpinning differences in survival in the solitary bee Megachile rotundata
.
J. Exp. Biol.
218
,
1060
-
1068
.
Trapnell
,
C.
,
Roberts
,
A.
,
Goff
,
L.
,
Pertea
,
G.
,
Kim
,
D.
,
Kelley
,
D. R.
,
Pimentel
,
H.
,
Salzberg
,
S. L.
,
Rinn
,
J. L.
and
Pachter
,
L.
(
2012
).
Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks
.
Nat. Protoc.
7
,
562
-
578
.
Vandesompele
,
J.
,
De Preter
,
K.
,
Pattyn
,
F.
,
Poppe
,
B.
,
Van Roy
,
N.
,
De Paepe
,
A.
and
Speleman
,
F.
(
2002
).
Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes
.
Genome Biol.
3
,
RESEARCH0034
.
Vasseur
,
D. A.
,
Delong
,
J. P.
,
Gilbert
,
B.
,
Greig
,
H. S.
,
Harley
,
C. D. G.
,
Mccann
,
K. S.
,
Savage
,
V.
,
Tunney
,
T. D.
,
Connor
,
M. I. O.
,
Delong
,
J. P.
, et al. 
(
2014
).
Increased temperature variation poses a greater risk to species than climate warming Increased temperature variation poses a greater risk to species than climate warming
.
Proc. Biol. Sci.
281
,
20132612
.
Whitfield
,
G. H.
and
Richards
,
K. W.
(
1992
).
Temperature dependent development and survival of immature stages of the alfalfa leafcutter bee, Megachile rotundata (Hymenoptera: Megachilidae)
.
Apidologie
23
,
11
-
23
.
Winston
,
P. W.
and
Bates
,
D. H.
(
1960
).
Saturated solutions for the control of humidity in biological research
.
Ecology
41
,
232
-
237
.
Yocum
,
G. D.
,
Žďárek
,
J.
,
Joplin
,
K. H.
,
Lee
,
R. E.
,
Smith
,
D. C.
,
Manter
,
K. D.
and
Denlinger
,
D. L.
(
1994
).
Alteration of the eclosion rhythm and eclosion behavior in the flesh fly, Sarcophaga crassipalpis, by low and high temperature stress
.
J. Insect Physiol.
40
,
13
-
21
.
Yocum
,
G. D.
,
Kemp
,
W. P.
,
Bosch
,
J.
and
Knoblett
,
J. N.
(
2006
).
Thermal history influences diapause development in the solitary bee Megachile rotundata
.
J. Insect Physiol.
52
,
1113
-
1120
.
Yocum
,
G. D.
,
Rinehart
,
J. P.
and
Kemp
,
W. P.
(
2012
).
Duration and frequency of a high temperature pulse affect survival of emergence-ready Megachile rotundata (Hymenoptera: Megachilidae) during low-temperature incubation
.
J. Econ. Entomol.
105
,
14
-
19
.
Zhao
,
Q.
,
Spradling
,
A. C.
,
Mattei
,
B.
,
An
,
H. J.
,
Weinstock
,
G. M.
,
Holt
,
R. A.
,
Zheng
,
L.
,
Harris
,
M.
,
Ke
,
Z.
,
Zaveri
,
J. S.
, et al. 
(
2000
).
The genome sequence of Drosophila melanogaster
.
Science.
287
,
2185
-
2195
.

Competing interests

The authors declare no competing or financial interests.

Supplementary information