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

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 false codling moth, Thaumatotibia leucotreta in response to high and low-temperature treatments

The false codling moth (FCM), Thaumatotibia leucotreta, is a major quarantine pest native to Africa. Physical postharvest phytosanitary measures such as cold and heat treatments are championed to control its spread to new regions. However, the molecular changes that T. leucotreta undergoes as it attempts to adjust to its surroundings during the treatments and withstand the extreme temperatures remain largely unknown. The current study employs RNA-seq using the next-generation Illumina HiSeq platform to produce transcriptome profiles for differential gene expression analysis of T. leucotreta larvae under thermal stress. The transcriptome assembly analysis revealed 226,067 transcripts, clustering into 127,018 unigenes. In comparison to the 25 °C treated group, 874, 91, 159, and 754 individual differentially expressed genes (DEGs) co-regulated at -10, 0, 40, and 50 °C, respectively were discovered. Annotation of the DEGs by gene ontology (GO) revealed several genes, previously implicated in low and high-temperature stresses, including heat shock proteins, cytochrome P450, cuticle proteins, odorant binding proteins, and immune system genes. Kyoto Encyclopedia of Genes and Genomics (KEGG) classification analysis revealed that substantive DEGs were those involved in metabolic pathways such as thiamine, purine, folate, and glycerolipid metabolism pathways. The RT-qPCR validation of several significantly up- and down-regulated DEGs showed congruence between RNA-seq and qPCR data. This baseline study lays a foundation for future research into the molecular mechanisms underlying T. leucotreta's cold/heat tolerance by providing a thorough differential gene expression analysis that has identified multiple genes that may be associated with the insect's ability to withstand cold and heat.

Reduced cold tolerance of viral-infected leafhoppers attenuates viral persistent epidemics

Most arthropod-borne viruses produce intermittent epidemics in infected plants. However, the underlying mechanisms of these epidemics are unclear. Here, we demonstrated that rice stripe mosaic virus (RSMV), a viral pathogen, significantly increases the mortality of its overwintering vector, the leafhopper species Recilia dorsalis. Cold-stress assays indicated that RSMV reduces the cold tolerance of leafhoppers, a process associated with the downregulation of leafhopper cuticular protein genes. An RSMV-derived small RNA (vsiR-t00355379) was found to facilitate the downregulation of a leafhopper endocuticle gene that is mainly expressed in the abdomen (named RdABD-5) and is conserved across dipteran species. The downregulation of RdABD-5 expression in R. dorsalis resulted in fewer and thinner endocuticle lamellae, leading to decreased cold tolerance. This effect was correlated with a reduced incidence rate of RSMV in early-planted rice plants. These findings contribute to our understanding of the mechanism by which viral pathogens reduce cold tolerance in arthropod vectors and suggest an approach to managing the fluctuating prevalence of arboviruses.

IMPORTANCE

Increasing arthropod vector dispersal rates have increased the susceptibility of crop to epidemic viral diseases. However, the incidence of some viral diseases fluctuates annually. In this study, we demonstrated that a rice virus reduces the cold tolerance of its leafhopper vector, Recilia dorsalis. This effect is linked to the virus-derived small RNA-mediated downregulation of a gene encoding a leafhopper abdominal endocuticle protein. Consequently, the altered structural composition of the abdominal endocuticle reduces the overwinter survival of leafhoppers, resulting in a lower incidence of RSMV infection in early-planted rice plants. Our findings illustrate the important roles of RNA interference in virus-vector insect-environment interactions and help explain the annual fluctuations of viral disease epidemics in rice fields.

Trehalose 6-phosphate synthase gene rdtps1 contributes to thermal acclimation in Rhyzopertha dominica

BACKGROUND

The lesser grain borer (Rhyzopertha dominica), a worldwide primary pest of stored grain, causes serious economic losses and threatens stored food safety. R. dominica can respond to changes in temperature, especially the adaptability to heat. In this study, transcriptome analysis of R. dominica exposed to different temperatures was performed to elucidate differences in gene expression and the underling molecular mechanism.

RESULTS

