Transcriptome response to temperature stress in the wolf spider Pardosa pseudoannulata (Araneae: Lycosidae)

Impact of high temperature and drought stress on the microbial community in wolf spiders

High temperatures and drought present significant abiotic challenges that can limit the survival of many arthropods, including wolf spiders, which are ectothermic and play a crucial role in controlling pest populations. However, the impact of these stress factors on the microbiota of spiders remains poorly understood. In this study, we utilized 16 S rRNA gene sequencing to explore the diversity and composition of bacterial communities within Pardosa pseudoannulata under conditions of high temperature and drought stress. We found that Firmicutes, Bacteroidetes, and Proteobacteria were the predominant bacterial phyla present. Analyses of alpha diversity indicated an increase in bacterial diversity under combined stress conditions, as reflected by various diversity indices such as Ace, Chao1, Shannon, and Simpson. Furthermore, co-occurrence network analysis highlighted intricate interactions among the microbial taxa (e.g., Enterobacter, Chitinophaga, and Eubacterium), revealing the adaptive complexity of the spider's microbiome to environmental stress. Functional prediction analysis suggested that combined stress conditions might enhance key metabolic pathways, particularly those related to oxidative phosphorylation and amino acid metabolism. Using Random Forest analysis, we determined that changes in three heat shock proteins were largely attributed to variations in bacterial communities, with Firmicutes being notably influential. Collectively, this in-depth analysis offers novel insights into the responses of microbial communities within spider microbiomes to combined abiotic stresses, providing valuable information for understanding extreme climate impacts and informing ecological management strategies.

Transcriptomics of differences in thermal plasticity associated with selection for an exaggerated male sexual trait

The information about the magnitude of differences in thermal plasticity both between and within populations, as well as identification of the underlying molecular mechanisms are key to understanding the evolution of thermal plasticity. In particular, genes underlying variation in the physiological response to temperature can provide raw material for selection acting on plastic traits. Using RNAseq, we investigate the transcriptional response to temperature in males and females from bulb mite populations selected for the increased frequency of one of two discrete male morphs (fighter- and scrambler-selected populations) that differ in relative fitness depending on temperature. We show that different mechanisms underlie the divergence in thermal response between fighter- and scrambler-selected populations at decreased vs. increased temperature. Temperature decrease to 18 °C was associated with higher transcriptomic plasticity of males with more elaborate armaments, as indicated by a significant selection-by-temperature interaction effect on the expression of 40 genes, 38 of which were upregulated in fighter-selected populations in response to temperature decrease. In response to 28 °C, no selection-by-temperature interaction in gene expression was detected. Hence, differences in phenotypic response to temperature increase likely depended on genes associated with their distinct morph-specific thermal tolerance. Selection of males also drove gene expression patterns in females. These patterns could be associated with temperature-dependent fitness differences between females from fighter- vs. scrambler-selected populations reported in previous studies. Our study shows that selection for divergent male sexually selected morphologies and behaviors has a potential to drive divergence in metabolic pathways underlying plastic response to temperature in both sexes.

Immunostimulation of Parasteatoda tepidariorum (Araneae: Theridiidae) in juvenile and adult stages. Immunity reactions to injury with foreign body and Bacillus subtilis infection

To assess the immune potential of spiders, in the present study juvenile and adult females of Parasteatoda tepidariorum were exposed to Bacillus subtilis infection, injury by a nylon monofilament and a combination of both. The expression level of selected immune-related genes: defensin 1 (PtDEF1), lysozyme 1 (PtLYS1), lysozyme C (PtLYSC), lysozyme M1 (PtLYSM1), autophagy-related protein 101 (PtATG101), dynamin (PtDYN) and heat shock proteins (HSP70) (PtHSPB, PtHSPB2A, PtHSPB2B), production of lysozyme and HSP70 proteins, and hemocytes viability were measured. The obtained results indicated expression of the lysozyme, autophagy-related protein and HSP70 genes in both ontogenetic stages of P. tepidariorum. It has been also shown that the simultaneous action of mechanical and biological factors causes higher level of lysozyme and HSP70, cell apoptosis intensity and lower level of hemocytes viability than in the case of exposure to a single immunostimulant. Moreover, mature females showed stronger early immune responses compared to juveniles.

