Cold tolerance and cold-induced modulation of gene expression in two Drosophila virilis group species with different distributions

Comparative transcriptome analysis of Liriomyza trifolii (Burgess) and Liriomyza sativae (Blanchard) (Diptera: Agromyzidae) in response to rapid cold hardening

The ability of insets to react efficiently to fluctuation in temperature is crucial for them to survive in variable surroundings. Rapid cold hardening (RCH) is a process that increase cold tolerance in most insect species. The molecular mechanisms of RCH remain largely unknown, and whether it is associated with transcriptional changes is unclear. In this study, we compared the transcriptomes of Liriomyza trifolii and L. sativae exposed to RCH to investigate the transcript abundance due to RCH in both species. RNA-seq revealed 93,166 assembled unigenes, and 34,303 of these were annotated in the L. trifolii and L. sativae transcriptome libraries. After a 4-h treatment at 1°C (RCH) compared with control, 268 and 606 unigenes were differentially expressed in L. trifolii and L. sativae, respectively. When comparing pupae exposed to 2h cold shock directly with pupae went through 4h acclimation prior to 2h cold shock, 60 and 399 unigenes were differentially expressed in L trifolii and L sativae, respectively. Genes that were commonly expressed in both L. trifolii and L. sativae, included cytochrome P450, cuticular protein, glucose dehydrogenase, solute carrier family 22 and cationic amino acid transporter. Additionally, several pathways including galactose metabolism and peroxisome were significantly enriched during RCH. Our results show that the transcriptional response is correlated with RCH in the pupal stage of the two Liriomyza species, but more transcriptional changes were identified in L sativae than in L. trifolii.

Divergence and introgression among the virilis group of Drosophila

Speciation with gene flow is now widely regarded as common. However, the frequency of introgression between recently diverged species and the evolutionary consequences of gene flow are still poorly understood. The virilis group of Drosophila contains 12 species that are geographically widespread and show varying levels of prezygotic and postzygotic isolation. Here, we use de novo genome assemblies and whole-genome sequencing data to resolve phylogenetic relationships and describe patterns of introgression and divergence across the group. We suggest that the virilis group consists of three, rather than the traditional two, subgroups. Some genes undergoing rapid sequence divergence across the group were involved in chemical communication and desiccation tolerance, and may be related to the evolution of sexual isolation and adaptation. We found evidence of pervasive phylogenetic discordance caused by ancient introgression events between distant lineages within the group, and more recent gene flow between closely related species. When assessing patterns of genome-wide divergence in species pairs across the group, we found no consistent genomic evidence of a disproportionate role for the X chromosome as has been found in other systems. Our results show how ancient and recent introgressions confuse phylogenetic reconstruction, but may play an important role during early radiation of a group.

Coping with global warming: Adult thermal thresholds in four pestiferous Anastrepha species determined under experimental laboratory conditions and development/survival times of immatures and adults under natural field conditions

Climate change, particularly global warming, is disturbing biological processes in unexpected ways and forcing us to re-study/reanalyze the effects of varying temperatures, among them extreme ones, on insect functional traits such as lifespan and fecundity/fertility. Here we experimentally tested, under both laboratory and field conditions, the effects of an extreme range of temperatures (5, 10, 15, 20, 30, 40, and 45 °C, and the naturally varying conditions experienced in the field), on survivorship/lifespan, fecundity, and fertility of four pestiferous fruit fly species exhibiting contrasting life histories and belonging to two phylogenetic groups within the genus Anastrepha: A. ludens, A. obliqua, A. striata, and A. serpentina. In the field, we also measured the length of the entire life cycle (egg to adult), and in one species (A. ludens), the effect on the latter of the host plant (mango and grapefruit). Under laboratory conditions, none of the adults, independent of species, could survive a single day when exposed to a constant temperature of 45 °C, but A. striata and A. serpentina females/males survived at the highly contrasting temperatures of 5 and 40 °C at least 7 days. Maximum longevity was achieved in all species at 15 °C (375, 225, 175 and 160 days in A. ludens, A. serpentina, A. striata and A. obliqua females, respectively). Anastrepha ludens layed many eggs until late in life (368 days) at 15 °C, but none eclosed. Eclosion was only observed in all species at 20 and 30 °C. Under natural conditions, flies lived ca. 100 days less than in the laboratory at 15 °C, likely due to the physiological cost of dealing with the highly varying environmental patterns over 24 h (minimum and maximum temperatures and relative humidity of ca. 10–40 °C, and 22–100%, respectively). In the case of A. ludens, the immature’s developmental time was shorter in mango, but adult survival was longer than in grapefruit. We discuss our results considering the physiological processes regulating the traits measured and tie them to the increasing problem of global warming and its hidden effects on the physiology of insects, as well as the ecological and pest management implications.

Transcriptome analysis of the response to low temperature acclimation in Calliptamus italicus eggs

Background

Calliptamus italicus is a dominant species in the desert and semi-desert grassland. It is widely distributed throughout many regions such as Asia, Europe, North Africa and the Mediterranean, and has enormous destructive potential for agriculture and animal husbandry. The C. italicus overwintering as eggs in the soil through diapause, and the cold tolerance of locust eggs is the key to their ability to survive the winter smoothly to maintain the population.

Results

Transcriptome analysis of C. italicus eggs was carried out in this paper in constant low temperature acclimation, natural low temperature acclimation and room temperature. The differentially expressed genes related to cold tolerance were screened out, the differences in expression patterns under different low temperature acclimation were analyzed, and the genes in the significantly up-regulated pathways may play an important role in cold tolerance. The results show that different domestication modes can induce C. italicus eggs to express a large number of genes to alleviate low temperature damage, but C. italicus eggs are more sensitive to changes in temperature. Compared with the control, there are 8689 DEGs at constant low temperature and 14,994 DEGs at natural low temperature. KEGG analysis showed that DEGs were mainly enriched in pathways related to metabolism and biological systems under constant low temperature, and were mainly enriched in pathways related to biological systems and environmental information processing under natural low temperature. In addition, RNAi technology was used to further verify the regulation of genes in the significantly enriched up-regulated pathways on C. italicus eggs, and it was confirmed that the hatching rate of C. italicus eggs at low temperature was significantly reduced after interference.

Conclusions

Transcriptome analysis of C. italicus eggs treated at different temperatures provided a theoretical basis for further understanding the adaptation mechanism of C. italicus eggs to low temperature. In addition, four potential RNAi target genes were verified in the eggs of C. italicus for the first time, providing new ideas for effective control of this species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08705-3.

Metabolomic and transcriptomic responses of ticks during recovery from cold shock reveal mechanisms of survival

