Upregulation of two actin genes and redistribution of actin during diapause and cold stress in the northern house mosquito, Culex pipiens

Transcription profiling reveals tissue-specific metabolic pathways in the fat body and ovary of the diapausing mosquito Culex pipiens

The northern house mosquito, Culex pipiens, employs diapause as an essential survival strategy during winter, inducing important phenotypic changes such as enhanced stress tolerance, lipid accumulation, and extended longevity. During diapause, the cessation of reproductive development represents another distinctive phenotypic change, underlining the need for adjusted modulation of gene expressions within the ovary. Although considerable advancements in screening gene expression profiles in diapausing and non-diapausing mosquitoes, there remains a gap in tissue-specific transcriptomic profiling that could elucidate the complicated formation of diverse diapause features in Cx. pipiens. Here, we filled this gap by utilizing RNA sequencing, providing a detailed examination of gene expression patterns in the fat body and ovary during diapause compared to non-diapause conditions. Functional annotation of upregulated genes identified associations with carbohydrate metabolism, stress tolerance, immunity, and epigenetic regulation. The validation of candidate genes using quantitative real-time PCR verified the differentially expressed genes identified in diapausing mosquitoes. Our findings contribute novel insights into potential regulators during diapause in Cx. pipiens, thereby opening possible avenues for developing innovative vector control strategies.

Variation in Thermal Sensitivity of Diapause Development among Individuals and over Time Predicts Life History Timing in a Univoltine Insect

AbstractPhysiological time is important for understanding the development and seasonal timing of ectothermic animals but has largely been applied to developmental processes that occur during spring and summer, such as morphogenesis. There is a substantial knowledge gap in the relationship between temperature and development during winter, a season that is increasingly impacted by climate change. Most temperate insects overwinter in diapause, a developmental process with little obvious morphological change. We used principles from the physiological time literature to measure and model the thermal sensitivity of diapause development rate in the apple maggot fly Rhagoletis pomonella, a univoltine fly whose diapause duration varies substantially within and among populations. We show that diapause duration can be predicted by modeling a relationship between temperature and development rate that is shifted toward lower temperatures compared with typical models of morphogenic, nondiapause development. However, incorporating interindividual variation and ontogenetic variation in the temperature-to-development rate relationship was critical for accurately predicting fly emergence, as diapause development proceeded more quickly at high temperatures later in diapause. We conclude that the conceptual framework may be flexibly applied to other insects and discuss possible mechanisms of diapause timers and implications for phenology with warming winters.

A rapid return to normal: temporal gene expression patterns following cold exposure in the bumble bee Bombus impatiens

Bumble bees are common in cooler climates and many species likely experience periodic exposure to very cold temperatures, but little is known about the temporal dynamics of cold response mechanisms following chill exposure, especially how persistent effects of cold exposure may facilitate tolerance of future events. To investigate molecular processes involved in the temporal response by bumble bees to acute cold exposure, we compared mRNA transcript abundance in Bombus impatiens workers exposed to 0°C for 75 min (inducing chill coma) and control bees maintained at a constant ambient temperature (28°C). We sequenced the 3' end of mRNA transcripts (TagSeq) to quantify gene expression in thoracic tissue of bees at several time points (0, 10, 30, 120 and 720 min) following cold exposure. Significant differences from control bees were only detectable within 30 min after the treatment, with most occurring at the 10 min recovery time point. Genes associated with gluconeogenesis and glycolysis were most notably upregulated, while genes related to lipid and purine metabolism were downregulated. The observed patterns of expression indicate a rapid recovery after chill coma, suggesting an acute differential transcriptional response during recovery from chill coma and return to baseline expression levels within an hour, with no long-term gene expression markers of this cold exposure. Our work highlights the functions and pathways important for acute cold recovery, provides an estimated time frame for recovery from cold exposure in bumble bees, and suggests that cold hardening may be less important for these heterothermic insects.

De novo genome assembly of the invasive mosquito species Aedes japonicus and Aedes koreicus

BACKGROUND

Recently, two invasive Aedes mosquito species, Ae. japonicus and Ae. koreicus, are circulating in several European countries posing potential health risks to humans and animals. Vector control is the main option to prevent mosquito-borne diseases, and an accurate genome sequence of these mosquitoes is essential to better understand their biology and to develop effective control strategies.

METHODS

A de novo genome assembly of Ae. japonicus (Ajap1) and Ae. koreicus (Akor1) has been produced based on a hybrid approach that combines Oxford Nanopore long-read and Illumina short-read data. Their quality was ascertained using various metrics. Masking of repetitive elements, gene prediction and functional annotation was performed.

RESULTS

Sequence analysis revealed a very high presence of repetitive DNA and, among others, thermal adaptation genes and insecticide-resistance genes. Through the RNA-seq analysis of larvae and adults of Ae. koreicus and Ae. japonicus exposed to different temperatures, we also identified genes showing a differential temperature-dependent activation.

CONCLUSIONS

The assembly of Akor1 and Ajap1 genomes constitutes the first updated collective knowledge of the genomes of both mosquito species, providing the possibility of understanding key mechanisms of their biology such as the ability to adapt to harsh climates and to develop insecticide-resistance mechanisms.

Multi-locus genomic signatures of local adaptation to snow across the landscape in California populations of a willow leaf beetle

Organisms living in mountains contend with extreme climatic conditions, including short growing seasons and long winters with extensive snow cover. Anthropogenic climate change is driving unprecedented, rapid warming of montane regions across the globe, resulting in reduced winter snowpack. Loss of snow as a thermal buffer may have serious consequences for animals overwintering in soil, yet little is known about how variability in snowpack acts as a selective agent in montane ecosystems. Here, we examine genomic variation in California populations of the leaf beetle Chrysomela aeneicollis, an emerging natural model system for understanding how organisms respond to climate change. We used a genotype-environment association approach to identify genomic signatures of local adaptation to microclimate in populations from three montane regions with variable snowpack and a coastal region with no snow. We found that both winter-associated environmental variation and geographical distance contribute to overall genomic variation across the landscape. We identified non-synonymous variation in novel candidate loci associated with cytoskeletal function, ion transport and membrane stability, cellular processes associated with cold tolerance in other insects. These findings provide intriguing evidence that variation in snowpack imposes selective gradients in montane ecosystems.

