Figure legends
FIGURE 1 Toxic and non-toxic communities of N. succinctus and N. variciferus collected from China coast.
FIGURE 2 Distribution of annotated species for both N. succinctus and N. variciferus which show similar annotated patterns. The Aplysia california species take accounted for the top proportion.
FIGURE 3 PCA plot for difference of expression level between the toxic and non-toxic communities (ZD & HD: non-toxic specimens, ZL & HL: toxic specimens). It indicated that the toxic and non-toxic communities generally show similar expression patterns.
FIGURE 4 DEGs patterns between toxic and non-toxic specimens ofN. succinctus and N. variciferus. HD inN.succinctus and ZD in N.variciferus (non-toxic samples) were set as control group and HL and ZL were set as treatment group. Examples: HD1-VS-HL3 means that the HL3 generated more up-regulated genes (blue) than down-regulated (orange) in comparison with HD1. ZD1-VS-ZL3 means that the ZL3 generated more up-regulated genes (orange) than down-regulated (blue) in comparison with ZD1.
FIGURE 5 Consistency of 2 -(∆∆Ct) and actual expression level by RT-PCR. Six potential DEGs related with TTX accumulation and resistance was performed.
FIGURE 6 Potential amino acids mutation site related with TTX resistance on sodium channel of N. succinctus. The site indicated in yellow is newly identified by our study.
FIGURE S1 Distribution of base content and quality for all specimens studied.
FIGURE S2 The top classification categories for the DEGs of both species.
TABLE S1 Functional annotation for both Nassariusspecies from 7 databases
Table 1 . The clean reads quality metrics after Quality control.