Comparative analysis of the transcriptional responses to low and high temperatures in three rice planthopper species

Chromosome-level genome of Osmia excavata (Hymenoptera: Megachilidae) provides insights into low-temperature tolerance of Osmia pollinator

The solitary bee Osmia excavata (Hymenoptera: Megachilidae) is a key pollinator managed on a large scale. It has been widely used for commercial pollination of fruit trees, vegetables, and other crops with high efficiency in increasing the crop seeding rate, yield, and seed quality in Northern hemisphere. Here, a high-quality chromosome-level genome of O. excavata was generated using PacBio sequencing along with Hi-C technology. The genome size was 207.02 Mb, of which 90.25% of assembled sequences were anchored to 16 chromosomes with a contig N50 of 9,485 kb. Approximately 186.83 Mb, accounting for 27.93% of the genome, was identified as repeat sequences. The genome comprises 12,259 protein-coding genes, 96.24% of which were functionally annotated. Comparative genomics analysis suggested that the common ancestor of O. excavata and Osmia bicornis (Hymenoptera: Megachilidae) lived 8.54 million years ago. Furthermore, cytochrome P450 family might be involved in the responses of O. excavata to low-temperature stress. Taken together, the chromosome-level genome assembly of O. excavata provides in-depth knowledge and will be a helpful resource for the pollination biology research.

Chemosensory protein 22 in Riptortus pedestris is involved in the recognition of three soybean volatiles

Riptortus pedestris (Hemiptera: Alydidae), a common agricultural pest, is the major causative agent of "soybean staygreen." However, the interactions between chemosensory proteins (CSPs) in R. pedestris and host plant volatiles have yet to be comprehensively studied. In this study, we performed real-time fluorescence quantitative polymerase chain reaction (PCR) to analyze the antennal expression of RpedCSP22 and subsequently analyzed the interactions between 21 soybean volatiles, five aggregation pheromones, and RpedCSP22 protein in vitro using a protein expression system, molecular docking, site-directed mutagenesis, and fluorescence competitive binding experiments. The RpedCSP22 protein showed binding affinity to three soybean volatiles (benzaldehyde, 4-ethylbenzaldehyde, and 1-octene-3-ol), with optimal binding observed under neutral pH conditions, and lost binding ability after site-directed mutagenesis. In subsequent RNA interference (RNAi) studies, gene silencing was more than 90 %, and in silenced insects, electroantennographic responses were reduced by more than 75 % compared to non-silenced insects. Moreover, Y-tube olfactory behavioral assessments revealed that the attraction of R. pedestris to the three soybean volatiles was significantly attenuated. These findings suggest that RpedCSP22 plays an important role in the recognition of host plant volatiles by R. pedestris andprovides a theoretical basis for the development of novel inhibitors targeting pest behavior.

Silencing NlFAR7 destroyed the pore canals and related structures of the brown planthopper

Fatty acyl-CoA reductase (FAR) is one of the key enzymes, which catalyses the conversion of fatty acyl-CoA to the corresponding alcohols. Among the FAR family members in the brown planthopper (Nilaparvata lugens), NlFAR7 plays a pivotal role in both the synthesis of cuticular hydrocarbons and the waterproofing of the cuticle. However, the precise mechanism by which NlFAR7 influences the formation of the cuticle structure in N. lugens remains unclear. Therefore, this paper aims to investigate the impact of NlFAR7 through RNA interference, transmission electron microscope, focused ion beam scanning electron microscopy (FIB-SEM) and lipidomics analysis. FIB-SEM is employed to reconstruct the three-dimensional (3D) architecture of the pore canals and related cuticle structures in N. lugens subjected to dsNlFAR7 and dsGFP treatments, enabling a comprehensive assessment of changes in the cuticle structures. The results reveal a reduction in the thickness of the cuticle and disruptions in the spiral structure of pore canals, accompanied by widened base and middle diameters. Furthermore, the lipidomics comparison analysis between dsNlFAR7- and dsGFP-treated N. lugens demonstrated that there were 25 metabolites involved in cuticular lipid layer synthesis, including 7 triacylglycerols (TGs), 5 phosphatidylcholines (PCs), 3 phosphatidylethanolamines (PEs) and 2 diacylglycerols (DGs) decreased, and 4 triacylglycerols (TGs) and 4 PEs increased. In conclusion, silencing NlFAR7 disrupts the synthesis of overall lipids and destroys the cuticular pore canals and related structures, thereby disrupting the secretion of cuticular lipids, thus affecting the cuticular waterproofing of N. lugens. These findings give significant attention with reference to further biochemical researches on the substrate specificity of FAR protein, and the molecular regulation mechanisms during N. lugens life cycle.

A Novel Approach to Comparative RNA-Seq Does Not Support a Conserved Set of Orthologs Underlying Animal Regeneration

Molecular studies of animal regeneration typically focus on conserved genes and signaling pathways that underlie morphogenesis. To date, a holistic analysis of gene expression across animals has not been attempted, as it presents a suite of problems related to differences in experimental design and gene homology. By combining orthology analyses with a novel statistical method for testing gene enrichment across large data sets, we are able to test whether tissue regeneration across animals shares transcriptional regulation. We applied this method to a meta-analysis of six publicly available RNA-Seq data sets from diverse examples of animal regeneration. We recovered 160 conserved orthologous gene clusters, which are enriched in structural genes as opposed to those regulating morphogenesis. A breakdown of gene presence/absence provides limited support for the conservation of pathways typically implicated in regeneration, such as Wnt signaling and cell pluripotency pathways. Such pathways are only conserved if we permit large amounts of paralog switching through evolution. Overall, our analysis does not support the hypothesis that a shared set of ancestral genes underlie regeneration mechanisms in animals. After applying the same method to heat shock studies and getting similar results, we raise broader questions about the ability of comparative RNA-Seq to reveal conserved gene pathways across deep evolutionary relationships.

Constitutively expressed small heat shock protein LsHsp21.5 not only enhances heat tolerance but also helps to maintain reproduction in female Laodelphax striatellus

Coping with stressful conditions and maintaining reproduction are two key biological processes that affect insect population dynamics. Small heat shock proteins (sHSPs) are involved in the stress response and the development of insects. The sHsp gene Laodelphax striatellus (Hemiptera: Delphacidae) sHsp 21.5 (LsHsp21.5) showed constitutive, stage- and organ-specific expression in L. striatellus, a pest that damages cultivated rice in east Asia. The expression of LsHsp21.5 was highest in the ovary, with 43.60, 12.99 and 1.45 time higher expression here than in the head, gut and female fat bodies, respectively. The expression of this gene was weakly affected by heat or cold shock. The gene provided in vitro protection against heat damage to malate dehydrogenase and in vivo protection against heat stress in Escherichia coli (Enterobacteriales: Enterobacteriaceae) BL21(DE3) and L. striatellus. Moreover, L. striatellus reproduction decreased by 1.85 times when the expression of LsHsp21.5 was inhibited by RNA interference. The expression of some genes related to reproduction, such as the homologous gene of chorion protein, also declined. These results suggest that LsHsp21.5 expression not only protects other proteins against stress but also helps maintain the stable expression of some reproduction-related genes under non-stressful conditions, with impacts on L. striatellus fecundity.

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

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

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

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