Isoform-sequencing generated 17,721,200 raw reads and yielded 20,416 full-length transcripts. A total of 18,880 (92.48%) transcripts were annotated. We extracted RNA from R. dominica reared at 5 °C (cold stress), 15 °C (cold stress), 27 °C (ambient temperature) and 40 °C (heat stress) for RNA-seq. Compared to those of control insects reared at 27 °C, 119, 342, and 875 differentially expressed genes (DEGs) were identified at 5 °C, 15 °C, and 40 °C, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that pathways associated with "fatty acid metabolism", "fatty acid biosynthesis", "AMPK signaling pathway", "neuroactive ligand receptor interaction", and "longevity regulating pathway-multiple species" were significantly enriched. The functional annotation revealed that the genes encoding heat shock proteins (HSPs), fatty acid synthase (FAS), phospholipases (PLA), trehalose transporter (TPST), trehalose 6-phosphate synthase (TPS), and vitellogenin (Vg) were most likely involved in temperature regulation, which was also validated by RT-qPCR. Seven candidate genes (rdhsp1, rdfas1, rdpla1, rdtpst1, rdtps1, rdvg1, and rdP450) were silenced in the RNA interference (RNAi) assay. RNAi of each candidate gene suggested that inhibiting rdtps1 expression significantly decreased the trehalose level and survival rate of R. dominica at 40 °C.

CONCLUSIONS

These results indicated that trehalose contributes to the high temperature resistance of R. dominica. Our study elucidates the molecular mechanisms underlying heat tolerance and provides a potential target for the pest management in R. dominica.

High and Low Temperatures Differentially Affect Survival, Reproduction, and Gene Transcription in Male and Female Moths of Spodoptera frugiperda

In this study, we found that both heat and cold stresses significantly affected the survival and reproduction of both sexes in Spodoptera frugiperda adults, with larvae showing relatively higher extreme temperature tolerance. Further transcriptomic analysis in adults found remarkable differences and similarities between sexes in terms of temperature stress responses. Metabolism-related processes were suppressed in heat stressed females, which did not occur to the same extend in males. Moreover, both heat and cold stress reduced immune activities in both sexes. Heat stress induced the upregulation of many heat shock proteins in both sexes, whereas the response to cold stress was insignificant. More cold tolerance-related genes, such as cuticle proteins, UDP-glucuronosyltransferase, and facilitated trehalose transporter Tret1, were found upregulated in males, whereas most of these genes were downregulated in females. Moreover, a large number of fatty acid-related genes, such as fatty acid synthases and desaturases, were differentially expressed under heat and cold stresses in both sexes. Heat stress in females induced the upregulation of a large number of zinc finger proteins and reproduction-related genes; whereas cold stress induced downregulation in genes linked to reproduction. In addition, TRPA1-like encoding genes (which have functions involved in detecting temperature changes) and sex peptide receptor-like genes were found to be differentially expressed in stressed moths. These results indicate sex-specific heat and cold stress responses and adaptive mechanisms and suggest sex-specific trade-offs between stress-resistant progresses and fundamental metabolic processes as well as between survival and reproduction.

Comparative Transcriptome Analysis of Galeruca daurica Reveals Cold Tolerance Mechanisms

Galeruca daurica (Joannis) is a pest species with serious outbreaks in the Inner Mongolian grasslands in recent years, and its larvae and eggs are extremely cold-tolerant. To gain a deeper understanding of the molecular mechanism of its cold-tolerant stress response, we performed de novo transcriptome assembly of G. daurica via RNA-Seq and compared the differentially expressed genes (DEGs) of first- and second-instar larvae grown and developed indoors and outdoors, respectively. The results show that cold tolerance in G. daurica is associated with changes in gene expression mainly involved in the glycolysis/gluconeogenesis pathway, the fatty acid biosynthesis pathway and the production of heat shock proteins (HSPs). Compared with the control group (indoor), the genes associated with gluconeogenesis, fatty acid biosynthesis and HSP production were up-regulated in the larvae grown and developed outdoors. While the changes in these genes were related to the physiological metabolism and growth of insects, it was hypothesized that the proteins encoded by these genes play an important role in cold tolerance in insects. In addition, we also investigated the expression of genes related to the metabolic pathway of HSPs, and the results show that the HSP-related genes were significantly up-regulated in the larvae of G. daurica grown and developed outdoors compared with the indoor control group. Finally, we chose to induce significant expression differences in the Hsp70 gene (Hsp70A1, Hsp70-2 and Hsp70-3) via RNAi to further illustrate the role of heat stress proteins in cold tolerance on G. daurica larvae. The results show that separate and mixed injections of dsHSP70A1, dsHsp70-2 and dsHsp70-3 significantly reduced expression levels of the target genes in G. daurica larvae. The super-cooling point (SCP) and the body fluid freezing point (FP) of the test larvae were determined after RNAi using the thermocouple method, and it was found that silencing the Hsp70 genes significantly increased the SCP and FP of G. daurica larvae, which validated the role of heat shock proteins in the cold resistance of G. daurica larvae. Our findings provide an important theoretical basis for further excavating the key genes and proteins in response to extremely cold environments and analyzing the molecular mechanism of cold adaptation in insects in harsh environments.

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.