Suppression of the gene encoding PDZ domain-containing protein decreases cold tolerance and overwintering survival of the mosquito, Culex pipiens (Culicidae: Diptera)

In diapausing mosquitoes, cold tolerance and prolonged lifespan are important features that are crucial for overwintering success. In the mosquito Culex pipiens, we suggest that PDZ domain-containing protein (PDZ) (post synaptic density protein [PSD95], drosophila disc large tumor suppressor [Dlg1], and zonula occludens-1 protein [zo-1]) domain-containing protein is involved with these diapause features for overwintering survival in Culex mosquitoes. The expression level of pdz was significantly higher in diapausing adult females in the early stage in comparison to their nondiapausing counterparts. Suppression of the gene that encodes PDZ by RNA interference significantly decreased actin accumulation in the midgut of early-stage adult diapausing females. Inhibition of pdz also significantly reduced the survivability of diapausing females which indicates that this protein could play a key role in preserving the midgut tissues during early diapause.

Protein regulation mechanism of cold tolerance in Haemaphysalis longicornis

Ticks are external parasitic arthropods that can transmit a variety of pathogens by sucking blood. Low-temperature tolerance is essential for ticks to survive during the cold winter. Exploring the protein regulation mechanism of low-temperature tolerance of Haemaphysalis longicornis could help to explain how ticks survive in winter. In this study, the quantitative proteomics of several tissues of H. longicornis exposed to low temperature were studied by data independent acquisition technology. Totals of 3 699, 3 422, and 1 958 proteins were identified in the salivary gland, midgut, and ovary, respectively. The proteins involved in energy metabolism, cell signal transduction, protein synthesis and repair, and cytoskeleton synthesis changed under low-temperature stress. The comprehensive analysis of the protein regulation of multiple tissues of female ticks exposed to low temperature showed that maintaining cell homeostasis, maintaining cell viability, and enhancing cell tolerance were the most important means for ticks to maintain vital signs under low temperature. The expression of proteins involved in and regulating the above cell activities was the key to the survival of ticks under low temperatures. Through the analysis of a large amount of data, we found that the expression levels of arylamine N-acetyltransferase, inositol polyphosphate multikinase, and dual-specificity phosphatase were up-regulated under low temperature. We speculated that they might have important significance in low-temperature tolerance. Then, we performed RNA interference on the mRNA of these 3 proteins, and the results showed that the ability of female ticks to tolerate low temperatures decreased significantly.

Molecular Mechanisms of Winter Survival

Winter provides many challenges for insects, including direct injury to tissues and energy drain due to low food availability. As a result, the geographic distribution of many species is tightly coupled to their ability to survive winter. In this review, we summarize molecular processes associated with winter survival, with a particular focus on coping with cold injury and energetic challenges. Anticipatory processes such as cold acclimation and diapause cause wholesale transcriptional reorganization that increases cold resistance and promotes cryoprotectant production and energy storage. Molecular responses to low temperature are also dynamic and include signaling events during and after a cold stressor to prevent and repair cold injury. In addition, we highlight mechanisms that are subject to selection as insects evolve to variable winter conditions. Based on current knowledge, despite common threads, molecular mechanisms of winter survival vary considerably across species, and taxonomic biases must be addressed to fully appreciate the mechanistic basis of winter survival across the insect phylogeny.

Transcriptional Regulation of Reproductive Diapause in the Convergent Lady Beetle, Hippodamia convergens

Simple Summary

Diapause is a dormant period typically controlled by daylength that ensures an insect’s survival through harsh environmental conditions. The convergent lady beetle, Hippodamia convergens, undergoes a reproductive diapause in winter, where female ovaries remain immature and no eggs are laid. This species is an important biological control agent, but during diapause, beetles are less likely to eat pest insects. Thus, knowledge of diapause mechanisms may facilitate manipulation thereof to improve biological control. Further, molecular studies of adult diapause and diapause in Coleoptera are relatively lacking. Here, we assembled and annotated a transcriptome for this species and quantified transcript expression changes during diapause. Female beetles were sampled at three times in diapause (early, mid, and late diapause), which allowed us to characterize the molecular processes occurring at distinct transitions throughout diapause. We found that transcripts involved in flight were consistently upregulated during diapause, which is consistent with dispersal flights at this stage, while transcripts involved in ovarian development were downregulated, which is consistent with the shutdown of reproduction in diapausing females. These findings identify key regulators of diapause in H. convergens and contribute to a growing body of literature on the molecular mechanisms of diapause across the insect phylogeny.

Abstract

Diapause is an alternate development program that synchronizes an insect’s life cycle with seasonally abundant resources and ensures survival in unfavorable conditions. The physiological basis of diapause has been well characterized, but the molecular mechanisms regulating it are still being elucidated. Here, we present a de novo transcriptome and quantify transcript expression during diapause in the convergent lady beetle Hippodamia convergens. H. convergens is used as an augmentative biocontrol agent, and adult females undergo reproductive diapause that is regulated by photoperiod. We sampled females at three stages (early, mid, and late diapause) and compared transcript expression to non-diapausing individuals. Based on principle component analysis, the transcriptomes of diapausing beetles were distinct from non-diapausing beetles, and the three diapausing points tended to cluster together. However, there were still classes of transcripts that differed in expression across distinct phases of diapause. In general, transcripts involved in muscle function and flight were upregulated during diapause, likely to support dispersal flights that occur during diapause, while transcripts involved in ovarian development were downregulated. This information could be used to improve biological control by manipulating diapause. Additionally, our data contribute to a growing understanding of the genetic regulation of diapause across diverse insects.

Transcriptome analysis uncovers different avenues for manipulating cold performance in Chrysomya megacephala (Diptera, Calliphoridae)