IMPORTANCE

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

A plant virus manipulates the long-winged morph of insect vectors

Wing dimorphism of insect vectors is a determining factor for viral long-distance dispersal and large-area epidemics. Although plant viruses affect the wing plasticity of insect vectors, the potential underlying molecular mechanisms have seldom been investigated. Here, we found that a planthopper-vectored rice virus, rice stripe virus (RSV), specifically induces a long-winged morph in male insects. The analysis of field populations demonstrated that the long-winged ratios of male insects are closely associated with RSV infection regardless of viral titers. A planthopper-specific and testis-highly expressed gene, Encounter, was fortuitously found to play a key role in the RSV-induced long-winged morph. Encounter resembles malate dehydrogenase in the sequence, but it does not have corresponding enzymatic activity. Encounter is upregulated to affect male wing dimorphism at early larval stages. Encounter is closely connected with the insulin/insulin-like growth factor signaling pathway as a downstream factor of Akt, of which the transcriptional level is activated in response to RSV infection, resulting in the elevated expression of Encounter. In addition, an RSV-derived small interfering RNA directly targets Encounter to enhance its expression. Our study reveals an unreported mechanism underlying the direct regulation by a plant virus of wing dimorphism in its insect vectors, providing the potential way for interrupting viral dispersal.

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.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Genomes of the cosmopolitan fruit pest Bactrocera dorsalis (Diptera: Tephritidae) reveal its global invasion history and thermal adaptation

INTRODUCTION

The oriental fruit fly Bactrocera dorsalis is one of the most destructive agricultural pests worldwide, with highly debated species delimitation, origin, and global spread routes.

OBJECTIVES

Our study intended to (i) resolve the taxonomic uncertainties between B. dorsalis and B. carambolae, (ii) reveal the population structure and global invasion routes of B. dorsalis across Asia, Africa, and Oceania, and (iii) identify genomic regions that are responsible for the thermal adaptation of B. dorsalis.

METHODS

Based on a high-quality chromosome-level reference genome assembly, we explored the population relationship using a genome-scale single nucleotide polymorphism dataset generated from the resequencing data of 487 B. dorsalis genomes and 25 B. carambolae genomes. Genome-wide association studies and silencing using RNA interference were used to identify and verify the candidate genes associated with extreme thermal stress.

RESULTS

We showed that B. dorsalis originates from the Southern India region with three independent invasion and spread routes worldwide: (i) from Northern India to Northern Southeast Asia, then to Southern Southeast Asia; (ii) from Northern India to Northern Southeast Asian, then to China and Hawaii; and (iii) from Southern India toward the African mainland, then to Madagascar, which is mainly facilitated by human activities including trade and immigration. Twenty-seven genes were identified by a genome-wide association study to be associated with 11 temperature bioclimatic variables. The Cyp6a9 gene may enhance the thermal adaptation of B. dorsalis and thus boost its invasion, which tended to be upregulated at a hardening temperature of 38 °C. Functional verification using RNA interference silencing against Cyp6a9, led to the specific decrease in Cyp6a9 expression, reducing the survival rate of dsRNA-feeding larvae exposed to extreme thermal stress of 45 °C after heat hardening treatments in B. dorsalis.

CONCLUSION

This study provides insights into the evolutionary history and genetic basis of temperature adaptation in B. dorsalis.

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.

Silencing a dehydration-responsive element-binding gene enhances the resistance of plants to a phloem-feeding herbivore

Herbivore-induced plant defence responses share common components with plant responses to abiotic stresses. However, whether abiotic stress-responsive factors influence the resistance of plants to herbivores by regulating these components remains largely unknown. Here, we cloned a dehydration-responsive element-binding gene in rice, OsDREB1A, and investigated its role in the resistance of rice to the phloem-feeding herbivore, brown planthopper (BPH, Nilaparvata lugens), under normal and low temperatures. We found that OsDREB1A localized to the nucleus, and its transcripts in rice were up-regulated in response to BPH infestation, low temperatures and treatment with methyl jasmonate or salicylic acid. Silencing OsDREB1A changed transcript levels of two defence-related WRKY and two PLD genes, enhanced levels of jasmonic acid (JA), JA-isoleucine and abscisic acid, and decreased the ethylene level in rice; these changes subsequently enhanced the resistance of plants to BPH, especially at 17°C, by decreasing the hatching rate and delaying the development of BPH eggs. Moreover, silencing OsDREB1A increased the growth of rice plants. These findings suggest that OsDREB1A, which positively regulates the resistance of rice to abiotic stresses, negatively regulates the resistance of rice to BPH.

Genomic adaptive potential to cold environments in the invasive red swamp crayfish

Biological invasion refers to the introduction, spread, and establishment of non-native species in a novel habitat. The ways in which invasive species successfully colonize new and different environments remain a fundamental topic of research in ecology and evolutionary biology. Here, we investigated the genomic and transcriptomic characteristics of the red swamp crayfish (Procambarus clarkii), a widespread invader in freshwater environments. Targeting a recently colonized population in Sapporo, Japan that appears to have acquired a high degree of cold tolerance, RNA-seq analysis revealed differentially expressed genes in response to cold exposure, and those involved in protease inhibitors and cuticle development were considered top candidates. We also found remarkable duplications for these gene families during evolution and their concerted expression patterns, suggesting functional amplification against low temperatures. Our study thus provides clues to the unique genetic characteristics of P. clarkii, possibly related to cold adaptation.

Heat shock protein 70 and Cathepsin B genes are involved in the thermal tolerance of Aphis gossypii

BACKGROUND

Elevated temperature can directly affect the insect pest population dynamics. Many experimental studies have indicated that high temperatures affect the biological and ecological characteristics of the widely distributed crop pest Aphis gossypii, but the molecular mechanisms underlying its response to heat stress remain unstudied. Here, we used transcriptomic analysis to explore the key genes and metabolic pathways involved in the regulation of thermotolerance in A. gossypii at 29 °C, 32 °C, and 35 °C.

RESULTS

The results of bioinformatics analysis show that few genes were consistently differentially expressed among the higher temperature treatments compared to 29 °C, and a moderate temperature increase of 3 °C can elicit gene expression changes that help A. gossypii adapt to warmer temperatures. Based on KEGG pathway enrichment analysis, we found that genes encoding four heat shock protein 70 s (Hsp70s) and nine cathepsin B (CathB) proteins were significantly upregulated at 35 °C compared with 32 °C. Genes related to glutathione production were also highly enriched between 32 °C and 29 °C. Silencing of two Hsp70s (ApHsp70A1-1 and ApHsp68) and two CathBs (ApCathB01 and ApCathB02) with RNA interference using a nanocarrier-based transdermal dsRNA delivery system significantly increased sensitivity of A. gossypii to high temperatures.

