Transcriptome analysis of the response of silkworm to drastic changes in ambient temperature

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.

Integrated transcriptomic and proteomic analyses reveal the molecular mechanism underlying the thermotolerant response of Spodoptera frugiperda

Spodoptera frugiperda (Lepidoptera: Noctuidae) is a highly destructive invasive pest with remarkable adaptability to extreme climatic conditions, posing a substantial global threat. Although the effects of temperature stress on the biological and ecological properties of S. frugiperda have been elucidated, the molecular mechanisms underlying its responses remain unclear. Herein, we combined transcriptomic and proteomic analyses to explore the key genes and proteins involved in thermotolerance regulation in S. frugiperda larvae at 42 °C. Overall, 1528 differentially expressed genes (DEGs) and 154 differentially expressed proteins (DEPs) were identified in S. frugiperda larvae under heat stress, including antioxidant enzymes, heat shock proteins (Hsps), cytochrome P450s, starch and sucrose metabolism genes, and insulin signaling pathway genes, indicating their involvement in heat tolerance regulation. Correlation analysis of DEGs and DEPs revealed that seven and eight had the same and opposite expression profiles, respectively. After nanocarrier-mediated RNA interference knockdown of SfHsp29, SfHsp20.4, SfCAT, and SfGST, the body weight and mortality of S. frugiperda larvae significantly decreased and increased under heat stress, respectively. This indicates that SfHsp29, SfHsp20.4, SfCAT, and SfGST play a crucial role in the thermotolerance of S. frugiperda larvae. These results provide insight into the mechanism of heat tolerance in S. frugiperda.

Chromosomal-Level Reference Genome for the Chinese Endemic Pygmy Grasshopper, Zhengitettix transpicula, Sheds Light on Tetrigidae Evolution and Advancing Conservation Efforts

The pygmy grasshopper, Zhengitettix transpicula, is a Chinese endemic species with an exceedingly limited distribution and fragile population structure, rendering it vulnerable to extinction. We present a high-continuity, chromosome-scale reference genome assembly to elucidate this species' distinctive biology and inform conservation. Employing an integrated sequencing approach, we achieved a 970.40 Mb assembly with 96.32% coverage across seven pseudo-chromosomes and impressive continuity (N50 > 220 Mb). Genome annotation achieves identification with 99.2% BUSCO completeness, supporting quality. Comparative analyses with 14 genomes from Orthoptera-facilitated phylogenomics and revealed 549 significantly expanded gene families in Z. transpicula associated with metabolism, stress response, and development. However, genomic analysis exposed remarkably low heterozygosity (0.02%), implying a severe genetic bottleneck from small, fragmented populations, characteristic of species vulnerable to extinction from environmental disruptions. Elucidating the genetic basis of population dynamics and specialization provides an imperative guideline for habitat conservation and restoration of this rare organism. Moreover, divergent evolution analysis of the CYP305m2 gene regulating locust aggregation highlighted potential structural and hence functional variations between Acrididae and Tetrigidae. Our chromosomal genomic characterization of Z. transpicula advances Orthopteran resources, establishing a framework for evolutionary developmental explorations and applied conservation genomics, reversing the trajectory of this unique grasshopper lineage towards oblivion.

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.

Influence of temperature variation on gene expression and cocoon production in Bombyx mori Linnaeus, 1758 (Lepidoptera: Bombycidae)

Silkworms (Bombyx mori) are lepidopterans of economic importance for global silk production. However, factors that directly affect the yield and quality of silkworm cocoon production, such as diseases and temperature fluctuations, cause great economic losses. Knowing how they respond to rearing temperature during the most critical stage of their life cycle (i.e., fifth instar) could provide information on their adaptation and improve silk production. In the current work, we analyzed transcriptional data from two groups of B. mori that were reared at 26 °C and 34 °C throughout the fifth instar. The silkworms and cocoons were weighed. In total, 3115 transcripts were differentially expressed (DE; including 1696 down-regulated and 1419 up-regulated) among the 29,157 sequences found by transcriptome assembly. We emphasize the genes associated with immunological response, transcription factors, silk biosynthesis, and heat shock proteins, among the DE transcripts in response to the temperature conditions. Silkworms reared at 34 °C presented a reduced mean body weight (-0.944 g in comparison to the 26 °C group), which had a direct impact on the weight of cocoons formed and the silk conversion rate. These changes were statistically significant when compared to silkworms reared at 26 °C. Mortality rates (6 and 9 %, at 26 °C and 34 °C, respectively) were similar to those obtained in breeding fields. The findings provide information on the biological processes involved in the temperature response mechanism of silkworms, as well as information that may be used in future climatization processes at rearing facilities and in breeding for improved thermotolerance.