The metabolism and detoxification effects of lead exposure on Aleurolyphus ovatus (Acari: Acaridae) via transcriptome analysis

Aleurolyphus ovatus Troupeau is one of the most predominant species of the Acaridae family worldwide. Recent reports have demonstrated that the accumulation of lead in stored grains and dietary items exceeds the required standards. However, the molecular mechanism of heavy metal stress on mites has not been reported. To understand the mechanism underlying the heavy metal response of A. ovatus, comparative transcriptome analysis was performed in this study using an Illumina high throughput mRNA sequencing (RNA-seq) platform. A. ovatus was fed on artificial diets containing two different concentrations of lead, namely, a low concentration of 12.5 mg/kg (LAO) and a high concentration of 100 mg/kg (HAO), while the mites in the control (NAO) group were not exposed to lead. A total of 44,362 unigenes, with an average length of 1547 bp, were identified. Of these, 996 unigenes were successfully annotated in seven functional databases. The number of differentially expressed genes (DEGs) in A. ovatus under different lead concentrations was compared. In NAO versus LAO group, including 310 up-regulated and 1580 down-regulated DEGs. In NAO versus HAO group, including 3928 up-regulated and 1761 down-regulated DEGs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment indicated that detoxification enzyme genes were significantly expressed in pathways, such as cytochrome P450 foreign body metabolism, glutathione metabolism and drug metabolism-cytochrome pathway. The results of gene annotation and quantitative real-time PCR showed that high concentration of lead significantly stimulated the expression of metabolic detoxification enzyme genes such as glutathione S transferase (GST) and superoxide dismutase (SOD), while low concentration inhibited their expression. This study will provide a basis for the molecular mechanism of A. ovatus in response to heavy metal lead stimulation in stored grain.

Integrated transcriptome and metabolome analysis reveals molecular responses of spider to single and combined high temperature and drought stress

High temperature and drought are abiotic stresses restricting many arthropods' survival and growth. Wolf spiders are poikilothermic arthropods that are vital in managing insects and pests. Nonetheless, investigating changes in spiders under temperature and drought stress are limited, especially at the molecular and gene expression levels. The study found that the combined effects of high temperature and drought stress significantly reduced survival rates and raised superoxide dismutase and malondialdehyde levels in the wolf spider Pardosa pseudoannulata. An integrated transcriptome and metabolome analysis showed that differentially expressed genes and metabolites were highly enriched in pathways involved in the proteolysis and oxidation-reduction process. The gene expression profiles displayed that heat shock protein (HSP) families (i.e., small heat shock protein, HSP70, HSP90, and HSP beta protein) were up-regulated under temperature and/or drought stresses. Additionally, a conjoint analysis revealed that under the combined stress, several important enzymes, including maltase-glucoamylase, glycerol-6-phosphate transporter, alanine-glyoxylate transaminase, and prostaglandin-H2 D-isomerase, were altered, affecting the metabolism of starch, sucrose, amino acids, and arachidonic acid. The protein interaction network further confirmed that under the combined stress, metabolic processes, peptide metabolic processes, and ATP generation from ADP were up-regulated, indicating that spiders could accelerate the metabolism of carbohydrates and proteins to combat stress and maintain homeostasis. Overall, this work showed that exposure to a combination of pressures might cause distinct defensive reactions in spiders and offered novel perspectives to research the molecular underpinnings of spider adaptation to a changing climate.

Temperature fluctuation alters optimal predator community composition for anticipated biological control

Alongside pesticides and specialist predators, natural communities of generalist beetle and spider predators play an important role in suppressing agricultural pests. However, the predation pressure of natural communities can be unpredictable. Overall predation pressure is influenced by a dense network of potential intraguild interactions, which are further shaped by species traits and environmental factors. Understanding how these different influences combine to impact pest control is especially important in the context of changing global temperatures. Recent empirical studies have demonstrated that the foraging behavior of arthropod predators is influenced by an interaction between temperature and predator body size. To explore the consequences of these findings for intraguild interactions and pest control, we expand a previously published model describing interactions between arthropod predators and a pest population. The model assumed that interaction strengths are influenced by body size and habitat preference. In our updated model, we incorporate the effect of temperature on predator foraging activity. We parameterize the model to match empirically observed predator community composition in 10 agricultural fields and use simulations to demonstrate how temperature-dependent behaviors change the expected efficiency of the natural predator community. Then, we use an optimization approach to identify the most efficient composition of natural predators for pest control. We then evaluate whether the most efficient predator compositions would change with increasing average daily temperature and its variability, as is expected under future temperature change. We find that optimal communities often include predators with complementary foraging activity and that in 2 fields, the optimal community changes drastically under future temperatures. We also note that at some temperatures, foraging activity reduces the negative effects of intraguild interference on pest consumption. This work allows us to assess the effect of climate change on the efficiency of natural predator communities to control pest populations and provide guidance for farmers to design pest management strategies tailored to different climate scenarios.