Ticks are blood-feeding ectoparasites but spend most of their life off-host where they may have to tolerate low winter temperatures. Rapid cold-hardening (RCH) is a process commonly used by arthropods, including ticks, to improve survival of acute low temperature exposure. However, little is known about the underlying mechanisms in ticks associated with RCH, cold shock, and recovery from these stresses. In the present study, we investigated the extent to which RCH influences gene expression and metabolism during recovery from cold stress in Dermacentor variabilis, the American dog tick, using a combined transcriptomics and metabolomics approach. Following recovery from RCH, 1,860 genes were differentially expressed in ticks, whereas only 99 genes responded during recovery to direct cold shock. Recovery from RCH resulted in an upregulation of various pathways associated with ion binding, transport, metabolism, and cellular structures seen in the response of other arthropods to cold. The accumulation of various metabolites, including several amino acids and betaine, corresponded to transcriptional shifts in the pathways associated with these molecules, suggesting congruent metabolome and transcriptome changes. Ticks receiving exogenous betaine and valine demonstrated enhanced cold tolerance, suggesting cryoprotective effects of these metabolites. Overall, many of the responses during recovery from cold shock in ticks were similar to those observed in other arthropods, but several adjustments may be distinct from other currently examined taxa.

Cold adaptation drives population genomic divergence in the ecological specialist, Drosophila montana

Detecting signatures of ecological adaptation in comparative genomics is challenging, but analysing population samples with characterised geographic distributions, such as clinal variation, can help identify genes showing covariation with important ecological variation. Here, we analysed patterns of geographic variation in the cold-adapted species Drosophila montana across phenotypes, genotypes and environmental conditions and tested for signatures of cold adaptation in population genomic divergence. We first derived the climatic variables associated with the geographic distribution of 24 populations across two continents to trace the scale of environmental variation experienced by the species, and measured variation in the cold tolerance of the flies of six populations from different geographic contexts. We then performed pooled whole genome sequencing of these six populations, and used Bayesian methods to identify SNPs where genetic differentiation is associated with both climatic variables and the population phenotypic measurements, while controlling for effects of demography and population structure. The top candidate SNPs were enriched on the X and fourth chromosomes, and they also lay near genes implicated in other studies of cold tolerance and population divergence in this species and its close relatives. We conclude that ecological adaptation has contributed to the divergence of D. montana populations throughout the genome and in particular on the X and fourth chromosomes, which also showed highest interpopulation F_ST . This study demonstrates that ecological selection can drive genomic divergence at different scales, from candidate genes to chromosome-wide effects.

Multiple paths to cold tolerance: the role of environmental cues, morphological traits and the circadian clock gene vrille

Background

Tracing the association between insect cold tolerance and latitudinally and locally varying environmental conditions, as well as key morphological traits and molecular mechanisms, is essential for understanding the processes involved in adaptation. We explored these issues in two closely-related species, Drosophila montana and Drosophila flavomontana, originating from diverse climatic locations across several latitudes on the coastal and mountainous regions of North America. We also investigated the association between sequence variation in one of the key circadian clock genes, vrille, and cold tolerance in both species. Finally, we studied the impact of vrille on fly cold tolerance and cold acclimation ability by silencing it with RNA interference in D. montana.

Results

We performed a principal component analysis (PCA) on variables representing bioclimatic conditions on the study sites and used latitude as a proxy of photoperiod. PC1 separated the mountainous continental sites from the coastal ones based on temperature variability and precipitation, while PC2 arranged the sites based on summer and annual mean temperatures. Cold tolerance tests showed D. montana to be more cold-tolerant than D. flavomontana and chill coma resistance (CT_min) of this species showed an association with PC2. Chill coma recovery time (CCRT) of both species improved towards northern latitudes, and in D. flavomontana this trait was also associated with PC1. D. flavomontana flies were darkest in the coast and in the northern mountainous populations, but coloration showed no linkage with cold tolerance. Body size decreased towards cold environments in both species, but only within D. montana populations largest flies showed fastest recovery from cold. Finally, both the sequence analysis and RNAi study on vrille suggested this gene to play an essential role in D. montana cold resistance and acclimation, but not in recovery time.

Conclusions

Our study demonstrates the complexity of insect cold tolerance and emphasizes the need to trace its association with multiple environmental variables and morphological traits to identify potential agents of natural selection. It also shows that a circadian clock gene vrille is essential both for short- and long-term cold acclimation, potentially elucidating the connection between circadian clock system and cold tolerance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01849-y.

Sex-specific responses to cold in a very cold-tolerant, northern Drosophila species

Organisms can plastically alter resource allocation in response to changing environmental factors. For example, in harsh conditions, organisms are expected to shift investment from reproduction toward survival; however, the factors and mechanisms that govern the magnitude of such shifts are relatively poorly studied. Here we compared the impact of cold on males and females of the highly cold-tolerant species Drosophila montana at the phenotypic and transcriptomic levels. Although both sexes showed similar changes in cold tolerance and gene expression in response to cold treatment, indicating that the majority of changes are concordant between the sexes, we identified a clear reduction in sexually dimorphic gene expression, suggesting that preparing for the colder season involves reducing investment in sex-specific traits. This reduction was larger in males than females, as expected if male sexual traits are more condition-dependent than female traits, as predicted by theory. Gene expression changes were primarily associated with shifts in metabolic profile, which likely play a role in increasing cold tolerance. Finally, we found that the expression of immune genes was reduced following cold treatment, suggesting that reduced investment in costly immune function may be important in helping flies survive colder periods.

Expression analysis of genes related to cold tolerance in Dendroctonus valens

Pine beetles are well known in North America for their widespread devastation of pine forests. However, Dendroctonus valens LeConte is an important invasive forest pest in China also. Adults and larvae of this bark beetle mainly winter at the trunks and roots of Pinus tabuliformis and Pinus sylvestris; larvae, in particular, result in pine weakness or even death. Since the species was introduced from the United States to Shanxi in 1998, its distribution has spread northward. In 2017, it invaded a large area at the junction of Liaoning, Inner Mongolia and Hebei provinces, showing strong cold tolerance. To identify genes relevant to cold tolerance and the process of overwintering, we sequenced the transcriptomes of wintering and non-wintering adult and larval D. valens using the Illumina HiSeq platform. Differential expression analysis methods for other non-model organisms were used to compare transcript abundances in adults and larvae at two time periods, followed by the identification of functions and metabolic pathways related to genes associated with cold tolerance. We detected 4,387 and 6,091 differentially expressed genes (DEGs) between sampling dates in larvae and adults, respectively, and 1,140 common DEGs, including genes encoding protein phosphatase, very long-chain fatty acids protein, cytochrome P450, and putative leucine-rich repeat-containing proteins. In a Gene Ontology (GO) enrichment analysis, 1,140 genes were assigned to 44 terms, with significant enrichment for cellulase activity, hydrolase activity, and carbohydrate metabolism. Kyoto Encyclopedia of Genes and Genomes (KEGG) classification and enrichment analyses showed that the lysosomal and purine metabolism pathways involved the most DEGs, the highly enriched terms included autophagy-animal, pentose and glucuronate interconversions and lysosomal processes. We identified 140 candidate genes associated with cold tolerance, including genes with established roles in this trait (e.g., genes encoding trehalose transporter, fructose-1,6-bisphosphatase, and trehalase). Our comparative transcriptome analysis of adult and larval D. valens in different conditions provides basic data for the discovery of key genes and molecular mechanisms underlying cold tolerance.

