Genome-wide identification and characterization of HSP gene superfamily in whitefly (Bemisia tabaci) and expression profiling analysis under temperature stress

Transcriptome-wide identification of the Hsp70 gene family in Pugionium cornutum and functional analysis of PcHsp70-5 under drought stress

MAIN CONCLUSION

The PcHsp70-5 enhances drought stress tolerance in transgenic Arabidopsis thaliana by upregulating stress tolerance genes and antioxidant enzyme activities. Heat shock proteins (HSPs) constitute a class of evolutionarily conserved proteins synthesized by organisms in response to various adverse environmental stimuli such as elevated temperatures, drought, hormonal fluctuations, high salt concentrations, and mechanical stress. However, research on HSPs has predominantly focused on model plants and crops, whereas their functions in desert plants have not been well investigated. This study analyzed the transcriptome of Pugionium cornutum and identified the complete ORFs of 25 genes of the PcHsp70 family genes. Their expression levels under drought stress were investigated using existing RNA-seq data. PcHsp70-5 genes exhibited high expression levels in both roots and leaves under drought stress. Consequently, the PcHsp70-5 genes were cloned and transformed into Arabidopsis thaliana for further analysis of their roles in drought stress response. Real-time fluorescence quantitative PCR (qRT-PCR) analysis demonstrated that both, drought stress and ABA, induced PcHsp70-5 expression. Under drought conditions, transgenic Arabidopsis plants exhibited markedly enhanced growth compared to wild-type plants, as evidenced by improved survival rates, root length, fresh weight, chlorophyll content, and reduced levels of malondialdehyde (MDA) and hydrogen peroxide (H2O2) in leaves, indicating that PcHsp70-5 overexpression mitigated growth inhibition and oxidative damage induced by drought stress. Subsequent research revealed that PcHsp70-5 overexpression significantly augmented the activities of superoxide dismutase (SOD), peroxidase (POD), catalase (CAT), and increased the proline content in transgenic Arabidopsis under drought conditions, alongside a significant increase in the expression levels of genes related to stress tolerance. This suggests that PcHsp70-5 enhances drought stress tolerance in transgenic Arabidopsis by upregulating stress tolerance genes and antioxidant enzyme activities.

Glycerol Metabolism is Important for the Low-Temperature Adaptation of a Global Quarantine Pest Anoplophora glabripennis Larvae

Anoplophora glabripennis is a critical global quarantine pest. Recently, its distribution has been extended to colder and higher-latitude regions. The adaptation to low temperatures is vital for the successful colonization of insects in new environments. However, the metabolic pathways of A. glabripennis larvae under cold stress remain undefined. This study analyzed the larval hemolymph under different low-temperature treatments using LC-MS/MS. The results showed that differential metabolites associated with sugar and lipid metabolism are pivotal in the larval chill coma process. Under low-temperature treatments, the glycerol content increased significantly compared with the control group. Cold stress significantly induced the expression of AglaGK2 and AglaGPDHs. After undergoing RNAi treatment for 48 h, larvae exposed to -20 °C for 1 h showed reduced recovery when injected with ds-AglaGK2 and ds-AglaGPDH1 compared to the control group, indicating that glycerol biosynthesis plays a role in the low-temperature adaptation of A. glabripennis larvae. Our results provide a theoretical basis for clarifying the molecular mechanism of A. glabripennis larvae in response to environmental stresses.

Dual Guardians of Immunity: FoRab10 and FoRab29 in Frankliniella occidentalis Confer Resistance to Tomato Spotted Wilt Orthotospovirus

Rab GTPase is critical for autophagy processes and is implicated in insect immunity against viruses. In this study, we aimed to investigate the role of FoRabs in the autophagic regulation of antiviral defense against tomato spotted wilt orthotospovirus (TSWV) in Frankliniella occidentalis. Transcriptome analysis revealed the downregulation of FoRabs in viruliferous nymph and adults of F. occidentalis in response to TSWV infection. Manipulation of autophagy levels with 3-MA and Rapa treatments resulted in a 5- to 15-fold increase and a 38-64% decrease in viral titers, respectively. Additionally, interference with FoRab10 in nymphs and FoRab29 in adults led to a 20-90% downregulation of autophagy-related genes, a decrease in ATG8-II (an autophagy marker protein), and an increase in the TSWV titers by 1.5- to 2.5-fold and 1.3- to 2.0-fold, respectively. In addition, the leaf disk and the living plant methods revealed increased transmission rates of 20.8-41.6 and 68.3-88.3%, respectively. In conclusion, FoRab10 and FoRab29 play a role in the autophagic regulation of the antiviral defense in F. occidentalis nymphs and adults against TSWV, respectively. These findings offer insights into the intricate immune mechanisms functional in F. occidentalis against TSWV, suggesting potential targeted strategies for F. occidentalis and TSWV management.

Comparative transcriptome analysis provides comprehensive insight into the molecular mechanisms of heat adaption in Plutella xylostella

Plutella xylostella is one of the most destructive pests for cruciferous vegetables, and is adaptability to different environmental stressors. However, we still know little about the molecular mechanisms of how P. xylostella adapt to thermal stress. Here, the comparative transcriptome analysis was conducted from the samples of control (27 °C, CK) and heat treatment (40 °C, 40 T) P. xylostella. The results showed 1253 genes were differentially expressed, with 624 and 629 genes up- and down-regulated respectively. The annotation analysis demonstrated that "Energy production and conversion", "Protein processing in endoplasmic reticulum", "Peroxisome" and "Tyrosine metabolism" pathways were significantly enriched. Additionally, we found the expression levels of heat shock protein genes (Hsps), cuticle related genes and mitochondrial genes were significantly up-regulated in 40 T insects, suggesting their vital roles in improving adaption to heat stress. Importantly, the SOD activity and MDA content of P. xylostella were both identified to be increased under high temperature stress, indicating the elevated antioxidant reactions might be involved in response to heat stress. In conclusion, the present study offered us an overview of gene expression changes after 40 °C treatments, and found some critical pathways and genes of P. xylostella might play the critical roles in resisting heat stress.

Genome-Wide Identification and Characterization of RdHSP Genes Related to High Temperature in Rhododendron delavayi

Heat shock proteins (HSPs) are molecular chaperones that play essential roles in plant development and in response to various environmental stresses. Understanding R. delavayi HSP genes is of great importance since R. delavayi is severely affected by heat stress. In the present study, a total of 76 RdHSP genes were identified in the R. delavayi genome, which were divided into five subfamilies based on molecular weight and domain composition. Analyses of the chromosome distribution, gene structure, and conserved motif of the RdHSP family genes were conducted using bioinformatics analysis methods. Gene duplication analysis showed that 15 and 8 RdHSP genes were obtained and retained from the WGD/segmental duplication and tandem duplication, respectively. Cis-element analysis revealed the importance of RdHSP genes in plant adaptations to the environment. Moreover, the expression patterns of RdHSP family genes were investigated in R. delavayi treated with high temperature based on our RNA-seq data, which were further verified by qRT-PCR. Further analysis revealed that nine candidate genes, including six RdHSP20 subfamily genes (RdHSP20.4, RdHSP20.8, RdHSP20.6, RdHSP20.3, RdHSP20.10, and RdHSP20.15) and three RdHSP70 subfamily genes (RdHSP70.15, RdHSP70.21, and RdHSP70.16), might be involved in enhancing the heat stress tolerance. The subcellular localization of two candidate RdHSP genes (RdHSP20.8 and RdHSP20.6) showed that two candidate RdHSPs were expressed and function in the chloroplast and nucleus, respectively. These results provide a basis for the functional characterization of HSP genes and investigations on the molecular mechanisms of heat stress response in R. delavayi.

Genome-wide identification and expression analysis of the heat shock protein gene superfamily in Hylurgus ligniperda

Hylurgus ligniperda belongs to Hylurgus Latreille, Curculionidae, Coleoptera. It primarily damages the base and roots of the trunk of pine plants. Short-term treatment at 42 °C can damage Hylurgus ligniperda; therefore, temperature is a vital factor limiting its spread. Heat shock proteins (HSPs) can protect, remove, and repair proteins to help H. ligniperda withstand high temperatures. However, information on HSP genes in H. ligniperda remains limited. In the study, we considered H. ligniperda as the focus of research and identified 56 HligHSP genes at the genome-wide level. These genes were mapped to the cytoplasm or nucleus. An identical subfamily exhibited a closely similar distribution of conserved domains. Combined with the transcriptome data collected in previous studies, we screened six candidate genes, namely HligsHSP-3, HligsHSP-4, HligHSP60-16, HligHSP70-3, HligHSP70-4, and HligHSP90-1, which are specifically expressed during different high-temperature treatments. A quantitative polymerase chain reaction was performed to measure the expression of these six HligHSPs in seven temperature treatment conditions. These genes may be involved in the heat resistance mechanism in adults. Our findings provided a foundation for further studying the heat resistance mechanism in H. ligniperda.

