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

Exploratory comparative transcriptomic analysis reveals potential gene targets associated with Cry1A.105 and Cry2Ab2 resistance in fall armyworm (Spodoptera frugiperda)

Genetically modified (GM) crops, expressing Bacillus thuringiensis (Bt) insecticidal toxins, have substantially transformed agriculture. Despite rapid adoption, their environmental and economic benefits face scrutiny due to unsustainable agricultural practices and the emergence of resistant pests like Spodoptera frugiperda, known as the fall armyworm (FAW). FAW's adaptation to Bt technology in corn and cotton compromises the long-term efficacy of Bt crops. To advance the understanding of the genetic foundations of resistance mechanisms, we conducted an exploratory comparative transcriptomic analysis of two divergent FAW populations. One population exhibited practical resistance to the Bt insecticidal proteins Cry1A.105 and Cry2Ab2, expressed in the genetically engineered MON-89Ø34 - 3 maize, while the other population remained susceptible to these proteins. Differential expression analysis supported that Cry1A.105 and Cry2Ab2 significantly affect the FAW physiology. A total of 247 and 254 differentially expressed genes were identified in the Cry-resistant and susceptible populations, respectively. By integrating our findings with established literature and databases, we underscored 53 gene targets potentially involved in FAW's resistance to Cry1A.105 and Cry2Ab2. In particular, we considered and discussed the potential roles of the differentially expressed genes encoding ABC transporters, G protein-coupled receptors, the P450 enzymatic system, and other Bt-related detoxification genes. Based on these findings, we emphasize the importance of exploratory transcriptomic analyses to uncover potential gene targets involved with Bt insecticidal proteins resistance, and to support the advantages of GM crops in the face of emerging challenges.

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

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

Comparative Transcriptome Analysis of Henosepilachna vigintioctomaculata Reveals Critical Pathways during Development

Henosepilachna vigintioctomaculata is distributed in several Asian countries. The larvae and adults often cause substantial economic losses to Solanaceae crops such as potato, tomato, eggplant, and Chinese boxthorn. Even though a chromosome-level genome has been documented, the expression profiles of genes involved in development are not determined. In this study, we constructed embryonic, larval, pupal, and adult transcriptomes, generated a comprehensive RNA-sequencing dataset including ~52 Gb of clean data, and identified 602,773,686 cleaned reads and 33,269 unigenes. A total of 18,192 unigenes were successfully annotated against NCBI nonredundant protein sequences, Swissprot, Eukaryotic Orthologous Groups, Gene Ontology (GO), or Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. There were 3580, 2040, 5160, 2496, 3008, and 3895 differentially expressed genes (DEGs) between adult/egg, egg/larval, larval/pupal, adult/pupal, egg/pupal, and adult/larval samples, respectively. GO and KEGG analyses of the DEGs highlighted several critical pathways associated with specific developing stages. This is the first comprehensive transcriptomic dataset encompassing all developmental stages in H. vigintioctomaculata. Our data may facilitate the exploitation of gene targets for pest control and can serve as a valuable gene resource for future molecular investigations.

Comparative transcriptome analysis of false codling moth, Thaumatotibia leucotreta in response to high and low-temperature treatments

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

Differential Antioxidant Enzyme Gene Expression and Functional Analysis of Pyridaben-Susceptible and -Resistant Strains of Tetranychus truncatus (Acari: Tetranychidae) under High Temperature Stress