Silencing the Myosin Regulatory Light Chain Gene sqh Reduces Cold Hardiness in Ophraella communa LeSage (Coleoptera: Chrysomelidae)

Ambrosia artemisiifolia is a noxious invasive alien weed, that is harmful to the environment and human health. Ophraella communa is a biocontrol agent for A. artemisiifolia, that was accidentally introduced to the Chinese mainland and has now spread throughout southern China. Recently, we found that upon artificial introduction, O. communa can survive in northern China as well. Therefore, it is necessary to study the cold hardiness of O. communa. Many genes have been identified to play a role in cold-tolerance regulation in insects, but the function of the gene encoding non-muscle myosin regulatory light chain (MRLC-sqh) remains unknown. To evaluate the role played by MRLC-sqh in the cold-tolerance response, we cloned and characterized MRLC-sqh from O. communa. Quantitative real-time PCR revealed that MRLC-sqh was expressed at high levels in the gut and pupae of O. communa. The expression of MRLC-sqh was shown to decrease after cold shock between 10 and 0 °C and ascend between 0 and -10 °C, but these did not show a positive association between MRLC-sqh expression and cold stress. Silencing of MRLC-sqh using dsMRLC-sqh increased the chill-coma recovery time of these beetles, suggesting that cold hardiness was reduced in its absence. These results suggest that the cold hardiness of O. communa may be partly regulated by MRLC-sqh. Our findings highlight the importance of motor proteins in mediating the cold response in insects.

Small heat shock protein Hsp67Bc plays a significant role in Drosophila melanogaster cold stress tolerance

Hsp67Bc in Drosophila melanogaster is a member of the small heat shock protein family, the main function of which is to prevent the aggregation of misfolded or damaged proteins. Hsp67Bc interacts with Starvin and Hsp23, which are known to be a part of the cold stress response in the fly during the recovery phase. In this study, we investigated the role of the Hsp67Bc gene in the cold stress response. We showed that in adult Drosophila, Hsp67Bc expression increases after cold stress and decreases after 1.5 h of recovery, indicating the involvement of Hsp67Bc in short-term stress recovery. We also implemented a deletion in the D. melanogaster Hsp67Bc gene using imprecise excision of a P-element, and analysed the cold tolerance of Hsp67Bc-null mutants at different developmental stages. We found that Hsp67Bc-null homozygous flies are viable and fertile but display varying cold stress tolerance throughout the stages of ontogenesis: the survival after cold stress is slightly impaired in late third instar larvae, unaffected in pupae, and notably affected in adult females. Moreover, the recovery from chill coma is delayed in Hsp67Bc-null adults of both sexes. In addition, the deletion in the Hsp67Bc gene caused more prominent up-regulation of Hsp70 following cold stress, suggesting the involvement of Hsp70 in compensation of the lack of the Hsp67Bc protein. Taken together, our results suggest that Hsp67Bc is involved in the recovery of flies from a comatose state and contributes to the protection of the fruit fly from cold stress.

Role of the Insect Neuroendocrine System in the Response to Cold Stress

Insects are the largest group of animals. They are capable of surviving in virtually all environments from arid deserts to the freezing permafrost of polar regions. This success is due to their great capacity to tolerate a range of environmental stresses, such as low temperature. Cold/freezing stress affects many physiological processes in insects, causing changes in main metabolic pathways, cellular dehydration, loss of neuromuscular function, and imbalance in water and ion homeostasis. The neuroendocrine system and its related signaling mediators, such as neuropeptides and biogenic amines, play central roles in the regulation of the various physiological and behavioral processes of insects and hence can also potentially impact thermal tolerance. In response to cold stress, various chemical signals are released either via direct intercellular contact or systemically. These are signals which regulate osmoregulation - capability peptides (CAPA), inotocin (ITC)-like peptides, ion transport peptide (ITP), diuretic hormones and calcitonin (CAL), substances related to the general response to various stress factors - tachykinin-related peptides (TRPs) or peptides responsible for the mobilization of body reserves. All these processes are potentially important in cold tolerance mechanisms. This review summarizes the current knowledge on the involvement of the neuroendocrine system in the cold stress response and the possible contributions of various signaling molecules in this process.

Rapid cold hardening: ecological relevance, physiological mechanisms and new perspectives

Rapid cold hardening (RCH) is a type of phenotypic plasticity that allows ectotherms to quickly enhance cold tolerance in response to brief chilling (lasting minutes to hours). In this Review, we summarize the current state of knowledge of this important phenotype and provide new directions for research. As one of the fastest adaptive responses to temperature known, RCH allows ectotherms to cope with sudden cold snaps and to optimize their performance during diurnal cooling cycles. RCH and similar phenotypes have been observed across a diversity of ectotherms, including crustaceans, terrestrial arthropods, amphibians, reptiles, and fish. In addition to its well-defined role in enhancing survival to extreme cold, RCH also protects against nonlethal cold injury by preserving essential functions following cold stress, such as locomotion, reproduction, and energy balance. The capacity for RCH varies across species and across genotypes of the same species, indicating that RCH can be shaped by selection and is likely favored in thermally variable environments. Mechanistically, RCH is distinct from other rapid stress responses in that it typically does not involve synthesis of new gene products; rather, the existing cellular machinery regulates RCH through post-translational signaling mechanisms. However, the protective mechanisms that enhance cold hardiness are largely unknown. We provide evidence that RCH can be induced by multiple triggers in addition to low temperature, and that rapidly induced tolerance and cross-tolerance to a variety of environmental stressors may be a general feature of stress responses that requires further investigation.

Evidence for a rapid cold hardening response in cultured Drosophila S2 cells

The ability to quickly respond to changes in environmental temperature is critical for organisms living in thermally variable environments. To cope with sudden drops in temperature, insects and other ectotherms are capable of rapid cold hardening (RCH), in which mild chilling significantly enhances cold tolerance within minutes. While the ecological significance of RCH is well established, the mechanisms underlying RCH are still poorly understood. Previous work has demonstrated that RCH is regulated at the cellular level by post-translational signaling mechanisms, and here we tested the hypothesis that cultured cells are capable of RCH. A 2 h cold shock at -8°C significantly reduced the metabolic viability of Drosophila S2 cells, but pre-treatment with RCH at 4°C for 2 h prevented this decrease in viability. Thus, S2 cells are capable of RCH in a similar manner to whole insects and provide a new system for investigating the cell biology of RCH.

