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

Influence of Cry1Ab protein on growth and development of a predatory spider, Pardosa pseudoannulata, from protective perspectives

The expression of Cry proteins in genetically modified rice varieties safeguards the crop from lepidopteran pests. These proteins have the potential to be transferred through the food chain to arthropods like planthoppers and predatory spiders, triggering defensive responses in these unintended organisms. Hence, we hypothesized that Cry protein might influence the growth and development of spiders by altering protective enzyme activities. The results showed that Cry1Ab protein could accumulate in tissues and subcellular organelles of Pardosa pseudoannulata from Nilaparvata lugens. Cry1Ab protein exposure prolonged the developmental duration in the 5th and 7th instar spiderlings but induced no alterations of other growth indicators, such as body length, median ocular area, and survival rate. In addition, Cry1Ab protein exerted no adverse impacts on several detoxifying enzymes (i.e., superoxide dismutase, catalase, glutathione peroxidase, and acetylcholine esterase) in muscle, midgut, ganglia, and hemolymph at subcellular components (i.e., microsome and cytoplasm). To further explore the effects of Cry1Ab protein on the spiderlings, we performed an integrated transcriptome analysis on spiderlings exposed to Cry1Ab protein. The results showed that Cry1Ab protein might prolong the development duration of P. pseudoannulata via the altered cuticle metabolism (e.g., chitin metabolic process and structural constituent of cuticle). In addition, the gene expression profile associated with detoxifying enzymes and three stress-responsive pathways (JAK/STAT, JNK/SAPK, and Hippo pathways) also displayed no significant alterations under Cry1Ab exposure. Collectively, this integrated analysis generates multidimensional insights to assess the effects of Cry1Ab protein on non-target spiders and demonstrates that Cry1Ab protein exerts no toxicity in P. pseudoannulata.

Response of the Pardosa astrigera bacterial community to Cry1B protein

While genetically modified (GM) crops bring economic benefits to human beings, their impact on non-target organisms has become an important part of environmental safety assessments. Symbiotic bacteria play an important role in eukaryotic biological functions and can adjust host communities to adapt to new environments. Therefore, this study examined the effects of Cry1B protein on the growth and development of non-target natural enemies of Pardosa astrigera (L. Koch) from the perspective of symbiotic bacteria. Cry1B protein had no significant effect on the health indicators of P. astrigera (adults and 2nd instar spiderlings). 16S rRNA sequencing results revealed that Cry1B protein did not change the symbiotic bacteria species composition of P. astrigera, but did reduce the number of OTU and species diversity. In 2nd instar spiderlings, neither the dominant phylum (Proteobacteria) nor the dominant genus (Acinetobacter) changed, but the relative abundance of Corynebacterium-1 decreased significantly; in adult spiders, the dominant bacteria genera of females and males were different. The dominant bacterial genera were Brevibacterium in females and Corynebacterium-1 in males, but Corynebacterium-1 was the dominant bacteria in both females and males feeding on Cry1B. The relative abundance of Wolbachia also increased significantly. In addition, bacteria in other genera varied significantly by sex. KEGG results showed that Cry1B protein only altered the significant enrichment of metabolic pathways in female spiders. In conclusion, the effects of Cry1B protein on symbiotic bacteria vary by growth and development stage and sex.

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.

No Effect of Bt-transgenic Rice on the Tritrophic Interaction of the Stored Rice, the Maize Weevil Sitophilus Zeamais and the Parasitoid Wasp Theocolax elegans

During Bt transgenic rice storage, Bt Cry1Ab/Cry1Ac fused protein is exposed to the maize weevil Sitophilus zeamais and the parasitoid wasp Theocolax elegans. We have carried out a long-term risk assessment for Bt rice to these non-target organisms in the storehouse. Effects of Bt rice on S. zeamais and T. elegans have been carefully detected in a laboratory experiment of over 5 years. The survival, development, fecundity, and longevity of the maize weevil were compared between Bt rice and non-Bt rice treatments for every 5 generations from generation 1 to 25. Moreover, the development, adult body size and sex ratio of T. elegans were compared between them parasitizing S. zeamais feeding on Bt rice or non-Bt rice. We found that although Bt Cry1Ab/Cry1Ac fused protein exists in the Bt rice grains and S. zeamais digestive tracts, Bt rice is not harmful to the maize weevil S. zeamais and its parasitoid T. elegans.