Tetranychus truncatus (Acari: Tetranychidae) has caused serious economic losses on some crops (soybean, corn, and cotton) in China, and has developed resistance to most acaricides. Our laboratory study found that T. truncatus was resistant to pyridaben and also adapted to high temperature (34-40 °C). High temperature stress may cause arthropods to produce a large amount of reactive oxygen species (ROS), causing oxidative damage. Antioxidant enzymes, as the main antioxidants, can reduce the damage caused by excessive ROS in arthropods. In order to study the adaptation mechanism of the pyridaben-resistant strain of T. truncatus to high temperature and the role of antioxidant enzyme genes under high temperature stress, four antioxidant enzyme genes, TtSOD, TtPOD3, TtPOD4, and TtGSTs2, were screened according to the transcriptome sequencing data of pyridaben-susceptible and -resistant strains in T. truncatus. Firstly, the phylogeny and structure analyses of these four genes were carried out. Then, real-time quantitative PCR (RT-qPCR) technology was used to analyze the gene expression patterns of antioxidant enzymes in two strains of T. truncatus at three different high temperature ranges (34 °C, 38 °C, and 42 °C). The results showed that the expression levels of four antioxidant enzyme genes of two strains of T. truncatus were induced by high temperature stress, and the expression levels of antioxidant enzyme genes were significantly different in each development state. The gene expression of antioxidant enzyme genes in resistant strains at the adult stage was significantly higher than that in susceptible strains. After the TtSOD and TtPOD4 genes of adult mites of the resistant strain were silenced by RNA interference (RNAi) technology, the mortality rate of mites with TtPOD4 gene silencing reached 41.11% after 96 h at 34 °C, which was significantly higher than that of the control and TtSOD gene silencing. It has been confirmed that the TtPOD4 gene plays a key role in the adaptation of pyridaben-resistant strain of T. truncatus to high temperature. It lays a theoretical foundation for revealing the thermal adaptation mechanism of T. truncatus.

Functional Characterization and Putative Regulatory Mechanism of an RNAi Efficiency-Related Nuclease (REase) in the Fall Armyworm, Spodoptera frugiperda

The lepidopteran-specific RNAi efficiency-related nuclease (REase) has been shown to contribute to double-strand RNA (dsRNA) degradation in several lepidopteran insects. However, little is known about its regulatory mechanism. In this study, we identified and characterized SfREase in Spodoptera frugiperda. The exposure of the third-instar larvae to dsEGFP and high temperature led to the upregulation of SfREase, whereas starvation treatment resulted in the downregulation of SfREase. Further experiments revealed that dsRNA degraded more slowly in the hemolymph or midgut fluid extracted from dsSfREase-injected or dsSfREase-ingested larvae compared with those from dsEGFP-treated larvae, and the recombinant SfREase degraded dsRNA in a concentration-dependent manner. Additionally, the knockdown of SfREase improved RNAi efficiency. Finally, both RNAi and dual-luciferase reporter assay in Sf9 cells revealed that SfREase is negatively regulated by FOXO. These data provide insights into the function and regulatory mechanism of REase and have applied implications for the development of an RNAi-based control strategy of S. frugiperda.

Acquired preferences for a novel food odor do not become stronger or stable after multiple generations of odor feeding in Bicyclus anynana butterfly larvae

Many herbivorous insects have specific host-plant preferences, and it is unclear how these preferences evolved. Previously, we found that Bicyclus anynana larvae can learn to prefer novel food odors from eating leaves with those odors and transmit those learned preferences to the next generation. It is uncertain whether such acquired odor preferences can increase across generations of repeated odor feeding and be maintained even in the absence of odor. In this study, we fed larvae with novel banana odor-coated leaves (odor-fed larvae) for five consecutive generations, without selection on behavioral choices, and measured how larval innate preferences changed over time. Then, we removed the odor stimulus from a larval subgroup, while the other group continued to be odor-fed. Our results show that larvae learned to prefer the novel odor within a generation of odor feeding and transmitted the learned preference to the next generation, as previously found. Odor-fed larvae preferred odor significantly more compared to control larvae across five generations of repeated odor or control feeding. However, this led neither to increased odor preference, nor its stabilization. This suggests that when butterfly larvae feed on a new host, a preference for that novel food plant may develop and be transmitted to the next generation, but this preference lasts for a single generation and disappears once the odor stimulus is removed.

Blue Light Exposure Caused Large-Scale Transcriptional Changes in the Abdomen and Reduced the Reproductive Fitness of the Fall Armyworm Spodoptera frugiperda