Cold acclimation triggers major transcriptional changes in Drosophila suzukii

BACKGROUND

Insects have the capacity to adjust their physiological mechanisms during their lifetime to promote cold tolerance and cope with sublethal thermal conditions, a phenomenon referred to as thermal acclimation. The spotted wing drosophila, Drosophila suzukii, is an invasive fruit pest that, like many other species, enhances its thermotolerance in response to thermal acclimation. However, little is known about the underlying mechanisms of this plastic response. Here, we promoted flies' cold tolerance by gradually increasing acclimation duration (i.e. pre-exposure from 2 h to 9 days at 10 °C), and then compared transcriptomic responses of cold hardy versus cold susceptible phenotypes using RNA sequencing.

RESULTS

Cold tolerance of D. suzukii increased with acclimation duration; the longer the acclimation, the higher the cold tolerance. Cold-tolerant flies that were acclimated for 9 days were selected for transcriptomic analyses. RNA sequencing revealed a total of 2908 differentially expressed genes: 1583 were up- and 1325 were downregulated in cold acclimated flies. Functional annotation revealed many enriched GO-terms among which ionic transport across membranes and signaling were highly represented in acclimated flies. Neuronal activity and carbohydrate metabolism were also enriched GO-terms in acclimated flies. Results also revealed many GO-terms related to oogenesis which were underrepresented in acclimated flies.

CONCLUSIONS

Involvement of a large cluster of genes related to ion transport in cold acclimated flies suggests adjustments in the capacity to maintain ion and water homeostasis. These processes are key mechanisms underlying cold tolerance in insects. Down regulation of genes related to oogenesis in cold acclimated females likely reflects that females were conditioned at 10 °C, a temperature that prevents oogenesis. Overall, these results help to understand the molecular underpinnings of cold tolerance acquisition in D. suzukii. These data are of importance considering that the invasive success of D. suzukii in diverse climatic regions relates to its high thermal plasticity.

The Desaturase Gene Family is Crucially Required for Fatty Acid Metabolism and Survival of the Brown Planthopper, Nilaparvata lugens

Desaturases are essentially required for unsaturated fatty acid (UFA) biosynthesis. We identified 10 genes encoding putative desaturases in the transcriptome database of the brown planthopper (BPH), Nilaparvata lugens. These include eight First Desaturase family genes, one cytochrome b5 fused desaturase gene (Nlug-Cytb5r) and one Sphingolipid Desaturase gene (Nlug-ifc). Transcript level profiling revealed significant variation in the expression patterns of these genes across tissues and developmental stages, which occur in a gene-specific manner. Interestingly, their expression was also modulated by the insect food source: the mRNA levels of Nlug-desatC and Nlug-Cytb5r were down-regulated, but the expression level of Nlug-desatA1-b and Nlug-desatA1-c were elevated in the BPH fed on the resistant rice variety Babawee as compared to the non-resistant variety Taichun Native 1 (TN1). Silencing Nlug-desatA1-b, Nlug-desatA1-c, or Nlug-Ifc reduced fatty acid composition and abundance in female BPH 1-d-old-adults compared to controls. Whereas, single knockdown of all ten desaturase genes significantly increased mortality of BPH nymphs compared with controls. Of the ten desaturase genes, knockdown of Nlug-desatA1-b and Nlug-desatA2 caused the highest mortality in BPH (91% and 97%, respectively). Our findings offer a base for expression and functional characterization of newly identified desaturase genes in BPH, and may contribute to RNA interference-based pest management strategies.

Inter and Intraspecific Genomic Divergence in Drosophila montana Shows Evidence for Cold Adaptation

The genomes of species that are ecological specialists will likely contain signatures of genomic adaptation to their niche. However, distinguishing genes related to ecological specialism from other sources of selection and more random changes is a challenge. Here, we describe the genome of Drosophila montana, which is the most extremely cold-adapted Drosophila species known. We use branch tests to identify genes showing accelerated divergence in contrasts between cold- and warm-adapted species and identify about 250 genes that show differences, possibly driven by a lower synonymous substitution rate in cold-adapted species. We also look for evidence of accelerated divergence between D. montana and D. virilis, a previously sequenced relative, but do not find strong evidence for divergent selection on coding sequence variation. Divergent genes are involved in a variety of functions, including cuticular and olfactory processes. Finally, we also resequenced three populations of D. montana from across its ecological and geographic range. Outlier loci were more likely to be found on the X chromosome and there was a greater than expected overlap between population outliers and those genes implicated in cold adaptation between Drosophila species, implying some continuity of selective process at these different evolutionary scales.

Effects of cold acclimation and dsRNA injections on Gs1l gene splicing in Drosophila montana

Alternative splicing, in which one gene produce multiple transcripts, may influence how adaptive genes respond to specific environments. A newly produced transcriptome of Drosophila montana shows the Gs1-like (Gs1l) gene to express multiple splice variants and to be down regulated in cold acclimated flies with increased cold tolerance. Gs1l’s effect on cold tolerance was further tested by injecting cold acclimated and non-acclimated flies from two distantly located northern and southern fly populations with double stranded RNA (dsRNA) targeting Gs1l. While both populations had similar cold acclimation responses, dsRNA injections only effected the northern population. The nature of splicing expression was then investigated in the northern population by confirming which Gs1l variants are present, by comparing the expression of different gene regions and by predicting the protein structures of splices using homology modelling. We find different splices of Gs1l not only appear to have independent impacts on cold acclimation but also elicit different effects in populations originating from two very different environments. Also, at the protein level, Gs1l appears homologous to the human HDHD1A protein and some splices might produce functionally different proteins though this needs to be verified in future studies by measuring the particular protein levels. Taken together, Gs1l appears to be an interesting new candidate to test how splicing influences adaptations.

Photosensitive Alternative Splicing of the Circadian Clock Gene timeless Is Population Specific in a Cold-Adapted Fly, Drosophila montana

To function properly, organisms must adjust their physiology, behavior and metabolism in response to a suite of varying environmental conditions. One of the central regulators of these changes is organisms’ internal circadian clock, and recent evidence has suggested that the clock genes are also important in the regulation of seasonal adjustments. In particular, thermosensitive splicing of the core clock gene timeless in a cosmopolitan fly, Drosophila melanogaster, has implicated this gene to be involved in thermal adaptation. To further investigate this link we examined the splicing of timeless in a northern malt fly species, Drosophila montana, which can withstand much colder climatic conditions than its southern relative. We studied northern and southern populations from two different continents (North America and Europe) to find out whether and how the splicing of this gene varies in response to different temperatures and day lengths. Interestingly, we found that the expression of timeless splice variants was sensitive to differences in light conditions, and while the flies of all study populations showed a change in the usage of splice variants in constant light compared to LD 22:2, the direction of the shift varied between populations. Overall, our findings suggest that the splicing of timeless in northern Drosophila montana flies is photosensitive, rather than thermosensitive and highlights the value of studying multiple species and populations in order to gain perspective on the generality of gene function changes in different kinds of environmental conditions.

