Cold hardening and transcriptional change in Drosophila melanogaster

Comparative Transcriptome Analysis Reveals Different Responses in Three Developmental Stages of Mythimna loreyi to Cold Stress

The loreyi leafworm Mythimna loreyi (Lepidoptera: Noctuidae) is a serious pest of agriculture that causes particular damage to Gramineae crops in Asia, Europe, Australia, Africa, and the Middle East. Low temperature is one of the important environmental factors that limits the survival, distribution, colonization, and abundance of M. loreyi. However, the metabolic synthesis pathways of cold-tolerant substances in M. loreyi and the key genes involved in the regulation under cold stress remain largely unknown. In this study, we sequenced the transcriptomes of three developmental stages (larvae, pupae, and adults) of M. loreyi to discover the molecular mechanisms of their responses to cold stress. In total, sequencing generated 120.64 GB of clean data from 18 samples, of which 19,459 genes and 1740 differentially expressed genes (DEGs) were identified. The enrichment analysis of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) revealed that many DEGs were mainly enriched in pathways associated with energy metabolism and hormone metabolism. Among these, genes encoding multiple metabolic enzymes, cuticle proteins (CPs), and heat shock proteins (HSPs) were differentially expressed. These results indicate that there are significant differences among the three developmental stages of M. loreyi exposed to cold stress and provide a basis for further studying the molecular mechanisms of cold tolerance in insects.

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.

Advances in Novel Animal Vitamin C Biosynthesis Pathways and the Role of Prokaryote-Based Inferences to Understand Their Origin

Vitamin C (VC) is an essential nutrient required for the optimal function and development of many organisms. VC has been studied for many decades, and still today, the characterization of its functions is a dynamic scientific field, mainly because of its commercial and therapeutic applications. In this review, we discuss, in a comparative way, the increasing evidence for alternative VC synthesis pathways in insects and nematodes, and the potential of myo-inositol as a possible substrate for this metabolic process in metazoans. Methodological approaches that may be useful for the future characterization of the VC synthesis pathways of Caenorhabditis elegans and Drosophila melanogaster are here discussed. We also summarize the current distribution of the eukaryote aldonolactone oxidoreductases gene lineages, while highlighting the added value of studies on prokaryote species that are likely able to synthesize VC for both the characterization of novel VC synthesis pathways and inferences on the complex evolutionary history of such pathways. Such work may help improve the industrial production of VC.

Gliotransmission of D-serine promotes thirst-directed behaviors in Drosophila

Thirst emerges from a range of cellular changes that ultimately motivate an animal to consume water. Although thirst-responsive neuronal signals have been reported, the full complement of brain responses is unclear. Here, we identify molecular and cellular adaptations in the brain using single-cell sequencing of water-deprived Drosophila. Water deficiency primarily altered the glial transcriptome. Screening the regulated genes revealed astrocytic expression of the astray-encoded phosphoserine phosphatase to bi-directionally regulate water consumption. Astray synthesizes the gliotransmitter D-serine, and vesicular release from astrocytes is required for drinking. Moreover, dietary D-serine rescues aay-dependent drinking deficits while facilitating water consumption and expression of water-seeking memory. D-serine action requires binding to neuronal NMDA-type glutamate receptors. Fly astrocytes contribute processes to tripartite synapses, and the proportion of astrocytes that are themselves activated by glutamate increases with water deprivation. We propose that thirst elevates astrocytic D-serine release, which awakens quiescent glutamatergic circuits to enhance water procurement.

Differential Transcriptome Analysis Reveals Genes Related to Low- and High-Temperature Stress in the Fall Armyworm, Spodoptera frugiperda

The fall armyworm (FAW), Spodoptera frugiperda, is regarded as one of the world’s most harmful plant pests. This research examines the molecular response processes of FAW to low temperature (4°C) and high temperature (40°C) when gene expression is compared to controls (25°C). A total of 211,967 unigenes were collected, at least 14,338 of which were annotated with gene descriptions, gene ontology terms, and metabolic pathways. There were 50 Gene Ontology (GO) functional sub-groups and 21 EggNOG words as a result. Differentially expresses genes (DEGs) with log₂FC ≥ 2 were identified and compared at various temperatures. In comparison to the 25°C treated group, we discovered 199 and 1,248 individual DEGs co-regulated at 4 and 40°C, respectively. Comparing transcriptome profiles for differential gene expression revealed a number of DEGs, including cytochrome P450, odorant binding proteins (OBPs), and immune system genes previously implicated in cold and high temperature stresses. The enrichment pathways were identified using Kyoto Encyclopedia of Genes and Genomics (KEGG) analysis, and heatmaps of similar unigenes from both treatment groups (T4 and T40) were plotted. We used quantitative reverse transcription PCR (RT-qPCR) to confirm the RNA-seq data on 10 up- and down-regulated DEGs. These findings provide a foundation for future understanding of FAW adaptation mechanisms and the underlying basis underlying the response to low and high temperatures.

Mitochondria as a target and central hub of energy division during cold stress in insects

Temperature stress is one of the crucial factors determining geographical distribution of insect species. Most of them are active in moderate temperatures, however some are capable of surviving in extremely high as well as low temperatures, including freezing. The tolerance of cold stress is a result of various adaptation strategies, among others the mitochondria are an important player. They supply cells with the most prominent energy carrier—ATP, needed for their life processes, but also take part in many other processes like growth, aging, protection against stress injuries or cell death. Under cold stress, the mitochondria activity changes in various manner, partially to minimize the damages caused by the cold stress, partially because of the decline in mitochondrial homeostasis by chill injuries. In the response to low temperature, modifications in mitochondrial gene expression, mtDNA amount or phosphorylation efficiency can be observed. So far study also showed an increase or decrease in mitochondria number, their shape and mitochondrial membrane permeability. Some of the changes are a trigger for apoptosis induced via mitochondrial pathway, that protects the whole organism against chill injuries occurring on the cellular level. In many cases, the observed modifications are not unequivocal and depend strongly on many factors including cold acclimation, duration and severity of cold stress or environmental conditions. In the presented article, we summarize the current knowledge about insect response to cold stress focusing on the role of mitochondria in that process considering differences in results obtained in different experimental conditions, as well as depending on insect species. These differentiated observations clearly indicate that it is still much to explore.

Genomic knockout of hsp23 both decreases and increases fitness under opposing thermal extremes in Drosophila melanogaster

Under exposure to harmful environmental stresses, organisms exhibit a general stress response involving upregulation of the expression of heat shock proteins (HSPs) which is thought to be adaptive. Small heat shock proteins (sHSPs) are key components of this response, although shsp genes may have other essential roles in development. However, the upregulation of expression of a suite of genes under stress may not necessarily be evidence of an adaptive response to stress that involves those genes. To explore this issue, we used the CRISPR/Cas9 system to investigate pleiotropic effects of the hsp23 gene in Drosophila melanogaster. Transgenic flies carrying a pCFD5 plasmid containing sgRNAs were created to generate a complete knockout of the hsp23 gene. The transgenic line lacking hsp23 showed an increased hatch rate and no major fitness costs under an intermediate temperature used for culturing the flies. In addition, hsp23 knockout affected tolerance to hot and cold temperature extremes but in opposing directions; knockout flies had reduced tolerance to cold, but increased tolerance to heat. Despite this, hsp23 expression (in wild type flies) was increased under both hot and cold conditions. The hsp23 gene was required for heat hardening at the pupal stage, but not at the 1st-instar larval stage, even though the gene was upregulated in wild type controls at that life stage. The phenotypic effects of hsp23 were not compensated for by expression changes in other shsps. Our study shows that the fitness consequences of an hsp gene knockout depends on environmental conditions, with potential fitness benefits of gene loss even under conditions when the gene is normally upregulated.

Phylogenomic analyses of the genus Drosophila reveals genomic signals of climate adaptation