Temperature strongly impacts the rates of physiological and biochemical processes, which in turn can determine the survival and population size of insects. At low temperatures performance is limited, however, cold tolerance and performance at low temperature can be improved after short- or long-term acclimation in many insect species. To understand mechanisms underlying acclimation, we sequenced and compared the transcriptome of the blowfly Chrysomya megacephala under rapid cold hardening (RCH) and long-term cold acclimation (LCA) conditions. The RCH response was dominated by genes related to immune response, spliceosome, and protein processing in endoplasmic reticulum with up-regulation during recovery from RCH. In contrast, LCA was associated with genes related to carbohydrate metabolism and cytoskeleton branching and stabilizing. Meanwhile, mRNA levels of genes related to glycerophospholipid metabolism, and some heat shock proteins (Hsps) were collectively up-regulated by both RCH and LCA. There were more genes and pathway adjustments associated with LCA than RCH. Overall, the transcriptome data provide basic information of molecular mechanisms underpinning the RCH and LCA response. The partly independent molecular responses to RCH and LCA suggest that several avenues for manipulating cold performance exist and RCH might be more effective as it only triggers fewer genes and affects the general metabolisms less. These observations provide some appropriate methods to improve cold tolerance of C. megacephala, and hold promise for developing an extended use of mass-reared C. megacephala with better cold performance as a pollinator of crops at low temperatures.

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.

Functional implication of heat shock protein 70/90 and tubulin in cold stress of Dermacentor silvarum

Background

The tick Dermacentor silvarum Olenev (Acari: Ixodidae) is a vital vector tick species mainly distributed in the north of China and overwinters in the unfed adult stage. The knowledge of the mechanism that underlies its molecular adaptation against cold is limited. In the present study, genes of hsp70 and hsp90 cDNA, named Dshsp70 and Dshsp90, and tubulin were cloned and characterized from D. silvarum, and their functions in cold stress were further evaluated.

Methods

The genome of the heat shock proteins and tubulin of D. silvarum were sequenced and analyzed using bioinformatics methods. Each group of 20 ticks were injected in triplicate with Dshsp90-, Dshsp70-, and tubulin-derived dsRNA, whereas the control group was injected with GFP dsRNA. Then, the total RNA was extracted and cDNA was synthesized and subjected to RT-qPCR. After the confirmation of knockdown, the ticks were incubated for 24 h and were exposed to − 20 °C lethal temperature (LT50), and then the mortality was calculated.

Results

Results indicated that Dshsp70 and Dshsp90 contained an open reading frame of 345 and 2190 nucleotides that encoded 114 and 729 amino acid residues, respectively. The transcript Dshsp70 showed 90% similarity with that identified from Dermacentor variabilis, whereas Dshsp90 showed 85% similarity with that identified from Ixodes scapularis. Multiple sequence alignment indicates that the deduced amino acid sequences of D. silvarum Hsp90, Hsp70, and tubulin show very high sequence identity to their corresponding sequences in other species. Hsp90 and Hsp70 display highly conserved and signature amino acid sequences with well-conserved MEEVD motif at the C-terminal in Hsp90 and a variable C-terminal region with a V/IEEVD-motif in Hsp70 that bind to numerous co-chaperones. RNA interference revealed that the mortality of D. silvarum was significantly increased after injection of dsRNA of Dshsp70 (P = 0.0298) and tubulin (P = 0.0448), whereas no significant increases were observed after the interference of Dshsp90 (P = 0.0709).

Conclusions

The above results suggested that Dshsp70 and tubulin play an essential role in the low-temperature adaptation of ticks. The results of this study can contribute to the understanding of the survival and acclimatization of overwintering ticks.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-05056-y.

Diapause differentially modulates the transcriptomes of fat body and flight muscle in the Colorado potato beetle

Many temperate insects, such as the Colorado potato beetle, enter diapause in winter, during which they arrest their development, suppress their metabolic rate and have high stress tolerance. Diapause phenotypes can be transcriptionally regulated, however many studies to date report only whole animal gene expression rather than tissue-specific processes during diapause. We used RNA-seq to measure gene expression in fat body and flight muscle of diapausing and non-diapausing beetles. We used differential expression and GO enrichment analyses to evaluate longstanding hypotheses about the mechanisms that drive arrested development, changes in energy metabolism, and increased stress tolerance during diapause. We found evidence of G2/M cell cycle arrest, juvenile hormone catabolism, increased antioxidant metabolism, epigenetic modification, transposable element regulation, and cytoskeletal remodeling in both the fat body and flight muscle of diapausing beetles. Beetles differentially modulated the fat body and flight muscle transcriptomes during diapause with fat body playing a larger role in the hypoxia response and immunity, whereas flight muscle had higher abundance of transcripts related to the chaperone response and proteostasis. Our transcriptome provides evidence for distinct roles and responses of fat body and flight muscle during diapause in the Colorado potato beetle, and we provide testable hypotheses for biological processes that appear to drive diapause phenotypes in insects.

The effects of genetic drift and genomic selection on differentiation and local adaptation of the introduced populations of Aedes albopictus in southern Russia

Background

Asian tiger mosquito Aedes albopictus is an arbovirus vector that has spread from its native habitation areal in Southeast Asia throughout North and South Americas, Europe, and Africa. Ae. albopictus was first detected in the Southern Federal District of the Russian Federation in the subtropical town of Sochi in 2011. In subsequent years, this species has been described in the continental areas with more severe climate and lower winter temperatures.

Methods

Genomic analysis of pooled Ae. albopictus samples collected in the mosquito populations in the coastal and continental regions of the Krasnodar Krai was conducted to look for the genetic changes associated with the spread and potential cold adaptation in Ae. albopictus.

Results