Differential expression patterns of two delta-9-acyl-CoA desaturases in Thitarodes pui (Lepidoptera: Hepialidae) during different seasons and cold exposure

Thitarodes pui larvae have a limited distribution in the Tibetan Plateau and are the host of a parasitic fungus, Ophiocordyceps sinensis. Low temperature is a main environmental stress. However, understanding of T. pui cold adaptation mechanisms is insufficient. Delta-9-acyl-CoA desaturase (D9D) is closely correlated with cold adaptation for many organisms. To further understand the cold adaptation processes in T. pui larvae, two D9Ds, TpdesatA and TpdesatB were sequenced, and expression patterns were investigated during different seasons and cold exposure (under 0°C) in the laboratory. The full lengths of two cDNAs are 1,290 bp and 1,603 bp, and the ORFs encode a polypeptide of 348 and 359 amino acids, respectively. Four transmembrane domains, three conserved histidine residues and five hydrophobic regions exist in these two sequences. The expression level of TpdesatA is up-regulated in the long-term cold exposure and negatively correlated with temperature in seasonal patterns. TpdesatB responds to cold temperature in short-term cold exposure and positively corresponds temporarily in seasonal expression. Two D9Ds may have different substrate specificities, TpdesatA tends to use C16:0 and C18:0 as substrate while TpdesatB prefers C18:0. In conclusion, TpdesatA may play a very important role in T. pui cold tolerance and TpdesatB regulates function in short-term cold exposure and content change of fatty acids in the body.

Inducing Cold-Sensitivity in the Frigophilic Fly Drosophila montana by RNAi

Cold acclimation is a critical physiological adaptation for coping with seasonal cold. By increasing their cold tolerance individuals can remain active for longer at the onset of winter and can recover more quickly from a cold shock. In insects, despite many physiological studies, little is known about the genetic basis of cold acclimation. Recently, transcriptomic analyses in Drosophila virilis and D. montana revealed candidate genes for cold acclimation by identifying genes upregulated during exposure to cold. Here, we test the role of myo-inositol-1-phosphate synthase (Inos), in cold tolerance in D. montana using an RNAi approach. D. montana has a circumpolar distribution and overwinters as an adult in northern latitudes with extreme cold. We assessed cold tolerance of dsRNA knock-down flies using two metrics: chill-coma recovery time (CCRT) and mortality rate after cold acclimation. Injection of dsRNAInos did not alter CCRT, either overall or in interaction with the cold treatment, however it did induced cold-specific mortality, with high levels of mortality observed in injected flies acclimated at 5°C but not at 19°C. Overall, injection with dsRNAInos induced a temperature-sensitive mortality rate of over 60% in this normally cold-tolerant species. qPCR analysis confirmed that dsRNA injection successfully reduced gene expression of Inos. Thus, our results demonstrate the involvement of Inos in increasing cold tolerance in D. montana. The potential mechanisms involved by which Inos increases cold tolerance are also discussed.

Transcriptional Differences between Diapausing and Non-Diapausing D. montana Females Reared under the Same Photoperiod and Temperature

BACKGROUND

A wide range of insects living at higher latitudes enter diapause at the end of the warm season, which increases their chances of survival through harsh winter conditions. In this study we used RNA sequencing to identify genes involved in adult reproductive diapause in a northern fly species, Drosophila montana. Both diapausing and non-diapausing flies were reared under a critical day length and temperature, where about half of the emerging females enter diapause enabling us to eliminate the effects of varying environmental conditions on gene expression patterns of the two types of female flies.

RESULTS

RNA sequencing revealed large differences between gene expression patterns of diapausing and non-diapausing females, especially in genes involved with metabolism, fatty acid biosynthesis, and metal and nucleotide binding. Differently expressed genes included several gene groups, including myosin, actin and cytochromeP450 genes, which have been previously associated with diapause. This study also identified new candidate genes, including some involved in cuticular hydrocarbon synthesis or regulation (desat1 and desat2), and acyl-CoA Δ11-desaturase activity (CG9747), and few odorant-binding protein genes (e.g. Obp44A). Also, several transposable elements (TEs) showed differential expression between the two female groups motivating future research on their roles in diapause.

CONCLUSIONS

Our results demonstrate that the adult reproductive diapause in D. montana involves changes in the expression level of a variety of genes involved in key processes (e.g. metabolism and fatty acid biosynthesis) which help diapausing females to cope with overwintering. This is consistent with the view that diapause is a complex adaptive phenotype where not only sexual maturation is arrested, but also changes in adult physiology are required in order to survive over the winter.

Knockdown expression of Syndecan in the fat body impacts nutrient metabolism and the organismal response to environmental stresses in Drosophila melanogaster

The heparan sulfate proteoglycan syndecans are transmembrane proteins involved in multiple physiological processes, including cell-matrix adhesion and inflammation. Recent evidence from model systems and humans suggest that syndecans have a role in energy balance and nutrient metabolism regulation. However, much remains to be learned about the mechanisms through which syndecans influence these phenotypes. Previously, we reported that Drosophila melanogaster Syndecan (Sdc) mutants had reduced metabolic activity compared to controls. Here, we knocked down endogenous Sdc expression in the fat body (the functional equivalent of mammalian adipose tissue and liver) to investigate whether the effects on metabolism originate from this tissue. We found that knocking down Sdc in the fat body leads to flies with higher levels of glycogen and fat and that survive longer during starvation, likely due to their extra energy reserves and an increase in gluconeogenesis. However, compared to control flies, they are also more sensitive to environmental stresses (e.g. bacterial infection and cold) and have reduced metabolic activity under normal feeding conditions. Under the same conditions, fat-body Sdc reduction enhances expression of genes involved in glyceroneogenesis and gluconeogenesis and induces a drastic decrease in phosphorylation levels of AKT and extracellular signal regulated kinase 1/2 (ERK1/2). Altogether, these findings strongly suggest that Drosophila fat body Sdc is involved in a mechanism that shifts resources to different physiological functions according to nutritional status.

Characterization of Multiple Heat-Shock Protein Transcripts from Cydia pomonella: Their Response to Extreme Temperature and Insecticide Exposure