Many Drosophila species differ widely in their distributions and climate niches, making them excellent subjects for evolutionary genomic studies. Here, we have developed a database of high‐quality assemblies for 46 Drosophila species and one closely related Zaprionus. Fifteen of the genomes were newly sequenced, and 20 were improved with additional sequencing. New or improved annotations were generated for all 47 species, assisted by new transcriptomes for 19. Phylogenomic analyses of these data resolved several previously ambiguous relationships, especially in the melanogaster species group. However, it also revealed significant phylogenetic incongruence among genes, mainly in the form of incomplete lineage sorting in the subgenus Sophophora but also including asymmetric introgression in the subgenus Drosophila. Using the phylogeny as a framework and taking into account these incongruences, we then screened the data for genome‐wide signals of adaptation to different climatic niches. First, phylostratigraphy revealed relatively high rates of recent novel gene gain in three temperate pseudoobscura and five desert‐adapted cactophilic mulleri subgroup species. Second, we found differing ratios of nonsynonymous to synonymous substitutions in several hundred orthologues between climate generalists and specialists, with trends for significantly higher ratios for those in tropical and lower ratios for those in temperate‐continental specialists respectively than those in the climate generalists. Finally, resequencing natural populations of 13 species revealed tropics‐restricted species generally had smaller population sizes, lower genome diversity and more deleterious mutations than the more widespread species. We conclude that adaptation to different climates in the genus Drosophila has been associated with large‐scale and multifaceted genomic changes.

Acclimation, duration and intensity of cold exposure determine the rate of cold stress accumulation and mortality in Drosophila suzukii

The spotted wing drosophila (SWD), Drosophila suzukii, is a major invasive fruit pest. There is strong consensus that low temperature is among the main drivers of SWD population distribution, and the invasion success of SWD is also linked to its thermal plasticity. Most studies on ectotherm cold tolerance focus on exposure to a single stressful temperature but here we investigated how cold stress intensity affected survival duration across a broad range of low temperatures (-7 to +3 °C). The analysis of Lt₅₀ at different stressful temperatures (Thermal Death Time curve - TDT) is based on the suggestion that cold injury accumulation rate increases exponentially with the intensity of thermal stress. In accordance with the hypothesis, Lt₅₀ of SWD decreased exponentially with temperature. Further, comparison of TDT curves from flies acclimated to 15, 19 and 23 °C, respectively, showed an almost full compensation with acclimation such that the temperature required to induce mortality over a fixed time decreased almost 1 °C per °C lowering of acclimation temperature. Importantly, this change in cold tolerance with acclimation was uniform across the range of moderate to intense cold stress exposures examined. To understand if cold stress at moderate and intense exposures affects the same physiological systems we examined how physiological markers/symptoms of chill injury developed at different intensities of the cold stress. Specifically, hsp23 expression and extracellular [K⁺] were measured in flies exposed to different intensities of cold stress (-6, -2 and +2 °C) and at various time points corresponding to the same progression of injury (equivalent to 1/3, 2/3 or 3/3 of Lt₅₀). The different cold stress intensities all triggered hsp23 expression following 2 h of recovery, but patterns of expression differed. At the most intense cold stress (-6 and -2 °C) a gradual increase with time was found. In contrast, at +2 °C an initial increase was followed by a dissipating expression. A gradual perturbation of ion balance (hyperkalemia) was also found at all three cold stress intensities examined, with only slight dissimilarities between treatment temperatures. Despite some differences between the three cold intensities examined, the results generally support the hypothesis that intense and moderate cold stress induces the same physiological perturbation. This suggests that cold stress experienced during natural fluctuating conditions is additive and the results also illustrate that the rate of injury accumulation increases dramatically (exponentially) with decreasing temperature (increasing stress).

Comparison of gene expression in the red imported fire ant (Solenopsis invicta) under different temperature conditions

The red imported fire ant (RIFA), Solenopsis invicta Buren is native to South America and is known as a global problematic invasive species. This study focused on the molecular response of RIFA by comparing gene expression profiles after exposing ants to low (10 °C) and high (40 °C) temperature stress and comparing them to untreated controls (30 °C). A total of 99,085 unigenes (the clustered non-redundant transcripts that are filtered from the longest assembled contigs) were obtained, of which 19,154 were annotated with gene descriptions, gene ontology terms, and metabolic pathways. 86 gene ontology (GO) functional sub-groups and 23 EggNOG terms resulted. Differentially expressed genes (DEGs) with log₂FC ≥ 10 were screened and were compared at different temperatures. We found 203, 48, and 66 specific DEGs co-regulated at 10, 20, and 40 °C. Comparing transcriptome profiles for differential gene expression resulted in various DE genes, including cytochrome P450, NADH dehydrogenase subunit 1, cuticle protein and heat shock protein (HSP), which have previously been reported to be involved in cold and high temperature resistance. GO analysis revealed that antioxidant activity is up-regulated under high temperature stress. We verified the RNA-seq data by qPCR on 20 up- and down-regulated DEGs. These findings provide a basis for future understanding of the adaptation mechanisms of RIFA and the molecular mechanisms underlying the response to low and high temperatures.

Molecular Characterization and Transcriptional Expression Analysis of ABC Transporter H Subfamily Genes in the Oriental Fruit Fly

The oriental fruit fly, Bactrocera dorsalis Hendel (Diptera: Tephretidae), is a serious pest of fruits and vegetables and has developed high levels of insecticide resistance. ATP-binding cassette transporter genes (ABC transporters) are involved in mediating the energy-driven transport of many substances across membranes and are closely associated with development and insecticide detoxification. In this study, three ABC transporters in the H subfamily were identified, and the possible roles of these genes in B. dorsalis are discussed. Bioinformatics analysis revealed that those genes are conserved, typical of half-transporters. The expression profiles of BdABCH genes (BdABCHs) in the developmental stages, tissues, and following insecticide exposure, extreme temperature, warm- and cold-acclimated strain, starvation, and desiccation stress were determined by quantitative real-time PCR. Expression of BdABCHs can be detected in various tissues and in different developmental stages. They were most highly expressed in the hindgut and in newly emerged adults. The mRNA levels of BdABCHs in males (including most tissues and body segments) were higher than in females. The expression of BdABCH1 was significantly upregulated 3.8-fold in the cold-acclimated strain, and was significantly upregulated by 1.9-, 3.8- and 4.1-fold in the 0°C, starvation, and desiccation treatments, respectively. Treatment with malathion and avermectin at LD20 and LD30 concentrations produced no obvious changes in the levels of BdABCHs. BdABCHs may be involved in the transport of related hormones during eclosion, as well as water and inorganic salts. BdABCH1 also demonstrated that it is related to the ability to cope with adverse environments.

Parallel and Population-specific Gene Regulatory Evolution in Cold-Adapted Fly Populations

Changes in gene regulation at multiple levels may comprise an important share of the molecular changes underlying adaptive evolution in nature. However, few studies have assayed within- and between-population variation in gene regulatory traits at a transcriptomic scale, and therefore inferences about the characteristics of adaptive regulatory changes have been elusive. Here, we assess quantitative trait differentiation in gene expression levels and alternative splicing (intron usage) between three closely related pairs of natural populations of Drosophila melanogaster from contrasting thermal environments that reflect three separate instances of cold tolerance evolution. The cold-adapted populations were known to show population genetic evidence for parallel evolution at the SNP level, and here we find evidence for parallel expression evolution between them, with stronger parallelism at larval and adult stages than for pupae. We also implement a flexible method to estimate cis- vs trans-encoded contributions to expression or splicing differences at the adult stage. The apparent contributions of cis- vs trans-regulation to adaptive evolution vary substantially among population pairs. While two of three population pairs show a greater enrichment of cis-regulatory differences among adaptation candidates, trans-regulatory differences are more likely to be implicated in parallel expression changes between population pairs. Genes with significant cis-effects are enriched for signals of elevated genetic differentiation between cold- and warm-adapted populations, suggesting that they are potential targets of local adaptation. These findings expand our knowledge of adaptive gene regulatory evolution and our ability to make inferences about this important and widespread process.

Involvement of small heat shock proteins (sHsps) in developmental stages of fall armyworm, Spodoptera frugiperda and its expression pattern under abiotic stress condition

