Abstract: To explore the molecular mechanism of bisphenol A (BPA) toxicity effects on echiuran worm (Urechisunicinctus),echiuran worm with body weight of (17.3±1.5)g was exposed to 3 mg/L BPA (0.15% of anhydrous ethanol as co-solvent) and sea water (control group) and then the transcriptome sequencing was adopted for coelomocytes of the echiuran worm under an Illumina Hiseq 2500 platform. The results showed that 13.74 Gb of high-quality sequencing data were obtained with Q30≥93.02%, 108 262 transcripts and 30 882 Unigene were assembled, with average length of 1 380 bp and N50 of 2 549 bp. Eight data bases (NR, GO, COG, Swiss-Prot, KOG, KEGG, eggNOG4.5 and Pfam) were selected for the Unigene annotation, with 14 581 annotated Unigene, in which 7 242 Unigene were annotated to genes in KEGG database, further being classified into 284 metabolic pathways. Gene differential expression analysis indicated that 1 073 transcripts were differentially expressed in response to BPA stress, including 649 up-regulated genes and 424 down-regulated genes. GO functional classification showed that the differentially expressed genes were primarily involved in the biological process of catalytic activity, binding, antioxidant, transporter, signal transducer and electron carrier. The KEGG enrich analysis revealed that 163 genes were mapped into 151 pathways, and three most significantly enriched pathways were ribosome biogenesis in eukaryotes, glutathione metabolism, and metabolism of xenobiotics by cytochrome P450. The synthesis, immonology and detoxification metabolism of the coelomocytes elucidated the adverse effects of BPA on echiuran worm. The findings provide the reference data with the marine biotoxicity research of environmental endocrine disruptors(EEDs) and anti-stress mechanism of echiuran worm.