Transcriptomic data supporting the paper ‘Chronic Paracetamol Exposure Disrupts Behavior and Induces Noncanonical Apoptosis in a Freshwater Bivalve’

The visceral mass RNA-seq transcriptome libraries were prepared using Illumina TruSeqTM RNA sample preparation Kit (San Diego, CA). To identify DEGs (differential expression genes) between two different samples, the expression level of each transcript was calculated according to the transcripts per million reads (TPM) method. We considered the genes with RPKM ≥ 1 in the samples as tissue-expressed genes and the rest as tissue-unexpressed genes(Cheng et al., 2019).

Differential expression analysis was performed using edgeR. Transcripts with a FDR ≤ 0.05 & ∣log2FC∣ ≥ 1 were considered differentially expressed genes (DEGs). The GOseq package was used to analyze the biological functional annotation and GO enrichment of the gene sets.

Furthermore, a pathway analysis of the DEGs was performed using the KEGG and GO databases. Pathways with a corrected p-value < 0.05 were designated as significantly enriched pathways.