Fall armyworm (FAW), Spodoptera frugiperda a recent invasive pest in India is reported to cause significant damage by feeding voraciously on maize and other economic crops from tropical to temperate provinces. It is becoming an arduous challenge to control the pest as it can survive in a wide range of temperature conditions and is already said to develop resistance towards certain insecticides. The small Heat shock proteins (hereafter, sHsps) are known to play an important role in adaptation of insects under such stress conditions. Our present study involved characterization of the three sHsps genes (sHsp19.74, sHsp20.7 and sHsp19.07) which encoded proteins of about 175, 176 and 165 amino acids with a conserved α-crystalline domain. Phylogenetic analysis of deduced amino acid sequences of the three genes showed strong similarity with the other lepidopteran sHsps. The effect of different growth stages on the expression profile of these stress proteins has also been studied and the Quantitative real time PCR (qRT-PCR) analysis revealed that the transcript level of sHsp19.07 and sHsp20.7 were significantly upregulated under extreme heat (44 °C) and cold (5 °C) stress. However, sHsp19.74 responded only to heat treatment but not to the cold treatment. In addition, the expression profile of all three sHsps was significantly lower in the larval stage (5^th instar). Chlorantraniliprole treatment resulted in maximum expression of sHsp19.07 and sHsp20.7 after 12hr of exposure to the insecticide. Meanwhile, the same expression was observed after 8hr of exposure in case of sHsp19.74. These results proved that the sHsp genes of S. frugiperda were induced and modulated in response to abiotic stress, thus influencing the physiological function leading to survival of FAW in diversified climate in India.

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.

Transcriptome analysis of Drosophila melanogaster laboratory strains of different geographical origin after long-term laboratory maintenance

Positive selection may be the main factor of the between-population divergence in gene expression. Expression profiles of two Drosophila melanogaster laboratory strains of different geographical origin and long-term laboratory maintenance were analyzed using microchip arrays encompassing probes for 18,500 transcripts. The Russian strain D18 and the North American strain Canton-S were compared. A set of 223 known or putative genes demonstrated significant changes in expression levels between these strains. Differentially expressed genes (DEG) were enriched in response to DDT (p = .0014), proteolysis (p = 2.285E-5), transmembrane transport (p = 1.03E-4), carbohydrate metabolic process (p = .0317), protein homotetramerization (p = .0444), and antibacterial humoral response (p = 425E-4). The expression in subset of genes from different categories was verified by qRT-PCR. Analysis of transcript abundance between Canton-S and D18 strains allowed to select several genes to estimate their participation in latitude adaptation. Expression of selected genes was analyzed in five D. melanogaster lines of different geographic origins by qRT-PCR, and we found two candidate genes that may be associated with latitude adaptation in adult flies-smp-30 and Cda9. Quite possible that several alleles of these genes may be important for insect survival in the environments of global warming. It is interesting that the number of genes involved in local adaptation demonstrates expression level appropriate to their geographical origin even after decades of laboratory maintenance.

Integrating GWAS and Transcriptomics to Identify the Molecular Underpinnings of Thermal Stress Responses in Drosophila melanogaster

Thermal tolerance of an organism depends on both the ability to dynamically adjust to a thermal stress and preparatory developmental processes that enhance thermal resistance. However, the extent to which standing genetic variation in thermal tolerance alleles influence dynamic stress responses vs. preparatory processes is unknown. Here, using the model species Drosophila melanogaster, we used a combination of Genome Wide Association mapping (GWAS) and transcriptomic profiling to characterize whether genes associated with thermal tolerance are primarily involved in dynamic stress responses or preparatory processes that influence physiological condition at the time of thermal stress. To test our hypotheses, we measured the critical thermal minimum (CTmin) and critical thermal maximum (CTmax) of 100 lines of the Drosophila Genetic Reference Panel (DGRP) and used GWAS to identify loci that explain variation in thermal limits. We observed greater variation in lower thermal limits, with CTmin ranging from 1.81 to 8.60°C, while CTmax ranged from 38.74 to 40.64°C. We identified 151 and 99 distinct genes associated with CTmin and CTmax, respectively, and there was strong support that these genes are involved in both dynamic responses to thermal stress and preparatory processes that increase thermal resistance. Many of the genes identified by GWAS were involved in the direct transcriptional response to thermal stress (72/151 for cold; 59/99 for heat), and overall GWAS candidates were more likely to be differentially expressed than other genes. Further, several GWAS candidates were regulatory genes that may participate in the regulation of stress responses, and gene ontologies related to development and morphogenesis were enriched, suggesting many of these genes influence thermal tolerance through effects on development and physiological status. Overall, our results suggest that thermal tolerance alleles can influence both dynamic plastic responses to thermal stress and preparatory processes that improve thermal resistance. These results also have utility for directly comparing GWAS and transcriptomic approaches for identifying candidate genes associated with thermal tolerance.

Artificial Selection to a Nonlethal Cold Stress in Trogoderma variabile Shows Associations With Chronic Cold Stress and Body Size

Extreme temperature has been used as an alternative to chemical treatments for stored product pests for years. Resistance to heat or cold treatments has not been documented in stored product insects, but repeated use of ineffective treatments could lead to adaptive tolerance. Trogoderma variabile (Dermestidae) is a common pest of stored products, and the larval stage is highly resistant to cold and destructive. We artificially selected populations by inducing chill coma at four different cold temperature treatments: 3 and 5 h at -10°C and 3 and 5 h at 0°C. Recovery time was highly heritable after selection for seven generations for decreased recovery time (cold tolerance) and increased recovery time (cold susceptibility) at all time and temperature combinations. Three replicate populations for each time and temperature combination varied substantially, suggesting different mutations in each population were probably responsible for selected phenotypes. Body size decreased in populations selected for cold susceptibility compared with those selected for cold tolerance and survivorship to long-term cold stress increased in the cold-tolerant populations compared with the susceptible populations. After the cessation of the selection experiment, cold tolerance dissipated within four generations from the populations at -10°C, but was maintained in populations exposed to 0°C. Our results suggest that warehouse beetles can adapt to cold stress quickly, but in the absence of cold stress, the proportion of cold-tolerant/susceptible individuals is quickly reduced, suggesting that some of the mutations responsible for these phenotypes may be associated with fitness costs under normal conditions.

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.

Phenotypic plasticity, but not adaptive tracking, underlies seasonal variation in post-cold hardening freeze tolerance of Drosophila melanogaster

In temperate regions, an organism's ability to rapidly adapt to seasonally varying environments is essential for its survival. In response to seasonal changes in selection pressure caused by variation in temperature, humidity, and food availability, some organisms exhibit plastic changes in phenotype. In other cases, seasonal variation in selection pressure can rapidly increase the frequency of genotypes that offer survival or reproductive advantages under the current conditions. Little is known about the relative influences of plastic and genetic changes in short-lived organisms experiencing seasonal environmental fluctuations. Cold hardening is a seasonally relevant plastic response in which exposure to cool, but nonlethal, temperatures significantly increases the organism's ability to later survive at freezing temperatures. In the present study, we demonstrate seasonal variation in cold hardening in Drosophila melanogaster and test the extent to which plasticity and adaptive tracking underlie that seasonal variation. We measured the post-cold hardening freeze tolerance of flies from outdoor mesocosms over the summer, fall, and winter. We bred outdoor mesocosm-caught flies for two generations in the laboratory and matched each outdoor cohort to an indoor control cohort of similar genetic background. We cold hardened all flies under controlled laboratory conditions and then measured their post-cold hardening freeze tolerance. Comparing indoor and field-caught flies and their laboratory-reared G1 and G2 progeny allowed us to determine the roles of seasonal environmental plasticity, parental effects, and genetic changes on cold hardening. We also tested the relationship between cold hardening and other factors, including age, developmental density, food substrate, presence of antimicrobials, and supplementation with live yeast. We found strong plastic responses to a variety of field- and laboratory-based environmental effects, but no evidence of seasonally varying parental or genetic effects on cold hardening. We therefore conclude that seasonal variation in post-cold hardening freeze tolerance results from environmental influences and not genetic changes.

Identification of candidate ATP-binding cassette transporter gene family members in Diaphorina citri (Hemiptera: Psyllidae) via adult tissues transcriptome analysis