Effects of urea application on the reproduction of Pardosa pseudoannulata: Field and laboratory studies

As an important chemical fertilizer, urea can greatly increase crop yields, but it also has negative effects on beneficial arthropods in the agricultural field, such as spiders. Here, we reported that urea application reduced the reproductive performance in Pardosa pseudoannulata, a dominant species of spider in rice fields, which preys on a range of insect pests, based on both field and laboratory studies. In a field test, urea application significantly reduced the egg production of adult and subadult females collected from the urea-treated fields. A laboratory test was set up to further evaluate the impact of urea application on P. pseudoannulata reproduction. In consistent with field test results, the spiders treated by urea for 14 d and 28 d had lower reproduction ability than their control counterparts, with regard to the mating rate, egg production, and egg hatchability. The transcriptomic sequencing of individuals treated by urea for 28 d showed that urea application caused a number of differentially expressed transcripts with several downregulated unigenes related to basic enzymes and several upregulated unigenes involved in stress resistance. The knockdown of a metalloproteinase gene caused a significant decrease in egg production, and the silencing of a carboxylesterase gene significantly reduced both the egg production and egg hatchability. Taken together, the present study found that urea application reduced P. pseudoannulata reproduction ability and the negative impact partially resulted from the downregulation of certain basic enzyme genes. The study provided a fresh view of fertilizers on beneficial arthropods with great potential in the protection of P. pseudoannulata in fields.

Ecdysteroid responses to urban heat island conditions during development of the western black widow spider (Latrodectus hesperus)

The steroid hormone 20-hydroxyecdysone (20E) controls molting in arthropods. The timing of 20E production, and subsequent developmental transitions, is influenced by a variety of environmental factors including nutrition, photoperiod, and temperature, which is particularly relevant in the face of climate change. Environmental changes, combined with rapid urbanization, and the increasing prevalence of urban heat islands (UHI) have contributed to an overall decrease in biodiversity making it critical to understand how organisms respond to elevating global temperatures. Some arthropods, such as the Western black widow spider, Latrodectus hesperus, appear to thrive under UHI conditions, but the physiological mechanism underlying their success has not been explored. Here we examine the relationship between hemolymph 20E titers and spiderling development under non-urban desert (27°C), intermediate (30°C), and urban (33°C) temperatures. We found that a presumptive molt-inducing 20E peak observed in spiders at non-urban desert temperatures was reduced and delayed at higher temperatures. Intermolt 20E titers were also significantly altered in spiders reared under UHI temperatures. Despite the apparent success of black widows in urban environments, we noted that, coincident with the effects on 20E, there were numerous negative effects of elevated temperatures on spiderling development. The differential effects of temperature on pre-molt and intermolt 20E titers suggest distinct hormonal mechanisms underlying the physiological, developmental, and behavioral response to heat, allowing spiders to better cope with urban environments.

The effects of prey lipid on female mating and reproduction of a wolf spider

As predators, the macronutrients spiders extract from their prey play important roles in their mating and reproduction. Previous studies of macronutrients on spider mating and reproduction focus on protein, the potential impact of prey lipid content on spider mating and reproduction remains largely unexplored. Here, we tested the influence of prey varying in lipid content on female mating, sexual cannibalism, reproduction, and offspring fitness in the wolf spider Pardosa pseudoannulata. We acquired 2 groups of fruit fly Drosophila melanogaster that differed significantly in lipid but not protein content by supplementing cultural media with a high or low dose of sucrose on which the fruit flies were reared (HL: high lipid and LL: low lipid). Subadult (i.e., 1 molt before adult) female spiders that fed HL flies matured with significantly higher lipid content than those fed LL flies. We found that the mated females fed with HL flies significantly shortened pre-oviposition time and resulted in a significantly higher fecundity. However, there was no significant difference in female spiders varying in lipid content on other behaviors and traits, including the latency to courtship, courtship duration, mating, copulation duration, sexual cannibalism, offspring body size, and survival. Hence, our results suggest that the lipid content of prey may be a limiting factor for female reproduction, but not for other behavioral traits in the wolf spiders P. pseudoannulata.