Gene Expression Analysis Reveals Potential Regulatory Factors Response to Temperature Stress in Bemisia tabaci Mediterranean

Exposure to extreme temperatures can hinder the development of insects and even reduce their survival rate. However, the invasive species Bemisia tabaci exhibits an impressive response to different temperatures. This study aims to identify important transcriptional changes of B. tabaci occupying different temperature habitats by performing RNA sequencing on populations originating from three regions of China. The results showed that the gene expression of B. tabaci populations inhabiting regions with different temperatures was altered and identified 23 potential candidate genes that respond to temperature stress. Furthermore, three potential regulatory factors' (the glucuronidation pathway, alternative splicing, and changes in the chromatin structure) response to different environmental temperatures were identified. Among these, the glucuronidation pathway is a notable regulatory pathway. A total of 12 UDP-glucuronosyltransferase genes were found in the transcriptome database of B. tabaci obtained in this study. The results of DEGs analysis suggest that UDP-glucuronosyltransferases with a signal peptide may help B. tabaci resist temperature stress by sensing external signals, such as BtUGT2C1 and BtUGT2B13, which are particularly important in responding to temperature changes. These results will provide a valuable baseline for further research on the thermoregulatory mechanisms of B. tabaci that contributes to its ability to effectively colonize regions with considerable temperature differences.

Comparative transcriptome analysis of Callosobruchus chinensis (L.) (Coleoptera: Chrysomelidae-Bruchinae) after heat and cold stress exposure

Callosobruchus chinensis is regarded as one of the cosmopolitan pests of legume crops and can cause tremendous losses to a variety of beans. This study focused on comparative transcriptome analyses of C. chinensis exposed to 45 °C (heat stress), 27 °C (ambient temperature) and -3 °C (cold stress) for 3 h to investigate the gene differences and underlying molecular mechanisms. There were 402 and 111 differentially expressed genes (DEGs) identified in the heat and cold stress treatments, respectively. "cell process", "cell" and "binding" were the main enriched functions and biological processes revealed by gene ontology (GO) analysis. The clusters of orthologous genes (COG) showed that DEGs were assigned to the categories: "posttranslational modification, protein turnover, chaperones", "lipid transport and metabolism", and "general function prediction only". With respect to the Kyoto Encyclopedia of Genes and Genomes (KEGG), the "longevity regulating pathway-multiple species", "carbon metabolism", "peroxisome", "protein processing in endoplasmic", "glyoxylate and dicarboxylate metabolism" pathways were significantly enriched. The annotation and enrichment analysis revealed that genes encoding heat shock proteins (Hsps) and cuticular proteins were significantly upregulated under high and low-temperature stresses, respectively. In addition, some DEGs encoding "Protein lethal essential for life", "Reverse transcriptase", "DnaJ domain", "Cytochrome" and "Zinc finger protein" were also upregulated to varying degrees. Transcriptomic data were validated using qRT‒PCR, which confirmed that they were consistent. In this paper, the temperature tolerance of C. chinensis adults was evaluated and the results showed that female adults were more sensitive to heat and cold stress than males, and the upregulation of heat shock protein and epidermal protein was the largest in DEGs after heat and cold stress, respectively. These findings provide a reference for further understanding the biological characteristics of C. chinensis adults and the molecular mechanisms underlying the response to low and high temperatures.

Selective effects of a short transient environmental fluctuation on a natural population

Natural populations experience continuous and often transient changes of environmental conditions. These in turn may result in fluctuating selection pressures leading to variable demographic and evolutionary population responses. Rapid adaptation as short-term response to a sudden environmental change has in several cases been attributed to polygenic traits, but the underlying genomic dynamics and architecture are poorly understood. In this study, we took advantage of a natural experiment in an insect population of the non-biting midge Chironomus riparius by monitoring genome-wide allele frequencies before and after a cold snap event. Whole genome pooled sequencing of time series samples revealed 10 selected haplotypes carrying ancient polymorphisms, partially with signatures of balancing selection. By constantly cold exposing genetically variable individuals in the laboratory, we could demonstrate with whole genome resequencing (i) that among the survivors, the same alleles rose in frequency as in the wild, and (ii) that the identified variants additively predicted fitness (survival time) of its bearers. Finally, by simultaneously sequencing the genome and the transcriptome of cold exposed individuals we could tentatively link some of the selected SNPs to the cis- and trans-regulation of genes and pathways known to be involved in cold response of insects, such as cytochrome P450 and fatty acid metabolism. Altogether, our results shed light on the strength and speed of selection in natural populations and the genomic architecture of its underlying polygenic trait. Population genomic time series data thus appear as promising tool for measuring the selective tracking of fluctuating selection in natural populations.