The ATP-binding cassette (ABC) transporters exist in all living organisms and play major roles in various biological functions by transporting a wide variety of substrates across membranes. The functions of ABC transporters in drug resistance have been extensively studied in vertebrates; however, they are rarely characterized in agricultural pests. The Asian citrus psyllid, Diaphorina citri, is one of the most damaging pests of the Citrus genus because of its transmission of Huanglongbing, also known as Yellow Dragon disease. In this study, the next-generation sequencing technique was applied to research the ABC transporters of D. citri. Fifty-three ABC transporter genes were found in the RNA-Seq data, and among these ABC transporters, 4, 4, 5, 2, 1, 4, 18 and 15 ABC proteins belonged to the ABCA-ABCH subfamilies, respectively. Different expression profiles of 52 genes between imidacloprid-resistant and imidacloprid-susceptible strains were studied by qRT-PCR; 5 ABCGs and 4 ABCHs were significantly upregulated in the imidacloprid-resistant strain. In addition, five of the nine upregulated genes were widely expressed in adult tissues in spatial expression analysis. The results suggest that these genes may play key roles in this phenotype. In general, this study contributed to our current understanding of D. citri resistance to insecticides.

Comparative Analysis of Transcriptome Responses to Cold Stress in Galeruca daurica (Coleoptera: Chrysomelidae)

Galeruca daurica (Joannis) has become a new insect pest in the Inner Mongolia grasslands since 2009, and its larvae and eggs have strong cold tolerance. To get a deeper insight into its molecular mechanisms of cold stress responses, we performed de novo transcriptome assembly for G. daurica by RNA-Seq and compared the transcriptomes of its larvae exposed to five different temperature treatments (-10, -5, 0, 5, and 25°C for 1 h and then recovered at 25°C for 1 h), respectively. Compared with the control (25°C), the numbers of differentially expressed genes (DEGs) decreased from 1,821 to 882, with the temperature declining from 5 to -10°C. Moreover, we obtained 323 coregulated DEGs under different low temperatures. Under four low temperatures (-10, -5, 0, and 5°C), a large number of genes were commonly upregulated during recovery from cold stresses, including those related to cuticle protein, followed by cytochrome P450, clock protein, fatty acid synthase, and fatty acyl-CoA reductase; meanwhile, lots of genes encoding cuticle protein, RNA replication protein, RNA-directed DNA polymerase, and glucose dehydrogenase were commonly downregulated. Our findings provide important clues for further investigations of key genes and molecular mechanisms involved in the adaptation of G. daurica to harsh environments.

Comparative proteomics of stenotopic caddisfly Crunoecia irrorata identifies acclimation strategies to warming

Species' ecological preferences are often deduced from habitat characteristics thought to represent more or less optimal conditions for physiological functioning. Evolution has led to stenotopic and eurytopic species, the former having decreased niche breadths and lower tolerances to environmental variability. Species inhabiting freshwater springs are often described as being stenotopic specialists, adapted to the stable thermal conditions found in these habitats. Whether due to past local adaptation these species have evolved or have lost intra-generational adaptive mechanisms to cope with increasing thermal variability has, to our knowledge, never been investigated. By studying how the proteome of a stenotopic species changes as a result of increasing temperatures, we investigate if the absence or attenuation of molecular mechanisms is indicative of local adaptation to freshwater springs. An understanding of compensatory mechanisms is especially relevant as spring specialists will experience thermal conditions beyond their physiological limits due to climate change. In this study, the stenotopic species Crunoecia irrorata (Trichoptera: Lepidostomatidae, Curtis 1834) was acclimated to 10, 15 and 20°C for 168 hr. We constructed a homology-based database and via liquid chromatography-tandem mass spectrometry (LC-MS/MS)-based shotgun proteomics identified 1,358 proteins. Differentially abundant proteins and protein norms of reaction revealed candidate proteins and molecular mechanisms facilitating compensatory responses such as trehalose metabolism, tracheal system alteration and heat-shock protein regulation. A species-specific understanding of compensatory physiologies challenges the characterization of species as having narrow tolerances to environmental variability if that characterization is based on occurrences and habitat characteristics alone.

Transcriptome-Based Analysis Reveals a Crucial Role of BxGPCR17454 in Low Temperature Response of Pine Wood Nematode (Bursaphelenchus xylophilus)

Background: The causal agent of pine wilt disease is the pine wood nematode (PWN) (Bursaphelenchus xylophilus), whose ability to adapt different ecological niches is a crucial determinant of their invasion to colder regions. To discover the molecular mechanism of low temperature response mechanism, we attempted to study the molecular response patterns under low temperature from B. xylophilus with a comprehensive RNA sequencing analysis and validated the differentially expressed genes (DEGs) with quantitative real-time polymerase chain reaction (qRT-PCR). Bioinformatic software was utilized to isolate and identify the low-temperature-related BxGPCR genes. Transcript abundance of six low-temperature-related BxGPCR genes and function of one of the BxGPCR genes are studied by qRT-PCR and RNA interference. Results: The results showed that we detected 432 DEGs through RNA sequencing between low-temperature-treated and ambient-temperature-treated groups nematodes. The transcript level of 6 low-temperature-related BxGPCR genes increased at low temperature. And, the survival rates of BxGPCR17454 silenced B. xylophilus revealed a significant decrease at low temperature. Conclusion: in conclusion, this transcriptome-based study revealed a crucial role of BxGPCR17454 in low temperature response process of pine wood nematode. These discoveries would assist the development of management and methods for efficient control of this devastating pine tree pest.

New Zealand Tree and Giant Wētā (Orthoptera) Transcriptomics Reveal Divergent Selection Patterns in Metabolic Loci

Exposure to low temperatures requires an organism to overcome physiological challenges. New Zealand wētā belonging to the genera Hemideina and Deinacrida are found across a wide range of thermal environments and therefore subject to varying selective pressures. Here we assess the selection pressures across the wētā phylogeny, with a particular emphasis on identifying genes under positive or diversifying selection. We used RNA-seq to generate transcriptomes for all 18 Deinacrida and Hemideina species. A total of 755 orthologous genes were identified using a bidirectional best-hit approach, with the resulting gene set encompassing a diverse range of functional classes. Analysis of ortholog ratios of synonymous to nonsynonymous amino acid changes found 83 genes that are under positive selection for at least one codon. A wide variety of Gene Ontology terms, enzymes, and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways are represented among these genes. In particular, enzymes involved in oxidative phosphorylation, melanin synthesis, and free-radical scavenging are represented, consistent with physiological and metabolic changes that are associated with adaptation to alpine environments. Structural alignment of the transcripts with the most codons under positive selection revealed that the majority of sites are surface residues, and therefore have the potential to influence the thermostability of the enzyme, with the exception of prophenoloxidase where two residues near the active site are under selection. These proteins provide interesting candidates for further analysis of protein evolution.

Impact of heat stress on development and fertility of Drosophila suzukii Matsumura (Diptera: Drosophilidae)

Drosophila suzukii is a globally invasive fruit pest that costs millions in yield losses and increased pest management costs. Management practices for D. suzukii currently rely heavily on calendar-based applications of broad-spectrum insecticides, but decision-based applications are theoretically possible with refined population modeling and monitoring. Temperature conditions are strongly deterministic of insect growth rates, fecundity, fertility, and resulting population densities. Therefore, information about the effects of temperature can be incorporated into population modeling to accurately predict D. suzukii population densities in the field which is crucial to maximize pesticide application efficiency and improve sustainability. Here, we investigated the effects of chronic heat stress during all of juvenile development on egg-to-adult viability and fertility. We also investigated egg-to-adult viability under heat stress after heat shock of the maternal parent. We found that heat stress during development results in lower egg-to-adult viability, and reduced lifespan and fertility for surviving adults. However, heat-shock treatment of females prior to egg laying increased the egg-to-adult viability of their eggs under heat stress. Female flies that developed at 30 °C had smaller ovaries than the untreated group and male flies had less sperm in their testes, and no sperm in their seminal vesicles. We conclude that heat stress during development is likely to have negative effect on D. suzukii population dynamics in the field. However, the intensity of such negative impact will depend on the phenotypic state of their maternal parents.