The results of the phylogenetic analysis based on mitochondrial genomes corresponded well with the hypothesis that Ae. albopictus haplotype A1a2a1 was introduced into the region from a single source. Population analysis revealed the role of dispersal and genetic drift in the local adaptation of the Asian tiger mosquito. The absence of shared haplotypes between the samples and high fixation indices suggest that gene flow between samples was heavily restricted. Mitochondrial and genomic differentiation together with different distances between dispersal routes, natural and anthropogenic barriers and local effective population size reduction could lead to difficulties in local climatic adaptations due to reduced selection effectiveness. We have found genomic regions with selective sweep patterns which can be considered as having been affected by recent selection events. The genes located in these regions participate in neural protection, lipid conservation, and cuticle formation during diapause. These processes were shown to be important for cold adaptation in the previous transcriptomic and proteomic studies. However, the population history and relatively low coverage obtained in the present article could have negatively affect sweep detection.

Transcriptome analyses and weighted gene coexpression network analysis reveal key pathways and genes involved in the rapid cold resistance of the Chinese white wax scale insect

The Chinese white wax scale insect, Ericerus pela, is an important resource insect in China. The rapid response of E. pela to decreasing temperatures plays key roles in the population distribution. In this study, we analyzed the gene expression of E. pela treated with low temperature using transcriptome analyses and weighted gene coexpression network analysis (WGCNA). The results showed that the cold resistance of E. pela involved changes in the expression of many genes. The genes were mainly involved in alcohol formation activity, lipid metabolism, membrane and structure, and oxidoreductase activity. According to the WGCNA results, some pathways related to cold resistance were found in the genes in the modules, such as cytoskeleton proteins, cytoskeleton protein pathway, biosynthesis of unsaturated fatty acids, glycerophospholipid metabolism, ether lipid metabolism, and thermogenesis. Some of the hub genes were nonspecific lipid-transfer proteins, DnaJ homolog subfamily C member 13, paramyosin, tropomodulin, and tubulin beta chain. In particular, the hub genes of the tan module included the heat shock protein (hsp) 10, hsp 60, hsp 70, and hsp 90 genes. Thirty-five antifreeze protein (afp) genes were identified according to the annotation results. Three afp genes were further identified among the hub genes. Six of these genes were selected for heterogeneous protein expression. One of them was expressed successfully. The thermal hysteresis activity (THA) analyses showed that the THA was 1.73°C. These results showed that the cytoskeleton, lipid metabolism, thermogenesis, HSPs and AFPs may play important roles in the cold resistance of E. pela.

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.

Transcriptome and gene expression analysis of Rhynchophorus ferrugineus (Coleoptera: Curculionidae) during developmental stages

Background

Red palm weevil, Rhynchophorus ferrugineus Olivier, is one of the most destructive pests harming palm trees. However, genomic resources for R. ferrugineus are still lacking, limiting the ability to discover molecular and genetic means of pest control.

Methods

In this study, PacBio Iso-Seq and Illumina RNA-seq were used to generate transcriptome from three developmental stages of R. ferrugineus (pupa, 7th-instar larva, adult) to increase the understanding of the life cycle and molecular characteristics of the pest.

Results

Sequencing generated 625,983,256 clean reads, from which 63,801 full-length transcripts were assembled with N50 of 3,547 bp. Expression analyses revealed 8,583 differentially expressed genes (DEGs). Moreover, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that these DEGs were mainly related to the peroxisome pathway which associated with metabolic pathways, material transportation and organ tissue formation. In summary, this work provides a valuable basis for further research on the growth and development, gene expression and gene prediction, and pest control of R. ferrugineus.

Characterization of Gromphadorhina coquereliana hemolymph under cold stress

Low temperatures in nature occur together with desiccation conditions, causing changes in metabolic pathways and cellular dehydration, affecting hemolymph volume, water content and ion homeostasis. Although some research has been conducted on the effect of low temperature on Gromphadorhina coquereliana, showing that it can survive exposures to cold or even freezing, no one has studied the effect of cold on the hemolymph volume and the immune response of this cockroach. Here, we investigated the effect of low temperature (4 °C) on the abovementioned parameters, hemocyte morphology and total number. Cold stress affected hemocytes and the immune response, but not hemolymph volume. After stress, the number of circulating hemocytes decreased by 44.7%, but the ratio of apoptotic cells did not differ significantly between stressed and control individuals: 8.06% and 7.18%, respectively. The number of phagocyting hemocytes decreased by 16.66%, the hemocyte morphology drastically changed, and the F-actin cytoskeleton differed substantially in cold-stressed insects compared to control insects. Moreover, the surface area of the cells increased from 393.69 µm² in the control to 458.38 µm² in cold-treated animals. Together, our results show the links between cold stress and the cellular immune response, which probably results in the survival capability of this species.

Chilling induces unidirectional solute leak through the locust gut epithelia

Chill-susceptible insects, like the migratory locust, often die when exposed to low temperatures from an accumulation of tissue damage that is unrelated to freezing (chilling injury). Chilling injury is often associated with a loss of ion balance across the gut epithelia. It has recently been suggested that this imbalance is at least partly caused by a cold-induced disruption of epithelial barrier function. Here, we aimed to test this hypothesis in the migratory locust (Locustamigratoria). First, chill tolerance was quantified by exposing locusts to -2°C and recording chill coma recovery time and survival 24 h post-cold exposure. Longer exposure times significantly increased recovery time and caused injury and death. Ion-selective microelectrodes were also used to test for a loss of ion balance in the cold. We found a significant increase of haemolymph K+ and decrease of haemolymph Na+ concentration over time. Next, barrier failure along the gut was tested by monitoring the movement of an epithelial barrier marker (FITC-dextran) across the gut epithelia during exposure to -2°C. We found a significant increase in haemolymph FITC-dextran concentration over time in the cold when assayed in the mucosal to serosal direction. However, when tested in the serosal to mucosal direction, we saw minimal marker movement across the gut epithelia. This suggests that while cold-induced barrier disruption is present, it is apparently unidirectional. It is important to note that these data reveal only the phenomenon itself. The location of this leak as well as the underlying mechanisms remain unclear and require further investigation.

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.

TMT-Based Quantitative Proteomic Profiling of Overwintering Lissorhoptrus oryzophilus

