BACKGROUND
genetically modified  technology has reshaped the agricultural industry since its insertion in the late1990s  <cit> .
from  <dig> to  <dig>  the global acreage of gm crops has increased dramatically from  <dig>  to  <dig>  million hectares  <cit> .
the ecological benefits from rapid development and adoption of gm crops include a significant reduction in both insecticide and herbicide usage and greenhouse gas emissions  <cit> .
a meta-analysis in  <dig> showed a 37% reduction in synthetic pesticide use, 22% increase in crop yield, and 68% increase in farmer profits  <cit> .
besides yield and profit gains and environmental benefits, non-monetary incentives include time savings, ease of use, and more flexibility in planning  <cit> .
with limited arable land in china, gm technology provides a potential solution to improve agricultural productivity and sustainability.
currently, transgenic bacillum thuringiensis  cotton, resistant to lepidoptera pests, is the most successful commercial gm crop in china  <cit> .
in  <dig>  acreage of bt cotton has reached  <dig>  million hectares, representing 80% of total cotton area in china.
bt cotton increased yield by 10%, reduced insecticide use by 60% and generated additional us $ <dig> profit per hectare on average  <cit> .
even with the success of bt cotton, consumers still have doubts about gm crops, partially due the lack of knowledge regarding the ecological risks  <cit> .
bt rice is facing the same challenges for the public acceptance.
the community structure of a rice field is primarily composed of soil organisms, rice, insect herbivores, predators, and parasitoids.
while insect herbivores are exposed to bt toxins by direct feeding, other community members can access bt toxins through trophic interactions.
previous risk assessment studies showed no harmful effect of bt rice on diversity, dominant species and abundance of non-target arthropods among the arthropod community in the field  <cit> .
laboratory studies, on one hand, did not detect adverse impacts of bt rice on non-target arthropods.
for example, the developmental time, fecundity and survival rate of herbivorous insects nilaparvata lugens and sogatella furcitera were unaffected when exposed to cry1c, cry2a, and cry1ac proteins, respectively  <cit> .
no significant effects were found on life history traits for predators as well, including chrysoperia sinica, propylea japonica, cyrtorhinus lividipennis, and ummeliata insecticeps .
on the other hand, some reports show non-target organisms may be susceptible to bt toxins.
a significant longer developmental time of pirata subpiraticus was recorded when it prayed on bt rice fed cnaphalocrocis medinalis  <cit> .
significantly lower catalase activity was found in fosomia candida fed on bt rice in comparison to those fed on non-bt rice  <cit> .
due to the varing degradation of bt toxin protein in soils with different physicochemical properties  <cit> , researchers did not find consistent differences in soil microorganism communities between bt and non bt rice fields  <cit> .
for parasitoids, effects of bt rice is also inconsistent, which depends on the host species, target or non-target insects  <cit> .
as a whole, risk assessment of bt rice has been focusing on the organismal level impacts, suborganismal impacts are largely unknown.
the advent of genomics era, however, allows us to evaluate ecological risks of transgenic bt rice on non-target organism at the transcription and translational level.
the wolf spider pardosa pseudoannulata is one of the dominant predators in south china, playing a crucial role in maintaining the stability of the rice agroecosystem  <cit> .
in this study, we carried out a comparative transcriptome analysis of the 5th instar spiders fed on n. lugens maintained on bt- and non-bt rice, respectively.
developmental time from the 2nd to 8th instars was recorded to reveal the potential impacts of bt rice on p. pseudoannulata and to correlate the biological impacts with differentially expressed genes.
methods
plant materials and nilaparvata lugens preparation
transgenic shanyou  <dig> rice expressing cry1ab protein  and its non-transgenic parental wild type shanyou  <dig> rice  were obtained from the life science college, hunan normal university.
both rice varieties were grown under nylon nets  without insecticide application during the entire experimental period.
nilaparvata lugens were collected from farmland in the hunan academy of agricultural science and reared on non-transgenic parental wild type, shanyou  <dig>  allowing for natural colonization.
the newly moulted 2nd instar n. lugens nymphs were then transplanted to transgenic and control rice lines.
after 15-day feeding, n. lugens was collected and used as spider diets  <cit>
spider sample collection
female spiders with egg sacs were collected from the experimental farmland in the hunan academy of agricultural science.
pardosa pseudoannulata larvae were collected immediately after hatch and placed individually in a glass tube with a moist cotton ball separately .
spiders in the test and control group were fed daily with n. lugens consumed bt and non-bt rice, respectively.
all tubes were marked and maintained in an artificial climate chamber .
developmental time of each spiderling at each instar was recorded until sexual maturity was reached.
in this analysis,  <dig> spiders were raised for the developmental time recording .
observation was made twice a day at 9 am and 9 pm, respectively.
quantification of the bt toxin, cry1ab, in spiderling
an enzyme-linked immunosorbent assay  was conducted for cry1ab protein detection using a qualipate kit for cry1ab/cry1ac .
for each treatment, five 5th instar spiderlings were weighed as a group , homogenized in 1 ml pbs buffer and centrifuged for 20 min at  <dig>  g. the supernatant was used to determine cry1ab concentration.
spectrophotometric measurements for three technical replications were obtained using a microplate reader  at 450 nm.
purified cry1ab toxin  at concentrations of  <dig>   <dig>   <dig>   <dig>   <dig>  and 24 ng/l was used to generate a standard curve.
three biological replicates were performed in elisa assays.
rna isolation and illumina sequencing
a total of ten 5th instar spiderlings from the test and control groups were collected, respectively, on the ninth day after moulting and submitted to oebiotech enterprise  for rna extraction and sequencing.
total rna was extracted from each sample using trizol  according to the manufacturer’s instructions.
rna quality was assessed using a nanodrop nd- <dig> spectrophotometer  using a standard of <dig>  ≤ od260/od280 ≤  <dig>  and was further confirmed by agarose gel electrophoresis.
rna sequencing libraries were constructed and sequenced on flow cells using an illumina hiseq  <dig> platform.
clean reads were assembled using the de novo transcriptome assembler trinity after removing adaptor sequences, low quantity reads , and duplicate sequences  <cit> .
the libraries were established and unigenes of length greater than 200 bp were subjected to subsequent sequence annotation analysis.
all raw reads were deposited in the ncbi sequence read archive .
sequence annotation
all unigenes were compared to those available in the ncbi non-redundant protein  database and swiss-prot database using blastx with an e-value cutoff of 10- <dig>  the blast2go program and wego software were used to obtain go annotation for all unigenes  <cit> .
kegg  metabolic pathway annotation and cog  classification of unigenes were determined by blastx searching against kegg and cog databases  <cit> .
the best aligning results were used to determine potential function of the unigenes.
identification of differentially expressed genes
the fpkm  method was used for quantifying gene expression levels and was able to eliminate the influence of different gene lengths and sequencing levels in the calculation of gene expression  <cit> .
the degseq software package  was used to screen differentially expressed genes  based on a statistical analysis of negative binomial distribution and to quantify the gene expression levels with basemean values  <cit> .
a threshold for false discovery rate of < <dig>  and absolute value of log <dig>  ratio >  <dig> were used to determine significant differences in gene expression.
functional annotation of degs
all degs were searched against five public databases, swiss prot, nr, cog, go, and kegg.
the hypergeometric test was used to find significantly enriched go terms in degs based on go annotation.
the calculated p value then underwent bonferroni correction, using corrected p value ≤  <dig>  as a threshold.
go terms fulfilling this condition are defined as significantly enriched go terms in degs.
similarly, pathway enrichment analysis was conducted to identify significantly enriched metabolic pathways or signal transduction pathways in degs, using p value ≤  <dig>  as a threshold.
quantitative real-time pcr analysis
transcriptome results were verified using quantitative real-time pcr .
total rna was isolated from each sample with trizol  and subjected to dnase i treatment  according to the manufacturers’ protocols.
cdna was synthesized with a revertaid™ h minus first strand cdna synthesis kit  and qpcr was performed using the abi  <dig> ht system  with a reaction volume of 25 μl containing 1 μl of 1: <dig> diluted cdna in ddh2o,  <dig>  μl of  <dig> × sybr green master mix  and  <dig> nm of each of the primers.
the qpcr conditions were  <dig> °c for 3 min, followed by 40 cycles of 94 °c for 30 s for denaturation, 55 °c for 30 s for annealing and 72 °c for 30 s for extension.
the experiment was repeated three times, and expression levels of each gene were normalized to 18s ribosomal rna .
all primers were designed using the primer  <dig>  program   <cit> .
data analysis
data on developmental time of spiderlings and qpcr were analysed using a t-test with spss  <dig>  software.
significant differences at p <  <dig>  were designated with *, and data were presented as the mean ± se.
RESULTS
bt rice affects the developmental time of p. pseudoannulata
the cry1ab protein content of 5th instar p. pseudoannulata spiderlings fed on n. lugens maintained on cry1ab rice was  <dig>  ng/g, while no detectable level of cry1ab was observed in the spiderlings fed on n. lugens maintained on wild type rice.
for total developmental time, the bt spiderlings spent more time to reach the mature stage when compared to the controls .
the intermoult period of each instar in bt spiderlings was also longer than that of controls at each instar, except for the 2nd and 8th  .fig.
<dig> developmental time of p. pseudoannulata spiderling.
a the left column showed developmental time of spiderlings in the control group, and the right column showed that from the test group.
t-test was used to determine the significance .
b qpcr analysis.
technical replicates were performed for each of three biological replicates.
data were presented as mean ± se
illumina sequencing and de novo assembly
sequences of mrnas pooled from the whole body of spiderlings were analyzed using an illumina  <dig> platform and resulted in  <dig> , <dig> and  <dig> , <dig> raw reads for non-bt and bt spiderlings, respectively.
trinity software was used for de novo assembly according to standard parameters.
the assembly yielded  <dig>  total transcripts  with an average length of 612 bp, and the unigene dataset included  <dig>  sequences with an average length of 537 bp.
all unigenes were used for the annotation.
annotation of all assembled unigenes
a total of  <dig>  unigene sequences were subjected to blast searching against five public available databases, including nr, swiss-prot, cog, go and kegg, with a cutoff e < 10− <dig>  of these,  <dig>   could be matched to nr,  <dig>   to swiss-prot,  <dig>   to kegg,  <dig>   to cog and  <dig>   unigenes to go database.
the ontology covers three domains: cellular component, molecular function and biological process  <cit> .
based on similarity search,  <dig>  sequences among all uingenes were categorized into  <dig>  go terms consisting of three domains: cellular component category , molecular function  and biological process  and a few fell into more than one gene function group .
most sequences were annotated as “cellular process”, “single-organism process”, “metabolic process” and “binding”.
the results from go classification presented a portion of unigenes that had the potential to impact molting, a necessary process for p. pseudoannulata to grow.
specifically, many unigenes corresponded to seen chitin-related go terms  and  <dig> cuticle-related go terms , including  <dig> in the category of “chitin binding”,  <dig> in “chitin metabolic process” ,  <dig> in “cuticle chitin catabolic process”, and  <dig> in “structural constituent of cuticle”.fig.
<dig> go classification for all unigenes of p. pseudoannulata spiderling.
go categories shown on the x-axis were grouped into three main ontologies: biological process, cellular component and molecular function.
the right y-axis indicated the number of genes in each category while the left y-axis represented the percentage of genes in that category.
the blast2go program was used to obtain go annotation of all unigenes
fig.
<dig> go categories associated with chitin synthesis or cuticle formation.
a go categories related to chitin synthesis.
b go categories related to cuticle formation.
c degs associated with chitin synthesis and cuticle formation between bt- and non-bt spiderdlings.
a unigene may be placed in more than one group increasing the total number of genes
cog is a database where orthologous genes are classified.
twenty eight thousand six hundred forty-six unigenes among all unigenes were classified and divided into  <dig> cog categories, among which a cluster for “signal transduction mechanisms” represented the largest group  followed by “general function prediction only”  and “posttranslational modification, protein turnover, chaperones” .
the categories of “cell motility” , “nuclear structure”  and “coenzyme transport and metabolism”  had the fewest matching genes .
a group of  <dig> unigenes were assigned to the category of “chitinase”, which is responsible for dissolving the endocuticle during molting.fig.
<dig> histogram of cog classification.
a total of  <dig>  unigenes were assigned to  <dig> categories in the cog classification.
the right legend shows a description of the  <dig> function categories.
some of them can be placed in more than one cog
identification of degs and functional analysis
in total,  <dig> degs, including  <dig> down- and four up-regulated, were detected in the bt spiderlings, compared to the controls .
among all degs,  <dig> had a fold change ranging from five to ten and the remaining  <dig> had a fold change of more than ten .
for all degs,  <dig>  were assigned to nr,  <dig>  to swiss prot,  <dig>  to cog,  <dig>  to kegg and  <dig>  to go database.
interestingly, several degs were found to regulate biological processes associated with cuticle or chitin, including  <dig> down-regulated degs that were categorized to the go term “structural constituent of cuticle”,  <dig> to “chitin binding”,  <dig> to “chitin metabolic process” and  <dig> to “cuticle hydrocarbon biosynthetic process”, implying a potential effect of cry1ab on spiderling molting .fig.
<dig> distribution of differentially expressed unigenes between bt- and non-bt spiderlings.
the threshold was based on fdr <  <dig>  and absolute value of log <dig> ratio > <dig>
the clusters of degs were further characterized by go and kegg enrichment analyses, respectively  <cit> .
degs were significantly enriched to  <dig> go terms , including six belong to “cellular component” ,  <dig> to “molecular function”  and  <dig> to “biological process”  .
six of the  <dig> enriched go terms were related to chitin and cuticle, including three terms which belong to molecular function categories, “structural constituent of cuticle”, “chitin binding” and “chitin synthase activity”, and three to biological process, “chitin metabolic process”, “cuticle chitin biosynthetic process” and “cuticle hydrocarbon biosynthetic process” .
metabolic pathway enrichment analysis demonstrated that  <dig> degs were involved in  <dig> pathways .
among them,  <dig> functional pathways were significantly enriched , seven of which have been associated with amino acid metabolism, including alanine, asparatate and glutamate metabolism ; glycine, serine and theronine metabolism ; cysteine and methionine metabolism ; tyrosine metabolism , phenylalanine metabolism ; phenylalanine, tyrosine and tryptophan biosynthesis  and biosynthesis of unsaturated fatty acids .
apparently, various amino acid metabolic pathways in p. pseudoannulata spiderlings were affected in response to cry1abfig.
<dig> categories of significantly enriched go terms within degs between bt and non-bt spiderdlings.
numbers in brackets indicated the sum of annotated degs
real-time pcr assays
to validate rna-seq results, five degs were randomly selected for qpcr analysis.
these genes were homologous to structural constituents of cuticle, chitin binding and chitin metabolic process .
the expression profile of all five degs was consistent with rna-seq data .
discussion
bt toxins can be transferred via the food web and accumulate in organisms to different degrees  <cit> .
the level of bt toxin protein in predators mainly depends on expression patterns of bt-protein in plants, and the feeding behavior of the herbivore  <cit> .
our tritrophic bioassay indicated the accumulative cry1ab content in 5th instar spdierling was  <dig>  ng/g when p. pseudoannulata was preyed on n. lugens maintained on cry1ab rice.
although this protein level is slightly lower than those in ummeliata insecticeps   <cit> , it is still informative.
developmental time of bt spiderlings was significantly prolonged, which is consistent with pirata subpiraticus  <cit> .
however, spiderlings were able to recover from the effect of bt rice at a later instar.
similar to other arthropods, p. pseudoannulata must molt periodically to grow.
the formation of new cuticle is a vital step during molting of arthropods  <cit> .
we speculated that the delayed development of spiderlings may be due to the disruption of chitin synthesis  during molting.
comparative transcriptome analysis identified  <dig> degs between bt- and non-bt spiderlings .
furthermore, go annotation and enrichment analysis both suggested potential impacts of bt rice on the chitin synthesis and cuticle formation .
as with other arthropods, the exoskeleton of spider is made of cuticle, of which one of the primary component is chitin  <cit> .
the molting process in spiders involves activation of hypodermal cells, secretion of exuvial fluid and apolysis, activation of enzymes in the exuvial fluid, and secretion of the new cuticle  <cit> .
functional analysis of degs suggested a disruption of new cuticle formation during molting.
in addition, go and kegg enrichment analyses indicated 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 affected in response to cry1ab.
the advent of genomic era offers new transcriptome resources for the study of wolf spiders.
meng et al.
sequenced cephalothoraxes of p. pseudoannulata adults and identified genes involved in insecticide metabolism and detoxification, including p450s, gsts, aches, achrs, gaba receptors, and gluci  <cit> .
xiao et al.
carried out rna-seq analysis in p. pseudoannulata and revealed an array of genes responding to temperature stress  <cit> .
in this study, we focused on the genes corresponding to ingested bt toxins.
as a non-model animal without a reference genome, omics resources, such as transcriptomes, lay the foundation for future functional genomic research.
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.
additional files
additional file 1: table s <dig>  degs between bt and non-bt spiderlings.
in all  <dig> degs, including  <dig> down- and  <dig> up-regulated, were detected in the bt spiderlings, compared to the controls.
additional file 2: table s <dig>  selected degs for qpcr and their primers.
these genes were homologous to structural constituents of cuticle, chitin binding and chitin metabolic process.
abbreviations
bphsbrown planthoppers
bt
bacillus thuringiensis
cogcluster of orthologous groups of proteins
degdifferential expression genes
fpkmnumber of reads per kb of exon region per million mapped reads
gmgenetically modified
gogene ontology
keggkyoto encyclopedia of genes and genomes database
pbs bufferphosphate buffer solution
qpcrreal-time pcr
rna-seqrna sequencing

 genetically modified technology has reshaped the agricultural industry since its insertion in the late1990s. from to the global acreage of gm crops has increased dramatically from to million hectares. the ecological benefits from rapid development and adoption of gm crops include a significant reduction in both insecticide and herbicide usage and greenhouse gas emissions. a meta-analysis in showed a 37% reduction in synthetic pesticide use, 22% increase in crop yield, and 68% increase in farmer profits. besides yield and profit gains and environmental benefits, non-monetary incentives include time savings, ease of use, and more flexibility in planning. with limited arable land in china, gm technology provides a potential solution to improve agricultural productivity and sustainability. currently, transgenic bacillum thuringiensis cotton, resistant to lepidoptera pests, is the most successful commercial gm crop in china. in acreage of bt cotton has reached million hectares, representing 80% of total cotton area in china. bt cotton increased yield by 10%, reduced insecticide use by 60% and generated additional us $ profit per hectare on average. even with the success of bt cotton, consumers still have doubts about gm crops, partially due the lack of knowledge regarding the ecological risks. bt rice is facing the same challenges for the public acceptance. the community structure of a rice field is primarily composed of soil organisms, rice, insect herbivores, predators, and parasitoids. while insect herbivores are exposed to bt toxins by direct feeding, other community members can access bt toxins through trophic interactions. previous risk assessment studies showed no harmful effect of bt rice on diversity, dominant species and abundance of non-target arthropods among the arthropod community in the field. laboratory studies, on one hand, did not detect adverse impacts of bt rice on non-target arthropods. for example, the developmental time, fecundity and survival rate of herbivorous insects nilaparvata lugens and sogatella furcitera were unaffected when exposed to cry1c, cry2a, and cry1ac proteins, respectively. no significant effects were found on life history traits for predators as well, including chrysoperia sinica, propylea japonica, cyrtorhinus lividipennis, and ummeliata insecticeps. on the other hand, some reports show non-target organisms may be susceptible to bt toxins. a significant longer developmental time of pirata subpiraticus was recorded when it prayed on bt rice fed cnaphalocrocis medinalis. significantly lower catalase activity was found in fosomia candida fed on bt rice in comparison to those fed on non-bt rice. due to the varing degradation of bt toxin protein in soils with different physicochemical properties, researchers did not find consistent differences in soil microorganism communities between bt and non bt rice fields. for parasitoids, effects of bt rice is also inconsistent, which depends on the host species, target or non-target insects. as a whole, risk assessment of bt rice has been focusing on the organismal level impacts, suborganismal impacts are largely unknown. the advent of genomics era, however, allows us to evaluate ecological risks of transgenic bt rice on non-target organism at the transcription and translational level. the wolf spider pardosa pseudoannulata is one of the dominant predators in south china, playing a crucial role in maintaining the stability of the rice agroecosystem. in this study, we carried out a comparative transcriptome analysis of the 5th instar spiders fed on n. lugens maintained on bt- and non-bt rice, respectively. developmental time from the 2nd to 8th instars was recorded to reveal the potential impacts of bt rice on p. pseudoannulata and to correlate the biological impacts with differentially expressed genes. methods plant materials and nilaparvata lugens preparation transgenic shanyou rice expressing cry1ab protein and its non-transgenic parental wild type shanyou rice were obtained from the life science college, hunan normal university. both rice varieties were grown under nylon nets without insecticide application during the entire experimental period. nilaparvata lugens were collected from farmland in the hunan academy of agricultural science and reared on non-transgenic parental wild type, shanyou allowing for natural colonization. the newly moulted 2nd instar n. lugens nymphs were then transplanted to transgenic and control rice lines. after 15-day feeding, n. lugens was collected and used as spider diets spider sample collection female spiders with egg sacs were collected from the experimental farmland in the hunan academy of agricultural science. pardosa pseudoannulata larvae were collected immediately after hatch and placed individually in a glass tube with a moist cotton ball separately. spiders in the test and control group were fed daily with n. lugens consumed bt and non-bt rice, respectively. all tubes were marked and maintained in an artificial climate chamber. developmental time of each spiderling at each instar was recorded until sexual maturity was reached. in this analysis, spiders were raised for the developmental time recording. observation was made twice a day at 9 am and 9 pm, respectively. quantification of the bt toxin, cry1ab, in spiderling an enzyme-linked immunosorbent assay was conducted for cry1ab protein detection using a qualipate kit for cry1ab/cry1ac. for each treatment, five 5th instar spiderlings were weighed as a group, homogenized in 1 ml pbs buffer and centrifuged for 20 min at  g. the supernatant was used to determine cry1ab concentration. spectrophotometric measurements for three technical replications were obtained using a microplate reader at 450 nm. purified cry1ab toxin at concentrations of and 24 ng/l was used to generate a standard curve. three biological replicates were performed in elisa assays. rna isolation and illumina sequencing a total of ten 5th instar spiderlings from the test and control groups were collected, respectively, on the ninth day after moulting and submitted to oebiotech enterprise for rna extraction and sequencing. total rna was extracted from each sample using trizol according to the manufacturer’s instructions. rna quality was assessed using a nanodrop nd- spectrophotometer using a standard of  ≤ od260/od280 ≤  and was further confirmed by agarose gel electrophoresis. rna sequencing libraries were constructed and sequenced on flow cells using an illumina hiseq platform. clean reads were assembled using the de novo transcriptome assembler trinity after removing adaptor sequences, low quantity reads, and duplicate sequences. the libraries were established and unigenes of length greater than 200 bp were subjected to subsequent sequence annotation analysis. all raw reads were deposited in the ncbi sequence read archive. sequence annotation all unigenes were compared to those available in the ncbi non-redundant protein database and swiss-prot database using blastx with an e-value cutoff of 10- the blast2go program and wego software were used to obtain go annotation for all unigenes. kegg metabolic pathway annotation and cog classification of unigenes were determined by blastx searching against kegg and cog databases. the best aligning results were used to determine potential function of the unigenes. identification of differentially expressed genes the fpkm method was used for quantifying gene expression levels and was able to eliminate the influence of different gene lengths and sequencing levels in the calculation of gene expression. the degseq software package was used to screen differentially expressed genes based on a statistical analysis of negative binomial distribution and to quantify the gene expression levels with basemean values. a threshold for false discovery rate of < and absolute value of log ratio >  were used to determine significant differences in gene expression. functional annotation of degs all degs were searched against five public databases, swiss prot, nr, cog, go, and kegg. the hypergeometric test was used to find significantly enriched go terms in degs based on go annotation. the calculated p value then underwent bonferroni correction, using corrected p value ≤  as a threshold. go terms fulfilling this condition are defined as significantly enriched go terms in degs. similarly, pathway enrichment analysis was conducted to identify significantly enriched metabolic pathways or signal transduction pathways in degs, using p value ≤  as a threshold. quantitative real-time pcr analysis transcriptome results were verified using quantitative real-time pcr. total rna was isolated from each sample with trizol and subjected to dnase i treatment according to the manufacturers’ protocols. cdna was synthesized with a revertaid™ h minus first strand cdna synthesis kit and qpcr was performed using the abi ht system with a reaction volume of 25 μl containing 1 μl of 1: diluted cdna in ddh2o,  μl of × sybr green master mix and nm of each of the primers. the qpcr conditions were °c for 3 min, followed by 40 cycles of 94 °c for 30 s for denaturation, 55 °c for 30 s for annealing and 72 °c for 30 s for extension. the experiment was repeated three times, and expression levels of each gene were normalized to 18s ribosomal rna. all primers were designed using the primer program. data analysis data on developmental time of spiderlings and qpcr were analysed using a t-test with spss software. significant differences at p <  were designated with *, and data were presented as the mean ± se.  bt rice affects the developmental time of p. pseudoannulata the cry1ab protein content of 5th instar p. pseudoannulata spiderlings fed on n. lugens maintained on cry1ab rice was  ng/g, while no detectable level of cry1ab was observed in the spiderlings fed on n. lugens maintained on wild type rice. for total developmental time, the bt spiderlings spent more time to reach the mature stage when compared to the controls. the intermoult period of each instar in bt spiderlings was also longer than that of controls at each instar, except for the 2nd and 8th.fig. developmental time of p. pseudoannulata spiderling. a the left column showed developmental time of spiderlings in the control group, and the right column showed that from the test group. t-test was used to determine the significance. b qpcr analysis. technical replicates were performed for each of three biological replicates. data were presented as mean ± se illumina sequencing and de novo assembly sequences of mrnas pooled from the whole body of spiderlings were analyzed using an illumina platform and resulted in, and, raw reads for non-bt and bt spiderlings, respectively. trinity software was used for de novo assembly according to standard parameters. the assembly yielded total transcripts with an average length of 612 bp, and the unigene dataset included sequences with an average length of 537 bp. all unigenes were used for the annotation. annotation of all assembled unigenes a total of unigene sequences were subjected to blast searching against five public available databases, including nr, swiss-prot, cog, go and kegg, with a cutoff e < 10− of these, could be matched to nr, to swiss-prot, to kegg, to cog and unigenes to go database. the ontology covers three domains: cellular component, molecular function and biological process. based on similarity search, sequences among all uingenes were categorized into go terms consisting of three domains: cellular component category, molecular function and biological process and a few fell into more than one gene function group. most sequences were annotated as “cellular process”, “single-organism process”, “metabolic process” and “binding”. the results from go classification presented a portion of unigenes that had the potential to impact molting, a necessary process for p. pseudoannulata to grow. specifically, many unigenes corresponded to seen chitin-related go terms and cuticle-related go terms, including in the category of “chitin binding”, in “chitin metabolic process”, in “cuticle chitin catabolic process”, and in “structural constituent of cuticle”.fig. go classification for all unigenes of p. pseudoannulata spiderling. go categories shown on the x-axis were grouped into three main ontologies: biological process, cellular component and molecular function. the right y-axis indicated the number of genes in each category while the left y-axis represented the percentage of genes in that category. the blast2go program was used to obtain go annotation of all unigenes fig. go categories associated with chitin synthesis or cuticle formation. a go categories related to chitin synthesis. b go categories related to cuticle formation. c degs associated with chitin synthesis and cuticle formation between bt- and non-bt spiderdlings. a unigene may be placed in more than one group increasing the total number of genes cog is a database where orthologous genes are classified. twenty eight thousand six hundred forty-six unigenes among all unigenes were classified and divided into cog categories, among which a cluster for “signal transduction mechanisms” represented the largest group followed by “general function prediction only” and “posttranslational modification, protein turnover, chaperones”. the categories of “cell motility”, “nuclear structure” and “coenzyme transport and metabolism” had the fewest matching genes. a group of unigenes were assigned to the category of “chitinase”, which is responsible for dissolving the endocuticle during molting.fig. histogram of cog classification. a total of unigenes were assigned to categories in the cog classification. the right legend shows a description of the function categories. some of them can be placed in more than one cog identification of degs and functional analysis in total, degs, including down- and four up-regulated, were detected in the bt spiderlings, compared to the controls. among all degs, had a fold change ranging from five to ten and the remaining had a fold change of more than ten. for all degs, were assigned to nr, to swiss prot, to cog, to kegg and to go database. interestingly, several degs were found to regulate biological processes associated with cuticle or chitin, including down-regulated degs that were categorized to the go term “structural constituent of cuticle”, to “chitin binding”, to “chitin metabolic process” and to “cuticle hydrocarbon biosynthetic process”, implying a potential effect of cry1ab on spiderling molting.fig. distribution of differentially expressed unigenes between bt- and non-bt spiderlings. the threshold was based on fdr <  and absolute value of log ratio > the clusters of degs were further characterized by go and kegg enrichment analyses, respectively. degs were significantly enriched to go terms, including six belong to “cellular component”, to “molecular function” and to “biological process”. six of the enriched go terms were related to chitin and cuticle, including three terms which belong to molecular function categories, “structural constituent of cuticle”, “chitin binding” and “chitin synthase activity”, and three to biological process, “chitin metabolic process”, “cuticle chitin biosynthetic process” and “cuticle hydrocarbon biosynthetic process”. metabolic pathway enrichment analysis demonstrated that degs were involved in pathways. among them, functional pathways were significantly enriched, seven of which have been associated with amino acid metabolism, including alanine, asparatate and glutamate metabolism; glycine, serine and theronine metabolism; cysteine and methionine metabolism; tyrosine metabolism, phenylalanine metabolism; phenylalanine, tyrosine and tryptophan biosynthesis and biosynthesis of unsaturated fatty acids. apparently, various amino acid metabolic pathways in p. pseudoannulata spiderlings were affected in response to cry1abfig. categories of significantly enriched go terms within degs between bt and non-bt spiderdlings. numbers in brackets indicated the sum of annotated degs real-time pcr assays to validate rna-seq results, five degs were randomly selected for qpcr analysis. these genes were homologous to structural constituents of cuticle, chitin binding and chitin metabolic process. the expression profile of all five degs was consistent with rna-seq data. discussion bt toxins can be transferred via the food web and accumulate in organisms to different degrees. the level of bt toxin protein in predators mainly depends on expression patterns of bt-protein in plants, and the feeding behavior of the herbivore. our tritrophic bioassay indicated the accumulative cry1ab content in 5th instar spdierling was  ng/g when p. pseudoannulata was preyed on n. lugens maintained on cry1ab rice. although this protein level is slightly lower than those in ummeliata insecticeps, it is still informative. developmental time of bt spiderlings was significantly prolonged, which is consistent with pirata subpiraticus. however, spiderlings were able to recover from the effect of bt rice at a later instar. similar to other arthropods, p. pseudoannulata must molt periodically to grow. the formation of new cuticle is a vital step during molting of arthropods. we speculated that the delayed development of spiderlings may be due to the disruption of chitin synthesis during molting. comparative transcriptome analysis identified degs between bt- and non-bt spiderlings. furthermore, go annotation and enrichment analysis both suggested potential impacts of bt rice on the chitin synthesis and cuticle formation. as with other arthropods, the exoskeleton of spider is made of cuticle, of which one of the primary component is chitin. the molting process in spiders involves activation of hypodermal cells, secretion of exuvial fluid and apolysis, activation of enzymes in the exuvial fluid, and secretion of the new cuticle. functional analysis of degs suggested a disruption of new cuticle formation during molting. in addition, go and kegg enrichment analyses indicated 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 affected in response to cry1ab. the advent of genomic era offers new transcriptome resources for the study of wolf spiders. meng et al. sequenced cephalothoraxes of p. pseudoannulata adults and identified genes involved in insecticide metabolism and detoxification, including p450s, gsts, aches, achrs, gaba receptors, and gluci. xiao et al. carried out rna-seq analysis in p. pseudoannulata and revealed an array of genes responding to temperature stress. in this study, we focused on the genes corresponding to ingested bt toxins. as a non-model animal without a reference genome, omics resources, such as transcriptomes, lay the foundation for future functional genomic research.  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. additional files additional file 1: table s degs between bt and non-bt spiderlings. in all degs, including down- and up-regulated, were detected in the bt spiderlings, compared to the controls. additional file 2: table s selected degs for qpcr and their primers. these genes were homologous to structural constituents of cuticle, chitin binding and chitin metabolic process. abbreviations bphsbrown planthoppers bt bacillus thuringiensis cogcluster of orthologous groups of proteins degdifferential expression genes fpkmnumber of reads per kb of exon region per million mapped reads gmgenetically modified gogene ontology keggkyoto encyclopedia of genes and genomes database pbs bufferphosphate buffer solution qpcrreal-time pcr rna-seqrna sequencing 