Detection of ryanodine receptor G4911E and I4754M mutation sites and analysis of binding modes of diamide insecticides with RyR on Galeruca daurica (Coleoptera: Chrysomelidae)

In recent years, the leaf beetle Galeruca daurica has broken out in the northern grasslands of Inner Mongolia, its management still mainly depends on chemical control using traditional insecticides or with novel action. The study was aim to identify mutation locus associated with resistance to diamide insecticides in field population of G. daurica, to provide a reference for rational selection of insecticides and to avoid the rapid resistance development to diamide insecticides. We cloned the full length of the ryanodine receptor gene of G. daurica (GdRyR), constructed 3D model and transmembrane regions by homologous modeling based on deduced amino acid sequence. Two potential mutation loci (Gly4911Glu and Ile4754Met) and allelic mutation frequencies were detected in individuals of G. daurica. In addition, their binding patterns to two diamide insecticides (chlorantraniliprole, cyantraniliprole) were analyzed separately using a molecular docking method. The full-length cDNA sequence of GdRyR (GenBank accession number: OP828593) was obtained by splicing and assembling, which is 15,399 bp in length and encodes 5,133 amino acids. The amino acid similarity of GdRyR with that of other Coleopteran insects were 86.70%-91.33%, which possessed the typical structural characteristics. An individual resistance allelic mutation frequency test on fifty field leaf beetles has identified 12% and 32% heterozygous individuals at two potential mutation loci Gly4911Glu and Ile4754Met, respectively. The affinity of the I4754M mutant model of GdRyR for chlorantraniliprole and cyantraniliprole was not significantly different from that of the wild type, and all had non-covalent interactions such as hydrogen bonding, hydrophobic interactions and π-cation interactions. However, the G4911E mutant model showed reduced affinity and reduced mode of action with two diamide insecticides, thus affecting the binding stability of the ryanodine receptor to the diamide insecticides. In conclusion, the G4911E mutation in GdRyR may be a potential mechanism for the development of resistance to diamide insecticides on G. daurica and should be a key concern for resistance risk assessment and reasonable applications of diamide insecticides for control in future. Moreover, this study could provide a reference for ryanodine receptor structure-based insecticides design.

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.

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

Background

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

Results

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

Conclusions

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

Supplementary Information

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

Expression Profiling and Functional Analysis of Candidate Odorant Receptors in Galeruca daurica

Simple Summary

Odorant receptors (ORs) play an important role in the olfactory system in insects. However, there is no functional research on the odorant receptors of Galeruca daurica. In this study, 21 OR genes were identified from the transcriptome database of G. daurica adults. Most GdauORs were mainly expressed in antennae, and the expression levels of GdauORs in adults were affected by age. When GdauOR4, GdauOR15, and GdauORco were silenced by RNAi, the electrophysiological responses to host plant volatiles were significantly decreased.

Abstract

Galeruca daurica (Joannis) is an oligophagous pest in the grasslands of Inner Mongolia, China, which feed mainly on Allium spp. Odorant receptors (ORs) play an important role in the olfactory system in insects, and function together with olfactory co-receptor (ORco). In this study, 21 OR genes were identified from the transcriptome database of G. daurica adults, and named GdauOR1-20 and GdauORco. The expression profiles were examined by RT-qPCR and RNA interference (RNAi) and electroantennogram (EAG) experiments were conducted to further identify the olfactory functions of GdauOR4, GdauOR11, GdauOR15, and GdauORco. It was found that 15 GdauORs (OR1, OR3-6, OR8, OR11-13, OR15, OR17-20, and ORco) were mainly expressed in antennae, and the expression levels of GdauORs in adults were affected by age. When GdauOR4, GdauOR15, and GdauORco were silenced by RNAi, the electrophysiological responses to host plant volatiles were significantly decreased in G. daurica. This study lays a necessary foundation for clarifying the mechanism on finding host plants in G. daurica.

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.