Sex dependent transcriptome responses of the diamondback moth, Plutella xylostella L. to cold stress

Temperature has fundamental influences on the performance and distribution of insects. While considerable attention has been devoted to extreme conditions, particularly extreme cold conditions, few studies have investigated effects of mild cold conditions on insects. We examined the transcriptomic changes in mid-fourth instar larvae of both sexes reared at 10 °C and 25 °C to investigate sex-dependent responses of Plutella xylostella to mild cold stress. There were 624 differentially expressed genes (DEGs) in females, the majority of which (n = 386) were down-regulated. In males 3239 genes were differentially expressed and the majority (n = 2341) were up-regulated. Only 280 DEGs were common to both sexes. In females, there were no DEGs encoding heat shock or cold shock proteins, but six of these DEGs were found in males. These differences suggest that females and males might adopt some different strategies to cope with cold stress and/or that they were affected by rearing under cold conditions to different degrees and in different ways. In addition, DEGs encoding antimicrobial peptides, cytochrome P450 monooxygenases, fatty acid-related enzymes, cuticle proteins, myofilament, and hormone-related proteins were found in both sexes under cold stress. The transcriptome study reveals unexpected sex-dependent thermal responses and provides new information of how an insect that does not diapause copes with low temperatures.

From perplexing to predictive: are we ready to forecast insect disease susceptibility in a warming world?

Insects are critical to our ecosystems, but we do not fully understand their future in our warming world. Rising temperatures are affecting insect physiology in myriad ways, including changes to their immune systems and the ability to fight infection. Whether predicted changes in temperature will contribute to insect mortality or success, and the role of disease in their future survival, remains unclear. Although heat can enhance immunity by activating the integrated defense system (e.g. via the production of protective molecules such as heat-shock proteins) and accelerating enzyme activity, heat can also compromise the immune system through energetic-resource trade-offs and damage. The responses to heat are highly variable among species. The reasons for this variability are poorly known, and we are lagging in our understanding of how and why the immune system responds to changes in temperature. In this Commentary, we highlight the variation in insect immune responses to heat and the likely underlying mechanisms. We suggest that we are currently limited in our ability to predict the effects of rising temperatures on insect immunity and disease susceptibility, largely owing to incomplete information, coupled with a lack of tools for data integration. Moreover, existing data are concentrated on a relatively small number of insect Orders. We provide suggestions for a path towards making more accurate predictions, which will require studies with realistic temperature exposures and housing design, and a greater understanding of both the thermal biology of the immune system and connections between immunity and the physiological responses to heat.

Structural insights into the substrate binding of phosphomevalonate kinase from the silkworm, Bombyx mori

Phosphomevalonate kinase (PMK) is an important enzyme involved in the juvenile hormone (JH) biosynthesis pathway that catalyzes the phosphorylation of mevalonate 5-phosphate into mevalonate 5-diphosphate in the mevalonate pathway. Herein, we report the crystal structure of insect PMK from Bombyx mori (BmPMK) at a resolution of 1.60 Å. The overall structure of BmPMK adopts a compact α/β conformation with two parts: the core and lid regions. The interface between the core and lid regions forms a continuous and negatively charged groove to accommodate the substrates. Using computational simulation combined with site-directed mutagenesis and biochemical analysis, we define the binding mode of BmPMK with the cofactor and the substrate, which provides a structural basis for understanding the catalytic mechanism and the design of inhibitors of PMK. Moreover, BmPMK showed the optimal enzyme activity at pH 8.0, and the optimal temperature was 30 °C, using mevalonate 5-phosphate as the substrate. The expression profiles and kinetic analyses of BmPMK indicated that it plays critical role in the control of JH biosynthesis in silkworms. Collectively, these findings provide a better understanding of the structural and biochemical features of insect PMK.