CONCLUSION

A. gossypii is able to fine-tune its response across a range of temperatures, and Hsp70 and CathB genes are essential for adaption of A. gossypii to warmer temperatures. © 2023 Society of Chemical Industry.

Functional characterization of chemosensory proteins in response to artificial light treatment in Thalassodes immissaria

BACKGROUND

Chemosensory proteins (CSPs) play a vital role in the response to environmental stimuli in insects. However, the involvement of insect CSPs in the stress response to night-time environmental changes has not been examined.

RESULTS

In the current study, four TiCSP genes were first cloned from Thalassodes immissaria by transcriptome and RACE-PCR techniques. TiCSPs had typical characteristics of insect CSPs, including a highly conserved four-cysteine motif and olfactory-specific protein D (OS-D) or OS-D superfamily domains. TiCSP1-4 were clustered classified within different clades in a phylogenetic analysis and were differentially expressed at all developmental stages. Under night-time artificial light stress, the expression levels of TiCSP1 in males were significantly decreased at 24 h, and those of TiCSP2 were decreased in both adult sexes at 48 h. In a molecular docking analysis, TiCSPs showed relatively higher binding affinities with sex pheromone components than with host plant volatile molecules.

CONCLUSION

Taking the reduced expression levels of TiCSPs and binding affinities into account, TiCSP1 and TiCSP2 are involved in the stress response processes of T. immissaria under light treatment. Our study supplies basic data for the evaluation of the effects of light interference control technology - an emerging physical control measure on nontarget pests of lychee orchards. © 2023 Society of Chemical Industry.

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.

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

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

Comparative transcriptome analysis of Liriomyza trifolii (Burgess) and Liriomyza sativae (Blanchard) (Diptera: Agromyzidae) in response to rapid cold hardening

The ability of insets to react efficiently to fluctuation in temperature is crucial for them to survive in variable surroundings. Rapid cold hardening (RCH) is a process that increase cold tolerance in most insect species. The molecular mechanisms of RCH remain largely unknown, and whether it is associated with transcriptional changes is unclear. In this study, we compared the transcriptomes of Liriomyza trifolii and L. sativae exposed to RCH to investigate the transcript abundance due to RCH in both species. RNA-seq revealed 93,166 assembled unigenes, and 34,303 of these were annotated in the L. trifolii and L. sativae transcriptome libraries. After a 4-h treatment at 1°C (RCH) compared with control, 268 and 606 unigenes were differentially expressed in L. trifolii and L. sativae, respectively. When comparing pupae exposed to 2h cold shock directly with pupae went through 4h acclimation prior to 2h cold shock, 60 and 399 unigenes were differentially expressed in L trifolii and L sativae, respectively. Genes that were commonly expressed in both L. trifolii and L. sativae, included cytochrome P450, cuticular protein, glucose dehydrogenase, solute carrier family 22 and cationic amino acid transporter. Additionally, several pathways including galactose metabolism and peroxisome were significantly enriched during RCH. Our results show that the transcriptional response is correlated with RCH in the pupal stage of the two Liriomyza species, but more transcriptional changes were identified in L sativae than in L. trifolii.

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

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

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

The fall armyworm (FAW), Spodoptera frugiperda, is regarded as one of the world’s most harmful plant pests. This research examines the molecular response processes of FAW to low temperature (4°C) and high temperature (40°C) when gene expression is compared to controls (25°C). A total of 211,967 unigenes were collected, at least 14,338 of which were annotated with gene descriptions, gene ontology terms, and metabolic pathways. There were 50 Gene Ontology (GO) functional sub-groups and 21 EggNOG words as a result. Differentially expresses genes (DEGs) with log₂FC ≥ 2 were identified and compared at various temperatures. In comparison to the 25°C treated group, we discovered 199 and 1,248 individual DEGs co-regulated at 4 and 40°C, respectively. Comparing transcriptome profiles for differential gene expression revealed a number of DEGs, including cytochrome P450, odorant binding proteins (OBPs), and immune system genes previously implicated in cold and high temperature stresses. The enrichment pathways were identified using Kyoto Encyclopedia of Genes and Genomics (KEGG) analysis, and heatmaps of similar unigenes from both treatment groups (T4 and T40) were plotted. We used quantitative reverse transcription PCR (RT-qPCR) to confirm the RNA-seq data on 10 up- and down-regulated DEGs. These findings provide a foundation for future understanding of FAW adaptation mechanisms and the underlying basis underlying the response to low and high temperatures.

Large-scale genome-wide study reveals climate adaptive variability in a cosmopolitan pest

Understanding the genetic basis of climatic adaptation is essential for predicting species' responses to climate change. However, intraspecific variation of these responses arising from local adaptation remains ambiguous for most species. Here, we analyze genomic data from diamondback moth (Plutella xylostella) collected from 75 sites spanning six continents to reveal that climate-associated adaptive variation exhibits a roughly latitudinal pattern. By developing an eco-genetic index that combines genetic variation and physiological responses, we predict that most P. xylostella populations have high tolerance to projected future climates. Using genome editing, a key gene, PxCad, emerged from our analysis as functionally temperature responsive. Our results demonstrate that P. xylostella is largely capable of tolerating future climates in most of the world and will remain a global pest beyond 2050. This work improves our understanding of adaptive variation along environmental gradients, and advances pest forecasting by highlighting the genetic basis for local climate adaptation.

Genome-wide gene expression profiles of the pea aphid (Acyrthosiphon pisum) under cold temperatures provide insights into body color variation

Cold temperatures are one of the factors influencing color polymorphisms in Acyrthosiphon pisum, resulting in a change from a red to greenish color. Here we characterized gene expression profiles of A. pisum under different low temperatures (1°C, 4°C, 8°C, and 14°C) and durations (3, 6, 12, and 24 h). The number of differentially expressed genes (DEGs) increased as temperatures decreased and time increased, but only a small number of significant DEGs were identified. Genes involved in pigment metabolism were downregulated. An interaction network analysis for 506 common DEGs in comparisons among aphids exposed to 1°C for four durations indicated that a cytochrome P450 gene (CYP, LOC112935894) significantly downregulated may interact with a carotenoid metabolism gene (LOC100574964), similar to other genes encoding CYP, lycopene dehydrogenase and fatty acid synthase. We proposed that the body color shift in A. pisum responding to low temperatures may be regulated by CYPs.

Comparison of gene expression in the red imported fire ant (Solenopsis invicta) under different temperature conditions