Molecular response uncovers neurotoxicity of Pardosa pseudoannulata exposed to cadmium pressure

Cadmium (Cd) is a widely distributed heavy metal in south of China. Growing evidence indicates that systemic exposure to Cd, particularly the long-term exposure, may cause neurotoxic effects. Nevertheless, mechanisms underlying Cd neurotoxicity remain not completely understood. In this report, we investigated the neural alterations in the spider Pardosa pseudoannulata (Bösenberg and Strand, 1906) exposed to long-term Cd (LCd) and short-term Cd (SCd) pressure. Cd stress lowered foraging ability and prey consuming time in the spiders. In addition, enzymatic analysis results indicated that Cd exposure reduced the level of acetylcholinesterase at subcellular level. We then identified differentially expressed genes (DEGs) in the Cd exposed spiders using pairwise comparisons and found that a large number of DEGs were related to neurotransmitter receptors and ion transport and binding proteins. Notably, LCd exposure harbored more altered genes in ion transporter activity comparing with SCd exposure. From six K-means clusters, 53 putative transcriptional factors (TFs) belonging to 21 families were characterized, and ZBTB subfamily displayed the most distinctive alterations in the characterized genes, which is assumed to play a key role in the regulation of ion transmembrane process under Cd stress. A protein-to-protein interaction network constructed by the yielded DEGs also showed that ion and receptor binding activities were affected under long-term Cd exposure. Four key modules from the network indicated that Cd may further down-regulate energy metabolism pathway in spiders. Collectively, this comprehensive analysis provides multi-dimensional insights to understand the molecular response of spiders to Cd exposure.

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.

Transcriptome analysis of Curcuma wenyujin from Haikou and Wenzhou, and a comparison of the main constituents and related genes of Rhizoma Curcumae

For more than a thousand years, Rhizoma Curcumae (known as E zhu), a Chinese herbal medicine, has been used to eradicate blood stasis and relieve aches. The plant Curcuma wenyujin, which is grown primarily in Wenzhou, China, is considered the best source of Rhizoma Curcumae. In this study, we sought to ascertain differences in transcript profiles of C. wenyujin grown in traditional (Wenzhou) and recently established (Haikou) production areas based on Illumina and RNA (RNA-seq) sequencing. We also examined differences in the main components of the volatile oil terpene; curcumin, polysaccharide, and starch constituents and related genes in the corresponding pathways, in C. wenyujin cultivated in the two production areas. We accordingly found that the essential oil (2.05%), curcumin (1.46%), and polysaccharide (8.90%) content in Wenzhou rhizomes was higher than that in the rhizomes of plants from Haikou (1.60%, 0.91%, and 6.15%, respectively). In contrast, the starch content of Wenzhou rhizomes (17.0%) was lower than that of Haikou rhizomes (23.8%). Furthermore, we detected significant differences in the oil components of Haikou and Wenzhou rhizomes, with curzerene (32.34%), curdione (21.35%), and germacrene B (9.39%) being the primary components of the essential oil derived from Wenzhou rhizomes, and curzerene (20.13%), curdione (14.73%), and cineole (9.76%) being the main constituents in Haikou rhizomes. Transcriptome and qPCR analyses revealed considerable differences in gene expression between Wenzhou and Haikou rhizomes. The expression of terpene, curcumin, and polysaccharide pathway-related genes in Wenzhou rhizomes was significantly up-regulated, whereas the expression of starch-associated genes was significantly down-regulated, compared with those in Haikou rhizomes. Difference in the content of terpene, curcumin, polysaccharides, and starch in rhizomes from the two production areas could be explained in terms of differences in expression of the related genes.

Transcriptome responses to elevated CO2 level and Wolbachia-infection stress in Hylyphantes graminicola (Araneae: Linyphiidae)