In the present study, we found that blue light stress negatively affected the development periods, body weight, survival and reproduction of Spodoptera frugiperda, and it showed a dose-dependent reaction, as longer irradiation caused severer effects. Further transcriptome analysis found blue light stress induced fast and large-scale transcriptional changes in the head, thorax and, particularly, the abdomen of female S. frugiperda adults. A functional enrichment analysis indicated that shorter durations of blue light irradiation induced the upregulation of more stress response- and defense-related genes or pathways, such as abiotic stimuli detection and response, oxidative stress, ion channels and protein-kinase-based signal pathways. In the abdomen, however, different durations of blue-light-exposure treatments all induced the downregulation of a large number genes and pathways related to cellular processes, metabolism, catalysis and reproduction, which may be a trade-off between antistress defense and other processes or a strategy to escape stressful conditions. These results indicate irradiation duration- and tissue-specific blue light stress responses and consequences, as well as suggest that the stress that results in transcriptional alterations is associated with the stress that causes a fitness reduction in S. frugiperda females.

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

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

The role of cytochrome P450 3A2 and 4V2 in response to high-temperature stress in Tetranychus truncatus (Acari: Tetranychidae)

Temperature is an important factor influencing the physiological activities of agricultural pests. Therefore, understanding pest physiological activities and the molecular response to high-temperature stress is of paramount importance for pest management. Tetranychus truncatus Ehara (Acari: Tetranychidae) is a harmful organism, that may cause serious harm to crops such as corn and cotton in high-temperature environments. Cytochrome P450 (CYP450) is induced by high-temperature stress, and it plays an important role in the resistance of spider mites to high temperatures. Because of their role in high-temperature tolerance, the cytochrome P450 genes of the spider mite are attracting more and more attention. In this study, we identified and analyzed CYP450 genes in T. truncatus to investigate their potential roles in growth and development and the resistance to high-temperature stress. Based on phylogenetic and structural analyses, we identified 17 CYP450 genes in T. truncatus. RNA-seq and real-time quantitative PCR (RT-qPCR) revealed differential expression patterns of these genes at different developmental stages and levels of high-temperature stress resistance. The RNA interference results of selected CYP450 showed that when TtCYP3A2 and TtCYP4V2 were silenced by feeding on dsRNA, respectively, the high-temperature resistance of T. truncatus was decreased, which indicated that the expression of these two CYP450 genes in this species may be related to high-temperature tolerance. Our results provide potential evidence for the response of spider mites to high-temperature stress and help to improve the understanding of T. truncatus's ability to resist high-temperature stress.

Museum Genomics of an Agricultural Super-Pest, the Colorado Potato Beetle, Leptinotarsa decemlineata (Chrysomelidae), Provides Evidence of Adaptation from Standing Variation

Despite extensive research on agricultural pests, our knowledge about their evolutionary history is often limited. A mechanistic understanding of the demographic changes and modes of adaptation remains an important goal, as it improves our understanding of organismal responses to environmental change and our ability to sustainably manage pest populations. Emerging genomic datasets now allow for characterization of demographic and adaptive processes, but face limits when they are drawn from contemporary samples, especially in the context of strong demographic change, repeated selection, or adaptation involving modest shifts in allele frequency at many loci. Temporal sampling, however, can improve our ability to reconstruct evolutionary events. Here, we leverage museum samples to examine whether population genomic diversity and structure has changed over time, and to identify genomic regions that appear to be under selection. We focus on the Colorado potato beetle (CPB), Leptinotarsa decemlineata (Say 1824; Coleoptera: Chrysomelidae), which is widely regarded as a super-pest due to its rapid, and repeated, evolution to insecticides. By combining whole genome resequencing data from 78 museum samples with modern sampling, we demonstrate that CPB expanded rapidly in the 19th century, leading to a reduction in diversity and limited genetic structure from the Midwest to Northeast United States. Temporal genome scans provide extensive evidence for selection acting in resistant field populations in Wisconsin and New York, including numerous known insecticide resistance genes. We also validate these results by showing that known selective sweeps in modern populations are identified by our genome scan. Perhaps most importantly, temporal analysis indicates selection on standing genetic variation, as we find evidence for parallel evolution in the two geographical regions. Parallel evolution involves a range of phenotypic traits not previously identified as under selection in CPB, such as reproductive and morphological functional pathways that might be important for adaptation to agricultural habitats.

Integrated Biological Control Using a Mixture of Two Entomopathogenic Bacteria, Bacillus thuringiensis and Xenorhabdus hominickii, against Spodoptera exigua and Other Congeners