The red imported fire ant (RIFA), Solenopsis invicta Buren is native to South America and is known as a global problematic invasive species. This study focused on the molecular response of RIFA by comparing gene expression profiles after exposing ants to low (10 °C) and high (40 °C) temperature stress and comparing them to untreated controls (30 °C). A total of 99,085 unigenes (the clustered non-redundant transcripts that are filtered from the longest assembled contigs) were obtained, of which 19,154 were annotated with gene descriptions, gene ontology terms, and metabolic pathways. 86 gene ontology (GO) functional sub-groups and 23 EggNOG terms resulted. Differentially expressed genes (DEGs) with log₂FC ≥ 10 were screened and were compared at different temperatures. We found 203, 48, and 66 specific DEGs co-regulated at 10, 20, and 40 °C. Comparing transcriptome profiles for differential gene expression resulted in various DE genes, including cytochrome P450, NADH dehydrogenase subunit 1, cuticle protein and heat shock protein (HSP), which have previously been reported to be involved in cold and high temperature resistance. GO analysis revealed that antioxidant activity is up-regulated under high temperature stress. We verified the RNA-seq data by qPCR on 20 up- and down-regulated DEGs. These findings provide a basis for future understanding of the adaptation mechanisms of RIFA and the molecular mechanisms underlying the response to low and high temperatures.

Cuticular Hydrocarbon Plasticity in Three Rice Planthopper Species

Insect cuticular hydrocarbons (CHCs) are organic compounds of the surface lipid layer, which function as a barrier against water loss and xenobiotic penetration, while also serving as chemical signals. Plasticity of CHC profiles can vary depending upon numerous biological and environmental factors. Here, we investigated potential sources of variation in CHC profiles of Nilaparvata lugens, Laodelphax striatellus and Sogatella furcifera, which are considered to be the most important rice pests in Asia. CHC profiles were quantified by GC/MS, and factors associated with variations were explored by conducting principal component analysis (PCA). Transcriptomes were further compared under different environmental conditions. The results demonstrated that CHC profiles differ among three species and change with different developmental stages, sexes, temperature, humidity and host plants. Genes involved in cuticular lipid biosynthesis pathways are modulated, which might explain why CHC profiles vary among species under different environments. Our study illustrates some biological and ecological variations in modifying CHC profiles, and the underlying molecular regulation mechanisms of the planthoppers in coping with changes of environmental conditions, which is of great importance for identifying potential vulnerabilities relating to pest ecology and developing novel pest management strategies.

Impact of Temperature on Survival Rate, Fecundity, and Feeding Behavior of Two Aphids, Aphis gossypii and Acyrthosiphon gossypii, When Reared on Cotton

Aphid performance is sensitive to temperature changes. Previous studies found that Acyrthosiphon gossypii (Mordviiko) was more sensitive to high temperature than Aphis gossypii (Glover). However, the effects of high temperatures on the survival, fecundity, and feeding behavior of these two aphid adults are not clear. This study examined the effect of different temperatures (29 °C, 32 °C, and 35 °C) on the adult survival rate, fecundity, and feeding behavior of these two aphid species. Our results showed that the adverse effects of high temperatures (32 °C and 35 °C) on aphid adult survival and fecundity were greater for Ac. gossypii than Ap. gossypii. The electrical penetration graph (EPG) data showed that Ac. gossypii spent more time feeding on xylem than phloem under all temperature treatments, which contrasted with Ap. gossypii. The time of phloem ingestion by Ap. gossypii at 32 °C was significantly higher than at 29 °C, while for Ac. gossypii, this value significantly decreased when temperature increased. These feeding patterns indicate that Ac. gossypii obtains less nutrition from phloem in support of its development and fecundity. Data generated in this study will serve as the basis for predicting the effects of increased temperature on these two cotton aphids.

Transcriptome responses to heat and cold stress in prepupae of Trichogramma chilonis

Trichogramma is a useful species that is widely applied in biocontrol. Temperature profoundly affects the commercial application of T. chilonis. Different developmental transcriptomes of prepupae and pupae of T. chilonis under 10, 25, and 40°C were obtained from our previous study. In this study, transcriptomic analysis was further conducted to gain a clear understanding of the molecular changes in the prepupae of T. chilonis under different thermal conditions. A total of 37,295 unigenes were identified from 3 libraries of prepupae of T. chilonis, 17,293 of which were annotated. Differential expression analysis showed that 408 and 108 differentially expressed genes (DEGs) were identified after heat and cold treatment, respectively. Under heat stress, the pathway of protein processing in endoplasmic reticulum was found to be active. Most of the genes involved in this pathway were annotated as lethal (2) essential for life [l(2)efl] and heat shock protein genes (hsps), which were both highly upregulated. Nevertheless, most of the genes involved in another significantly enriched pathway of starch and sucrose metabolism were downregulated, including 1 alpha-glucosidase gene and 2 beta-glucuronidase genes. Under cold stress, no significantly enriched pathway was found, and the significantly enriched GO terms were related to the interaction with host and immune defenses. Together, these results provide us with a comprehensive view of the molecular mechanisms of T. chilonis in response to temperature stresses and will provide new insight into the mass rearing and utilization of T. chilonis.

Comparative transcriptome analysis of the newly discovered insect vector of the pine wood nematode in China, revealing putative genes related to host plant adaptation

BACKGROUND

In many insect species, the larvae/nymphs are unable to disperse far from the oviposition site selected by adults. The Sakhalin pine sawyer Monochamus saltuarius (Gebler) is the newly discovered insect vector of the pine wood nematode (Bursaphelenchus xylophilus) in China. Adult M. saltuarius prefers to oviposit on the host plant Pinus koraiensis, rather than P. tabuliformis. However, the genetic basis of adaptation of the larvae of M. saltuarius with weaken dispersal ability to host environments selected by the adult is not well understood.

RESULTS

In this study, the free amino and fatty acid composition and content of the host plants of M. saltuarius larvae, i.e., P. koraiensis and P. tabuliformis were investigated. Compared with P. koraiensis, P. tabuliformis had a substantially higher content of various free amino acids, while the opposite trend was detected for fatty acid content. The transcriptional profiles of larval populations feeding on P. koraiensis and P. tabuliformis were compared using PacBio Sequel II sequencing combined with Illumina sequencing. The results showed that genes relating to digestion, fatty acid synthesis, detoxification, oxidation-reduction, and stress response, as well as nutrients and energy sensing ability, were differentially expressed, possibly reflecting adaptive changes of M. saltuarius in response to different host diets. Additionally, genes coding for cuticle structure were differentially expressed, indicating that cuticle may be a potential target for plant defense. Differential regulation of genes related to the antibacterial and immune response were also observed, suggesting that larvae of M. saltuarius may have evolved adaptations to cope with bacterial challenges in their host environments.

CONCLUSIONS

The present study provides comprehensive transcriptome resource of M. saltuarius relating to host plant adaptation. Results from this study help to illustrate the fundamental relationship between transcriptional plasticity and adaptation mechanisms of insect herbivores to host plants.

HSP70/DNAJ Family of Genes in the Brown Planthopper, Nilaparvata lugens: Diversity and Function