Transcriptome Analysis Reveals Antioxidant Defense Mechanisms in the Silkworm Bombyx mori after Exposure to Lead

Lead (Pb) is a major source of heavy metal contamination, and poses a threat to biodiversity and human health. Elevated levels of Pb can hinder insect growth and development, leading to apoptosis via mechanisms like oxidative damage. The midgut of silkworms is the main organ exposed to heavy metals. As an economically important lepidopteran model insect in China, heavy metal-induced stress on silkworms causes considerable losses in sericulture, thereby causing substantial economic damage. This study aimed to investigate Pb-induced detoxification-related genes in the midgut of silkworms using high-throughput sequencing methods to achieve a deeper comprehension of the genes' reactions to lead exposure. This study identified 11,567 unigenes and 14,978 transcripts. A total of 1265 differentially expressed genes (DEGs) were screened, comprising 907 upregulated and 358 downregulated genes. Subsequently, Gene Ontology (GO) classification analysis revealed that the 1265 DEGs were distributed across biological processes, cellular components, and molecular functions. This suggests that the silkworm midgut may affect various organelle functions and biological processes, providing crucial clues for further exploration of DEG function. Additionally, the expression levels of 12 selected detoxification-related DEGs were validated using qRT-PCR, which confirmed the reliability of the RNA-seq results. This study not only provides new insights into the detoxification defense mechanisms of silkworms after Pb exposure, but also establishes a valuable foundation for further investigation into the molecular detoxification mechanisms in silkworms.

Thermal acclimation uncovers a simple genetic basis of adaptation to high temperature in a cosmopolitan pest

Understanding a population's fitness heterogeneity and genetic basis of thermal adaptation is essential for predicting the responses to global warming. We examined the thermotolerance and genetic adaptation of Plutella xylostella to exposure to hot temperatures. The population fitness parameters of the hot-acclimated DBM strains varied in the thermal environments. Using genome scanning and transcription profiling, we find a number of genes potentially involved in thermal adaptation of DBM. Editing two ABCG transporter genes, PxWhite and PxABCG, confirmed their role in altering cuticle permeability and influencing thermal responses. Our results demonstrate that SNP mutations in genes and changes in gene expression can allow DBM to rapidly adapt to thermal environment. ABCG transporter genes play an important role in thermal adaptation of DBM. This work improves our understanding of genetic adaptation mechanisms of insects to thermal stress and our capacity to predict the effects of rising global temperatures on ectotherms.

HSP gene superfamily in Aspongopus chinensis Dallas: unravelling identification, characterisation and expression patterns during diapause and non-diapause stages

Aspongopus chinensis Dallas 1851, an insect of important economic value, faces challenges in artificial breeding due to mandatory diapause and limited access to wild resources. Heat shock proteins (Hsps) are thought to influence diapause in insects, but little is known about their role in A. chinensis during diapause. This study used genomic methods to identify 25 Hsp genes in A. chinensis, including two Hsp90, 14 Hsp70, four Hsp60 and five small Hsp genes, were located on seven chromosomes, respectively. The gene structures among the same families are relatively conserved. Meanwhile, the motif compositions and secondary structures of A. chinensis Hsps (AcHsps) were predicted. RNA-seq data and fluorescence quantitative PCR analysis showed that there were differences in the expression patterns of AcHsps in diapause and non-diapause stages, and AcHsp70-5 was significantly differentially expressed in both analysis, which was enriched in the pathway of response to hormone. All the results showed that Hsps play an important role in the diapause mechanism of A. chinensis. Our observations highlight the molecular evolution of the Hsp gene and their effect on diapause in A. chinensis.

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.

The influence of high-temperature frequency variation on the life-history traits of pyridaben-sensitive and -resistant strains of Tetranychus truncatus

With a generally warming global climate, the number of Tetranychus truncatus specimens in the Hexi region in China has been increasing. As ectotherms, the growth and development of T. truncatus are greatly affected by changes in environmental temperature. The effect of heatwaves on organisms depends on a delicate balance between damage and repair periods. Therefore, we simulated nine patterns of periodically recurring changes in the frequency of high-temperature days using an intraday gradual temperature change model to study and compare the effects on the development and reproduction of pyridaben-sensitive and -resistant strains of T. truncatus. The results showed that the influence of the frequency of high-temperature days on developmental stages, longevity and fecundity was different between the two strains. The egg and immature stages of the sensitive strain were all affected by hot days, whereas the adult stage was less affected by the frequency. The egg stage of the resistant strain was less affected; it was mainly affected in the immature and adult stages. Under the moderate condition of increasing the proportion of days at normal temperature, the longevity of the resistant strain gradually increased and reached a maximum at a 1:3 frequency, and then it decreased with the increase in high-temperature days. The longevity of the sensitive strain was less affected by frequency, and there was no significant difference between most treatment and control groups. In addition, both sensitive and resistant strains were able to complete growth and development under all nine frequencies of high-temperature days, but the reproductive rate was lower than it was at normal temperatures, indicating that both strains of T. truncatus adapted to high temperatures at the expense of reduced reproduction rates. This lays a key theoretical foundation for predicting the occurrence of agricultural pest populations under the background of climate warming and developing appropriate control strategies.

Characteristics of the Accessible Chromatin Landscape and Transcriptome under Different Temperature Stresses in Bemisia tabaci

Bemisia tabaci is an important invasive pest with worldwide distribution and strong temperature tolerance. Previous studies have shown that temperature tolerance varies significantly between the different invasive populations. Several key factors involved in epigenetic regulation have been identified and verified in B. tabaci; therefore, epigenetic adaptation mechanisms may also exist. This study aimed to detect changes in the chromatin accessibility landscape and genome-wide transcriptome under different temperature stresses in B. tabaci. Assay for transposase-accessible chromatin with high-throughput sequencing and RNA-seq analyses indicated that transcriptional activity of the genes strongly correlates with chromatin accessibility. Chromatin transcription-activated gene expression regulation is dominant during high-temperature stress in B. tabaci, mainly through the transcriptional repression of genes related to low-temperature stress resistance. Furthermore, B. tabaci resists low-temperature stress by regulating enzyme activities and withstands high-temperature stress by regulating metabolism and synthesis of organic substances, both achieved by altering chromatin accessibility. In summary, this study provides a theoretical basis for exploring changes in gene expression and chromatin accessibility under different temperature stresses, offering a new approach to unravelling regulatory mechanisms underlying the onset of molecular regulation in response to various temperature stress conditions.

Catalase promotes whitefly adaptation to high temperature by eliminating reactive oxygen species

Thermal stress usually leads to excessive production of reactive oxygen species (ROS) in all aerobic organisms. Catalases (CAT) are the key antioxidant enzymes, which act as the first line of defense against ROS in the antioxidant pathway. The highly invasive and widely distributed whitefly Bemisia tabaci MED damages plants by feeding as well as by transmitting many plant viruses. Previous studies have shown that strong adaptability to high temperature helps explain the spread of MED around the world. However, the mechanism underlying high temperature adaptation of this pest is not well understood. In this study, 6 CAT genes were identified from the MED genome and transcriptome dataset, among which BtCAT1, BtCAT2, and BtCAT3 were found to be highly expressed in adults. The expression of BtCAT1, BtCAT2, or BtCAT3 increased with induction temperature and induction time. The MED was exposed with mean high temperature (30 °C or 35 °C) and a short-term extremely high temperature (39 °C or 41 °C) after the silencing of BtCAT1, BtCAT2, or BtCAT3 to significantly increased ROS levels by at least 0.5 times and significantly decreased survival rate and fecundity of MED adults. The ROS level in the treated specimens gradually returned to a normal level after 24 h at 25 °C, but the survival rate still declined significantly. Taken together, our results demonstrate that CAT could help B. tabaci adapt to long-term mean high temperatures and short-term extremely high temperatures by eliminating excessive ROS.

Characterization and expression of heat shock and immune genes in natural populations of Prodiamesa olivacea (Diptera) exposed to thermal stress