Larvae of Drosophila melanogaster exhibit transcriptional activation of immune response pathways and antimicrobial peptides during recovery from supercooling stress

The biochemical and molecular mechanisms underlying insect cold acclimation prior to cold stress are relatively well explored, but the mechanisms linked to recovery and repair after cold stress have received much less attention. Here we focus on recovery from cold stress in the larvae of the vinegar fly (Drosophila melanogaster) that were exposed to two physiologically distinct cold stress situations: supercooling (S, survival > 95%) and freezing (F, survival < 10%), both at -5 °C. We analysed the metabolic and transcriptomic responses to cold stress via GC-MS/LC-MS and whole-genome microarrays, respectively. Both stresses (S and F) caused metabolic perturbations which were transient in supercooled larvae but deeper and irreversible in frozen larvae. Differential gene expression analysis revealed a clear disparity in responses to supercooling and freezing (less than 10% of DE genes overlapped between S and F larvae). Using GO term enrichment analysis and KEGG pathway mapping, we identified the stimulation of immune response pathways as a strong candidate mechanism for coping with supercooling. Supercooling caused complex transcriptional activation of innate immunity potential: from Lysozyme-mediated degradation of bacterial cell walls, recognition of pathogen signals, through phagocytosis and lysosomal degradation, Toll and Imd signaling, to upregulation of genes coding for different antimicrobial peptides. The transcriptomic response to freezing was instead dominated by degradation of macromolecules and death-related processes such as autophagy and apoptosis. Of the 45 upregulated DE genes overlapping in responses to supercooling and freezing, 26 were broadly ascribable to defense and repair functions.

Novel microRNAs are associated with population divergence in transcriptional response to thermal stress in an intertidal copepod

The role of gene expression in adaptation to differing thermal environments has been assayed extensively. Yet, in most natural systems, analyses of gene expression reveal only one level of the complexity of regulatory machineries. MicroRNAs (miRNAs) are small noncoding RNAs which are key components of many gene regulatory networks, and they play important roles in a variety of cellular pathways by modulating post-transcriptional quantities of mRNA available for protein synthesis. The characterization of miRNA loci and their regulatory dynamics in nonmodel systems are still largely understudied. In this study, we examine the role of miRNAs in response to high thermal stress in the intertidal copepod Tigriopus californicus. Allopatric populations of this species show varying levels of local adaptation with respect to thermal regimes, and previous studies showed divergence in gene expression between populations from very different thermal environments. We examined the transcriptional response to temperature stress in two populations separated by only 8 km by utilizing RNA-seq to quantify both mRNA and miRNA levels. Using the currently available genome sequence, we first describe the repertoire of miRNAs in T. californicus and assess the degree to which transcriptional response to temperature stress is governed by miRNA activity. The two populations showed large differences in the number of genes involved, the magnitude of change in commonly used genes and in the number of miRNAs involved in transcriptional modulation during stress. Our results suggest that an increased level of regulatory network complexity may underlie improved survivorship under thermal stress in one of the populations.

Genomic signatures of experimental adaptive radiation in Drosophila

Abiotic environmental factors play a fundamental role in determining the distribution, abundance and adaptive diversification of species. Empowered by new technologies enabling rapid and increasingly accurate examination of genomic variation in populations, researchers may gain new insights into the genomic background of adaptive radiation and stress resistance. We investigated genomic variation across generations of large-scale experimental selection regimes originating from a single founder population of Drosophila melanogaster, diverging in response to ecologically relevant environmental stressors: heat shock, heat knock down, cold shock, desiccation and starvation. When compared to the founder population, and to parallel unselected controls, there were more than 100,000 single nucleotide polymorphisms (SNPs) displaying consistent allelic changes in response to selective pressures across generations. These SNPs were found in both coding and noncoding sequences, with the highest density in promoter regions, and involved a broad range of functionalities, including molecular chaperoning by heat-shock proteins. The SNP patterns were highly stressor-specific despite considerable variation among line replicates within each selection regime, as reflected by a principal component analysis, and co-occurred with selective sweep regions. Only ~15% of SNPs with putatively adaptive changes were shared by at least two selective regimes, while less than 1% of SNPs diverged in opposite directions. Divergent stressors driving evolution in the experimental system of adaptive radiation left distinct genomic signatures, most pronounced in starvation and heat-shock selection regimes.

Functional characterization of the Hyles euphorbiae hawkmoth transcriptome reveals strong expression of phorbol ester detoxification and seasonal cold hardiness genes

Background

The European spurge hawkmoth, Hyles euphorbiae (Lepidoptera, Sphingidae), has been intensively studied as a model organism for insect chemical ecology, cold hardiness and evolution of species delineation. To understand species isolation mechanisms at a molecular level, this study aims at determining genetic factors underlying two adaptive ecological trait candidates, phorbol ester (TPA) detoxification and seasonal cold acclimation.

Method

A draft transcriptome of H. euphorbiae was generated using Illumina sequencing, providing the first genomic resource for the hawkmoth subfamily Macroglossinae. RNA expression levels in tissues of experimental TPA feeding larvae and cooled pupae was compared to levels in control larvae and pupae using 26 bp RNA sequence tag libraries (DeepSuperSAGE). Differential gene expression was assessed by homology searches of the tags in the transcriptome.

Results

In total, 389 and 605 differentially expressed transcripts for detoxification and cold hardiness, respectively, could be identified and annotated with proteins. The majority (22 of 28) of differentially expressed detox transcripts of the four 'drug metabolism' enzyme groups (cytochrome P450 (CYP), carboxylesterases (CES), glutathione S-transferases (GST) and lipases) are up-regulated. Triacylglycerol lipase was significantly over proportionally annotated among up-regulated detox transcripts. We record several up-regulated lipases, GSTe2, two CESs, CYP9A21, CYP6BD6 and CYP9A17 as candidate genes for further H. euphorbiae TPA detoxification analyses. Differential gene expression of the cold acclimation treatment is marked by metabolic depression with enriched Gene Ontology terms among down-regulated transcripts almost exclusively comprising metabolism, aerobic respiration and dissimilative functions. Down-regulated transcripts include energy expensive respiratory proteins like NADH dehydrogenase, cytochrome oxidase and ATP synthase. Gene expression patterns show shifts in carbohydrate metabolism towards cryoprotectant production. The Glycolysis enzymes, G1Pase, A1e, Gpi and an Akr isoform are up-regulated. Glycerol, an osmolyte which lowers the body liquid supercooling point, appears to be the predominant polyol cryoprotectant in H. euphorbiae diapause pupae. Several protein candidates involved in glucose, glycerol, myo-inositol and potentially sorbitol and trehalose synthesis were identified.

Conclusions

A majority of differently expressed transcripts unique for either detoxification or cold hardiness indicates highly specialized functional adaptation which may have evolved from general cell metabolism and stress response.The transcriptome and extracted candidate biomarkers provide a basis for further gene expression studies of physiological processes and adaptive traits in H. euphorbiae.

A transcriptomics assessment of oxygen-temperature interactions reveals novel candidate genes underlying variation in thermal tolerance and survival