Heat shock 70kDa proteins (HSP70s) and their cochaperones DNAJs are ubiquitous molecular chaperones, which function as the “HSP70/DNAJ machinery” in a myriad of biological processes. At present, a number of HSP70s have been classified in many species, but studies on DNAJs, especially in insects, are lacking. Here, we first systematically identified and characterized the HSP70 and DNAJ family members in the brown planthopper (BPH), Nilaparvata lugens, a destructive rice pest in Asia. A total of nine HSP70 and 31 DNAJ genes were identified in the BPH genome. Sequence and phylogenetic analyses revealed the high diversity of the NlDNAJ family. Additionally, spatio-temporal expression analysis showed that most NlHSP70 and NlDNAJ genes were highly expressed in the adult stage and gonads. Furthermore, RNA interference (RNAi) revealed that seven NlHSP70s and 10 NlDNAJs play indispensable roles in the nymphal development, oogenesis, and female fertility of N. lugens under physiological growth conditions; in addition, one HSP70 (NlHSP68) was found to be important in the thermal tolerance of eggs. Together, our results in this study shed more light on the biological roles of HSP70/DNAJ in regulating life cycle, coping with environmental stresses, and mediating the interactions within, or between, the two gene families in insects.

Chromosomal-level genomes of three rice planthoppers provide new insights into sex chromosome evolution

The brown planthopper Nilaparvata lugens, white-backed planthopper Sogatella furcifera, and small brown planthopper Laodelphax striatellus are three major insect pests of rice. They are genetically close; however, they differ in several ecological traits such as host range, migration capacity, and in their sex chromosomes. Though the draft genome of these three planthoppers have been previously released, the quality of genome assemblies need to be improved. The absence of chromosome-level genome resources has hindered in-depth research of these three species. Here, we performed a de novo genome assembly for N. lugens to increase its genome assembly quality with PacBio and Illumina platforms, increasing the contig N50 to 589.46 Kb. Then, with the new N. lugens genome and previously reported S. furcifera and L. striatellus genome assemblies, we generated chromosome-level scaffold assemblies of these three planthopper species using HiC scaffolding technique. The scaffold N50s significantly increased to 77.63 Mb, 43.36 Mb and 29.24 Mb for N. lugens, S. furcifera and L. striatellus, respectively. To identify sex chromosomes of these three planthopper species, we carried out genome re-sequencing of males and females and successfully determined the X and Y chromosomes for N. lugens, and X chromosome for S. furcifera and L. striatellus. The gene content of the sex chromosomes showed high diversity among these three planthoppers suggesting the rapid evolution of sex-linked genes, and all chromosomes showed high synteny. The chromosome-level genome assemblies of three planthoppers would provide a valuable resource for a broad range of future research in molecular ecology, and subsequently benefits development of modern pest control strategies.

Comparative transcriptome analysis of three invasive leafminer flies provides insights into interspecific competition

Liriomyza spp. (Diptera: Agromyzidae) represent a group of economically-significant highly polyphagous pests of plants grown in field and greenhouse conditions. Liriomyza spp. share similar biological and morphological characteristics, and complex interspecific interactions have been documented among these species in various geographical regions. Where the displacement of one of these species by the other has been studied, no unique mechanisms have been identified as causing it. The impact of competitive factors (such as, insecticide tolerance, thermotolerance, and adaptability to cropping systems) may be unique to specific geographic regions of Liriomyza spp., but more research is needed to confirm these hypotheses. In this study, RNA-seq was used to determine the transcriptomes of three closely-related leafminers, e.g. L. sativae, L. trifolii, and L. huidobrensis. Over 20 Gb of clean reads were generated and assembled into unique transcriptomes, and 38,747 unigenes were annotated in different databases. In pairwise comparisons, L. trifolii and L. sativae had more up-regulated genes than L. huidobrensis. With respect to common differentially-expressed genes (Co-DEGs), the three leafminers exhibited distinct groups of highly-expressed gene clusters. When genes related to competitive factors were compared, expression patterns in L. trifolii and L. sativae were more closely related to each other than to L. huidobrensis. The data suggest that DEGs involved in competitive factors may play a key role in competition and displacement of leafminers. The divergent genes identified in this study will be valuable in revealing possible mechanisms of invasion, displacement and interspecific competition in Liriomyza spp.

Echolocation call frequency variation in horseshoe bats: molecular basis revealed by comparative transcriptomics

Recently diverged taxa with contrasting phenotypes offer opportunities for unravelling the genetic basis of phenotypic variation in nature. Horseshoe bats are a speciose group that exhibit a derived form of high-duty cycle echolocation in which the inner ear is finely tuned to echoes of the narrowband call frequency. Here, by focusing on three recently diverged subspecies of the intermediate horseshoe bat (Rhinolophus affinis) that display divergent echolocation call frequencies, we aim to identify candidate loci putatively involved in hearing frequency variation. We used de novo transcriptome sequencing of two mainland taxa (himalayanus and macrurus) and one island taxon (hainanus) to compare expression profiles of thousands of genes. By comparing taxa with divergent call frequencies (around 15 kHz difference), we identified 252 differentially expressed genes, of which six have been shown to be involved in hearing or deafness in human/mouse. To obtain further validation of these results, we applied quantitative reverse transcription-PCR to the candidate gene FBXL15 and found a broad association between the level of expression and call frequency across taxa. The genes identified here represent strong candidate loci associated with hearing frequency variation in bats.

A novel negative thermotaxis behavior in rice planthoppers is regulated by TRPA1 channel

BACKGROUND

In the field, we observed that many white-backed planthoppers (Sogatella furcifera, WBPHs) stayed in the top region of rice plants exposed to direct sunshine. It was known that WBPHs frequently took flight when the ground temperature was about 25 °C, then climbed to and flew in a dense layer corresponding to an air temperature of about 16 °C in the sky. Its migration height was higher than the top of the surface temperature inversion. It is still unclear whether WBPHs prefer warm or cold regions, and therefore we studied the thermal responses of WBPHs and other insects using a simulated system.

RESULTS

We found that WBPHs preferred a cold region to a warm one, unexpectedly below their comfort temperature zone. After comparative analysis with other insect species, such as small brown planthoppers, brown planthoppers, Trialeurodes vaporariorum (stinkbugs, a predator of planthoppers) and Bemisia tabaci (whitefly), only three planthoppers showed cold preference behavior. RNA interference experiments revealed that this behavior of WBPHs can be regulated by the transient receptor potential (TRP) channels TRPA1 channel. Furthermore, podocarpic acid, an agonist of TRPA1, weakened the cold preference, whereas A-967079, an antagonist of TRPA1, had the opposite effect.

CONCLUSION

We reported a novel cold preference (negative thermotaxis) in rice planthoppers that was regulated in WBPHs by the TRPA1 channel. Cold preference of rice planthoppers is probably related to its choice behavior of the special migratory temperature layer. Our results expanded a new perspective to develop novel strategies for behavioral manipulation and management of rice planthoppers. © 2020 Society of Chemical Industry.