Hylyphantes graminicola is a resident spider species found in maize and cotton fields and is an important biological control agent of various pests. Previous studies have demonstrated that stress from elevated CO₂ and Wolbachia infection can strongly affect spider species. Thus, based on CO₂ levels (400 ppm, current atmospheric CO₂ concentration and 800 ppm, high CO₂ concentration) and Wolbachia status (Wolbachia-infected, W⁺ and Wolbachia-uninfected, W^- ), we divided H. graminicola individuals into four treatment groups: W^- 400 ppm, W^- 800 ppm, W⁺ 400 ppm, and W⁺ 800 ppm. To investigate the effects of elevated CO₂ levels (W^- 400 vs W^- 800), Wolbachia infection (W^- 400 vs W⁺ 400), and the interactions between these two factors (W^- 400 vs W⁺ 800), high-throughput transcriptome sequencing was employed to characterize the de novo transcriptome of the spiders and identify stress-related differentially expressed genes (DEGs). De novo assembly of complementary DNA sequences generated 86 688 unigenes, 23 938 of which were annotated in public databases. A total of 84, 21, and 157 DEGs were found among W^- 400 vs W^- 800, W^- 400 vs W⁺ 400, and W^- 400 vs W⁺ 800, respectively. Functional enrichment analysis revealed that metabolic processes, signaling, and catalytic activity were significantly affected by elevated CO₂ levels and Wolbachia infection. Our findings suggest that the impact of elevated CO₂ levels and Wolbachia infection on the H. graminicola transcriptome was, to a large extent, on genes involved in metabolic processes. This study is the first description of transcriptome changes in response to elevated CO₂ levels and Wolbachia infection in spiders.

Genetic Constraints, Transcriptome Plasticity, and the Evolutionary Response to Climate Change

In situ adaptation to climate change will be critical for the persistence of many ectotherm species due to their relative lack of dispersal capacity. Climate change is causing increases in both the mean and the variance of environmental temperature, each of which may act as agents of selection on different traits. Importantly, these traits may not be heritable or have the capacity to evolve independently from one another. When genetic constraints prevent the "baseline" values of thermal performance traits from evolving rapidly, phenotypic plasticity driven by gene expression might become critical. We review the literature for evidence that thermal performance traits in ectotherms are heritable and have genetic architectures that permit their unconstrained evolution. Next, we examine the relationship between gene expression and both the magnitude and duration of thermal stress. Finally, we identify genes that are likely to be important for adaptation to a changing climate and determine whether they show patterns consistent with thermal adaptation. Although few studies have measured narrow-sense heritabilities of thermal performance traits, current evidence suggests that the end points of thermal reaction norms (tolerance limits) are moderately heritable and have the potential to evolve rapidly. By contrast, performance at intermediate temperatures has substantially lower evolutionary potential. Moreover, evolution in many species appears to be constrained by genetic correlations such that populations can adapt to either increases in mean temperature or temperature variability, but not both. Finally, many species have the capacity for plastic expression of the transcriptome in response to temperature shifts, with the number of differentially expressed genes increasing with the magnitude, but not the duration, of thermal stress. We use these observations to develop a conceptual model that describes the likely trajectory of genome evolution in response to changes in environmental temperature. Our results indicate that extreme weather events, rather than gradual increases in mean temperature, are more likely to drive genetic and phenotypic change in wild ectotherms.

Cooler temperatures slow the repair of DNA damage in tadpoles exposed to ultraviolet radiation: Implications for amphibian declines at high altitude

Ultraviolet B radiation (UVBR) damages the DNA of exposed cells, causing dimers to form between adjacent pyrimidine nucleotides. These dimers block DNA replication, causing mutations and apoptosis. Most organisms utilize biochemical or biophysical DNA repair strategies to restore DNA structure; however, as with most biological reactions, these processes are likely to be thermally sensitive. Tadpoles exposed to elevated UVBR at low environmental temperatures have significantly higher rates of mortality and developmental deformities compared with tadpoles exposed to the same levels of UVBR at higher environmental temperatures. We hypothesized that low environmental temperatures impair the primary enzymatic (photolyase) DNA repair pathway in amphibians, leading to the accumulation of DNA damage. To test this hypothesis, we compared DNA repair rates and photolyase gene expression patterns in Limnodynastes peronii. Tadpoles were acutely exposed to UVBR for 1 hr at either 20 or 30°C, and we measured DNA damage and photolyase expression levels at intervals following this exposure. Temperature had a significant effect on the rate of DNA repair, with repair at 30°C occurring twice as fast as repair at 20°C. Photolyase gene expression (6-4 PP and CPD) was significantly upregulated by UVBR exposure, with expression levels increasing within 6 hr of UVBR exposure. CPD expression levels were not significantly affected by temperature, but 6-4 PP expression was significantly higher in tadpoles in the 30°C treatment within 12 hr of UVBR exposure. These data support the hypothesis that DNA repair rates are thermally sensitive in tadpoles and may explain why enigmatic amphibian declines are higher in montane regions where UVBR levels are naturally elevated and environmental temperatures are lower.