While single stress responses are fairly well researched, multiple, interactive stress responses are not-despite the obvious importance thereof. Here, using D. melanogaster, we investigated the effects of simultaneous exposures to low O₂ (hypoxia) and varying thermal conditions on mortality rates, estimates of thermal tolerance and the transcriptome. We used combinations of 21 (normoxia), 10 or 5kPa O₂ with control (23°C), cold (4°C) or hot (31°C) temperature exposures before assaying chill coma recovery time (CCRT) and heat knock down time (HKDT) as measures of cold and heat tolerance respectively. We found that mortality was significantly affected by temperature, oxygen partial pressure (PO₂) and the interaction between the two. Cold treatments resulted in low mortality (<5%), regardless of PO₂ treatment; while hot treatments resulted in higher mortality (∼20%), especially at 5kPa O₂ which was lethal for most flies (∼80%). Both CCRT and HKDT were significantly affected by temperature, but not PO₂, of the treatments, and the interaction of temperature and PO₂ was non-significant. Hot treatments led to significantly longer CCRT, and shorter HKDT in comparison to cold treatments. Global gene expression profiling provided the first transcriptome level response to the combined stress of PO₂ and temperature, showing that stressful treatments resulted in higher mortality and induced transcripts that were associated with protein kinases, catabolic processes (proteases, hydrolases, peptidases) and membrane function. Several genes and pathways that may be responsible for the protective effects of combined PO₂ and cold treatments were identified. We found that urate oxidase was upregulated in all three cold treatments, regardless of the PO₂. Small heat shock proteins Hsp22 and Hsp23 were upregulated after both 10 and 21kPa O₂-hot treatments. Collectively, the data from PO₂-hot treatments suggests that hypoxia does exacerbate heat stress, through an as yet unidentified mechanism. Hsp70B and an unannotated transcript (CG6733) were significantly differentially expressed after 5kPa O₂-cold and 10kPa O₂-hot treatments relative to their controls. Downregulation of these transcripts was correlated with reduced thermal tolerance (longer CCRT and shorter HKDT), suggesting that these genes may be important candidates for future research.

Identification and expression analysis of multiple small heat shock protein genes in spruce budworm, Choristoneura fumiferana (L.)

Fifteen small heat shock protein (sHSP) genes were identified from spruce budworm, Choristoneura fumiferana (L.), an important native forest pest in North America. The transcript levels of each CfHSP were measured under non-stress conditions in all life stages from egg to adult and in five different larval tissues. CfHSP transcript levels showed variation during development, with highest levels in adults and lowest in eggs. Most CfHSP transcripts are highly expressed in larval fat body and Malpighian tubules; two CfHSPs display extremely high expression in the head and epidermis. Upon heat stress, nine CfHSP genes are significantly upregulated, increasing by 50- to 2500-fold depending on developmental stage and tissue type. Upon starvation, eight CfHSPs are upregulated or downregulated, whereas six others retain constant expression. These results suggest that CfHSPs have important and multiple roles in spruce budworm development and in response to heat stress and starvation.

Differential transcriptome analysis reveals genes related to cold tolerance in seabuckthorn carpenter moth, Eogystia hippophaecolus

Seabuckthorn carpenter moth, Eogystia hippophaecolus (Lepidoptera: Cossidae), is an important pest of sea buckthorn (Hippophae rhamnoides), which is a shrub that has significant ecological and economic value in China. E. hippophaecolus is highly cold tolerant, but limited studies have been conducted to elucidate the molecular mechanisms underlying its cold resistance. Here we sequenced the E. hippophaecolus transcriptome using RNA-Seq technology and performed de novo assembly from the short paired-end reads. We investigated the larval response to cold stress by comparing gene expression profiles between treatments. We obtained 118,034 unigenes, of which 22,161 were annotated with gene descriptions, conserved domains, gene ontology terms, and metabolic pathways. These resulted in 57 GO terms and 193 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. By comparing transcriptome profiles for differential gene expression, we identified many differentially expressed proteins and genes, including heat shock proteins and cuticular proteins which have previously been reported to be involved in cold resistance of insects. This study provides a global transcriptome analysis and an assessment of differential gene expression in E. hippophaecolus under cold stress. We found seven differential expressed genes in common between developmental stages, which were verified with qPCR. Our findings facilitate future genomic studies aimed at improving our understanding of the molecular mechanisms underlying the response of insects to low temperatures.

Metabolic insights into the cold survival strategy and overwintering of the common cutworm, Spodoptera litura (Fabricius) (Lepidoptera: Noctuidae)

The common cutworm, Spodoptera litura (Fabricius) (Lepidoptera: Noctuidae), is a destructive pest in Asia. Although overwintering in the field has not been reported for this species, their larvae are capable of long-term survival in fluctuating temperatures, i.e., 5°C (12h) plus 13°C (12h), if food is available. With an increase in climate change due to global warming and the widespread use of greenhouses, further understanding of their cold survival strategy is needed to predict and control their population in the future. In this study, metabolomics was performed to analyze the metabolic features of S. litura larvae exposed to two typical low temperatures: 15°C and 4°C, at which the development, locomotion and feeding activities are maintained or halted, respectively. The results showed that the strategies that regulate lipid and amino acid metabolism were similar at 15°C and 4°C. Cold exposure induced a metabolic shift of energy from carbohydrate to lipid and decreased free amino acids level. Biosynthesis likely contributed to the decrease in amino acids levels even at 4°C, a non-feeding temperature, suggesting an insufficient suppression of anabolism. This explains why food and high temperature pulses are necessary for their long-term cold survival. Glycometabolism was different between 15°C and 4°C. Carbohydrates were used rapidly at 15°C, while trehalose accumulated at 4°C. Interestingly, abundant trehalose and serine are prominent features of Spodoptera exigua larvae, an overwintering species, when compared to S. litura larvae. Exposure to 4°C also induced up-regulation of carbohydrase and protease in the guts of S. litura. Therefore, it is likely that concurrence of food supplement and fluctuating temperatures could facilitate the cold survival of S. litura larvae. We also found that exposure to 4°C could activate the mevalonate pathway in S. litura larvae, which might be related to glycometabolism at 4°C. Overall, our study describes systematically the responses of a cold susceptible insect, S. litura, to low temperatures and explains how fluctuating temperatures facilitate their long-term cold survival indicating the possibility for overwintering of S. litura larvae with global warming and agricultural reforms.

Transcriptomic analysis to uncover genes affecting cold resistance in the Chinese honey bee (Apis cerana cerana)

The biological activity and geographical distribution of honey bees is strongly temperature-dependent, due to their ectothermic physiology. In China, the endemic Apis cerana cerana exhibits stronger cold hardiness than Western honey bees, making the former species important pollinators of winter-flowering plants. Although studies have examined behavioral and physiological mechanisms underlying cold resistance in bees, data are scarce regarding the exact molecular mechanisms. Here, we investigated gene expression in A. c. cerana under two temperature treatments, using transcriptomic analysis to identify differentially expressed genes (DEGs) and relevant biological processes, respectively. Across the temperature treatments, 501 DEGs were identified. A gene ontology analysis showed that DEGs were enriched in pathways related to sugar and amino acid biosynthesis and metabolism, as well as calcium ion channel activity. Additionally, heat shock proteins, zinc finger proteins, and serine/threonine-protein kinases were differentially expressed between the two treatments. The results of this study provide a general digital expression profile of thermoregulation genes responding to cold hardiness in A. c. cerana. Our data should prove valuable for future research on cold tolerance mechanisms in insects, and may be beneficial in breeding efforts to improve bee hardiness.

Effects of cold-acclimation on gene expression in Fall field cricket (Gryllus pennsylvanicus) ionoregulatory tissues

Background

Cold tolerance is a key determinant of temperate insect distribution and performance. Chill-susceptible insects lose ion and water homeostasis during cold exposure, but prior cold acclimation improves both cold tolerance and defense of homeostasis. The mechanisms underlying these processes are mostly unknown; cold acclimation is thought to enhance ion transport in the cold and/or prevent leak of water and ions. To identify candidate mechanisms of cold tolerance plasticity we generated transcriptomes of ionoregulatory tissues (hindgut and Malpighian tubules) from Gryllus pennsylvanicus crickets and compared gene expression in warm- and cold-acclimated individuals.

Results

We assembled a G. pennsylvanicus transcriptome de novo from 286 million 50-bp reads, yielding 70,037 contigs (~44% of which had putative BLAST identities). We compared the transcriptomes of warm- and cold-acclimated hindguts and Malpighian tubules. Cold acclimation led to a ≥ 2-fold change in the expression of 1493 hindgut genes (733 downregulated, 760 upregulated) and 2008 Malpighian tubule genes (1009 downregulated, 999 upregulated). Cold-acclimated crickets had altered expression of genes putatively associated with ion and water balance, including: a downregulation of V-ATPase and carbonic anhydrase in the Malpighian tubules and an upregulation of Na⁺-K⁺ ATPase in the hindgut. We also observed acclimation-related shifts in the expression of cytoskeletal genes in the hindgut, including actin and actin-anchoring/stabilizing proteins, tubulin, α-actinin, and genes involved in adherens junctions organization. In both tissues, cold acclimation led to differential expression of genes encoding cytochrome P450s, glutathione-S-transferases, apoptosis factors, DNA repair, and heat shock proteins.