Comparative analysis of diet-associated responses in two rice planthopper species

Background

Host adaptation is the primary determinant of insect diversification. However, knowledge of different host ranges in closely related species remains scarce. The brown planthopper (Nilaparvata lugens, BPH) and the small brown planthopper (Laodelphax striatellus, SBPH) are the most destructive insect pests within the family Delphacidae. These two species differ in their host range (SBPH can well colonize rice and wheat plants, whereas BPH survives on only rice plants), but the underlying mechanism of this difference remains unknown. High-throughput sequencing provides a powerful approach for analyzing the association between changes in gene expression and the physiological responses of insects. Therefore, gut transcriptomes were performed to elucidate the genes associated with host adaptation in planthoppers. The comparative analysis of planthopper responses to different diets will improve our knowledge of host adaptation regarding herbivorous insects.

Results

In the present study, we analyzed the change in gene expression of SBPHs that were transferred from rice plants to wheat plants over the short term (rSBPH vs tSBPH) or were colonized on wheat plants over the long term (rSBPH vs wSBPH). The results showed that the majority of differentially expressed genes in SBPH showed similar changes in expression for short-term transfer and long-term colonization. Based on a comparative analysis of BPH and SBPH after transfer, the genes associated with sugar transporters and heat-shock proteins showed similar variation. However, most of the genes were differentially regulated between the two species. The detoxification-related genes were upregulated in SBPH after transfer from the rice plants to the wheat plants, but these genes were downregulated in BPH under the same conditions. In contrast, ribosomal-related genes were downregulated in SBPH after transfer, but these genes were upregulated in BPH under the same conditions.

Conclusions

The results of this study provide evidence that host plants played a dominant role in shaping gene expression and that the low fitness of BPH on wheat plants might be determined within 24 h after transfer. This study deepens our understanding of different host ranges for the two planthopper species, which may provide a potential strategy for pest management.

Ten fatty acyl-CoA reductase family genes were essential for the survival of the destructive rice pest, Nilaparvata lugens

BACKGROUND

Fatty alcohols are the precursors of sex pheromone components, wax esters and hydrocarbons in insects. Fatty acyl-CoA reductases (FARs) are important enzymes required for the reduction of fatty alcohol and thereby contribute to the production of cuticular hydrocarbon (CHC).

RESULTS

Based on bioinformatics analyses we identified 17 FAR genes in the brown planthopper, Nilaparvata lugens. RNA interference against these genes demonstrated that ten NlFAR genes were essential for the survival of N. lugens. For instance, knockdown of NlFAR5, 6, 11 or 15 was lethal and caused a slender body shape, while the old cuticles of the respective animals remained attached to the abdomen or failed to split open from the nota. Knockdown of NlFAR9 resulted in a phenotype, with a smooth body surface and a decrease in CHC amounts. Similarly, CHC deficiency in N. lugens resulted in increased adhesion of water droplets and secreted honeydew to the insect surface and the inability of N. lugens to survive in paddy fields with varying humidity. Knockdown of NlFAR1, 4, 5, 6, 8, 9, 11 and 13 additionally resulted in female adult infertility.

CONCLUSION

The present study illustrates the structural and functional differences of FAR family genes and provides potential targets for RNA interference-based rice planthopper management. © 2020 Society of Chemical Industry.

Temperature-dependent oviposition and nymph performance reveal distinct thermal niches of coexisting planthoppers with similar thresholds for development

The brown planthopper (Nilapavata lugens: BPH) and whitebacked planthopper (Sogatella furcifera: WBPH) co-occur as the principal pests of rice in Asia. A review of previous studies suggests that the two species have similar temperature tolerances and similar temperature thresholds for development. However, the distribution and seasonality of WBPH suggest that its temperature optima for performance (survival, oviposition and growth) may be lower than for BPH. We compared adult longevity, oviposition, nymph survival and development success, as well as nymph biomass in both species across a gradient of constant temperatures from 15°C-40°C, at 5°C intervals. The most suitable temperatures for oviposition, nymph biomass and development success were 5-10°C lower for WBPH than for BPH. Furthermore, compared to BPH, WBPH demonstrated clear differences in oviposition on different rice subspecies and on rice at different growth stages at 25°C and 30°C, but not at other temperatures. The results suggest that aspects of herbivore performance within tolerable temperature ranges, which are not often included in temperature models, may be more useful than thermal tolerances or development thresholds in predicting the effects of global warming on pest damage to crops.

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.

Comparative transcriptome analysis of differentially expressed genes in Bradysia odoriphaga Yang et Zhang (Diptera: Sciaridae) at different acute stress temperatures

The gnat, Bradysia odoriphaga Yang et Zhang, is an important underground pest in Asia. B. odoriphaga differ in heat and cold tolerance and exhibit quite different developmental strategies. To understand the underlying mechanisms, we sequenced and compared the transcriptome of B. odoriphaga under 40 °C (a stressful high temperature), 25 °C, and 4 °C (a stressful low temperature) for 1 h. We found that metabolism- and ribosome-related genes were modulated. In high temperature (40 °C), heat shock protein (HSP) genes, detoxication genes, metabolism genes, protein turnover genes, and stress signal transduction genes were differentially expressed. In low temperature (4 °C), genes related with heat shock protein (HSP) and detoxication were differentially expressed. Our study increases our understanding of the complex molecular mechanisms involved in the responses of B. odoriphaga to acute temperature stress and provides a potential strategy for pest management.

Effect of glycogen synthase and glycogen phosphorylase knockdown on the expression of glycogen- and insulin-related genes in the rice brown planthopper Nilaparvata lugens

Nilaparvata lugens is a serious threat to rice growth. Glycogen metabolism is one of the important physiological processes of insects, which is mainly regulated by glycogen synthase (GS) and glycogen phosphorylase (GP). In the present study, trehalose content was significantly reduced at 72 h after NlGP and NlGS knockdown, whereas glucose content was significantly increased at both 48 h and 72 h after GS knockdown. RNAi combined with RNA-Seq was used to identify NlGP- and NlGS-related pathways and genes in N. lugens. A total of 593 genes were up-regulated and 5969 genes were down-regulated after NlGP and NlGS knockdown, respectively. Moreover, the NlGS-knockdown group was mapped to 10,967 pathways, whereas the NlGP-knockdown group was mapped to 7948 pathways, and the greatest differences between the groups were associated with carbohydrate, lipid, amino acid and energy metabolism. Meanwhile, 1800, 1217, and 1211 transcripts in the NlGP-knockdown group and 2511, 1666, and 1727 transcripts in the NlGS-knockdown group were involved in bioprocess, cellular ingredients and molecular function, respectively. Almost all these genes were down-regulated by either NlGP or NlGS knockdown, with significant down-regulation of the 6-trehalose phosphate synthase (TPS), trehalase (TRE), GS, GP, phosphoacetylglucosamine mutase (PGM, n = 2), Insulin receptors (InRs) and insulin-like peptides (Ilps) genes. These results have demonstrated that RNAi-mediated NlGP and NlGS knockdown could lead to content of trehalose and glucose out of balance, but have no obvious effect on glycogen content, and have suggested that GS plays more complex role in other metabolism pathway of N. lugens.