The economically important fruit pest Cydia pomonella (L.) exhibits a strong adaptability and stress tolerance to environmental stresses. Heat-shock proteins (HSPs) play key roles in insects in coping with environmental stresses. However, little is known about the spatiotemporal expression patterns of HSPs and their response to stresses in C. pomonella. In this study, a thermal treatment-recovery test was performed, and the expression profiles of a novel isolated HSP, named CpHSP40, and six CpHSPs were determined. Third-instar larvae were able to recover from cold shock (0 °C) and heat shock (40 °C). Escherichia coli BL21 (DE3) cells harboring recombinant pET-28a (+)-CpHSP40 plasmid showed significant temperature tolerance. CpHSPs were developmentally and tissue-specifically expressed. The responses of CpHSPs to 0 and 40 °C (with or without recovery) and insecticide exposure were varied. All of these indicated that the expression of HSPs plays a role in the development and in environmental adaptation in C. pomonella.

Preparing for Winter: The Transcriptomic Response Associated with Different Day Lengths in Drosophila montana

At northern latitudes, the most robust cue for assessing the onset of winter is the shortening of day lengths. Many species use day length as a cue to increase their cold tolerance and/or enter into diapause, but little is known about changes in gene expression that occur under different day lengths. We investigate the gene expression changes associated with differences in light/dark cycles in Drosophila montana, a northerly distributed species with a strong adult photoperiodic reproductive diapause. To examine gene expression changes induced by light both prior to and during diapause, we used both nondiapausing and diapausing flies. We found that the majority of genes that are differentially expressed between different day lengths in nondiapausing and diapausing flies differ. However, the biological processes involved were broadly similar. These included neuron development and metabolism, which are largely consistent with an increase in cold tolerance previously observed to occur in these flies. We also found that many genes associated with reproduction change in expression level between different day lengths, suggesting that D. montana use changes in day length to cue changes in reproduction both before and after entering into diapause. Finally, we also identified several interesting candidate genes for light-induced changes including Lsp2, para, and Ih.

Seasonal cues induce phenotypic plasticity of Drosophila suzukii to enhance winter survival

BACKGROUND

As global climate change and exponential human population growth intensifies pressure on agricultural systems, the need to effectively manage invasive insect pests is becoming increasingly important to global food security. Drosophila suzukii is an invasive pest that drastically expanded its global range in a very short time since 2008, spreading to most areas in North America and many countries in Europe and South America. Preliminary ecological modeling predicted a more restricted distribution and, for this reason, the invasion of D. suzukii to northern temperate regions is especially unexpected. Investigating D. suzukii phenology and seasonal adaptations can lead to a better understanding of the mechanisms through which insects express phenotypic plasticity, which likely enables invasive species to successfully colonize a wide range of environments.

RESULTS

We describe seasonal phenotypic plasticity in field populations of D. suzukii. Specifically, we observed a trend of higher proportions of flies with the winter morph phenotype, characterized by darker pigmentation and longer wing length, as summer progresses to winter. A laboratory-simulated winter photoperiod and temperature (12:12 L:D and 10 °C) were sufficient to induce the winter morph phenotype in D. suzukii. This winter morph is associated with increased survival at 1 °C when compared to the summer morph, thus explaining the ability of D. suzukii to survive cold winters. We then used RNA sequencing to identify gene expression differences underlying seasonal differences in D. suzukii physiology. Winter morph gene expression is consistent with known mechanisms of cold-hardening such as adjustments to ion transport and up-regulation of carbohydrate metabolism. In addition, transcripts involved in oogenesis and DNA replication were down-regulated in the winter morph, providing the first molecular evidence of a reproductive diapause in D. suzukii.

CONCLUSIONS

To date, D. suzukii cold resistance studies suggest that this species cannot overwinter in northern locations, e.g. Canada, even though they are established pests in these regions. Combining physiological investigations with RNA sequencing, we present potential mechanisms by which D. suzukii can overwinter in these regions. This work may contribute to more accurate population models that incorporate seasonal variation in physiological parameters, leading to development of better management strategies.

Thermal reactionomes reveal divergent responses to thermal extremes in warm and cool-climate ant species

Background

The distributions of species and their responses to climate change are in part determined by their thermal tolerances. However, little is known about how thermal tolerance evolves. To test whether evolutionary extension of thermal limits is accomplished through enhanced cellular stress response (enhanced response), constitutively elevated expression of protective genes (genetic assimilation) or a shift from damage resistance to passive mechanisms of thermal stability (tolerance), we conducted an analysis of the reactionome: the reaction norm for all genes in an organism’s transcriptome measured across an experimental gradient. We characterized thermal reactionomes of two common ant species in the eastern U.S, the northern cool-climate Aphaenogaster picea and the southern warm-climate Aphaenogaster carolinensis, across 12 temperatures that spanned their entire thermal breadth.

Results

We found that at least 2 % of all genes changed expression with temperature. The majority of upregulation was specific to exposure to low temperatures. The cool-adapted A. picea induced expression of more genes in response to extreme temperatures than did A. carolinensis, consistent with the enhanced response hypothesis. In contrast, under high temperatures the warm-adapted A. carolinensis downregulated many of the genes upregulated in A. picea, and required more extreme temperatures to induce down-regulation in gene expression, consistent with the tolerance hypothesis. We found no evidence for a trade-off between constitutive and inducible gene expression as predicted by the genetic assimilation hypothesis.

Conclusions

These results suggest that increases in upper thermal limits may require an evolutionary shift in response mechanism away from damage repair toward tolerance and prevention.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2466-z) contains supplementary material, which is available to authorized users.

Parallel molecular routes to cold adaptation in eight genera of New Zealand stick insects

The acquisition of physiological strategies to tolerate novel thermal conditions allows organisms to exploit new environments. As a result, thermal tolerance is a key determinant of the global distribution of biodiversity, yet the constraints on its evolution are not well understood. Here we investigate parallel evolution of cold tolerance in New Zealand stick insects, an endemic radiation containing three montane-occurring species. Using a phylogeny constructed from 274 orthologous genes, we show that stick insects have independently colonized montane environments at least twice. We compare supercooling point and survival of internal ice formation among ten species from eight genera, and identify both freeze tolerance and freeze avoidance in separate montane lineages. Freeze tolerance is also verified in both lowland and montane populations of a single, geographically widespread, species. Transcriptome sequencing following cold shock identifies a set of structural cuticular genes that are both differentially regulated and under positive sequence selection in each species. However, while cuticular proteins in general are associated with cold shock across the phylogeny, the specific genes at play differ among species. Thus, while processes related to cuticular structure are consistently associated with adaptation for cold, this may not be the consequence of shared ancestral genetic constraints.

Seasonal gene expression kinetics between diapause phases in Drosophila virilis group species and overwintering differences between diapausing and non-diapausing females