Evolutionary genomics can improve prediction of species' responses to climate change

Global climate change (GCC) increasingly threatens biodiversity through the loss of species, and the transformation of entire ecosystems. Many species are challenged by the pace of GCC because they might not be able to respond fast enough to changing biotic and abiotic conditions. Species can respond either by shifting their range, or by persisting in their local habitat. If populations persist, they can tolerate climatic changes through phenotypic plasticity, or genetically adapt to changing conditions depending on their genetic variability and census population size to allow for de novo mutations. Otherwise, populations will experience demographic collapses and species may go extinct. Current approaches to predicting species responses to GCC begin to combine ecological and evolutionary information for species distribution modelling. Including an evolutionary dimension will substantially improve species distribution projections which have not accounted for key processes such as dispersal, adaptive genetic change, demography, or species interactions. However, eco-evolutionary models require new data and methods for the estimation of a species' adaptive potential, which have so far only been available for a small number of model species. To represent global biodiversity, we need to devise large-scale data collection strategies to define the ecology and evolutionary potential of a broad range of species, especially of keystone species of ecosystems. We also need standardized and replicable modelling approaches that integrate these new data to account for eco-evolutionary processes when predicting the impact of GCC on species' survival. Here, we discuss different genomic approaches that can be used to investigate and predict species responses to GCC. This can serve as guidance for researchers looking for the appropriate experimental setup for their particular system. We furthermore highlight future directions for moving forward in the field and allocating available resources more effectively, to implement mitigation measures before species go extinct and ecosystems lose important functions.

Comparative Transcriptome Analysis of Megacopta cribraria (Hemiptera: Plataspidae) in Response to High-Temperature Stress

Kudzu bug, Megacopta cribraria (Fabricius), is mainly distributed in southern China and has been considered an invasive species in the southeastern United States. Megacopta cribraria is a soybean pest with high-temperature resistance, but the molecular mechanisms underlying its thermal adaptation are largely unknown. Here, we performed comparative transcriptome analysis to unravel the molecular response of M. cribraria toward high-temperature stress. Following RNA-seq, we identified 93,959 assembled unigenes, 14,073 of which were annotated in M. cribraria transcriptome libraries. In addition, 127 differentially expressed unigenes (DEGs) were detected, 88 of them were significantly upregulated, whereas the remaining 39 genes were significantly downregulated. Functional classification revealed that the pathways of metabolic process, cellular processes, and single-organism processes were considered to be significantly enriched. In the COG classification, DEGs were mainly localized into O: post-translational modification, protein turnover, chaperone. Moreover, protein processing in endoplasmic reticulum and linoleic acid metabolism were significantly enriched among the 38 KEGG pathways. Further gene annotation analysis indicated that nine heat shock protein-related genes were significantly upregulated. Finally, five HSP DEGs were selected for real-time quantitative polymerase chain reaction validation and demonstrated a similar upregulation trend with RNA-seq expression profiles. Taken altogether, these findings provide new insights into the molecular mechanisms of thermal adaptation in M. cribraria.

Effects of adult temperature on gene expression in a butterfly: identifying pathways associated with thermal acclimation

BACKGROUND

Phenotypic plasticity is a pervasive property of all organisms and considered to be of key importance for dealing with environmental variation. Plastic responses to temperature, which is one of the most important ecological factors, have received much attention over recent decades. A recurrent pattern of temperature-induced adaptive plasticity includes increased heat tolerance after exposure to warmer temperatures and increased cold tolerance after exposure to cooler temperatures. However, the mechanisms underlying these plastic responses are hitherto not well understood. Therefore, we here investigate effects of adult acclimation on gene expression in the tropical butterfly Bicyclus anynana, using an RNAseq approach.