Effects of abiotic stress on the expression of Hsp70 genes in Sogatella furcifera (Horváth)

Sogatella furcifera (Horváth), a prominent rice pest in Asia, is a typical R-strategic and highly adaptable insect. Heat shock proteins (Hsps) are highly conserved molecular chaperones regulating responses to various abiotic stresses; however, limited information is available regarding their role in responding to abiotic stress in S. furcifera. This study aimed to investigate the effect of abiotic stresses on the expression of Hsp70 genes in the S. furcifera. Five Hsp70 genes were isolated from S. furcifera, and the expression patterns at different developmental stages and temperatures, upon treatment with different insecticides and ultraviolet A (UV-A) stress, were analyzed. Hsp70 genes were expressed at different developmental stages. Hsp70-2, Hsp70-5, and Hsp70-6 were significantly upregulated upon heat shock at 40 °C for 30 min. Hsp70-3 and Hsp70-4 were significantly upregulated upon heat shock at 30 °C for 30 min. Under UV-A stress, Hsp70-3, Hsp70-4, Hsp70-5, and Hsp70-6 were significantly upregulated. Conversely, Hsp70-2 was significantly downregulated under UV-A stress. The five Hsp70 genes were significantly downregulated in 3rd-instar nymphs on exposure to thiamethoxam, buprofezin, and avermectin at LC10 and LC25 concentrations. Hence, Hsp70 genes significantly contribute to the tolerance of S. furcifera to temperature and UV-A stress; however, they are not involved in the response to insecticides.

Comparison of upper sublethal and lethal temperatures in three species of rice planthoppers

Temperature is an important environmental factor for ectotherms' fitness and survival. The upper sublethal and lethal temperatures were compared between adults of three closely related destructive planthopper species, the small brown planthopper (Laodelphax striatellus, SBPH), the brown planthopper (Nilaparvata lugens, BPH), and the white-backed planthopper (Sogatella furcifera, WBPH) in the absence and presence of the host plant (Oryza sativa, var. Taichong1). Values of the critical thermal maxima (CT_max) were higher in SBPH than in both BPH and WBPH and higher in BPH than in WBPH, and values of the heat coma temperatures (HCT) were higher in both BPH and SBPH than in WBPH. CT_max and HCT values were higher in the presence than in the absence of plant material. Between sexes, females generally showed higher CT_max and HCT than males. The upper lethal temperatures (ULT₅₀) measured in the absence of plant material were not significantly different among the planthopper species. The planthoppers also exhibited different behaviors in an increasing temperature regime, with fewer insects dropping-off from the plant in SBPH than in BPH and WBPH. These results indicate that SBPH and BPH are more heat tolerant than WBPH. The findings highlight the biological divergence in closely related planthopper species and the importance of performing the heat tolerance measurement in an ecologically relevant setting, which serves to predict seasonal and spatial occurrence patterns of the destructive planthopper species.

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

Galeruca daurica (Joannis) has become a new insect pest in the Inner Mongolia grasslands since 2009, and its larvae and eggs have strong cold tolerance. To get a deeper insight into its molecular mechanisms of cold stress responses, we performed de novo transcriptome assembly for G. daurica by RNA-Seq and compared the transcriptomes of its larvae exposed to five different temperature treatments (-10, -5, 0, 5, and 25°C for 1 h and then recovered at 25°C for 1 h), respectively. Compared with the control (25°C), the numbers of differentially expressed genes (DEGs) decreased from 1,821 to 882, with the temperature declining from 5 to -10°C. Moreover, we obtained 323 coregulated DEGs under different low temperatures. Under four low temperatures (-10, -5, 0, and 5°C), a large number of genes were commonly upregulated during recovery from cold stresses, including those related to cuticle protein, followed by cytochrome P450, clock protein, fatty acid synthase, and fatty acyl-CoA reductase; meanwhile, lots of genes encoding cuticle protein, RNA replication protein, RNA-directed DNA polymerase, and glucose dehydrogenase were commonly downregulated. Our findings provide important clues for further investigations of key genes and molecular mechanisms involved in the adaptation of G. daurica to harsh environments.

Comparison between transcriptomic responses to short-term stress exposures of a common Holarctic and endemic Lake Baikal amphipods

Background

Lake Baikal is one of the oldest freshwater lakes and has constituted a stable environment for millions of years, in stark contrast to small, transient bodies of water in its immediate vicinity. A highly diverse endemic endemic amphipod fauna is found in one, but not the other habitat. We ask here whether differences in stress response can explain the immiscibility barrier between Lake Baikal and non-Baikal faunas. To this end, we conducted exposure experiments to increased temperature and the toxic heavy metal cadmium as stressors.

Results

Here we obtained high-quality de novo transcriptome assemblies, covering mutiple conditions, of three amphipod species, and compared their transcriptomic stress responses. Two of these species, Eulimnogammarus verrucosus and E. cyaneus, are endemic to Lake Baikal, while the Holarctic Gammarus lacustris is a potential invader.

Conclusions

Both Baikal species possess intact stress response systems and respond to elevated temperature with relatively similar changes in their expression profiles. G. lacustris reacts less strongly to the same stressors, possibly because its transcriptome is already perturbed by acclimation conditions.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6024-3) contains supplementary material, which is available to authorized users.

Mitochondrial variation in small brown planthoppers linked to multiple traits and probably reflecting a complex evolutionary trajectory

While it has been proposed in several taxa that the mitochondrial genome is associated with adaptive evolution to different climatic conditions, making links between mitochondrial haplotypes and organismal phenotypes remains a challenge. Mitonuclear discordance occurs in the small brown planthopper (SBPH), Laodelphax striatellus, with one mitochondrial haplogroup (HGI) more common in the cold climate region of China relative to another form (HGII) despite strong nuclear gene flow, providing a promising model to investigate climatic adaptation of mitochondrial genomes. We hypothesized that cold adaptation through HGI may be involved, and considered mitogenome evolution, population genetic analyses, and bioassays to test this hypothesis. In contrast to our hypothesis, chill-coma recovery tests and population genetic tests of selection both pointed to HGII being involved in cold adaptation. Phylogenetic analyses revealed that HGII is nested within HGI, and has three nonsynonymous changes in ND2, ND5 and CYTB in comparison to HGI. These molecular changes likely increased mtDNA copy number, cold tolerance and fecundity of SBPH, particularly through a function-altering amino acid change involving M114T in ND2. Nuclear background also influenced fecundity and chill recovery (i.e., mitonuclear epistasis) and protein modelling indicates possible nuclear interactions for the two nonsynonymous changes in ND2 and CYTB. The high occurrence frequency of HGI in the cold climate region of China remains unexplained, but several possible reasons are discussed. Overall, our study points to a link between mtDNA variation and organismal-level evolution and suggests a possible role of mitonuclear interactions in maintaining mtDNA diversity.