This paper characterizes the heat stress response (HSR) and explores the impact of temperatures on the immune response of larvae from two chironomid species, Prodiamesa olivacea and Chironomus riparius. Genes involved in crucial metabolic pathways were de novo identified in P. olivacea: Hsp27, Hsp60, Hsp70, Hsc70, Cdc37, and HSF for the heat stress response (HSR) and TOLL, PGRP, C-type lectin, and JAK/hopscotch for the immune system response (ISR). Quantitative real-time PCR was used to evaluate the expression levels of the selected genes in short-term treatments (up to 120') at high temperatures (35 °C and 39 °C). Exposing P. olivacea to elevated temperatures resulted in HSR induction with increased expression of specific heat shock genes, suggesting the potential of HSPs as early indicators of acute thermal stress. Surprisingly, we found that heat shock represses multiple immune genes, revealing the antagonist relation between the heat shock response and the innate immune response in P. olivacea. Our results also showed species-dependent gene responses, with more significant effects in P. olivacea, for most of the biomarkers studied, demonstrating a higher sensitivity in this species to environmental stress conditions than that of C. riparius. This work shows a multi-species approach that enables a deeper understanding of the effects of heat stress at the molecular level in aquatic dipterans.

Comparative evolutionary analyses of eight whitefly Bemisia tabaci sensu lato genomes: cryptic species, agricultural pests and plant-virus vectors

BACKGROUND

The group of > 40 cryptic whitefly species called Bemisia tabaci sensu lato are amongst the world's worst agricultural pests and plant-virus vectors. Outbreaks of B. tabaci s.l. and the associated plant-virus diseases continue to contribute to global food insecurity and social instability, particularly in sub-Saharan Africa and Asia. Published B. tabaci s.l. genomes have limited use for studying African cassava B. tabaci SSA1 species, due to the high genetic divergences between them. Genomic annotations presented here were performed using the 'Ensembl gene annotation system', to ensure that comparative analyses and conclusions reflect biological differences, as opposed to arising from different methodologies underpinning transcript model identification.

RESULTS

We present here six new B. tabaci s.l. genomes from Africa and Asia, and two re-annotated previously published genomes, to provide evolutionary insights into these globally distributed pests. Genome sizes ranged between 616-658 Mb and exhibited some of the highest coverage of transposable elements reported within Arthropoda. Many fewer total protein coding genes (PCG) were recovered compared to the previously published B. tabaci s.l. genomes and structural annotations generated via the uniform methodology strongly supported a repertoire of between 12.8-13.2 × 10³ PCG. An integrative systematics approach incorporating phylogenomic analysis of nuclear and mitochondrial markers supported a monophyletic Aleyrodidae and the basal positioning of B. tabaci Uganda-1 to the sub-Saharan group of species. Reciprocal cross-mating data and the co-cladogenesis pattern of the primary obligate endosymbiont 'Candidatus Portiera aleyrodidarum' from 11 Bemisia genomes further supported the phylogenetic reconstruction to show that African cassava B. tabaci populations consist of just three biological species. We include comparative analyses of gene families related to detoxification, sugar metabolism, vector competency and evaluate the presence and function of horizontally transferred genes, essential for understanding the evolution and unique biology of constituent B. tabaci. s.l species.

CONCLUSIONS

These genomic resources have provided new and critical insights into the genetics underlying B. tabaci s.l. biology. They also provide a rich foundation for post-genomic research, including the selection of candidate gene-targets for innovative whitefly and virus-control strategies.

Characterization of HSP70 family in watermelon (Citrullus lanatus): identification, structure, evolution, and potential function in response to ABA, cold and drought stress

Watermelon (Citrullus lanatus) as a crop with important economic value, is widely cultivated around the world. The heat shock protein 70 (HSP70) family in plant is indispensable under stress conditions. However, no comprehensive analysis of watermelon HSP70 family is reported to date. In this study, 12 ClHSP70 genes were identified from watermelon, which were unevenly located in 7 out of 11 chromosomes and divided into three subfamilies. ClHSP70 proteins were predicted to be localized primarily in cytoplasm, chloroplast, and endoplasmic reticulum. Two pairs of segmental repeats and 1 pair of tandem repeats existed in ClHSP70 genes, and ClHSP70s underwent strong purification selection. There were many abscisic acid (ABA) and abiotic stress response elements in ClHSP70 promoters. Additionally, the transcriptional levels of ClHSP70s in roots, stems, true leaves, and cotyledons were also analyzed. Some of ClHSP70 genes were also strongly induced by ABA. Furthermore, ClHSP70s also had different degrees of response to drought and cold stress. The above data indicate that ClHSP70s may be participated in growth and development, signal transduction and abiotic stress response, laying a foundation for further analysis of the function of ClHSP70s in biological processes.

Chromosome-level genome assembly reveals potential epigenetic mechanisms of the thermal tolerance in the oriental fruit fly, Bactrocera dorsalis

The oriental fruit fly, Bactrocera dorsalis (Hendel), has very strong ecological adaptability and phenotypic plasticity. Here, the genome of B. dorsalis was assembled into 549.45 Mb sequences with a contig N50 length of 12.81 Mb. Among, 95.67 % assembled genome sequences were anchored on six chromosomes with an N50 length of 94.63 Mb. According to the basic characteristics of the sex chromosomes of Tephritidae, the X chromosome of B. dorsalis was identified. Significant gene expansions were detected in several important gene families related to adaptability. In particular, we annotated 50 histone modification enzymes (HMEs) in this genome. A comparative transcriptome analysis indicated that 12 HME genes were differentially expressed in two thermo-tolerant strains (heat and cold). Interestingly, four and seven of the 12 HME genes responded to heat shock or cold hardening, respectively. These evidences suggested that the histone modification as an epigenetic modification may be involved in the thermal tolerance of B. dorsalis, but with different regulation mechanisms in thermal acclimation and hardening. The high quality genome of B. dorsalis provides an invaluable resource for further functional genomic study. Moreover, comparative genomic analysis will shed insights on revealing the mechanisms of adaptive evolution in this fly.

Transcription dynamics of heat-shock proteins (Hsps) and endosymbiont titres in response to thermal stress in whitefly, Bemisia tabaci (Asia-I)

The sweet potato whitefly, Bemisia tabaci (Gennadius), is one of the several species complexes of whitefly that are currently significant agricultural pests. Bemisia tabaci infests more than 600 plant species and thrives under a wide range of temperature conditions. In addition to the direct damage caused by sucking plant sap, it vectors several plant viruses. Heat-shock proteins play a pivotal role in enabling the insect to extend its geographical location, survival, and reproduction under different stress conditions. B. tabaci harbours several endosymbionts under the genera Portiera, Rickettsia, Hamiltonella, Wolbachia, Arsenophonus, Cardinium, and Fritschea that directly or indirectly affect its fitness. By accelerating cuticle biosynthesis and sclerotisation, symbiotic microbes can reduce or enhance tolerance to extreme temperatures and detoxify heavy metals. Thus, symbionts or microbial communities can expand or constrain the abiotic niche space of their host and affect its ability to adapt to changing conditions. The present study delineates the effect of thermal stress on the expression of heat-shock genes and endosymbionts in B. tabaci. Studies of the expression level of heat-shock proteins with the help of quantitative real-time polymerase chain reaction (qRT-PCR) showed that heat- and cold-shock treatment fuels the increased expression of heat-shock proteins (Hsp40 and Hsp70). However, Hsp90 was not induced by a heat- and cold-shock treatment. A significant decrease in the relative titre of secondary endosymbionts, such as Rickettsia, Arsenophonus, and Wolbachia, were recorded in B. tabaci upon heat treatment. However, the titre of the primary symbiont, C. Portiera, was relatively unaffected by both cold and heat treatments. These results are indicative of the fact that Hsp genes and endosymbionts in B. tabaci are modulated in response to thermal stress, and this might be responsible for the adaptation of whitefly under changing climatic scenario.

Heat-shock protein 70-a hub gene-underwent adaptive evolution involved in whitefly-wild tomato interaction

BACKGROUND

The whitefly Bemisia tabaci causes severe damage to cultivated tomato plants, but actively avoids the wild tomato Solanum habrochaites. Moreover, the mortality of whitefly increases significantly after feeding with the wild tomato. However, additional experiments are warranted to more carefully elucidate the specific molecular elements underlying the interaction between whitefly and wild tomato.

RESULTS

Our results showed that S. habrochaites significantly increases the mortality of whitefly adults and decreases both their fertility and fecundity. In addition, the expression of stress-response genes in whitefly after exposure to S. habrochaites was analyzed using RNA sequencing. Weighted gene co-expression network analysis was conducted to identify the hub genes to determine their potential associations with the mortality of whitefly. These results suggested that the expression of heat-shock protein (HSP), multicopper oxidase, and 2-Oxo-4-hydroxy-4-carboxy-5-ureidoimidazoline (OHCU) decarboxylase genes were induced in whitefly. To validate the gene associations with whitefly mortality, a high-throughput in vivo model system and RNAi-based gene silencing were used. The results revealed that the RNAi-mediated depletion of the HSP gene, which belongs to the HSP70 subfamily, increased the mortality of whitefly. Furthermore, the selection pressure analysis showed that a total of five amino acid sites of positive selection were identified, three of which were located in the nucleotide-binding domain and the other two in the substrate-binding domain.

CONCLUSIONS

This is the first report on the potential implication of HSPs in whitefly-wild plant interactions. This study could more precisely identify the molecular mechanisms of whitefly in response to wild tomatoes. © 2022 Society of Chemical Industry.

Molecular characterization of six heat shock protein 70 genes from Arma chinensis and their expression patterns in response to temperature stress

Arma chinensis is an important predatory enemy of many agricultural and forest pests. Heat shock protein 70 (Hsp70) plays an essential role in insect adaptation to various stress factors. To explore the functions of Hsp70s in relation to thermal tolerance of A. chinensis, full-length cDNAs of six Hsp70 genes (AcHsp70Ba, AcHsp70-4, AcHsp68a, AcHsp68b, AcHsp70-2, and AcHsc70-4) were cloned. Their open reading frames (ORFs) were 1902, 2454, 1884, 1905, 1872, and 1947 bp, respectively. Developmental expression profiles showed that AcHsp70Ba, AcHsp70-4, and AcHsc70-4 were extremely highly expressed in adult stages. AcHsp68a and AcHsp70-2 showed the highest level of expression in nymph stages, and AcHsp68b was mainly expressed in male adults. Tissue distribution analysis demonstrated that the AcHsp70s were ubiquitously expressed but showing gene-specific and sex-driven patterns of expression. High temperature induced the expression of the six AcHsp70s. Among them, AcHsp70Ba, AcHsp70-4, AcHsp68a, and AcHsc70-4 were significantly induced at 38 °C for 6 h, while all six AcHsp70s were significantly induced at 38 °C for 24 h. There were differences in responses of the six AcHsp70s to low-temperature stress. The expressions of AcHsp70-4, AcHsp68a, and AcHsp68b in male adults were significantly repressed at 4 °C for 6 h, whereas AcHsp70Ba and AcHsp70-2 were significantly induced. The levels of AcHsp70Ba, AcHsp68b, and AcHsp70-2 in female adults were significantly repressed at 4 °C for 24 h, whereas AcHsc70-4 was significantly induced. These results suggested that AcHsp70s play important roles in various developmental stages and tissue function, and contribute to the tolerance of A. chinensis to extreme temperatures.

A heat shock 70kDa protein MaltHSP70-2 contributes to thermal resistance in Monochamus alternatus (Coleoptera: Cerambycidae): quantification, localization, and functional analysis

Background

Heat Shock Proteins 70 (HSP70s) in insects act on a diverse range of substrates to assist with overcoming extreme high temperatures. MaltHSP70-2, a member of HSP70s, has been characterized to involve in the thermotolerance of Monochamus alternatus in vitro, while quantification and localization of MaltHSP70-2 in various tissues and its functional analysis in vivo remain unclear.

Results

In this study, temporal expression of MaltHSP70-2 indicated a long-last inductive effect on MaltHSP70-2 expression maintained 48 hours after heat shock. MaltHSP70-2 showed a global response to heat exposure which occurring in various tissues of both males and females. Particularly in the reproductive tissues, we further performed the quantification and localization of MaltHSP70-2 protein using Western Blot and Immunohistochemistry, suggesting that enriched MaltHSP70-2 in the testis (specifically in the primary spermatocyte) must be indispensable to protect the reproductive activities (e.g., spermatogenesis) against high temperatures. Furthermore, silencing MaltHSP70-2 markedly influenced the expression of other HSP genes and thermotolerance of adults in bioassays, which implied a possible interaction of MaltHSP70-2 with other HSP genes and its role in thermal resistance of M. alternatus adults.

Conclusions

These findings shed new insights into thermo-resistant mechanism of M. alternatus to cope with global warming from the perspective of HSP70s functions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08858-1.

A chromosome-level genome assembly provides insights into the environmental adaptability and outbreaks of Chlorops oryzae

Chlorops oryzae is a pest of rice that has caused severe damage to crops in major rice-growing areas in recent years. We generated a 447.60 Mb high-quality chromosome-level genome with contig and scaffold N50 values of 1.17 Mb and 117.57 Mb, respectively. Hi-C analysis anchored 93.22% scaffolds to 4 chromosomes. The relatively high expression level of Heat Shock Proteins (HSPs) and antioxidant genes in response to thermal stress suggests these genes may play a role in the environmental adaptability of C. oryzae. The identification of multiple pathways that regulate reproductive development (juvenile hormone, 20-hydroxyecdsone, and insulin signaling pathways) provides evidence that these pathways also play an important role in vitellogenesis and thus insect population maintenance. These findings identify possible reasons for the increased frequency of outbreaks of C. oryzae in recent years. Our chromosome-level genome assembly may provide a basis for further genetic studies of C. oryzae, and promote the development of novel, sustainable strategies to control this pest.

A chromosome-level genome assembly for the rice pest, Chlorops oryzae, pinpoints molecular pathways that might contribute toward increased outbreaks for this important crop pest.

Genome-Wide Identification and Analysis of the Heat-Shock Protein Gene Superfamily in Bemisia tabaci and Expression Pattern Analysis under Heat Shock

Simple Summary

Bemisia tabaci MED is an invasive pest that had caused considerable economic damage in the past decades. Its successful colonization is closely related to heat-shock proteins (HSPs), which are related to heat resistance. In this study, 33 BtaHsps were identified based on the sequenced genome of B. tabaci MED belonging to six HSP families, among which 22 BtaHsps were newly identified. Analysis of the secondary structure and evolutionary relationship showed that they were all closely related. In addition, BtaHsp90A3 of the HSP90 family was screened by analyzing the expression level changes of these genes under 42 °C heat shock and RNAi was performed on the BtaHsp90A3. The results showed that the silencing of BtaHsp90A3 is closely related to the heat resistance of B. tabaci MED. Taken together, this study conducted an in-depth identification of BtaHsps that clarifies their evolutionary relationships and their response to thermal stress in B. tabaci MED.

Abstract

The thermal tolerance of Bemisia tabaci MED, an invasive whitefly species with worldwide distribution, plays an important role in its ecological adaptation during the invasion process. Heat-shock proteins (HSPs) are closely related to heat resistance. In this study, 33 Hsps (BtaHsps) were identified based on sequenced genome of B. tabaci MED belonging to six HSP families, among which 22 Hsps were newly identified. The secondary structures of a further 22 BtaHsps were also predicted. The results of RT-qPCR showed that heat shock could affect the expression of 14 of the 22 Hsps newly identified in this study. Among them, the expression level of six Hsps increased under 42 °C treatment. As the unstudied gene, BtaHsp90A3 had the highest increase rate. Therefore, BtaHsp90A3 was chosen for the RNAi test, and silencing BtaHsp90A3 by RNAi decreased the survival rate of adult B. tabaci at 42 °C. The results indicated that only a few Hsps were involved in the thermal tolerance of host whitefly although many Hsps would response under heat stress. This study conducted a more in-depth and comprehensive identification that demonstrates the evolutionary relationship of BtaHsps and illustrates the response of BtaHsps under the influence of thermal stress in B. tabaci MED.

Characterization of Three Heat Shock Protein Genes in Pieris melete and Their Expression Patterns in Response to Temperature Stress and Pupal Diapause

Simple Summary

Pieris melete, a major pest of crucifers, undergoes obligatory diapause as pupae to survive unfavorable temperature extremes during hot summers and cold winters. Heat shock proteins 70 (Hsp70) participate in this process; however, little is known about the underlying changes in Hsp70 expression both during the summer and winter diapause. The study aimed to investigate expression patterns of Hsp70s (PmHsc70/PmHsp70a, b) in response to diapause and short-term temperature stresses. The results showed that the expression of PmHsc70 and PmHsp70b were upregulated both in summer and winter diapause. Heat shock significantly induced up-regulation of the three genes in both summer and winter diapause. In non-diapause pupae, none of the genes responded to cold or heat stress. Further, it was found that 39 °C for 30 min was the most sensitive heat stress condition for PmHsc70 expressions in summer diapause and all three genes’ expressions in winter diapause. During summer diapause, the expression of the genes was up-regulated in response to high temperature acclimation at 31 °C. Meanwhile, only PmHsp70a and PmHsp70b were up-regulated when acclimated to a low temperature of 4 °C in winter diapause. In conclusion, the current results indicate that PmHsp70s plays a crucial role during both summer and winter diapause, in response to temperature stresses; and our findings may contribute to the increasing knowledge on seasonal diapause adaption.

Abstract

Heat shock protein 70 genes participate in obligatory pupal diapause in Pieris melete to survive unfavorable conditions. In this study, three full-length cDNAs of PmHsc70, PmHsp70a and PmHsp70b were identified, and their expression patterns in response to diapause and short-term temperature stresses were investigated. Summer and winter diapause were induced in the pupae and non-diapause individuals were used as a control. The pupae from each diapause group were subjected to either hot or cold conditions and the expression levels of the HSP genes were measured. Our results showed that up-regulation of PmHsc70 and PmHsp70b were detected both in summer and winter diapause, but not for PmHsp70a. Under cold stress, PmHsp70a and PmHsp70b were upregulated in summer and winter diapause, while heat shock significantly induced upregulation of all three genes. In non-diapause pupae, none of the genes responded to cold or heat stress. Furthermore, we found that incubation at 39 °C for 30 min was the most sensitive heat stress condition for PmHsc70 expression in summer diapause. On the other hand, the same temperature was effective for PmHsc70, PmHsp70a, and PmHsp70b expression in winter diapause. During summer diapause, expression of all three genes was upregulated in response to high-temperature acclimation at 31 °C, but only PmHsp70a and PmHsp70b were upregulated when acclimated to a low temperature of 4 °C in winter diapause. These results suggest that the PmHsc70, PmHsp70a, and PmHsp70b respond differently to pupal diapause and temperature stress, and that PmHsc70 is more sensitive to heat shock than to cold stress.

Effects of heat stress on virus transmission and virus-mediated apoptosis in whitefly Bemisia tabaci

Tomato yellow leaf curl virus (TYLCV), a plant DNA virus of the genus Begomovirus, is transmitted by whiteflies of the Bemisia tabaci species complex in a persistent manner. Our previous study indicated that activation of the apoptosis pathway in whiteflies could facilitate TYLCV accumulation and transmission. Considering that temperature change can influence the spread of insect-borne plant viruses, we focused on plant virus induced-apoptosis to investigate the underlying mechanism of temperature regulation on plant virus transmission via an insect vector. We found that heat stress (40°C) on whiteflies could facilitate TYLCV accumulation and increase transmission to tomato plants. Despite upregulation of caspase-1 and caspase-3 gene expression, heat stress failed to induce an increase in the activation of cleaved caspase-3 and DNA fragmentation in TYLCV-infected whiteflies. However, our data failed to determine the role of heat stress in apoptosis modulation of insect-plant virus interplay while still providing clues to understand insect vectors and their transmitted plant viruses.

Functional characterization of four Hsp70 genes involved in high-temperature tolerance in Aphis aurantii (Hemiptera: Aphididae)

The tea aphid, Aphis aurantii (Boyer de Fonscolombe), is a serious pest that can infest many economically important plants. Tea aphids damage plants by directly sucking phloem sap, transmitting viruses, and secreting honeydew to cause sooty mold. At present, tea aphids has become one of the most important pests in tropical and subtropical tea plants. The heat shock protein 70 (Hsp70) is a key protein involved in heat stress tolerance. In this study, we cloned four Hsp70 genes that are highly expressed in tea aphids after heat shock. Bioinformatic analysis of the deduced amino acid sequences showed that these four AaHsp70s had a close genetic relationship to Hsp70 in Hemiptera insects and shared a conserved ATPase domain. After incubation at low (14 °C) or high (36 °C) temperature, the expression of four AaHsp70s was significantly up-regulated compared to the control (25 °C); however, the up-regulation of the AaHsp70s in the low-temperature treatment was far less than that of the high-temperature treatment. The ATPase activity of the four purified recombinant AaHsp70 proteins after high-temperature treatment was significantly increased compared to the control. In addition, these proteins effectively improved the heat tolerance of Escherichia coli in vivo. Our data indicate that AaHsp701, AaHsp702, AaHsp703, AaHsp704 play important roles in response to the high-temperature tolerance in tea aphids.

Comparative transcriptome analysis of Tamarixia radiata (Hymenoptera: Eulophidae) reveals differentially expressed genes upon heat shock

Temperature is a key parameter that affects insect population, abundance, and distribution in tropical and subtropical regions. Tamarixia radiata Waterson (Hymenoptera: Eulophidae) is a species-specific ectoparasitoid widely used as a biological control agent for the major citrus pest Diphornia citri Kuwayama (Hemiptera: Liviidea). To date, T. radiata response to high temperature at the molecular level still is unclear. In this study, we conducted a comparative analysis of the transcriptomes of T. radiata exposed at 25 °C and 38 °C for 15 min. A total of 51,072 unigenes were obtained, 22,413 annotated with a mean length of 1054 bp. Differential expression analysis showed that 502 genes were identified, including 476 genes significantly up-regulated and 26 genes down-regulated after heat stress exposure. The Gene Ontology analysis showed that most enriched DEGs are categorized into "cellular process", "metabolic process" and "DNA binding." In addition, "Lysosome," "Longevity regulating pathway-multiple species," and "starch and sucrose metabolism" were highly enriched in Kyoto Encyclopedia of Genes and Genomes pathways. Transcriptome analyses showed that heat stress significantly induced the transcription of the molecular chaperone, immune response, stress signaling transduction, and oxidation resistance, including highly expressed heat shock proteins, ATPases, and detoxifying enzymes. Furthermore, the expression patterns of thirteen genes including heat shock proteins (HSP), glutathione S-transferase (GST) and cytochrome P450 were consistent with the transcriptome results obtained through qRT-PCR. Together, our results provided a comprehensive study of the molecular response of T. radiata to heat stress and provides new insight for the future functional validation of heat resistance-related genes.

Genome-wide identification and expression analysis of the Hsp gene superfamily in Asian long-horned beetle (Anoplophora glabripennis)

The environmental adaptability of insects has been a key focus of ecological research. As molecular chaperones, Heat shock proteins (Hsps) play an important role in insect responses to environmental stress. Anoplophora glabripennis is a destructive pest of broad-leaved trees such as poplars. The ability to adapt to low temperature is an important factor for successful colonization of A. glabripennis in new diffusion area. However, the roles of Hsp in the stress responses in A. glabripennis have not been established. In this study, we identified 47 Hsp genes, including 3 Hsp90, 14 Hsp70, 9 Hsp60, and 21 sHsp genes. The Hsp gene family expanded substantially in A. glabripennis. The differences in expression patterns may be related to the type and intensity of stress. Larval overwintering transcriptomes showed that 13 Hsp genes were not induced during overwintering and 21 Hsp genes were involved in the regulation of life activities under non-stress conditions. In a quantitative RT-PCR analysis, AglaHsp90-2 responded more quickly under gradient cooling treatments; AglaHsp90-2 and AglaHsp90-3 were sensitive to treatment at 0 °C for 6 h under instantaneous cooling. Our results provide a theoretical basis for clarifying the molecular mechanism of Hsp genes in A. glabripennis in responsing to environmental stresses.

Genome-Wide Identification and Expression Analysis of Heat Shock Protein 70 (HSP70) Gene Family in Pumpkin (Cucurbita moschata) Rootstock under Drought Stress Suggested the Potential Role of these Chaperones in Stress Tolerance

Heat shock protein 70s (HSP70s) are highly conserved proteins that are involved in stress responses. These chaperones play pivotal roles in protein folding, removing the extra amounts of oxidized proteins, preventing protein denaturation, and improving the antioxidant system activities. This conserved family has been characterized in several crops under drought stress conditions. However, there is no study on HSP70s in pumpkin (Cucurbita moschata). Therefore, we performed a comprehensive analysis of this gene family, including phylogenetic relationship, motif and gene structure analysis, gene duplication, collinearity, and promoter analysis. In this research, we found 21 HSP70s that were classified into five groups (from A to E). These genes were mostly localized in the cytoplasm, chloroplast, mitochondria, nucleus, and endoplasmic reticulum (ER). We could observe more similarity in closely linked subfamilies in terms of motifs, the number of introns/exons, and the corresponding cellular compartments. According to the collinearity analysis, gene duplication had occurred as a result of purifying selection. The results showed that the occurrence of gene duplication for all nine gene pairs was due to segmental duplication (SD). Synteny analysis revealed a closer relationship between pumpkin and cucumber than pumpkin and Arabidopsis. Promoter analysis showed the presence of various cis-regulatory elements in the up-stream region of the HSP70 genes, such as hormones and stress-responsive elements, indicating a potential role of this gene family in stress tolerance. We furtherly performed the gene expression analysis of the HSP70s in pumpkin under progressive drought stress. Pumpkin is widely used as a rootstock to improve stress tolerance, as well as fruit quality of cucumber scion. Since stress-responsive mobile molecules translocate through vascular tissue from roots to the whole plant body, we used the xylem of grafted materials to study the expression patterns of the HSP70 (potentially mobile) gene family. The results indicated that all CmoHSP70s had very low expression levels at 4 days after stress (DAS). However, the genes showed different expression patterns by progressing he drought period. For example, the expression of CmoHSP70-4 (in subgroup E) and CmoHSP70-14 (in subgroup C) sharply increased at 6 and 11 DAS, respectively. However, the expression of all genes belonging to subgroup A did not change significantly in response to drought stress. These findings indicated the diverse roles of this gene family under drought stress and provided valuable information for further investigation on the function of this gene family, especially under stressful conditions.

A heat shock protein protects against oxidative stress induced by lambda-cyhalothrin in the green peach aphid Myzus persicae

Lambda-cyhalothrin (LCT) is a pyrethroid insecticide widely used to control insect pests. Insect exposure to LCT may cause abnormal accumulation of reactive oxygen species (ROS) and result in oxidative damage. Heat shock proteins (HSPs) may help protect against oxidative stress. However, little is known about the role of HSPs in response to LCT in the green peach aphid, Myzus persicae. This insect is an important agricultural pest causing severe yield losses in crops. In this study, we characterized a cDNA sequence (MpHsp70) encoding a member of the HSP70 family in M. persicae. MpHsp70 encoded a 623 amino acid protein putatively localized in the cytosol. The highest expression level of MpHsp70 occurred in fourth-instar nymphs. Treatment of M. persicae with LCT resulted in oxidative stress and significantly increased H₂O₂ and malondialdehyde levels. This led to an elevated transcription level of MpHsp70. Injection of H₂O₂ into M. persicae also upregulated the MpHsp70 expression level, suggesting that MpHsp70 is responsive to ROS, particularly H₂O₂, induced by LCT. Recombinant MpHSP70 protein was expressed in Escherichia coli. E. coli cells overexpressing MpHSP70 exhibited significant tolerance to H₂O₂ and the ROS generators, cumene hydroperoxide and paraquat. This indicated that MpHSP70 protects against oxidative stress. Furthermore, knockdown of MpHsp70 by RNA interference resulted in increased susceptibility in apterous adults of M. persicae to LCT. These findings indicate that MpHsp70 plays an important role in defense against LCT-induced oxidative stress and insecticide susceptibility in M. persicae.

Heat Shock Protein 70 Family in Response to Multiple Abiotic Stresses in the Silkworm

Simple Summary

Heat shock protein 70 family is widely distributed in all the organisms, which plays important roles in protein folding and preventing protein denaturation. Heat or cold stress response has been studied in some insects, but there is a lack of systematic investigation on the response of the same species to multiple stressors. Here, we performed genome-wide identification of heat shock protein 70 family in the silkworm, Bombyx mori. Using the silkworm as a model, the transcription profiles of all the genes against heat, cold, and pesticides were studied. Our findings would provide insights into the functional diversification of heat shock proteins 70 in insects.

Abstract

The 70 kDa heat shock proteins play important roles in protecting organisms against environmental stresses, which are divided into stress-inducible forms (HSP70s) and heat shock cognates (HSC70s). In this study, heat shock protein 70 family was identified in the whole genome of the silkworm. Based on the known nomenclature and phylogenetic analysis, four HSP70s and five HSC70s were classified. Relatively, heat shock cognates were more conservative and were constitutively expressed in various tissues of the silkworm larvae. Under thermal (37 °C and 42 °C) and cold (2 °C) stresses, the expressions of HSP70–1, HSP70–2, and HSP70–3 were up-regulated, and the highest induction reached 4147.3, 607.1, and 1987.3 times, respectively. Interestingly, HSC70–1, HSC70–4, and HSC70–5 also showed slight induced expressions in the fat body and/or midgut under thermal stresses. In addition, the expression of HSP70–1 was induced by dichlorvos and phoxim insecticides, while most HSC70 genes were inhibited. The results suggested that stress-inducible forms play more important roles in adaptation to various stresses than HSC70s.

Evolution and genomic organization of the insect sHSP gene cluster and coordinate regulation in phenotypic plasticity

Background

Conserved syntenic gene complexes are rare in Arthropods and likely only retained due to functional constraint. Numerous sHSPs have been identified in the genomes of insects, some of which are located clustered in close proximity. Previous phylogenetic analyses of these clustered sHSP have been limited to a small number of holometabolous insect species and have not determined the pattern of evolution of the clustered sHSP genes (sHSP-C) in insect or Arthropod lineages.

Results

Using eight genomes from representative insect orders and three non-insect arthropod genomes we have identified that a syntenic cluster of sHSPs (sHSP-C) is a hallmark of most Arthropod genomes. Using 11 genomes from Hymenopteran species our phylogenetic analyses have refined the evolution of the sHSP-C in Hymenoptera and found that the sHSP-C is order-specific with evidence of birth-and-death evolution in the hymenopteran lineage. Finally we have shown that the honeybee sHSP-C is co-ordinately expressed and is marked by genomic features, including H3K27me3 histone marks consistent with coordinate regulation, during honeybee ovary activation.

Conclusions

The syntenic sHSP-C is present in most insect genomes, and its conserved coordinate expression and regulation implies that it is an integral genomic component of environmental response in arthropods.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01885-8.

Characterization, expression profiling, and thermal tolerance analysis of heat shock protein 70 in pine sawyer beetle, Monochamus alternatus hope (Coleoptera: Cerambycidae)

Monochamus alternatus Hope (Coleoptera: Cerambycidae) warrants attention as a dominant transmission vector of the pinewood nematode, and it exhibits tolerance to high temperature. Heat shock protein 70 (HSP70) family members, including inducible HSP70 and heat shock cognate protein 70 (HSC70), are major contributors to the molecular chaperone networks of insects under heat stress. In this regard, we specifically cloned and characterized three MaltHSP70s and three MaltHSC70s. Bioinformatics analysis on the deduced amino acid sequences showed these genes, having close genetic relationships with HSP70s of Coleopteran species, collectively shared conserved signature structures and ATPase domains. Subcellular localization prediction revealed the HSP70s of M. alternatus were located not only in the cytoplasm and endoplasmic reticulum but also in the nucleus and mitochondria. The transcript levels of MaltHSP70s and MaltHSC70s in each state were significantly upregulated by exposure to 35-50°C for early 3 h, while MaltHSP70s reached a peak after exposure to 45°C for 2-3 h in contrast to less-upregulated MaltHSC70s. In terms of MaltHSP70s, the expression threshold in females was lower than that in males. Also, both fat bodies and Malpighian tubules were the tissues most sensitive to heat stress in M. alternatus larvae. Lastly, the ATPase activity of recombinant MaltHSP70-2 in vitro remained stable at 25-40°C, and this recombinant availably enhanced the thermotolerance of Escherichia coli. Overall, our findings unraveled HSP70s might be the intrinsic mediators of the strong heat tolerance of M. alternatus due to their stabilized structure and bioactivity.

Functional Analysis of Alkaline Phosphatase in Whitefly Bemisia tabaci (Middle East Asia Minor 1 and Mediterranean) on Different Host Plants

Alkaline phosphatases (ALPs: EC 3.1.3.1) are ubiquitous enzymes and play crucial roles in the fundamental phosphate uptake and secretory processes. Although insects are regarded as the most diverse group of organisms, the current understanding of ALP roles in insects is limited. As one type of destructive agricultural pest, whitefly Bemisia tabaci, a phloem feeder and invasive species, can cause extensive crop damage through feeding and transmitting plant diseases. In this study, we retrieved five ALP genes in MEAM1 whitefly, nine ALP genes in MED whitefly via comparative genomics approaches. Compared with nine other insects, whiteflies' ALP gene family members did not undergo significant expansion during insect evolution, and whiteflies' ALP genes were dispersed. Moreover, whiteflies' ALP gene family was conserved among insects and emerged before speciation via phylogenetic analysis. Whiteflies' ALP gene expression profiles presented that most ALP genes have different expression patterns after feeding on cotton or tobacco plants. Female/male MED whiteflies possessed higher ALP activities on both cotton and tobacco plants irrespective of sex, relative to MEAM1 whiteflies. Meanwhile, adult MED whiteflies possessed higher ALP activity in both whole insect and salivary samples, relative to MEAM1 whiteflies. We also found that both MED and MEAM1 whiteflies could upregulate ALP activities after feeding on cotton compared with feeding on tobacco plants. These findings demonstrated the functions of whiteflies ALPs and will assist the further study of the genomic evolution of insect ALPs.

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.

A plant reovirus hijacks endoplasmic reticulum-associated degradation machinery to promote efficient viral transmission by its planthopper vector under high temperature conditions

In the field, many insect-borne crop viral diseases are more suitable for maintenance and spread in hot-temperature areas, but the mechanism remains poorly understood. The epidemic of a planthopper (Sogatella furcifera)-transmitted rice reovirus (southern rice black-streaked dwarf virus, SRBSDV) is geographically restricted to southern China and northern Vietnam with year-round hot temperatures. Here, we reported that two factors of endoplasmic reticulum-associated degradation (ERAD) machinery, the heat shock protein DnaJB11 and ER membrane protein BAP31, were activated by viral infection to mediate the adaptation of S. furcifera to high temperatures. Infection and transmission efficiencies of SRBSDV by S. furcifera increased with the elevated temperatures. We observed that high temperature (35°C) was beneficial for the assembly of virus-containing tubular structures formed by nonstructural protein P7-1 of SRBSDV, which facilitates efficient viral transmission by S. furcifera. Both DnaJB11 and BAP31 competed to directly bind to the tubule protein P7-1 of SRBSDV; however, DnaJB11 promoted whereas BAP31 inhibited P7-1 tubule assembly at the ER membrane. Furthermore, the binding affinity of DnaJB11 with P7-1 was stronger than that of BAP31 with P7-1. We also revealed that BAP31 negatively regulated DnaJB11 expression through their direct interaction. High temperatures could significantly upregulate DnaJB11 expression but inhibit BAP31 expression, thereby strongly facilitating the assembly of abundant P7-1 tubules. Taken together, we showed that a new temperature-dependent protein quality control pathway in the ERAD machinery has evolved for strong activation of DnaJB11 for benefiting P7-1 tubules assembly to support efficient transmission of SRBSDV in high temperatures. We thus deduced that ERAD machinery has been hitchhiked by insect-borne crop viruses to enhance their transmission in tropical climates.

Evaluation of Reference Genes and Age Estimation of Forensically Useful Aldrichina grahami (Diptera: Calliphoridae) During Intrapuparial Period

The minimum postmortem interval (PMImin) could be evaluated from the developmental stage of forensically important insects colonize a corpse, such as blow flies (Diptera: Calliphoridae). Unlike larvae, the developmental stage of which is well established according to their morphology, estimating the age of pupae is proven to be challenging. Recently, several studies reported the regulation of special genes during the development of blow fly pupae. However, gene regulation in Aldrichina grahami during the intrapuparial period remains to be studied. Therefore, we set out to investigate the mRNA levels of heat shock protein 23 (Hsp23), heat shock protein 24 (Hsp24), and 1_16 during the metamorphosis of A. grahami pupae. First, we examined seven candidate reference genes (ribosomal protein 49 (RP49), 18S ribosomal RNA (18S rRNA), 28S ribosomal RNA (28S rRNA), beta-tubulin at 56D (β-tubulin), Ribosomal protein L23 (RPL23), glutathione S-transferase (GST1), and Actin. Three widely used algorithms (NormFinder, BestKeeper, and geNorm) were applied to evaluate the mRNA levels of reference gene candidates in puparium at three stable temperatures (15, 22, and 27°C). Next, mRNA expression of Hsp23, Hsp24, and 1_16 during A. grahami metamorphosis was examined. We demonstrated that mRNA expression levels of Hsp23, Hsp24, and 1_16 showed time-specific regulation. In summary, our study identified three gene markers for the intrapuparial period of A. grahami and might provide a potential application in PMImin estimation.

Induction of Heat Shock Protein Genes is the Hallmark of Egg Heat Tolerance in Agasicles hygrophila (Coleoptera: Chrysomelidae)

Insects are ecotothermic organisms. Their development, survival, reproduction as well as distribution and abundance are affected by temperature. Heat shock protein (HSP) gene expression is closely associated with temperature variation and influences the adaptation of organisms to adverse environments. The beetle Agasicles hygrophila has successfully been used for biological control of the invasive plant alligator weed (Alternanthera philoxeroides). As A. hygrophila populations are substantially inhibited by high temperatures in the summer, increasing global temperatures may limit the efficacy of this control agent. We previously established that A. hygrophila eggs have low tolerance to heat and this factored into the decreased numbers of A. hygrophila beetles at temperatures of 37.5°C and above. Here, we identified 26 HSP genes in A. hygrophila and examined the relationship between the transcript levels of these genes and heat tolerance. The temperature at which the expression of these 21 HSP genes peaked (Tpeak) was 37.5°C, which is in line with the limit of the high temperatures that A. hygrophila eggs tolerate. Therefore, we speculate that the Tpeak of HSP gene expression in eggs indicates the upper limit of temperatures that A. hygrophila eggs tolerate. This study identifies HSP genes as potential robust biomarkers and emphasizes that determining species' heat tolerance in their natural habitats remains an important consideration for biocontrol. HSP gene expression data provide information about a species' heat tolerance and may be used to predict its geographical distribution.

Preferential gene retention increases the robustness of cold regulation in Brassicaceae and other plants after polyploidization

Cold stress profoundly affects plant growth and development and is a key factor affecting the geographic distribution and evolution of plants. Plants have evolved adaptive mechanisms to cope with cold stress. Here, through the genomic analysis of Arabidopsis, three Brassica species and 17 other representative species, we found that both cold-related genes (CRGs) and their collinearity were preferentially retained after polyploidization followed by genome instability, while genome-wide gene sets exhibited a variety of other expansion mechanisms. The cold-related regulatory network was increased in Brassicaceae genomes, which were recursively affected by polyploidization. By combining our findings regarding the selective retention of CRGs from this ecological genomics study with the available knowledge of cold-induced chromosome doubling, we hypothesize that cold stress may have contributed to the success of polyploid plants through both increasing polyploidization and selectively maintaining CRGs during evolution. This hypothesis requires further biological and ecological exploration to obtain solid supporting evidence, which will potentially contribute to understanding the generation of polyploids and to the field of ecological genomics.

Dynamic transcriptome profiling exploring cold tolerance in forensically important blow fly, Aldrichina grahami (Diptera: Calliphoridae)

Background

Aldrichina grahami (Diptera: Calliphoridae) is a forensically important fly, which has been widely applied to practical legal investigations. Unlike other necrophagous flies, A. grahami exhibits cold tolerance which helps to maintain its activity during low-temperature months, when other species are usually not active. Hence, A. grahami is considered an important forensic insect especially in cold seasons. In this study, we aim to explore the molecular mechanisms of cold tolerance of A. grahami through transcriptome.

Results

We collected eggs and larvae (first-instar, second-instar and third-instar) at three different temperatures (4 °C, 12 °C and 20 °C) and performed RNA-seq analyses. The differentially expressed genes (DEGs) associated with the cold-tolerance were screened out. The Venn analysis of DEGs from egg to third-instar larvae at three different temperatures showed there were 9 common genes. Candidate biological processes and genes were identified which refer to growth, and development of different temperatures, especially the chitin and cuticle metabolic process. The series-clusters showed crucial and unique trends when the temperature changed. Moreover, by comparing the results of growth and developmental transcriptomes from different temperatures, we found that DEGs belonging to the family of larval cuticle proteins (LCP), pupal cuticle protein (CUP), and heat shock proteins (HSP) have certain differences.

Conclusions

This study identified functional genes and showed differences in the expression pattern of diverse temperatures. The DEGs series-clusters with increasing or decreasing trends were analyzed which may play an important role in cold-tolerance. Moreover, the findings in LCP, CUP and HSP showed more possible modulations in a cold environment. This work will provide valuable information for the future investigation of the molecular mechanism of cold tolerance in A. grahami.

Comparative Transcriptome Analysis of the Heat Stress Response in Monochamus alternatus Hope (Coleoptera: Cerambycidae)

Temperature is a critical factor of insect population abundance and distribution. Monochamus alternatus Hope (Coleoptera: Cerambycidae) is a significant concern since it is transmitted vector of the pinewood nematode posing enormous economic and environmental losses. This pest shows tolerance to heat stress, especially extremely high temperatures. Exposing for 6, 12, 24, 48, or 96 h, the 50% median lethal temperatures (Ltem₅₀) for fourth-instar larvae were 47.5, 45.5, 43.9, 43.4, and 42.3°C, respectively. A total of 63,360 unigenes were obtained from complementary DNA libraries of M. alternatus fourth-instar larvae (kept at 25°C and exposed to 40°C for 3 h) and annotated with six databases. Five hundred sixty-one genes were significantly upregulated, and 245 genes were downregulated after heat stress. The Gene Ontology enrichment analysis showed that most different expression genes are categorized into “protein folding” and “unfold protein binding” terms. In addition, “Longevity regulating pathway-multiple species,” “Antigen processing and presentation” as well as “MAPK signaling pathway” were significantly enriched Kyoto Encyclopedia of Genes and Genomes pathways. Further analysis of different expression genes showed that metabolism processes were suppressed, while ubiquitin proteolytic system, heat shock proteins, immune response, superoxide dismutase, cytochrome P450s, and aldehyde dehydrogenase were induced after heat shock. The stress signaling transduction pathways such as MAPK, Hippo, and JAK-STAT might be central convergence points in M. alternatus heat tolerance mechanism. The expression levels from quantitative real-time PCR of 13 randomly selected genes were consistent with the transcriptome results. These results showed that M. alternatus possessed strong heat tolerance and genes related to protein activity, immune response, and signal transduction composed of a complicated heat tolerance mechanism of M. alternatus. This research provided new insights into the mechanisms of thermal tolerance in other insects and aided in exploring the function of heat resistance-related genes.

Transcription Analysis of the Stress and Immune Response Genes to Temperature Stress in Ostrinia furnacalis

Ostrinia furnacalis is one of the most important pests on maize. O. furnacalis larvae are frequently exposed to the temperature challenges such as high temperature in summer and cold temperature in winter in the natural environment. High and low temperature stress, like any abiotic stress, impairs the physiology and development of insects. Up to now, there is limited information about gene regulation and signaling pathways related to the high and cold stress response in O. furnacalis. High-throughput sequencing of transcriptome provides a new approach for detecting stress and immune response genes under high and low temperature stresses in O. furnacalis. In the present study, O. furnacalis larvae were treated with the temperature at 8 and 40°C, and the responses of O. furnacalis larvae to the temperature stress were investigated through RNA-sequencing and further confirmation. The results showed that immune responses were up-regulated in larvae by the cold stress at 8°C while some stress response genes, such as HSP family, GST-2, Bax inhibitor and P450, were significantly increased at 40°C. Furthermore, quantitative real time polymerase chain reaction were performed to quantify the expression levels of immune related genes, such as PGRP-LB, antimicrobial peptides, lysozyme, serine protease and stress response genes such as small HSPs and HSP90, and the expression levels of these genes were similar to the RNA-seq results. In addition, the iron storage protein Ferritin was found to be involved in the response to temperature stress, and the changes of total iron concentration in the hemolymph were, in general, consistent with the expression levels of Ferritin. Taken together, our results suggested that the stress response genes were involved in the defense against the heat stress at 40°C, and the immune responses triggered by cold stress might provide protection for larvae from cold stress at 8°C. More interestingly, our results showed that during the responses to temperature stress, the total iron concentration in hemolymph regulated by Ferritin increased, which might help O. furnacalis in surviving the low and high temperature stress.

Ambient temperature-mediated enzymic activities and intestinal microflora in Lymantria dispar larvae

To understand how ambient temperature affect the gypsy moth larvae, and provide a theoretical basis for pest control in different environments. Fourth instar gypsy moth larvae were incubating for 3 hr at 15℃, 20℃, 25℃, 30℃, 35℃, and 40℃, respectively. Afterward, digestive and antioxidant enzyme activities, total antioxidant capacity, and intestinal microflora community were analyzed to reveal how the caterpillars respond to ambient temperature stress. Results showed that both digestive and antioxidant enzymes were regulated by the ambient temperature. The optimum incubation temperatures of protease, amylase, trehalase, and lipase in gypsy moth larvae were 30℃, 25℃, and 20℃, respectively. When the incubation temperature was deviated optimum temperatures, digestive enzyme activities would be downregulated depending on the extent of temperature stress. In addition, glutathione S-transferase, peroxidase, catalase, and polyphenol oxidase would be activated under a sufferable temperature stress, but superoxide dismutase and carboxylesterase (CarE) would be inhibited. In addition, results showed that the top two abundant phyla were Proteobacteria and Firmicutes. The phylum Firmicutes abundance was decreased and phylum Proteobacteria abundance was increased by ambient temperature stress. Moreover, it suggested that gypsy moth caterpillars at different ambient temperature mainly differed from each other by Escherichia-Shigella and Bifidobacterium in control, Acinetobacter in T15, and Lactobacillus in T40, respectively.

Conservation of transcriptional elements in the obligate symbiont of the whitefly Bemisia tabaci

Background

Bacterial symbiosis is widespread in arthropods, especially in insects. Some of the symbionts undergo a long-term co-evolution with the host, resulting in massive genome decay. One particular consequence of genome decay is thought to be the elimination of transcriptional elements within both the coding region and intergenic sequences. In the whitefly Bemisia tabaci species complex, the obligate symbiont Candidatus Portiera aleyrodidarum is of vital importance in nutrient provision, and yet little is known about the regulatory capacities of it.

Methods

Portiera genomes of two whitefly species in China were sequenced and assembled. Gene content of these two Portiera genomes was predicted, and then subjected to Kyoto Encyclopedia of Genes and Genomes pathway analysis. Together with two other Portiera genomes from whitefly species available previously, four Portiera genomes were utilized to investigate regulatory capacities of Portiera, focusing on transcriptional elements, including genes related with transcription and functional elements within the intergenic spacers.

Results

Comparative analyses of the four Portiera genomes of whitefly B. tabaci indicate that the obligate symbionts Portiera is similar in different species of whiteflies, in terms of general genome features and possible functions in the biosynthesis of essential amino acids. The screening of transcriptional factors suggests compromised ability of Portiera to regulate the essential amino acid biosynthesis pathways. Meanwhile, thermal tolerance ability of Portiera is indicated with the detection of a σ³² factor, as well as two predicted σ³² binding sites. Within intergenic spacers, functional elements are predicted, including 37 Shine-Dalgarno sequences and 34 putative small RNAs.

Genomewide identification and analysis of heat-shock proteins 70/110 to reveal their potential functions in Chinese soft-shelled turtle Pelodiscus sinensis

Heat-shock proteins 70/110 (Hsp70/110) are vital molecular chaperones and stress proteins whose expression and production are generally induced by extreme temperatures or external stresses. The Hsp70/110 family is largely conserved in diverse animals. Although many reports have studied and elaborated on the characteristics of Hsp70/110 in various species, the systematic identification and analysis of Hsp70/110 are still poor in turtles. In this study, a genomewide search was performed, and 18 candidate PsHSP70/110 family genes were identified in Chinese soft-shelled turtle, Pelodiscus sinensis. These PsHSP70/110 proteins contained the conserved "heat shock protein 70" domain. Phylogenetic analysis of PsHSP70/110 and their homologs revealed evolutionary conservation of Hsp70/110 across different species. Tissue-specific expression analysis showed that these PsHSP70/110 genes were differentially expressed in different tissues of P. sinensis. Furthermore, to examine the putative biological functions of PsHSP70/110, the dynamic expression of PsHSP70/110 genes was analyzed in the testis of P. sinensis during seasonal spermatogenesis following germ cell apoptosis. Notably, genes such as PsHSPA1B-L, PsHSPA2, and PsHSPA8 were significantly upregulated in P. sinensis testes along with a seasonal decrease in apoptosis. Protein interaction prediction revealed that PsHSPA1B-L, PsHSPA2, and PsHSPA8 may interact with each other and participate in the MAPK signaling pathway. Moreover, immunohistochemical analysis showed that PsHSPA1B-L, PsHSPA2, and PsHSPA8 protein expression was associated with seasonal temperature variation. The expression profiling and interaction relationships of the PsHSPA1B-L, PsHSPA2, and PsHSPA8 proteins implied their potential roles in inhibiting the apoptosis of germ cells in P. sinensis. These results provide insights into PsHSP70/110 functions and will serve as a rich resource for further investigation of HSP70/110 family genes in P. sinensis and other turtles.

Identification and Expression Analysis of Four Small Heat Shock Protein Genes in Cigarette Beetle, Lasioderma serricorne (Fabricius)

Small heat shock proteins (sHsps) are molecular chaperones that play crucial roles in the stress adaption of insects. In this study, we identified and characterized four sHsp genes (LsHsp19.4, 20.2, 20.3, and 22.2) from the cigarette beetle, Lasioderma serricorne (Fabricius). The four cDNAs encoded proteins of 169, 180, 181, and 194 amino acids with molecular weights of 19.4, 20.2, 20.3, and 22.2 kDa, respectively. The four LsHsp sequences possessed a typical sHsp domain structure. Quantitative real-time PCR analyses revealed that LsHsp19.4 and 20.3 transcripts were most abundant in pupae, whereas the transcript levels of LsHsp20.2 and 22.2 were highest in adults. Transcripts of three LsHsp genes were highly expressed in the larval fat body, whereas LsHsp20.2 displayed an extremely high expression level in the gut. Expression of the four LsHsp genes was dramatically upregulated in larvae exposed to 20-hydroxyecdysone. The majority of the LsHsp genes were significantly upregulated in response to heat and cold treatments, while LsHsp19.4 was insensitive to cold stress. The four genes were upregulated when challenged by immune triggers (peptidoglycan isolated from Staphylococcus aureus and from Escherichia coli 0111:B4). Exposure to CO₂ increased LsHsp20.2 and 20.3 transcript levels, but the LsHsp19.4 transcript level declined. The results suggest that different LsHsp genes play important and distinct regulatory roles in L. serricorne development and in response to diverse stresses.