Conclusions

This is the first G. pennsylvanicus transcriptome, and our tissue-specific approach yielded new candidate mechanisms of cold tolerance plasticity. Cold acclimation may reduce loss of hemolymph volume in the cold by 1) decreasing primary urine production via reduced expression of carbonic anhydrase and V-ATPase in the Malpighian tubules and 2) by increasing Na⁺ (and therefore water) reabsorption across the hindgut via increase in Na⁺-K⁺ ATPase expression. Cold acclimation may reduce chilling injury by remodeling and stabilizing the hindgut epithelial cytoskeleton and cell-to-cell junctions, and by increasing the expression of genes involved in DNA repair, detoxification, and protein chaperones.

Electronic supplementary material

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

Evolutionary adaptation to environmental stressors: a common response at the proteomic level

Mechanistic trade-offs between traits under selection can shape and constrain evolutionary adaptation to environmental stressors. However, our knowledge of the quantitative and qualitative overlap in the molecular machinery among stress tolerance traits is highly restricted by the challenges of comparing and interpreting data between separate studies and laboratories, as well as to extrapolating between different levels of biological organization. We investigated the expression of the constitutive proteome (833 proteins) of 35 Drosophila melanogaster replicate populations artificially selected for increased resistance to six different environmental stressors. The evolved proteomes were significantly differentiated from replicated control lines. A targeted analysis of the constitutive proteomes revealed a regime-specific selection response among heat-shock proteins, which provides evidence that selection also adjusts the constitutive expression of these molecular chaperones. Although the selection response in some proteins was regime specific, the results were dominated by evidence for a "common stress response." With the exception of high temperature survival, we found no evidence for negative correlations between environmental stress resistance traits, meaning that evolutionary adaptation is not constrained by mechanistic trade-offs in regulation of functional important proteins. Instead, standing genetic variation and genetic trade-offs outside regulatory domains likely constrain the evolutionary responses in natural populations.

CRISPR-induced null alleles show that Frost protects Drosophila melanogaster reproduction after cold exposure

The ability to survive and reproduce after cold exposure is important in all kingdoms of life. However, even in a sophisticated genetic model system like Drosophila melanogaster, few genes have been identified as functioning in cold tolerance. The accumulation of the Frost (Fst) gene transcript increases after cold exposure, making it a good candidate for a gene that has a role in cold tolerance. However, despite extensive RNAi knockdown analysis, no role in cold tolerance has been assigned to Fst. CRISPR is an effective technique for completely knocking down genes, and less likely to produce off-target effects than GAL4-UAS RNAi systems. We have used CRISPR-mediated homologous recombination to generate Fst null alleles, and these Fst alleles uncovered a requirement for FST protein in maintaining female fecundity following cold exposure. However, FST does not have a direct role in survival following cold exposure. FST mRNA accumulates in the Malpighian tubules, and the FST protein is a highly disordered protein with a putative signal peptide for export from the cell. Future work is needed to determine whether FST is exported from the Malpighian tubules and directly interacts with female reproductive tissues post-cold exposure, or if it is required for other repair/recovery functions that indirectly alter energy allocation to reproduction.

Multiple Independent Retroelement Insertions in the Promoter of a Stress Response Gene Have Variable Molecular and Functional Effects in Drosophila

Promoters are structurally and functionally diverse gene regulatory regions. The presence or absence of sequence motifs and the spacing between the motifs defines the properties of promoters. Recent alternative promoter usage analyses in Drosophila melanogaster revealed that transposable elements significantly contribute to promote diversity. In this work, we analyzed in detail one of the transposable element insertions, named FBti0019985, that has been co-opted to drive expression of CG18446, a candidate stress response gene. We analyzed strains from different natural populations and we found that besides FBti0019985, there are another eight independent transposable elements inserted in the proximal promoter region of CG18446. All nine insertions are solo-LTRs that belong to the roo family. We analyzed the sequence of the nine roo insertions and we investigated whether the different insertions were functionally equivalent by performing 5'-RACE, gene expression, and cold-stress survival experiments. We found that different insertions have different molecular and functional consequences. The exact position where the transposable elements are inserted matters, as they all showed highly conserved sequences but only two of the analyzed insertions provided alternative transcription start sites, and only the FBti0019985 insertion consistently affects CG18446 expression. The phenotypic consequences of the different insertions also vary: only FBti0019985 was associated with cold-stress tolerance. Interestingly, the only previous report of transposable elements inserting repeatedly and independently in a promoter region in D. melanogaster, were also located upstream of a stress response gene. Our results suggest that functional validation of individual structural variants is needed to resolve the complexity of insertion clusters.

Canalization of gene expression is a major signature of regulatory cold adaptation in temperate Drosophila melanogaster

BACKGROUND

Transcriptome analysis may provide means to investigate the underlying genetic causes of shared and divergent phenotypes in different populations and help to identify potential targets of adaptive evolution. Applying RNA sequencing to whole male Drosophila melanogaster from the ancestral tropical African environment and a very recently colonized cold-temperate European environment at both standard laboratory conditions and following a cold shock, we seek to uncover the transcriptional basis of cold adaptation.

RESULTS

In both the ancestral and the derived populations, the predominant characteristic of the cold shock response is the swift and massive upregulation of heat shock proteins and other chaperones. Although we find ~25 % of the genome to be differentially expressed following a cold shock, only relatively few genes (n = 16) are up- or down-regulated in a population-specific way. Intriguingly, 14 of these 16 genes show a greater degree of differential expression in the African population. Likewise, there is an excess of genes with particularly strong cold-induced changes in expression in Africa on a genome-wide scale.

CONCLUSIONS

The analysis of the transcriptional cold shock response most prominently reveals an upregulation of components of a general stress response, which is conserved over many taxa and triggered by a plethora of stressors. Despite the overall response being fairly similar in both populations, there is a definite excess of genes with a strong cold-induced fold-change in Africa. This is consistent with a detrimental deregulation or an overshooting stress response. Thus, the canalization of European gene expression might be responsible for the increased cold tolerance of European flies.

Thermal fluctuations affect the transcriptome through mechanisms independent of average temperature

Terrestrial ectotherms are challenged by variation in both mean and variance of temperature. Phenotypic plasticity (thermal acclimation) might mitigate adverse effects, however, we lack a fundamental understanding of the molecular mechanisms of thermal acclimation and how they are affected by fluctuating temperature. Here we investigated the effect of thermal acclimation in Drosophila melanogaster on critical thermal maxima (CTmax) and associated global gene expression profiles as induced by two constant and two ecologically relevant (non-stressful) diurnally fluctuating temperature regimes. Both mean and fluctuation of temperature contributed to thermal acclimation and affected the transcriptome. The transcriptomic response to mean temperatures comprised modification of a major part of the transcriptome, while the response to fluctuations affected a much smaller set of genes, which was highly independent of both the response to a change in mean temperature and to the classic heat shock response. Although the independent transcriptional effects caused by fluctuations were relatively small, they are likely to contribute to our understanding of thermal adaptation. We provide evidence that environmental sensing, particularly phototransduction, is a central mechanism underlying the regulation of thermal acclimation to fluctuating temperatures. Thus, genes and pathways involved in phototransduction are likely of importance in fluctuating climates.

Cold acclimation wholly reorganizes the Drosophila melanogaster transcriptome and metabolome

