Weighted gene coexpression network analysis reveals ESR1, FLNA and Furin as hub genes for DEHP-induced prepubertal testicular injury

Di-(2-ethylhexyl) phthalate (DEHP) is an environmental endocrine disruptor that accumulates in organisms in various ways and induces male reproductive system disorders. In this study, we established a testicular injury model by gavage with different concentrations of DEHP. The testes were then collected for RNA sequencing (RNA-seq), and the results were analyzed by bioinformatics and verified by experiments. Our research results show that different concentrations of DEHP interfere with testicular development differently. Weighted gene coexpression network analysis (WGCNA) generated sixteen modules and identified the turquoise module as key. Then, estrogen receptor 1 (ESR1), filamin A (Flna) and Furin were identified as hub genes. qPCR and immunohistochemistry results revealed that all three hub genes were upregulated. We detected the locations of these genes by immunohistochemistry. ESR1 was mainly located in Leydig cells; Flna immunostaining is observed in the Leydig and some germ cells and Furin staining was seen in almost all types of testicular cells. Furthermore, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses showed enrichment mainly in MAPK signaling pathways, p53 signaling pathways, HIF-1 signaling pathways, protein processing in the endoplasmic reticulum, apoptosis, the cell cycle, RNA degradation, etc. This is the first study using WGCNA to investigate the mechanism of DEHP-induced injury in the prepubertal testis, providing new research angles to further understand the mechanism of DEHP-induced injury in the prepubertal testis.