RESULTS

We show that several antioxidant markers (e.g. peroxidase, cytochrome P450) were up-regulated at a higher temperature compared with a lower adult temperature, which might play an important role in the acclamatory responses subsequently providing increased heat tolerance. Furthermore, several metabolic pathways were up-regulated at the higher temperature, likely reflecting increased metabolic rates. In contrast, we found no evidence for a decisive role of the heat shock response.

CONCLUSIONS

Although the important role of antioxidant defence mechanisms in alleviating detrimental effects of oxidative stress is firmly established, we speculate that its potentially important role in mediating heat tolerance and survival under stress has been underestimated thus far and thus deserves more attention.

Metatranscriptome analysis of the intestinal microorganisms in Pardosa pseudoannulata in response to cadmium stress

Cadmium (Cd) generates a variety of physiological and ecological toxicity to spiders. However, little is known about the effects of Cd on symbiotic bacteria of spiders. Metatranscriptomics is increasing our knowledge of microorganisms in environment. To better understand the impact of Cd on the symbiotic bacteria of spiders, we generated and compared the metatranscriptomes of the intestinal microorganisms of Pardosa pseudoannulata with and without Cd stress. The community structure of intestinal microorganisms in P. pseudoannulata was composed of 4 kingdoms, namely bacteria, viruses, eukaryotes and archaea, including 46 phyla, 97 classes, 184 orders, 339 families, 470 genera, and 598 species. The abundance of eukaryotes, bacteria and viruses was decreased by 0.14%, 1.22% and 2.52% respectively while the archaea was increased by 99.16% when under Cd stress. We identified 1519 differentially expressed genes (DEGs), including 770 up-regulated and 749 down-regulated genes. The results of KEGG annotation revealed that the expression of genes that are involved in the carbon metabolism, protein and amino acid metabolism and synthesis, glucose metabolism, oxidative phosphorylation, and glutathione metabolism were influenced by Cd. Collectively, these findings showed that Cd significantly impacted the community structure and expression of related functional genes of intestinal microorganisms in P. pseudoannulata.

Transcriptome assembly and expression profiling of the molecular responses to cadmium toxicity in cerebral ganglia of wolf spider Pardosa pseudoannulata (Araneae: Lycosidae)

Cadmium (Cd) is a heavy metal that can cause irreversible toxicity to animals, and is an environmental pollutant in farmlands. Spiders are considered to be an excellent model for investigating the impacts of heavy metals on the environment. To date, the changes at the molecular level in the cerebral ganglia of spiders are poorly understood. Cd exposure leads to strong damage in the nervous system, such as apoptosis and necrosis of nerve cells, therefore we conducted a transcriptomic analysis of Pardosa pseudoannulata cerebral ganglia under Cd stress to profile differential gene expression (DGE). We obtained a total of 123,328 assembled unigenes, and 1441 Cd stress-associated DEGs between the Cd-treated and control groups. Expression profile analysis demonstrated that many genes involved in calcium signaling, cGMP-PKG signaling, tyrosine metabolism, phototransduction-fly, melanogenesis and isoquinoline alkaloid biosynthesis were up-regulated under Cd stress, whereas oxidative phosphorylation-related, nervous disease-associated, non-alcoholic fatty liver disease-associated, and ribosomal-associated genes were down-regulated. Here, we provide a comprehensive set of DEGs influenced by Cd stress, and heavy metal stress, and provide new information for elucidating the neurotoxic mechanisms of Cd stress in spiders.

Comparative transcriptome analysis of Glyphodes pyloalis Walker (Lepidoptera: Pyralidae) reveals novel insights into heat stress tolerance in insects

BACKGROUND

Heat tolerance is a key parameter that affects insect distribution and abundance. Glyphodes pyloalis Walker (Lepidoptera: Pyralidae) is a devastating pest of mulberry in the main mulberry-growing regions and can cause tremendous losses to sericulture by directly feeding on mulberry leaves and transmitting viruses to Bombyx mori. Moreover, G. pyloalis shows a prominent capacity for adaptation to daily and seasonal temperature fluctuations and can survive several hours under high temperature. To date, the molecular mechanism underlying the outstanding adaptability of this pest to high temperature remains unclear.