Comparative Transcriptome Analyses Reveal Potential Mechanisms of Enhanced Drought Tolerance in Transgenic Salvia Miltiorrhiza Plants Expressing AtDREB1A from Arabidopsis

In our previous study, drought-resistant transgenic plants of Salvia miltiorrhiza were produced via overexpression of the transcription factor AtDREB1A. To unravel the molecular mechanisms underpinning elevated drought tolerance in transgenic plants, in the present study we compared the global transcriptional profiles of wild-type (WT) and AtDREB1A-expressing transgenic plants using RNA-sequencing (RNA-seq). Using cluster analysis, we identified 3904 differentially expressed genes (DEGs). Compared with WT plants, 423 unigenes were up-regulated in pRD29A::AtDREB1A-31 before drought treatment, while 936 were down-regulated and 1580 and 1313 unigenes were up- and down-regulated after six days of drought. COG analysis revealed that the 'signal transduction mechanisms' category was highly enriched among these DEGs both before and after drought stress. Based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation, DEGs associated with "ribosome", "plant hormone signal transduction", photosynthesis", "plant-pathogen interaction", "glycolysis/gluconeogenesis" and "carbon fixation" are hypothesized to perform major functions in drought resistance in AtDREB1A-expressing transgenic plants. Furthermore, the number of DEGs associated with different transcription factors increased significantly after drought stress, especially the AP2/ERF, bZIP and MYB protein families. Taken together, this study substantially expands the transcriptomic information for S. miltiorrhiza and provides valuable clues for elucidating the mechanism of AtDREB1A-mediated drought tolerance in transgenic plants.

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.

Rice Field Spiders in China: A Review of the Literature

Many laboratory and field studies have been conducted on rice field spiders in China. There are 375 species, 108 genera, and 22 families of rice field spiders distributed within the major rice growing areas and 17 dominant species. The biological and ecological characteristics of 17 rice field spider species have been reported in detail. The biology and ecology of these species show significant differences among regions, farmland habitats, and agricultural practices. Future research should focus on rice field habitat diversity, enhancing the insecticide resistance of dominant spider populations, implementing large-scale breeding of spiders and augmentative release, breeding more leaf dominant species, conducting biosafety assessment of spiders in transgenic crops.

De novo characterization of venom apparatus transcriptome of Pardosa pseudoannulata and analysis of its gene expression in response to Bt protein

Background

Pardosa pseudoannulata is a prevailing spider species, and has been regarded as an important bio-control agent of insect pests in farmland of China. However, the available genomic and transcriptomic databases of P. pseudoannulata and their venom are limited, which severely hampers functional genomic analysis of P. pseudoannulata. Recently high-throughput sequencing technology has been proved to be an efficient tool for profiling the transcriptome of relevant non-target organisms exposed to Bacillus thuringiensis (Bt) protein through food webs.

Results

In this study, the transcriptome of the venom apparatus was analyzed. A total of 113,358 non-redundant unigenes were yielded, among which 34,041 unigenes with complete or various length encoding regions were assigned biological function annotations and annotated with gene ontology and karyotic orthologous group terms. In addition, 3726 unigenes involved in response to stimulus and 720 unigenes associated with immune-response pathways were identified. Furthermore, we investigated transcriptomic changes in the venom apparatus using tag-based DGE technique. A total of 1724 differentially expressed genes (DEGs) were detected, while 75 and 372 DEGs were functionally annotated with KEGG pathways and GO terms, respectively. qPCR analyses were performed to verify the DEGs directly or indirectly related to immune and stress responses, including genes encoding heat shock protein, toll-like receptor, GST and NADH dehydrogenase.

Conclusion

This is the first study conducted to specifically investigate the venom apparatus of P. pseudoannulata in response to Bt protein exposure through tritrophic chain. A substantial fraction of transcript sequences was generated by high-throughput sequencing of the venom apparatus of P. pseudoannulata. Then a comparative transcriptome analysis showing a large number of candidate genes involved in immune response were identified by the tag-based DGE technology. This transcriptome dataset will provide a comprehensive sequence resource for furture molecular genetic research of the venom apparatus of P. pseudoannulata.

Electronic supplementary material

The online version of this article (10.1186/s12896-017-0392-z) contains supplementary material, which is available to authorized users.