Transcriptome-Based Analysis Reveals a Crucial Role of BxGPCR17454 in Low Temperature Response of Pine Wood Nematode (Bursaphelenchus xylophilus)

Background: The causal agent of pine wilt disease is the pine wood nematode (PWN) (Bursaphelenchus xylophilus), whose ability to adapt different ecological niches is a crucial determinant of their invasion to colder regions. To discover the molecular mechanism of low temperature response mechanism, we attempted to study the molecular response patterns under low temperature from B. xylophilus with a comprehensive RNA sequencing analysis and validated the differentially expressed genes (DEGs) with quantitative real-time polymerase chain reaction (qRT-PCR). Bioinformatic software was utilized to isolate and identify the low-temperature-related BxGPCR genes. Transcript abundance of six low-temperature-related BxGPCR genes and function of one of the BxGPCR genes are studied by qRT-PCR and RNA interference. Results: The results showed that we detected 432 DEGs through RNA sequencing between low-temperature-treated and ambient-temperature-treated groups nematodes. The transcript level of 6 low-temperature-related BxGPCR genes increased at low temperature. And, the survival rates of BxGPCR17454 silenced B. xylophilus revealed a significant decrease at low temperature. Conclusion: in conclusion, this transcriptome-based study revealed a crucial role of BxGPCR17454 in low temperature response process of pine wood nematode. These discoveries would assist the development of management and methods for efficient control of this devastating pine tree pest.

A transcriptional and functional analysis of heat hardening in two invasive fruit fly species, Bactrocera dorsalis and Bactrocera correcta

Many insects have the capacity to increase their resistance to high temperatures by undergoing heat hardening at nonlethal temperatures. Although this response is well established, its molecular underpinnings have only been investigated in a few species where it seems to relate at least partly to the expression of heat shock protein (Hsp) genes. Here, we studied the mechanism of hardening and associated transcription responses in larvae of two invasive fruit fly species in China, Bactrocera dorsalis and Bactrocera correcta. Both species showed hardening which increased resistance to 45°C, although the more widespread B. dorsalis hardened better at higher temperatures compared to B. correcta which hardened better at lower temperatures. Transcriptional analyses highlighted expression changes in a number of genes representing different biochemical pathways, but these changes and pathways were inconsistent between the two species. Overall B. dorsalis showed expression changes in more genes than B. correcta. Hsp genes tended to be upregulated at a hardening temperature of 38°C in both species, while at 35°C many Hsp genes tended to be upregulated in B. correcta but not B. dorsalis. One candidate gene (the small heat shock protein gene, Hsp23) with a particularly high level of upregulation was investigated functionally using RNA interference (RNAi). We found that RNAi may be more efficient in B. dorsalis, in which suppression of the expression of this gene removed the hardening response, whereas in B. correcta RNAi did not decrease the hardening response. The different patterns of gene expression in these two species at the two hardening temperatures highlight the diverse mechanisms underlying hardening even in closely related species. These results may provide target genes for future control efforts against such pest species.

The influence of temperature on the toxicity of insecticides to Nilaparvata lugens (Stål)

The toxicity of insecticides is associated with a variety of factors including temperature, and global warming is bound to lead to the outbreak of pests; therefore, it is important to study the influence of temperature on insecticide toxicity and pest control. In this study, the influence of temperature on the toxicity of insecticides to Nilaparvata lugens (BPH) was determined. The results showed that the sensitivity of BPH to cycloxaprid (LC₅₀ = 42.5-0.388 mg/L), nitenpyram (LC₅₀ = 3.49-0.187 mg/L), triflumezopyrim (LC₅₀ = 0.354-0.0533 mg/L) and chlorpyrifos (LC₅₀ = 36.3-7.41 mg/L) increased significantly when the temperature changed from 18 °C to 36 °C. BPH sensitivity to etofenprox (LC₅₀ = 9.04-54.2 mg/L) was also affected by temperature. Additionally, the feeding amount and the activities of three detoxification enzymes [cytochrome P450 (P450), glutathione S-transferase (GST) and carboxylesterase (CarE)] of BPH at different temperatures were also measured. The feeding amounts were positively correlated with temperature increases while the activities of P450 and GST were significantly inhibited. The correlation analysis showed that changes in P450 activity (but not GST activity) were closely related to the sensitivity of BPH to cycloxaprid, nitenpyram, chlorpyrifos, and etofenprox according to the variation in temperatures. This study provides a theoretical basis for the rational use of chemical pesticides under the global warming trend and provides a reference for the integrated management of BPH in the field.

The fatty acid elongase gene family in the brown planthopper, Nilaparvata lugens

The cuticular hydrocarbon (CHC) biosynthetic pathways branches off from the synthesis of fatty acids. Fatty acid elongases (ELOs) are enzymes catalyzing the synthesis of long-chain fatty acids and thereby contribute to the diversification of CHCs. Based on bioinformatics analyses we identified 20 ELO genes in the brown planthopper, Nilaparvata lugens. RNA interference against these genes demonstrated that 9 NlELO genes were essential for the survival of N. lugens nymphs and adults. Indeed, knockdown of NlELOs 1, 3, 4, 7, 8, 9, 10, 12 and 18 caused lethal phenotypes with a thin and wizened body and reduced lipids in the fat body. Surface analysis by scanning electron microscopy and CHC quantification indicated that knockdown of NlELOs 2, 3, 8 and 16 additionally resulted in a smooth body surface and a decrease in CHC amounts. Therefore, we speculate that long-chain CHCs are needed for CHC attachment to the cuticle surface. CHC deficiency, in turn, resulted in increased adhesion of water droplets and secreted honeydew to the animal surface and the inability of N. lugens to survive in paddy fields with varying humidity. Our present study provides an initial comprehensive analysis of ELO gene functions in an insect, and may serve to better understand the biology of CHCs.

Comparative transcriptome analysis reveals differentially expressed genes in the Asian citrus psyllid (Diaphorina citri) upon heat shock

Temperature is an important factor influencing insect distribution. In the tropical and subtropical regions, insects always suffer the extreme temperature. Therefore, appropriate molecular response to temperature change is crucial for their survival. To understand how Diaphorina citri responds to high temperature at the molecular level, we conducted a comparative analysis of the transcriptomes of D. citri under room temperature and 40 °C heat shock treatment. The RNA sequencing analysis identified a total of 451 differentially expressed genes upon heat stress, including 167 down-regulated genes and 284 up-regulated genes. Thermal stress mainly significantly induced the transcription of molecular chaperone, protein biosynthesis and oxidation resistance, including members of the heat shock protein families, ATPases, and detoxifying enzymes. This study provides a preliminary insight into the transcriptional response of D. citri to heat stress and provides a foundation for the future functional validation of genes involved in thermotolerance in this important insect pest in tropical and subtropical regions.