Global Profiling of Genes Expressed in the Silk Glands of Philippine-Reared Mulberry Silkworms (Bombyx mori)

RNA sequencing was used to assemble transcriptome data for Philippine-reared silkworm and compare gene expression profiles of strains reared in high- and low-temperature environments. RNA was isolated from the silk glands of fifth instar larvae and mRNA-enriched libraries were sequenced using Illumina NextSeq 500. Transcriptome reads were assembled using reference-based and de novo assemblers, and assemblies were evaluated using different metrics for transcriptome quality, including the read mapping rate, E90N50, RSEM-eval, and the presence of single-copy orthologs. All transcriptome assemblies were able to reconstruct >40,000 transcripts. Differential expression analysis found 476 differentially expressed genes (DEGs; 222 upregulated, 254 downregulated) in strains reared in different temperatures. Among the top DEGs were myrosinase, heat shock proteins, serine protease inhibitors, dehydrogenases, and regulators of the juvenile hormone. Validation of some of the top DEGs using qPCR supported the findings of the in silico analysis. GO term enrichment analysis reveals an overrepresentation of GO terms related to nucleotide metabolism and biosynthesis, lipid and carbohydrate metabolic processes, regulation of transcription, nucleotide binding, protein binding, metal binding, catalytic activity, oxidoreductase activity, and hydrolase activity. The data provided here will serve as a resource for improving local strains and increasing silk production of Philippine-reared B. mori strains.

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.

Heat shock protein 19.9 (Hsp19.9) from Bombyx mori is involved in host protection against viral infection

Adverse environmental conditions cause serious economic losses in sericulture; Bombyx mori nucleopolyhedrovirus (BmNPV) is the primary biotic stress and high temperature is the major abiotic stress in this industry. B. mori heat shock protein 19.9 (Bmhsp19.9) overexpression was previously demonstrated to protect transgenic silkworm H19.9 against extreme temperature. This study analyzed the role of Bmhsp19.9 in H19.9A and H19.9B silkworm lines and BmE cells infected with BmNPV at regular and high temperatures. qPCR results showed that Bmhsp19.9 expression was upregulated in BmE cells and silkworm after BmNPV challenge. Bmhsp19.9 overexpression significantly inhibited BmNPV proliferation in BmE cells. The viral DNA content was significantly decreased in transgenic H19.9 silkworm compared to the control. These results suggested that Bmhsp19.9 was involved in antiviral immunity against BmNPV. Furthermore, Bmhsp19.9 overexpression protected BmE cells against BmNPV under high temperature shock. This indicates that Bmhsp19.9 is a promising candidate for improving silkworm resistance to biotic and abiotic stresses, thereby reducing sericulture losses.

Toxicological evaluation of transgenic silkworms

Safety of transgenic silkworms must be evaluated before their commercial application. We assessed subacute toxicity using a 28-day feeding study in rats. Eighty rats were evenly allocated into four groups, with each group containing 10 male and 10 female rats. Rats of three groups were fed dried transgenic silkworm H19.9A pupae with overexpressed endogenous Bmhsp19.9, transgenic silkworm A4SOR pupae with overexpressed exogenous SOR, or normal silkworm pupae at a dose of 3.0 g/kg/day, respectively. The fourth group served as a normal feeding control. The body weight, feed consumption, hematology response variables, serum biochemical parameters, organ weights, gross necropsy, and histopathologic of animals were evaluated. No mortality, adverse effects, or major differences in the evaluated parameters were observed in the groups fed transgenic pupae in comparison with the control, suggesting that transgenic silkworms are toxicologically equivalent to normal silkworms and are safe for consumption in rats.

Comparative transcriptome analysis of the rice leaf folder (Cnaphalocrocis medinalis) to heat acclimation

Background

The rice leaf folder Cnaphalocrocis medinalis Güenée is a serious insect pest of rice in Asia. This pest occurs in summer, and it is sensitive to high temperature. However, the larvae exhibit heat acclimation/adaptation. To understand the underlying mechanisms, we established a heat-acclimated strain via multigenerational selection at 39 °C. After heat shock at 41 °C for 1 h, the transcriptomes of the heat-acclimated (S-39) and unacclimated (S-27) larvae were sequenced, using the unacclimated larvae without exposure to 41 °C as the control.