Most northern insect species experience a period of developmental arrest, diapause, which enables them to survive over the winter and postpone reproduction until favorable conditions. We studied the timing of reproductive diapause and its long-term effects on the cold tolerance of Drosophila montana, D. littoralis and D. ezoana females in seasonally varying environmental conditions. At the same time we traced expression levels of 219 genes in D. montana using a custom-made microarray. We show that the seasonal switch to reproductive diapause occurs over a short time period, and that overwintering in reproductive diapause has long-lasting effects on cold tolerance. Some genes, such as Hsc70, Jon25Bi and period, were upregulated throughout the diapause, while others, including regucalcin, couch potato and Thor, were upregulated only at its specific phases. Some of the expression patterns induced during the sensitive stage, when the females either enter diapause or not, remained induced regardless of the later conditions. qPCR analyses confirmed the findings of the microarray analysis in D. montana and revealed similar gene expression changes in D. littoralis and D. ezoana. The present study helps to achieve a better understanding of the genetic regulation of diapause and of the plasticity of seasonal responses in general.

The Role of Inducible Hsp70, and Other Heat Shock Proteins, in Adaptive Complex of Cold Tolerance of the Fruit Fly (Drosophila melanogaster)

Background

The ubiquitous occurrence of inducible Heat Shock Proteins (Hsps) up-regulation in response to cold-acclimation and/or to cold shock, including massive increase of Hsp70 mRNA levels, often led to hasty interpretations of its role in the repair of cold injury expressed as protein denaturation or misfolding. So far, direct functional analyses in Drosophila melanogaster and other insects brought either limited or no support for such interpretations. In this paper, we analyze the cold tolerance and the expression levels of 24 different mRNA transcripts of the Hsps complex and related genes in response to cold in two strains of D. melanogaster: the wild-type and the Hsp70^- null mutant lacking all six copies of Hsp70 gene.

Principal Findings

We found that larvae of both strains show similar patterns of Hsps complex gene expression in response to long-term cold-acclimation and during recovery from chronic cold exposures or acute cold shocks. No transcriptional compensation for missing Hsp70 gene was seen in Hsp70^- strain. The cold-induced Hsps gene expression is most probably regulated by alternative splice variants C and D of the Heat Shock Factor. The cold tolerance in Hsp70^- null mutants was clearly impaired only when the larvae were exposed to severe acute cold shock. No differences in mortality were found between two strains when the larvae were exposed to relatively mild doses of cold, either chronic exposures to 0°C or acute cold shocks at temperatures down to -4°C.

Conclusions

The up-regulated expression of a complex of inducible Hsps genes, and Hsp70 mRNA in particular, is tightly associated with cold-acclimation and cold exposure in D. melanogaster. Genetic elimination of Hsp70 up-regulation response has no effect on survival of chronic exposures to 0°C or mild acute cold shocks, while it negatively affects survival after severe acute cold shocks at temperaures below -8°C.

Localization of quantitative trait loci for diapause and other photoperiodically regulated life history traits important in adaptation to seasonally varying environments

Seasonally changing environments at high latitudes present great challenges for the reproduction and survival of insects, and photoperiodic cues play an important role in helping them to synchronize their life cycle with prevalent and forthcoming conditions. We have mapped quantitative trait loci (QTL) responsible for the photoperiodic regulation of four life history traits, female reproductive diapause, cold tolerance, egg-to-eclosion development time and juvenile body weight in Drosophila montana strains from different latitudes in Canada and Finland. The F2 progeny of the cross was reared under a single photoperiod (LD cycle 16:8), which the flies from the Canadian population interpret as early summer and the flies from the Finnish population as late summer. The analysis revealed a unique QTL for diapause induction on the X chromosome and several QTL for this and the other measured traits on the 4th chromosome. Flies' cold tolerance, egg-to-eclosion development time and juvenile body weight had several QTL also on the 2nd, 3rd and 5th chromosome, some of the peaks overlapping with each other. These results suggest that while the downstream output of females' photoperiodic diapause response is partly under a different genetic control from that of the other traits in the given day length, all traits also share some QTL, possibly involving genes with pleiotropic effects and/or multiple tightly linked genes. Nonoverlapping QTL detected for some of the traits also suggest that the traits are potentially capable of independent evolution, even though this may be restricted by epistatic interactions and/or correlations and trade-offs between the traits.

Insect capa neuropeptides impact desiccation and cold tolerance

The success of insects is linked to their impressive tolerance to environmental stress, but little is known about how such responses are mediated by the neuroendocrine system. Here we show that the capability (capa) neuropeptide gene is a desiccation- and cold stress-responsive gene in diverse dipteran species. Using targeted in vivo gene silencing, physiological manipulations, stress-tolerance assays, and rationally designed neuropeptide analogs, we demonstrate that the Drosophila melanogaster capa neuropeptide gene and its encoded peptides alter desiccation and cold tolerance. Knockdown of the capa gene increases desiccation tolerance but lengthens chill coma recovery time, and injection of capa peptide analogs can reverse both phenotypes. Immunohistochemical staining suggests that capa accumulates in the capa-expressing Va neurons during desiccation and nonlethal cold stress but is not released until recovery from each stress. Our results also suggest that regulation of cellular ion and water homeostasis mediated by capa peptide signaling in the insect Malpighian (renal) tubules is a key physiological mechanism during recovery from desiccation and cold stress. This work augments our understanding of how stress tolerance is mediated by neuroendocrine signaling and illustrates the use of rationally designed peptide analogs as agents for disrupting protective stress tolerance.

How consistent are the transcriptome changes associated with cold acclimation in two species of the Drosophila virilis group?

For many organisms the ability to cold acclimate with the onset of seasonal cold has major implications for their fitness. In insects, where this ability is widespread, the physiological changes associated with increased cold tolerance have been well studied. Despite this, little work has been done to trace changes in gene expression during cold acclimation that lead to an increase in cold tolerance. We used an RNA-Seq approach to investigate this in two species of the Drosophila virilis group. We found that the majority of genes that are differentially expressed during cold acclimation differ between the two species. Despite this, the biological processes associated with the differentially expressed genes were broadly similar in the two species. These included: metabolism, cell membrane composition, and circadian rhythms, which are largely consistent with previous work on cold acclimation/cold tolerance. In addition, we also found evidence of the involvement of the rhodopsin pathway in cold acclimation, a pathway that has been recently linked to thermotaxis. Interestingly, we found no evidence of differential expression of stress genes implying that long-term cold acclimation and short-term stress response may have a different physiological basis.

Analysis of survival, gene expression and behavior following chill-coma in the medfly Ceratitis capitata: effects of population heterogeneity and age