Adaptations to low temperature play a critical role in restricting the geographical distribution of insects. Decreasing day lengths and temperatures trigger seasonal cold adaptations in insects. These adaptions include changes in expression at the miRNA, mRNA and protein levels. The rice water weevil (RWW), Lissorhoptrus oryzophilus, introduced from the Mississippi River, is a globally invasive pest of wetland rice that can survive at the northern border of China. To investigate the changes in expression at the protein level in overwintering female RWW adults, 6-plex tandem mass tags (TMTs) were used in overwintering and summer adults. By using a proteome database available for Curculionidae, 1077 proteins were quantified, 183 of which differed significantly between the overwintering and summer samples. To further understand these differentially expressed proteins (DEPs), bioinformatics analyses such as gene ontology (GO) enrichment analyses were performed. DEPs associated with the terms binding, structural molecule activity, catalytic activity, multicellular organismal process, extracellular region, chitin binding, metabolic process, intracellular part and organic cyclic compound binding were altered by selection during winter. The changes in the expression of these proteins suggest that the proteins are important for RWW survival in winter.

Transcriptional analysis of insect extreme freeze tolerance

Few invertebrates can survive cryopreservation in liquid nitrogen, and the mechanisms by which some species do survive are underexplored, despite high application potential. Here, we turn to the drosophilid Chymomyza costata to strengthen our fundamental understanding of extreme freeze tolerance and gain insights about potential avenues for cryopreservation of biological materials. We first use RNAseq to generate transcriptomes of three C. costata larval phenotypic variants: those warm-acclimated in early or late diapause (weak capacity to survive cryopreservation), and those undergoing cold acclimation after diapause entry (extremely freeze tolerant, surviving cryopreservation). We identify mRNA transcripts representing genes and processes that accompany the physiological transition to extreme freeze tolerance and relate cryopreservation survival to the transcriptional profiles of select candidate genes using extended sampling of phenotypic variants. Enhanced capacity for protein folding, refolding and processing appears to be a central theme of extreme freeze tolerance and may allow cold-acclimated larvae to repair or eliminate proteins damaged by freezing (thus mitigating the toxicity of denatured proteins, endoplasmic reticulum stress and subsequent apoptosis). We also find a number of candidate genes (including both known and potentially novel, unannotated sequences) whose expression profiles tightly mirror the change in extreme freeze tolerance status among phenotypic variants.

Understanding the regulation of overwintering diapause molecular mechanisms in Culex pipiens pallens through comparative proteomics

To reveal overwintering dormancy (diapause) mechanisms of Culex pipiens pallens (L.), global protein expression differences at three separate time points represent nondiapause, diapause preparation and overwintering diapause phases of Cx. pipiens pallens were compared using iTRAQ. Cx. pipiens pallens females accumulate more lipid droplets during diapause preparation and overwintering diapause maintenance than during the nondiapause phase. A total of 1030 proteins were identified, among which 1020 were quantified and compared. Gene Ontology, Kyoto Encyclopedia of Genes and Genomes (KEGG), Domain and Clusters of Orthologous Groups (COG) analyses revealed key groups of proteins, pathways and domains differentially regulated during diapause preparation and overwintering diapause maintenance phases in this mosquito, including major shifts in energy production and conversion, fatty acid metabolism, the citrate (TCA) cycle, and the cytoskeletal reorganization pathway. Our results provide novel insight into the molecular bases of diapause in mosquitoes and corroborate previously reported diapause-associated features in invertebrates. More interestingly, the phototransduction pathway exists in Cx. pipiens pallens, in particular, actin, rather than other proteins, appears to have substantial role in diapause regulation. In addition, the differential changes in calmodulin protein expression in each stage implicate its important regulatory role of the Cx. pipiens pallens biological clock. Finally, 24 proteins were selected for verification of differential expression using a parallel reaction monitoring strategy. The findings of this study provide a unique opportunity to explore the molecular modifications underlying diapause in mosquitoes and might therefore enable the future design and development of novel genetic tools for improving management strategies in mosquitoes.

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.

Transcriptome Analysis Reveals Candidate Genes for Cold Tolerance in Drosophila ananassae

Coping with daily and seasonal temperature fluctuations is a key adaptive process for species to colonize temperate regions all over the globe. Over the past 18,000 years, the tropical species Drosophila ananassae expanded its home range from tropical regions in Southeast Asia to more temperate regions. Phenotypic assays of chill coma recovery time (CCRT) together with previously published population genetic data suggest that only a small number of genes underlie improved cold hardiness in the cold-adapted populations. We used high-throughput RNA sequencing to analyze differential gene expression before and after exposure to a cold shock in cold-tolerant lines (those with fast chill coma recovery, CCR) and cold-sensitive lines (slow CCR) from a population originating from Bangkok, Thailand (the ancestral species range). We identified two candidate genes with a significant interaction between cold tolerance and cold shock treatment: GF14647 and GF15058. Further, our data suggest that selection for increased cold tolerance did not operate through the increased activity of heat shock proteins, but more likely through the stabilization of the actin cytoskeleton and a delayed onset of apoptosis.

Insect fat body cell morphology and response to cold stress is modulated by acclimation

Mechanistic understanding about the nature of cellular cryoinjury and mechanisms by which some animals survive freezing while others do not is currently lacking. Here, we exploited the broadly manipulable freeze tolerance of larval malt flies (Chymomyza costata) to uncover cell and tissue morphological changes associated with freeze mortality. Diapause induction, cold acclimation and dietary proline supplementation generate malt fly variants ranging from weakly to extremely freeze tolerant. Using confocal microscopy and immunostaining of the fat body, Malpighian tubules and anterior midgut, we described tissue and cytoskeletal (F-actin and α-tubulin) morphologies among these variants after exposure to various cold stresses (from chilling at -5°C to extreme freezing at -196°C), and upon recovery from cold exposure. Fat body tissue appeared to be the most susceptible to cryoinjury: freezing caused coalescence of lipid droplets, loss of α-tubulin structure and apparent aggregation of F-actin. A combination of diapause and cold acclimation substantially lowered the temperature at which these morphological disruptions occurred. Larvae that recovered from a freezing challenge repaired F-actin aggregation but not lipid droplet coalescence or α-tubulin structure. Our observations indicate that lipid coalescence and damage to α-tubulin are non-lethal forms of freeze injury, and suggest that repair or removal (rather than protection) of actin proteins is a potential mechanism of acquired freeze tolerance.