Results

Five generations of selection at 39 °C led larvae to acclimate to this heat stress. Exposure to 41 °C induced 1160 differentially expressed genes (DEGs) between the heat-acclimated and unacclimated larvae. Both the heat-acclimated and unacclimated larvae responded to heat stress via upregulating genes related to sensory organ development and structural constituent of eye lens, whereas the unacclimated larvae also upregulated genes related to structural constituent of cuticle. Compared to unacclimated larvae, heat-acclimated larvae downregulated oxidoreductase activity-related genes when encountering heat shock. Both the acclimated and unacclimated larvae adjusted the longevity regulating, protein processing in endoplasmic reticulum, antigen processing and presentation, MAPK and estrogen signaling pathway to responsed to heat stress. Additionally, the unacclimated larvae also adjusted the spliceosome pathway, whereas the heat-acclimated larvae adjusted the biosynthesis of unsaturated fatty acids pathway when encountering heat stress. Although the heat-acclimated and unacclimated larvae upregulated expression of heat shock protein genes under heat stress including HSP70, HSP27 and CRYAB, their biosynthesis, metabolism and detoxification-related genes expressed differentially.

Conclusions

The rice leaf folder larvae could acclimate to a high temperature via multigenerational heat selection. The heat-acclimated larvae induced more DEGs to response to heat shock than the unacclimated larvae. The changes in transcript level of genes were related to heat acclimation of larvae, especially these genes in sensory organ development, structural constituent of eye lens, and oxidoreductase activity. The DEGs between heat-acclimated and unacclimated larvae after heat shock were enriched in the biosynthesis and metabolism pathways. These results are helpful to understand the molecular mechanism underlying heat acclimation of insects.

Transcriptome analysis of the immune response of silkworm at the early stage of Bombyx mori bidensovirus infection

Bombyx mori bidensovirus (BmBDV) infects silkworm midgut and causes chronic flacherie disease; however, the interaction between BmBDV and silkworm is unclear. Twenty-four hours after BmBDV infection, the midgut was extracted for RNA-seq to analyze the factors associated with BmBDV-invasion and the early antiviral immune response in silkworms. The total reads from each sample were more than 16100000 and the number of expressed genes exceeded 8200. There were 334 upregulated and 272 downregulated differentially expressed genes (DEGs). Gene ontology analysis of DEGs showed that structural constituents of cuticle, antioxidant, and immune system processes were upregulated. Further analysis revealed BmBDV-mediated induction of BmorCPR23 and BmorCPR44, suggesting possible involvement in viral invasion. Antioxidant genes that protect host cells from virus-induced oxidative stress, were significantly upregulated after BmBDV infection. Several genes related to peroxisomes, apoptosis, and autophagy-which may be involved in antiviral immunity-were induced by BmBDV. These results provide insights into the mechanism of BmBDV infection and host defense.

Distinct Functions of Bombyx mori Peptidoglycan Recognition Protein 2 in Immune Responses to Bacteria and Viruses

Peptidoglycan recognition protein (PGRP) is an important pattern recognition receptor in innate immunity that is vital for bacterial recognition and defense in insects. Few studies report the role of PGRP in viral infection. Here we cloned two forms of PGRP from the model lepidopteran Bombyx mori: BmPGRP2-1 is a transmembrane protein, whereas BmPGRP2-2 is an intracellular protein. BmPGRP2-1 bound to diaminopimelic acid (DAP)-type peptidoglycan (PGN) to activate the canonical immune deficiency (Imd) pathway. BmPGRP2-2 knockdown reduced B. mori nucleopolyhedrovirus (BmNPV) multiplication and mortality in cell lines and in silkworm larvae, while its overexpression increased viral replication. Transcriptome and quantitative PCR (qPCR) results confirmed that BmPGRP2 negatively regulated phosphatase and tensin homolog (PTEN). BmPGRP2-2 expression was induced by BmNPV, and the protein suppressed PTEN-phosphoinositide 3-kinase (PI3K)/Akt signaling to inhibit cell apoptosis, suggesting that BmNPV modulates BmPGRP2-2-PTEN-PI3K/Akt signaling to evade host antiviral defense. These results demonstrate that the two forms of BmPGRP2 have different functions in host responses to bacteria and viruses.