Insect immunity defends against the virulence of various entomopathogens, including Bacillus thuringiensis (Bt). This study tested a hypothesis that any suppression of immune responses enhances Bt virulence. In a previous study, the entomopathogenic bacterium, Xenorhabdus hominickii (Xh), was shown to produce secondary metabolites to suppress insect immune responses. Indeed, the addition of Xh culture broth (XhE) significantly enhanced the insecticidal activity of Bt against S. exigua. To analyze the virulence enhanced by the addition of Xh metabolites, four bacterial secondary metabolites were individually added to the Bt treatment. Each metabolite significantly enhanced the Bt insecticidal activity, along with significant suppression of the induced immune responses. A bacterial mixture was prepared by adding freeze-dried XhE to Bt spores, and the optimal mixture ratio to kill the insects was determined. The formulated bacterial mixture was applied to S. exigua larvae infesting Welsh onions in a greenhouse and showed enhanced control efficacy compared to Bt alone. The bacterial mixture was also effective in controlling other Spodopteran species such as S. litura and S. frugiperda but not other insect genera or orders. This suggests that Bt+XhE can effectively control Spodoptera-associated pests by suppressing the immune defenses.

Suppression of PBAN receptor expression reduces fecundity in the fall armyworm, Spodoptera frugiperda

The fall armyworm, Spodoptera frugiperda, native to the tropical and subtropical areas of the American continent is one of the world's most destructive insect pests. In most insects, sex pheromone production is initiated following the activation of a pheromone-biosynthesis-activating neuropeptide (PBAN) receptor, which belongs to G protein-coupled receptor. We explored expression level of S. frugiperda PBAN receptor (Sf-PBANr) gene and validated the physiological function by assessing the fecundity of adult females subjected to its specific RNA interference (RNAi). Sf-PBANr was predicted from a transcriptome of S. frugiperda. Reverse-transcription polymerase chain reaction assay showed its expression in all developmental stages of S. frugiperda. Specific suppression of Sf-PBANr by RNAi in either sex significantly reduced the total number of laid eggs per adult female. Matings between both RNAi-treated males and female resulted in 63.3% reduction in fecundity. In contrast, the RNAi effect was less 47.5%-49.5% at the matings from single-parent RNAi treatment. These results suggest that the Sf-PBANr is associated with female of S. frugiperda.

Satellitome Analysis and Transposable Elements Comparison in Geographically Distant Populations of Spodoptera frugiperda

Spodoptera frugiperda (fall armyworm) is a member of the superfamily Noctuoidea that accounts for more than a third of all Lepidoptera and includes a considerable number of agricultural and forest pest species. Spodoptera frugiperda is a polyphagous species that is a significant agricultural pest worldwide, emphasizing its economic importance. Spodoptera frugiperda’s genome size, assembly, phylogenetic classification, and transcriptome analysis have all been previously described. However, the different studies reported different compositions of repeated DNA sequences that occupied the whole assembled genome, and the Spodoptera frugiperda genome also lacks the comprehensive study of dynamic satellite DNA. We conducted a comparative analysis of repetitive DNA across geographically distant populations of Spodoptera frugiperda, particularly satellite DNA, using publicly accessible raw genome data from eight different geographical regions. Our results showed that most transposable elements (TEs) were commonly shared across all geographically distant samples, except for the Maverick and PIF/Harbinger elements, which have divergent repeat copies. The TEs age analysis revealed that most TEs families consist of young copies 1–15 million years old; however, PIF/Harbinger has some older/degenerated copies of 30–35 million years old. A total of seven satellite DNA families were discovered, accounting for approximately 0.65% of the entire genome of the Spodoptera frugiperda fall armyworm. The repeat profiling analysis of satellite DNA families revealed differential read depth coverage or copy numbers. The satellite DNA families range in size from the lowest 108 bp SfrSat06-108 families to the largest (1824 bp) SfrSat07-1824 family. We did not observe a statistically significant correlation between monomer length and K2P divergence, copy number, or abundance of each satellite family. Our findings suggest that the satellite DNA families identified in Spodoptera frugiperda account for a considerable proportion of the genome’s repetitive fraction. The satellite DNA families’ repeat profiling revealed a point mutation along the reference sequences. Limited TEs differentiation exists among geographically distant populations of Spodoptera frugiperda.