Genomic signatures of local adaptation to the degree of environmental predictability in rotifers

Environmental fluctuations are ubiquitous and thus essential for the study of adaptation. Despite this, genome evolution in response to environmental fluctuations —and more specifically to the degree of environmental predictability– is still unknown. Saline lakes in the Mediterranean region are remarkably diverse in their ecological conditions, which can lead to divergent local adaptation patterns in the inhabiting aquatic organisms. The facultatively sexual rotifer Brachionus plicatilis shows diverging local adaptation in its life-history traits in relation to estimated environmental predictability in its habitats. Here, we used an integrative approach —combining environmental, phenotypic and genomic data for the same populations– to understand the genomic basis of this diverging adaptation. Firstly, a novel draft genome for B. plicatilis was assembled. Then, genome-wide polymorphisms were studied using genotyping by sequencing on 270 clones from nine populations in eastern Spain. As a result, 4,543 high-quality SNPs were identified and genotyped. More than 90 SNPs were found to be putatively under selection with signatures of diversifying and balancing selection. Over 140 SNPs were correlated with environmental or phenotypic variables revealing signatures of local adaptation, including environmental predictability. Putative functions were associated to most of these SNPs, since they were located within annotated genes. Our results reveal associations between genomic variation and the degree of environmental predictability, providing genomic evidence of adaptation to local conditions in natural rotifer populations.

How crickets become freeze tolerant: The transcriptomic underpinnings of acclimation in Gryllus veletis

Some ectotherms can survive internal ice formation. In temperate regions, freeze tolerance is often induced by decreasing temperature and/or photoperiod during autumn. However, we have limited understanding of how seasonal changes in physiology contribute to freeze tolerance, and how these changes are regulated. During a six week autumn-like acclimation, late-instar juveniles of the spring field cricket Gryllus veletis (Orthoptera: Gryllidae) become freeze tolerant, which is correlated with accumulation of low molecular weight cryoprotectants, elevation of the temperature at which freezing begins, and metabolic rate suppression. We used RNA-Seq to assemble a de novo transcriptome of this emerging laboratory model for freeze tolerance research. We then focused on gene expression during acclimation in fat body tissue due to its role in cryoprotectant production and regulation of energetics. Acclimated G. veletis differentially expressed >3000 transcripts in fat body. This differential expression may contribute to metabolic suppression in acclimated G. veletis, but we did not detect changes in expression that would support cryoprotectant accumulation or enhanced control of ice formation, suggesting that these latter processes are regulated post-transcriptionally. Acclimated G. veletis differentially regulated transcripts that likely coordinate additional freeze tolerance mechanisms, including upregulation of enzymes that may promote membrane and cytoskeletal remodelling, cryoprotectant transporters, cytoprotective proteins, and antioxidants. Thus, while accumulation of cryoprotectants and controlling ice formation are commonly associated with insect freeze tolerance, our results support the hypothesis that many other systems contribute to surviving internal ice formation. Together, this information suggests new avenues for understanding the mechanisms underlying insect freeze tolerance.

Mechanisms underpinning the beneficial effects of fluctuating thermal regimes in insect cold tolerance

Insects exposed to low temperature often have high mortality or exhibit sublethal effects. A growing number of recent studies have shown beneficial effects of exposing insects to recurrent brief warm pulses during low-temperature stress (fluctuating thermal regime, FTR). The physiological underpinnings of the beneficial effects of FTR on cold survival have been extensively studied over the past few years. Profiling with various ‘-omics’ techniques has provided supporting evidence for different physiological responses between insects exposed to FTR and constant low temperature. Evidence from transcriptomic, metabolomic and lipidomic studies points to a system-wide loss of homeostasis at low temperature that can be counterbalanced by repair mechanisms under FTR. Although there has been considerable progress in understanding the physiological mechanisms underlying the beneficial effects of FTR, here we discuss how many areas still lack clarity, such as the precise role(s) of heat shock proteins, compatible solutes or the identification of regulators and key players involved in the observed homeostatic responses. FTR can be particularly beneficial in applied settings, such as for model insects used in research, integrated pest management and pollination services. We also explain how the application of FTR techniques in large-scale facilities may require overcoming some logistical and technical constraints. FTR definitively enhances survival at low temperature in insects, but before it can be widely used, we suggest that the possible fitness and energy costs of FTR must be explored more thoroughly. Although FTR is not ecologically relevant, similar processes may operate in settings where temperatures fluctuate naturally.

Diapause in a tropical oil-collecting bee: molecular basis unveiled by RNA-Seq

BACKGROUND

Diapause is a natural phenomenon characterized by an arrest in development that ensures the survival of organisms under extreme environmental conditions. The process has been well documented in arthropods. However, its molecular basis has been mainly studied in species from temperate zones, leaving a knowledge gap of this phenomenon in tropical species. In the present study, the Neotropical and solitary bee Tetrapedia diversipes was employed as a model for investigating diapause in species from tropical zones. Being a bivoltine insect, Tetrapedia diversipes produce two generations of offspring per year. The first generation, normally born during the wet season, develops faster than individuals from the second generation, born after the dry season. Furthermore, it has been shown that the development of the progeny, of the second generation, is halted at the 5th larval instar, and remains in larval diapause during the dry season. Towards the goal of gaining a better understanding of the diapause phenomenon we compared the global gene expression pattern, in larvae, from both reproductive generations and during diapause. The results demonstrate that there are similarities in the observed gene expression patterns to those already described for temperate climate models, and also identify diapause-related genes that have not been previously reported in the literature.

RESULTS

