Analysis of digital gene expression profiling in the gonad of male silkworms (Bombyx mori) under fluoride stress

Fluorine is an essential element, but excessive fluoride can cause serious effects on the respiratory, digestive, and reproductive systems. Fluorine has been suggested to cause reproductive toxicity in vertebrates, but its potential to reproductively affect invertebrates remains unknown. In the present study, the lepidopteran model insect Bombyx mori was used to assess the reproductive toxicity of NaF. The underlying molecular mechanisms were explored by RNA sequencing, and we investigated the testes transcriptomic profile of B. mori treated with NaF via a digital gene expression (DGE) analysis. Among 520 candidate genes, 297 and 223 were identified as significantly upregulated or downregulated, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were carried out on all genes to determine their biological functions and associated processes. The results indicated that numerous differentially expressed genes are involved in the stress response, detoxification, antibacterial, transport, oxidative phosphorylation, and ribosome. The reliability of the data was confirmed by a quantitative real-time polymerase chain reaction (qRT-PCR) analysis. The changed Glutathione S-transferase (GST) activity and glutathione (GSH) content in the NaF-treated groups were increased and reduced respectively. This study reveals that using RNA-sequencing for the transcriptome profiling of B. mori testes can lead to better comprehension of the male reproductive toxicity effects of NaF. Furthermore, we expect that these results will aid future molecular studies on the reproductive toxicity of NaF in other species.