Comparative RNA-Seq Analyses of Solenopsis japonica (Hymenoptera: Formicidae) Reveal Gene in Response to Cold Stress

Solenopsis japonica, as a fire ant species, shows some predatory behavior towards earthworms and woodlice, and preys on the larvae of other ant species by tunneling into a neighboring colony's brood chamber. This study focused on the molecular response process and gene expression profiles of S. japonica to low (9 °C)-temperature stress in comparison with normal temperature (25 °C) conditions. A total of 89,657 unigenes (the clustered non-redundant transcripts that are filtered from the longest assembled contigs) were obtained, of which 32,782 were annotated in the NR (nonredundant protein) database with gene ontology (GO) terms, gene descriptions, and metabolic pathways. The results were 81 GO subgroups and 18 EggNOG (evolutionary genealogy of genes: Non-supervised Orthologous Groups) keywords. Differentially expressed genes (DEGs) with log₂fold change (FC) > 1 and log₂FC < -1 with p-value ≤ 0.05 were screened for cold stress temperature. We found 215 unigenes up-regulated and 115 unigenes down-regulated. Comparing transcriptome profiles for differential gene expression resulted in various DE proteins and genes, including fatty acid synthases and lipid metabolism, which have previously been reported to be involved in cold resistance. We verified the RNA-seq data by qPCR on 20 up- and down-regulated DEGs. These findings facilitate the basis for the future understanding of the adaptation mechanisms of S. japonica and the molecular mechanisms underlying the response to low temperatures.

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.

Effect of a single cold stress exposure on the reproductive behavior of male crickets

Cold stress is an important abiotic factor that can impact insect physiology, behavior, and overall fitness. Upon exposure to cold temperature, many insects enter a reversible state of immobility called chill coma. If the cold stress is brief and mild enough, insects can recover and regain full mobility upon return to warmer temperatures. However, the long-term impact of sublethal cold stress on insect behavior has been understudied. Here, sexually naïve adult male Acheta domesticus crickets were exposed to a single 0 °C cold stress for 6 h. One week later, the ability of these males to mate with a female was examined. For mating trials, a cold stressed male cricket was paired with a non-cold stressed, control female. Control pairs were comprised of a non-cold stressed control male and control female. Cold exposed males were less successful at mating than control males because most did not carry a spermatophore at the time of their mating trials. However, when these cold stressed males were allowed 1 h of chemosensory contact with a female, most produced a spermatophore. Males that produced spermatophores were given the opportunity to mate once with a female, and stressed males that successfully mated sired as many offspring as did control males. However, our results support that a single cold stress exposure can negatively impact the reproductive fitness of male crickets since it reduced their capacity to carry spermatophores and, as a consequence, to attract females.

Functional Characterization of Olfactory Proteins Involved in Chemoreception of Galeruca daurica

Odorant-binding proteins (OBPs) and chemosensory proteins (CSPs) play a fundamental role in insect olfaction. Galeruca daurica (Joannis) is a new pest with outbreak status in the Inner Mongolia grasslands, northern China. In this study, six olfactory protein genes (GdauOBP1, GdauOBP6, GdauOBP10, GdauOBP15, GdauCSP4, and GdauCSP5) were cloned by RACE and expressed by constructing a prokaryotic expression system. Their binding affinities to 13 compounds from host volatiles (Allium mongolicum) were determined by fluorescence-binding assay. In order to further explore the olfactory functions of GdauOBP15 and GdauCSP5, RNA interference (RNAi) and electroantennogram (EAG) experiments were conducted. Ligand-binding assays showed that the binding properties of the six recombinant proteins to the tested volatiles were different. GdauOBP6, GdauOBP15, GdauCSP4, and GdauCSP5 could bind several tested ligands of host plants. It was suspected that GdauOBP6, GdauOBP15, GdauCSP4, and GdauCSP5 were related to the host location in G. daurica. We also found that there were different EAG responses between males and females when the GdauOBP15 and GdauCSP5 genes were silenced by RNAi. The EAG response of G. daurica females to 2-hexenal was significantly decreased in dsRNA-OBP15-injected treatment compared to the control, and the dsRNA-CSP5-treated females significantly reduced EAG response to eight tested host volatiles (1,3-dithiane, 2-hexenal, methyl benzoate, dimethyl trisulfide, myrcene, hexanal, 1,3,5-cycloheptatriene, and p-xylene). However, the EAG response had no significant difference in males. Both GdauOBP15 and GdauCSP5 may have different functions between males and females in G. daurica and may play more important roles in females searching for host plants.

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.