The RNA-Seq analysis identified 2275 differentially expressed transcripts, of which 1167 were annotated. Of these genes, during diapause, 352 were upregulated and 815 were downregulated. According to their biological functions, these genes were categorized into the following groups: cellular detoxification, cytoskeleton, cuticle, sterol and lipid metabolism, cell cycle, heat shock proteins, immune response, circadian clock, and epigenetic control.

CONCLUSION

Many of the identified genes have already been described as being related to diapause; however, new genes were discovered, for the first time, in this study. Among those, we highlight: Niemann-Pick type C1, NPC2 and Acyl-CoA binding protein homolog (all involved in ecdysteroid synthesis); RhoBTB2 and SASH1 (associated with cell cycle regulation) and Histone acetyltransferase KAT7 (related to epigenetic transcriptional regulation). The results presented here add important findings to the understanding of diapause in tropical species, thus increasing the comprehension of diapause-related molecular mechanisms.

Recovery from supercooling, freezing, and cryopreservation stress in larvae of the drosophilid fly, Chymomyza costata

Physiological adjustments accompanying insect cold acclimation prior to cold stress have been relatively well explored. In contrast, recovery from cold stress received much less attention. Here we report on recovery of drosophilid fly larvae (Chymomyza costata) from three different levels of cold stress: supercooling to −10 °C, freezing at −30 °C, and cryopreservation at −196 °C. Analysis of larval CO₂ production suggested that recovery from all three cold stresses requires access to additional energy reserves to support cold-injury repair processes. Metabolomic profiling (targeting 41 metabolites using mass spectrometry) and custom microarray analysis (targeting 1,124 candidate mRNA sequences) indicated that additional energy was needed to: clear by-products of anaerobic metabolism, deal with oxidative stress, re-fold partially denatured proteins, and remove damaged proteins, complexes and/or organelles. Metabolomic and transcriptomic recovery profiles were closely similar in supercooled and frozen larvae, most of which successfully repaired the cold injury and metamorphosed into adults. In contrast, the majority of cryopreseved larvae failed to proceed in ontogenesis, showed specific metabolic perturbations suggesting impaired mitochondrial function, and failed to up-regulate a set of 116 specific genes potentially linked to repair of cold injury.

Effects of cold acclimation on rectal macromorphology, ultrastructure, and cytoskeletal stability in Gryllus pennsylvanicus crickets

Cold-acclimated insects maintain ion and water balance in the cold, potentially by reducing permeability or increasing diffusion distance across ionoregulatory epithelia such as the rectum. We explored whether cold acclimation induces structural modifications that minimize water and ion diffusion across the rectum and maintain rectal cell integrity. We investigated rectal structure and cytoskeletal stability in chill-susceptible adult Gryllus pennsylvanicus crickets acclimated for one week to either warm (25 °C) or cold (12 °C) conditions. After acclimation, we used light and transmission electron microscopy to examine rectal macromorphology and rectal pad paracellular ultrastructure. We also used fluorescence microscopy and a filamentous-actin (F-actin) specific phalloidin stain to compare the polymerization state of the actin cytoskeleton for each of the acclimation groups before and after a cold shock (1 h at -4 °C). Cold acclimation did not alter rectal pad cell density, or the thickness of the rectal pads, muscle, or cuticle. The tortuosity and width of the rectal pad paracellular channels also did not differ between warm- and cold-acclimated crickets. Rectal pad cells had clear basal and apical regions with differing densities of F-actin. Cold shock reduced the density of F-actin in warm-acclimated crickets, whereas cold-acclimated crickets appeared to have unchanged (basal) or enhanced (apical) F-actin density after cold shock. This suggests that while cold acclimation does not modify rectal permeability through structural modifications to increase diffusion distance for water and ions, cold-acclimated crickets have a modified cytoskeleton that resists the depolymerising effects of cold shock.

Comparative analysis of the transcriptome of the overwintering desert beetle Microdera punctipennis

The cold tolerance mechanisms of insect have been studied extensively on the model species Drosophila and a few other species at the transcriptional level. However studies on insects that inherit strong cold tolerance are limited. Cold hardy Tenebrionid beetle Microdera punctipennis is endemic to Gurbantonggut Desert, northwest of China. However, its genomic information is lacking. To investigate the overwintering mechanisms of M. punctipennis adult, RNA-seq was performed on the winter adults and the control adults that were kept in laboratory at 30 °C. A total of 175,247 unigenes were acquired with an average length of 645 bp. By using DESeq package, we identified 3367 unigenes that were up-regulated and 7988 down-regulated in the winter adults compared with the controls. To further our understanding of these differentially expressed genes (DEGs), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed. Pathway analysis showed that the "ECM-receptor interaction", "PI3K-Akt signaling pathway", "Estrogen signaling pathway", "Tight junction", and "Regulation of actin cytoskeleton", etc. might play important roles in M. punctipennis overwintering. The DEGs results from the RNA-Seq were confirmed partially by qRT-PCR for 13 DEGs, which showed high consistence with a Pearson's correlation coefficient of 0.851. Overall, the sequence data will provide basic information for subsequent bioinformatical analysis and mining of the genes responsible for cold tolerance in M. punctipennis, as well as for understanding the molecular mechanisms of desert beetle overwintering.

Multiple Trait Covariance Association Test Identifies Gene Ontology Categories Associated with Chill Coma Recovery Time in Drosophila melanogaster

The genomic best linear unbiased prediction (GBLUP) model has proven to be useful for prediction of complex traits as well as estimation of population genetic parameters. Improved inference and prediction accuracy of GBLUP may be achieved by identifying genomic regions enriched for causal genetic variants. We aimed at searching for patterns in GBLUP-derived single-marker statistics, by including them in genetic marker set tests, that could reveal associations between a set of genetic markers (genomic feature) and a complex trait. GBLUP-derived set tests proved to be powerful for detecting genomic features, here defined by gene ontology (GO) terms, enriched for causal variants affecting a quantitative trait in a population with low degree of relatedness. Different set test approaches were compared using simulated data illustrating the impact of trait- and genomic feature-specific factors on detection power. We extended the most powerful single trait set test, covariance association test (CVAT), to a multiple trait setting. The multiple trait CVAT (MT-CVAT) identified functionally relevant GO categories associated with the quantitative trait, chill coma recovery time, in the unrelated, sequenced inbred lines of the Drosophila melanogaster Genetic Reference Panel.