The medfly Ceratitis capitata is an agricultural pest distributed worldwide thanks, in part, to its phenotypic plasticity of thermal tolerance. Cold exposure has been shown to reduce C. capitata survival, which may affect its distribution in areas with subfreezing temperatures. When insects are increasingly cooled, they attain a critical thermal threshold and enter a chill-coma state characterized by cessation of movement. It is not clear how a rapid cold exposure affects the physiological state of medflies, and how this is influenced by age and population heterogeneity. In order to approach these questions, C. capitata single-sex laboratory populations of 15 and 30 days old were subjected to a chill-coma recovery assay, and separated according to their recovery time in three subgroups: Fast-Subgroups, Intermediate-Subgroups, and Slow-Subgroups. Thereafter, we analyzed their survival, behavioral, and gene expression outputs. In female and old male populations, we found that flies with the slowest recovery time had a reduced life expectancy, a higher initial mortality rate, and a worse climbing performance compared with flies that recovered faster. Therefore, we were able to separate subgroups that developed chilling-injury from subgroups that had a reversible full recovery after cold exposure. The gene expression analysis of the heat shock protein genes hsp70 and hsp83 showed no clear association with the parameters studied. Interestingly, thorax expression levels of the Cu/Zn superoxide dismutase gene were elevated during the recovery phase in the Fast-Subgroups, but remained constant in the Slow-Subgroups that developed chilling-injury. On the other hand, none of the young male subgroups seemed to have suffered irreversible damage. Thus, we concluded that depending on age and population heterogeneity, chill-coma recovery time points out significant differences on individual cold tolerance. Moreover, the inability to properly induce the antioxidant defense system to counteract the oxidative damage caused by cold seems to contribute to the development of chilling-injury.

Anoxic stress and rapid cold hardening enhance cold tolerance of the migratory locust

Anoxia and rapid cold hardening (RCH) can increase the cold tolerance of many animals. However, mechanisms underlying these two kinds of stresses remain unclear. In this study, we aimed to explore the relationship of acclimation to cold stress with acclimation to anoxic stress in the migratory locust, Locusta migratoria. RCH at 0°C for 3h promoted the survival of cold stress-exposed locusts. Anoxic hypercapnia (CO2 anoxic treatment) for 40 min exerted an effect similar to that of RCH. Anoxic hypercapnia within 1h can all promote the cold hardiness of locusts. We investigated the transcript levels of six heat shock protein (Hsp) genes, namely, Hsp20.5, Hsp20.6, Hsp20.7, Hsp40, Hsp70, and Hsp90. Four genes, namely, Hsp90, Hsp40, Hsp20.5, and Hsp20.7, showed differential responses to RCH and anoxic hypercapnia treatments. Under cold stress, locusts exposed to the two regimens showed different responses for Hsp90, Hsp20.5, and Hsp20.7. However, the varied responses disappeared after recovery from cold stress. Compared with the control group, the transcript levels of six Hsp genes were generally downregulated in locusts subjected to anoxic hypercapnia or/and RCH. These results indicate that anoxic stress and RCH have different mechanisms of regulating the transcription of Hsp family members even if the two treatments exerted similar effects on cold tolerance of the migratory locust. However, Hsps may not play a major role in the promotion of cold hardiness by the two treatments.

Response to cold acclimation in diapause pupae of Hyles euphorbiae (Lepidoptera: Sphingidae): candidate biomarker identification using proteomics

The distribution range of Hyles euphorbiae covers distinct climates across the Palaearctic. Previous investigations showed a correlation between mitochondrial DNA identity of populations and climatic conditions related to winter; however, the lack of biomarkers hampers investigations to test whether geographically distinct populations do show specific molecular responses to low temperatures or whether they possess specific genetic identity at loci functionally related to cold response. The present study was designed to identify candidate protein biomarkers and biological processes that are associated with cold acclimation of overwintering H. euphorbiae diapause pupae. Specimens taken from a single central European population were gradually cooled from 20 °C to -2 °C over 36 days and 12 differentially abundant proteins were identified. In addition, DeepSuperSAGE sequencing technology was applied to study differentially regulated genes. There was incongruence between differentially abundant proteins and differentially expressed genes, but functional characteristics of regulated proteins and analyses of gene ontology term enrichment among differentially regulated genes pointed to activation of the same biological processes, e.g. oxidative stress response. As proteins represent biologically active molecules, candidate biomarkers derived from proteomics are considered well suited to explore intraspecific patterns of local adaptation to different climates.

The rapid cold hardening response of Drosophila melanogaster: complex regulation across different levels of biological organization

Rapid cold hardening (RCH) is a form of thermal acclimation that allows ectotherms to fine-tune their physiological state to match rapid changes in thermal environment. Despite progress in recent years, there is still a considerable uncertainty regarding the physiological basis of RCH in insects. Here we investigated the physiological response of adult Drosophila melanogaster to a gradual reduction of temperature from 25 to 0°C followed by 1h at 0°C. As expected, this RCH treatment promoted cold tolerance, and so we hypothesized that this change could be detected at the proteomic level. Using 2D-DIGE, we found that only a few proteins significantly changed in abundance, and of these, we identified a set of four proteins of particular interest. These were identified as two different variants of glycogen phosphorylase (GlyP) of which three spots were up-regulated and another was down regulated. In subsequent experiments, we quantified upstream events by measuring the GlyP mRNA amount, but we found no marked effect of RCH. We also examined downstream events by measuring GlyP activity and the level of free sugars. We found no effect of RCH on GlyP activity. On the other hand, screening of whole animal sugar contents revealed a small increase in glucose levels following RCH while trehalose content was unaltered. This study highlights a complex regulation of GlyP in relation to RCH where we found associations between the cold tolerance, the protein abundance and the metabolite concentrations but no changes in mRNA expression and enzyme activity. These data stress the necessity of combining the hypothesis-generating power of an 'Omics' approach with subsequent targeted validations across several levels of the biological organization. We discuss reasons why different biological linked levels do not necessarily change stoichiometrically.

Sequence variation in couch potato and its effects on life-history traits in a northern malt fly, Drosophila montana

Couch potato (cpo) has previously been connected to reproductive diapause in several insect species including Drosophila melanogaster, where it has been suggested to provide a link between the insulin signalling pathway and the hormonal control of diapause. In the first part of the study we sequenced nearly 3.6 kb of this gene in a northern Drosophila species (Drosophila montana) with a robust photoperiodically determined diapause and found several types of polymorphisms along the sequenced area. We also found variation among five Drosophila virilis group species in the length of the 5th exon of cpo and in the site of the stop codon at the end of this exon. The second part of the study was targeted on a deletion of six amino acids located in the last section of exon 5, which in D. melanogaster, is translated only in one short transcript lacking the following exons. The studied deletion appeared to be extremely rare in the wild D. montana population where it was found, but its frequency rapidly increased during laboratory culture. qPCR analyses showed the expression level of the deletion allele to be significantly downregulated in both the diapausing and non-diapausing females compared to the wild type allele. At the phenotypic level, the deletion and the decreased expression of cpo transcript involving it did not have direct effect on the incidence of female reproductive diapause, but it was associated with a reduction in development time under diapause-inducing conditions. This suggests that while the cpo transcript containing the prolonged version of the 5th exon with a stop codon is clearly associated with fly development time, the exons with RNA domains included in other transcripts of the gene may be more directly related to diapause regulation.