Cold tolerance is a key determinant of insect distribution and abundance, and thermal acclimation can strongly influence organismal stress tolerance phenotypes, particularly in small ectotherms like Drosophila. However, there is limited understanding of the molecular and biochemical mechanisms that confer such impressive plasticity. Here, we use high-throughput mRNA sequencing (RNA-seq) and liquid chromatography – mass spectrometry (LC-MS) to compare the transcriptomes and metabolomes of D. melanogaster acclimated as adults to warm (rearing) (21.5 °C) or cold conditions (6 °C). Cold acclimation improved cold tolerance and led to extensive biological reorganization: almost one third of the transcriptome and nearly half of the metabolome were differentially regulated. There was overlap in the metabolic pathways identified via transcriptomics and metabolomics, with proline and glutathione metabolism being the most strongly-supported metabolic pathways associated with increased cold tolerance. We discuss several new targets in the study of insect cold tolerance (e.g. dopamine signaling and Na⁺-driven transport), but many previously identified candidate genes and pathways (e.g. heat shock proteins, Ca²⁺ signaling, and ROS detoxification) were also identified in the present study, and our results are thus consistent with and extend the current understanding of the mechanisms of insect chilling tolerance.

Reproductive arrest and stress resistance in winter-acclimated Drosophila suzukii

Overwintering insects must survive the multiple-stress environment of winter, which includes low temperatures, reduced food and water availability, and cold-active pathogens. Many insects overwinter in diapause, a developmental arrest associated with high stress tolerance. Drosophila suzukii (Diptera: Drosophilidae), spotted wing drosophila, is an invasive agricultural pest worldwide. Its ability to overwinter and therefore establish in temperate regions could have severe implications for fruit crop industries. We demonstrate here that laboratory populations of Canadian D. suzukii larvae reared under short-day, low temperature, conditions develop into dark 'winter morph' adults similar to those reported globally from field captures, and observed by us in southern Ontario, Canada. These winter-acclimated adults have delayed reproductive maturity, enhanced cold tolerance, and can remain active at low temperatures, although they do not have the increased desiccation tolerance or survival of fungal pathogen challenges that might be expected from a more heavily melanised cuticle. Winter-acclimated female D. suzukii have underdeveloped ovaries and altered transcript levels of several genes associated with reproduction and stress. While superficially indicative of reproductive diapause, the delayed reproductive maturity of winter-acclimated D. suzukii appears to be temperature-dependent, not regulated by photoperiod, and is thus unlikely to be 'true' diapause. The traits of this 'winter morph', however, likely facilitate overwintering in southern Canada, and have probably contributed to the global success of this fly as an invasive species.

PROTEOMICS ANALYSIS OF OVEREXPRESSED PLASMA PROTEINS IN RESPONSE TO COLD ACCLIMATION IN Ostrinia furnacalis

Many insects in temperate regions overwinter in diapause. In these insects, one of the metabolic adaptations to cold stress is the synthesis of responsive proteins. Using proteomic analysis, an investigation aimed to a better understanding of the molecular adaptation mechanisms to cold stress was carried out in Ostrinia furnacalis larva. Proteins were extracted from the larval hemolymph collected from both control and overwintering larva. By polyethylene glycol precipitation, approximately 560 protein spots were separated and visualized on two-dimensional (2D) gels after silver staining. Eighteen protein spots were found to be upregulated in overwinter larval plasma in different patterns. As an initial work, 13 of these proteins were identified using MALDI TOF/TOF MS. The differentially overexpressed proteins include heat shock 70 kDa cognate protein, small heat shock protein (sHSP), putative aliphatic nitrilase, arginine kinase, phosphoglyceromutase, triosephosphateisomerase, and glutathione transferase. Alterations in the levels of these proteins were further confirmed by qPCR. This study is the first analysis of differentially expressed plasma proteins in O. furnacalis diapause larvae under extremely low temperature conditions and gives new insights into the acclimation mechanisms responsive to cold stress. Our results also support the idea that energy metabolism, alanine and proline metabolism, and antioxidative reaction act in the cold acclimation of O. furnacalis diapause larvae.

Expression profiling of heat shock genes in a scuttle fly Megaselia scalaris (Diptera, Phoridae)

The Phoridae are a family of small, hump-backed flies, dominating in post-fire areas. Some of these flies are probably able to survive a fire as an egg, larva, or pupa, and may be adapted to the fire-altered environment at the genomic level. In this study, we describe the influence of short-term temperature treatment on the expression of seven heat shock protein genes in the third-instar larvae and imagoes of a scuttle fly Megaselia scalaris-one of the cosmopolitan and polyphagous phorids. In terms of the response to temperature treatment, these genes tested against tubulin as a reference split into three classes. The first class consists of hsp22 (larvae), hsp23 (larvae), and hsp26 (both larvae and imagoes), and is upregulated at the lowest temperature (33°C). The second class consists of hsp22 (imagoes), hsp23 (imagoes), hsp40 (larvae and imagoes), and hsp70 (larvae and imagoes), and is upregulated or induced at 37°C. Expression of genes of the third class (hsp27 and hsp83-larvae and imagoes) increased after treatment at 41°C temperature. Expression of the first two classes of genes occurred at a temperature lower than the LT₅₀ of larvae and imagoes. The fact that there is a gap between the temperature upregulating hsp genes and the temperature leading to the loss of viability suggests that not only the level of hsp gene expression but also the temperature at which gene expression increased is important in an adaptation of M. scalaris to harsh environment. J. Exp. Zool. 323A: 704-713, 2015. © 2015 Wiley Periodicals, Inc.

The novel ABC transporter ABCH1 is a potential target for RNAi-based insect pest control and resistance management

Insect pests cause serious crop damage and develop high-level resistance to chemical insecticides and Bacillus thuringiensis (Bt) insecticidal Cry toxins. A new promising approach for controlling them and overcoming this resistance is RNA interference (RNAi). The RNAi-based insect control strategy depends on the selection of suitable target genes. In this study, we cloned and characterized a novel ABC transporter gene PxABCH1 in diamondback moth, Plutella xylostella (L.). Phylogenetic analysis showed that PxABCH1 is closely related to ABCA and ABCG subfamily members. Spatial-temporal expression detection revealed that PxABCH1 was expressed in all tissues and developmental stages, and highest expressed in head and male adult. Midgut sequence variation and expression analyses of PxABCH1 in all the susceptible and Bt-resistant P. xylostella strains and the functional analysis by sublethal RNAi demonstrated that Cry1Ac resistance was independent of this gene. Silencing of PxABCH1 by a relatively high dose of dsRNA dramatically reduced its expression and resulted in larval and pupal lethal phenotypes in both susceptible and Cry1Ac-resistant P. xylostella strains. To our knowledge, this study provides the first insight into ABCH1 in lepidopterans and reveals it as an excellent target for RNAi-based insect pest control and resistance management.

Constraints, independence, and evolution of thermal plasticity: probing genetic architecture of long- and short-term thermal acclimation

Seasonal and daily thermal variation can limit species distributions because of physiological tolerances. Low temperatures are particularly challenging for ectotherms, which use both basal thermotolerance and acclimation, an adaptive plastic response, to mitigate thermal stress. Both basal thermotolerance and acclimation are thought to be important for local adaptation and persistence in the face of climate change. However, the evolutionary independence of basal and plastic tolerances remains unclear. Acclimation can occur over longer (seasonal) or shorter (hours to days) time scales, and the degree of mechanistic overlap is unresolved. Using a midlatitude population of Drosophila melanogaster, we show substantial heritable variation in both short- and long-term acclimation. Rapid cold hardening (short-term plasticity) and developmental acclimation (long-term plasticity) are positively correlated, suggesting shared mechanisms. However, there are independent components of these traits, because developmentally acclimated flies respond positively to short-term acclimation. A strong negative correlation between basal cold tolerance and developmental acclimation suggests that basal cold tolerance may constrain developmental acclimation, whereas a weaker negative correlation between basal cold tolerance and short-term acclimation suggests less constraint. Using genome-wide association mapping, we show the genetic architecture of rapid cold hardening and developmental acclimation responses are nonoverlapping at the SNP and corresponding gene level. However, genes associated with each trait share functional similarities, including genes involved in apoptosis and autophagy, cytoskeletal and membrane structural components, and ion binding and transport. These results indicate substantial opportunity for short-term and long-term acclimation responses to evolve separately from each other and for short-term acclimation to evolve separately from basal thermotolerance.