Large scale phosphoprotein profiling to explore Drosophila cold acclimation regulatory mechanisms

The regulatory mechanisms involved in the acquisition of thermal tolerance are unknown in insects. Reversible phosphorylation is a widespread post-translational modification that can rapidly alter proteins function(s). Here, we conducted a large-scale comparative screening of phosphorylation networks in adult Drosophila flies that were cold-acclimated versus control. Using a modified SIMAC method followed by a multiple MS analysis strategy, we identified a large collection of phosphopeptides (about 1600) and phosphoproteins (about 500) in both groups, with good enrichment efficacy (80%). The saturation curves from the four biological replicates revealed that the phosphoproteome was rather well covered under our experimental conditions. Acclimation evoked a strong phosphoproteomic signal characterized by large sets of unique and differential phosphoproteins. These were involved in several major GO superclusters of which cytoskeleton organization, positive regulation of transport, cell cycle, and RNA processing were particularly enriched. Data suggest that phosphoproteomic changes in response to acclimation were mainly localized within cytoskeletal network, and particularly within microtubule associated complexes. This study opens up novel research avenues for exploring the complex regulatory networks that lead to acquired thermal tolerance.

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.

Transcriptome and proteome analyses reveal complex mechanisms of reproductive diapause in the two-spotted spider mite, Tetranychus urticae

Although a variety of factors underlying diapause have been identified in arthropods and other organisms, the molecular mechanisms regulating diapause are still largely unknown. Here, to better understand this process, we examined diapause-associated genes in the two-spotted spider mite, Tetranychus urticae, by comparing the transcriptomes and proteomes of early diapausing and reproductive adult females. Amongst genes underlying diapause revealed by the transcriptomic and proteomic data sets, we described the noticeable change in Ca²⁺ -associated genes, including 65 Ca²⁺ -binding protein genes and 23 Ca²⁺ transporter genes, indicating that Ca²⁺ signalling has a substantial role in diapause regulation. Other interesting changes in diapause included up-regulation of (1) glutamate receptors that may be involved in synaptic plasticity changes, (2) genes involved in cytoskeletal reorganization including genes encoding each of the components of thick and thin filaments, tubulin and members of integrin signalling and (3) genes involved in anaerobic energy metabolism, which reflects a shift to anaerobic energy metabolism in early diapausing mites.

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.

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.

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.

Insight into the possible mechanism of the summer diapause of Delia antiqua (Diptera: Anthomyiidae) through digital gene expression analysis

The onion fly, Delia antiqua, is a major underground agricultural pest that can enter pupal diapause in the summer and winter seasons. However, little is known about its molecular regulation due to the lack of genomic resources. To gain insight into the possible mechanism of summer diapause (SD), high-throughput RNA-Seq data were generated from non-diapause (ND) and SD (initial, maintenance and quiescence phase) pupae. Three pair-wise comparisons were performed and identified, 1380, 1471 and 435, and were significantly regulated transcripts. Further analysis revealed that the enrichment of several functional terms related to juvenile hormone regulation, cell cycle, carbon hydrate and lipid metabolism, innate immune and stress responses, various signalling transductions, ubiquitin-dependent proteosome, and variation in cuticular and cytoskeleton components were found between ND and SD and between different phases of SD. Global characterization of transcriptome profiling between SD and ND contributes to the in-depth elucidation of the molecular mechanism of SD. Our results also offer insights into the evolution of insect diapause and support the importance of using the onion fly as a model to compare the molecular regulation events of summer and winter diapauses.

Drosophila americana Diapausing Females Show Features Typical of Young Flies

Diapause is a period of arrested development which is controlled physiologically, preprogrammed environmentally and characterized by metabolic depression that can occur during any stage of insect development. Nevertheless, in the genus Drosophila, diapause is almost always associated with the cessation of ovarian development and reproductive activity in adult females. In this work, we show that, in D. americana (a temperate species of the virilis group), diapause is a genetically determined delay in ovarian development that is triggered by temperature and/or photoperiod. Moreover, we show that in this species diapause incidence increases with latitude, ranging from 13% in the southernmost to 91% in the northernmost range of the distribution. When exposed to diapause inducing conditions, both diapausing and non-diapausing females show a 10% increase in lifespan, that is further increased by 18.6% in diapausing females, although senescence is far from being negligible. ActinD1 expression levels suggest that diapausing females are biologically much younger than their chronological age, and that the fly as a whole, rather than the ovarian development alone, which is phenotypically more evident, is delayed by diapause. Therefore, diapause candidate genes that show expression levels that are compatible with flies younger than their chronological age may not necessarily play a role in reproductive diapause and in adaptation to seasonally varying environmental conditions.

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.

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.

Identification of FOXO targets that generate diverse features of the diapause phenotype in the mosquito Culex pipiens

Insulin and juvenile hormone signaling direct entry of the mosquito Culex pipiens into its overwintering adult diapause, and these two critical signaling pathways appear to do so by converging on the regulation of forkhead transcription factor (FOXO). Diapause is a complex phenotype, and FOXO emerges as a prime candidate for activating many of the diverse physiological pathways that generate the diapause phenotype. Here, we used ChIP sequencing to identify direct targets of FOXO. The nearest gene in a 10-kb region surrounding a predicted binding site was extracted for each binding site, resulting in a dataset containing genes potentially regulated by FOXO. By selecting candidate genes based on their functional relevance to diapause, we identified five gene categories of potential interest, including stress tolerance, metabolic pathways, lifespan extension, cell cycle and growth regulation, and circadian rhythms. Twelve targets were prioritized for further analysis, 10 of which were validated by ChIP-quantitative PCR and quantitative real-time PCR. These 10 genes activated by FOXO are highly up-regulated during diapause and are thus strong candidates for implementation of the diapause syndrome.