RESULTS

In this study, we performed comparative transcriptome analyses on G. pyloalis exposed to 25 and 40 °C for 4 h. We obtained 34,034 unigenes and identified 1275 and 1222 genes significantly upregulated or downregulated, respectively, by heat stress. Data from the transcriptome analyses indicated that some processes involved in heat tolerance are conserved, such as high expression of heat shock protein (HSP) genes and partial repression of metabolism progress. In addition, vitamin digestion and absorption pathways and detoxification pathways identified here provided new insights for the investigation of the molecular mechanisms of heat stress tolerance. Furthermore, transcriptome analysis indicated that immune and phosphatidylinositol signaling system have a close relationship with heat tolerance. In addition, the expression patterns of ten randomly selected genes, such as HSP and cytochrome P450, were consistent with the transcriptome results obtained through quantitative real-time PCR.

CONCLUSIONS

Comparisons among transcriptome results revealed the upregulation of HSPs and genes involved in redox homeostasis, detoxication, and immune progress. However, many metabolism progresses, such as glycolysis/gluconeogenesis and fatty acid biosynthesis, were partially repressed. The results reflected that the heat tolerance of G. pyloalis is a fairly complicated process and related to a broad range of physiological regulations. Our study can improve understanding on the mechanisms of insect thermal tolerance.

Transcriptomic response of wolf spider, Pardosa pseudoannulata, to transgenic rice expressing Bacillus thuringiensis Cry1Ab protein

BACKGROUND

Bacillum thuringiensis (Bt) toxin produced in Cry1-expressing genetically modified rice (Bt rice) is highly effective to control lepidopteran pests, which reduces the needs for synthetic insecticides. Non-target organisms can be exposed to Bt toxins through direct feeding or trophic interactions in the field. The wolf spider Pardosa pseudoannulata, one of the dominant predators in South China, plays a crucial role in the rice agroecosystem. In this study, we investigated transcriptome responses of the 5th instar spiders fed on preys maintained on Bt- and non-Bt rice.

RESULTS

Comparative transcriptome analysis resulted in 136 differentially expressed genes (DEGs) between spiderlings preying upon N. lugens fed on Bt- and non-Bt rice (Bt- and non-Bt spiderlings). Functional analysis indicated a potential impact of Bt toxin on the formation of new cuticles during molting. GO and KEGG enrichment analyses suggested that GO terms associated with chitin or cuticle, including "chitin binding", "chitin metabolic process", "chitin synthase activity", "cuticle chitin biosynthetic process", "cuticle hydrocarbon biosynthetic process", and "structural constituent of cuticle", and an array of amino acid metabolic pathways, including "alanine, asparatate and glutamate metabolism", "glycine, serine and theronine metabolism", "cysteine and methionine metabolism", "tyrosine metabolism", "phenylalanine metabolism and phenylalanine", and "tyrosine and tryptophan biosynthesis" were significantly influenced in response to Cry1Ab.

CONCLUSIONS

The Cry1Ab may have a negative impact on the formation of new cuticles during molting, which is contributed to the delayed development of spiderlings. To validate these transcriptomic responses, further examination at the translational level will be warranted.

Transcriptome Profiling Analysis of Wolf Spider Pardosa pseudoannulata (Araneae: Lycosidae) after Cadmium Exposure

Pardosa pseudoannulata is one of the most common wandering spiders in agricultural fields and a potentially good bioindicator for heavy metal contamination. However, little is known about the mechanisms by which spiders respond to heavy metals at the molecular level. In the present study, high-throughput transcriptome sequencing was employed to characterize the de novo transcriptome of the spiders and to identify differentially expressed genes (DEGs) after cadmium exposure. We obtained 60,489 assembled unigenes, 18,773 of which were annotated in the public databases. A total of 2939 and 2491 DEGs were detected between the libraries of two Cd-treated groups and the control. Functional enrichment analysis revealed that metabolism processes and digestive system function were predominately enriched in response to Cd stress. At the cellular and molecular levels, significantly enriched pathways in lysosomes and phagosomes as well as replication, recombination and repair demonstrated that oxidative damage resulted from Cd exposure. Based on the selected DEGs, certain critical genes involved in defence and detoxification were analysed. These results may elucidate the molecular mechanisms underlying spiders' responses to heavy metal stress.