pardosa pseudoannulatacry1abdevelopmentrna-seqchitincuticlehttp://dx.doi.org/ /501100001809national natural science foundation of china3107194331272339wang zhi scientific research key fund of hunanprovincial science and technology department10a054wang zhi agricultural science and technology innovation program of chinacaas-astip-2016-ibfcpeng yuande issue-copyright-statement© the author 2017 the wolf spider pardosa pseudoannulata, one of the dominant predators in south china, plays a crucial role in the rice agroecosystem. to validate these transcriptomic responses, further examination at the translational level will be warranted. in this study, we investigated transcriptome responses of the 5th instar spiders fed on preys maintained on bt- and non-bt rice. non-target organisms can be exposed to bt toxins through direct feeding or trophic interactions in the field. functional analysis indicated a potential impact of bt toxin on the formation of new cuticles during molting. comparative transcriptome analysis resulted in differentially expressed genes between spiderlings preying upon n. lugens fed on bt- and non-bt rice. the cry1ab may have a negative impact on the formation of new cuticles during molting, which is contributed to the delayed development of spiderlings. 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.
 bacillum thuringiensis toxin produced in cry1-expressing genetically modified 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.  comparative transcriptome analysis resulted in differentially expressed genes between spiderlings preying upon n. lugens fed on bt- and non-bt rice. 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.  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. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords pardosa pseudoannulatacry1abdevelopmentrna-seqchitincuticlehttp://dx.doi.org/ /501100001809national natural science foundation of china3107194331272339wang zhi scientific research key fund of hunanprovincial science and technology department10a054wang zhi agricultural science and technology innovation program of chinacaas-astip-2016-ibfcpeng yuande issue-copyright-statement© the author 2017 