Transcriptome analysis reveals gene expression changes of the fat body of silkworm (Bombyx mori L.) in response to selenium treatment

A comparative transcriptome analysis was conducted to investigate the gene expression changes in the fat body of silkworm after treatment with different concentrations (50 μM and 200 μM) of selenium (Se). 912 differential expression genes (DEGs) (371 up-regulated and 541 down-regulated) and 1420 DEGs (1078 up-regulated and 342 down-regulated) were identified in silkworm fat body treated with 50 μM and 200 μM of Se, respectively. In case of 50 μM group, DEGs were mainly enriched in the peroxisome pathway and fatty acid metabolism pathway, and later were associated with antioxidant defense and nutrition regulation. After 200 μM Se-treatment, DEGs were mainly located in the glycerolipid metabolism and arachidonic acid metabolism pathways, which further encoded detoxification related genes. Furthermore, 32 candidate DEGs from these pathways had been selected to confirm the RNA-seq data. Among these DEGs, 14 genes were up-regulated in the 50 μM Se-treated group (only three genes in the 200 μM Se-treated group) which were involved in lipid metabolism and antioxidant defense, and 13 up-regulated genes (only two genes were up-regulated in the 50 μM Se-treated group) were involved in detoxification of the 200 μM Se-treated group. These changes showed that lower concentration of Se could regulate the nutrition and promote antioxidation pathways; whereas, high levels of Se promoted the detoxification of silkworm. These findings can be helpful to understand the possible mechanisms of Se action and detoxification in silkworm and other insects.