Identification of key genes in hepatocellular carcinoma associated with exposure to TCDD and α-endosulfan by WGCNA

2,3,7,8-tet-rachlorodibenzo-p-dioxin (TCDD) and α-endosulfan are two typical persistent organic pollutants (POPs), both of which accumulate in the liver and have potential carcinogenic hepatic effects. The underlying molecular mechanisms of pathogenesis of hepatocellular carcinoma (HCC) remain elusive when exposure to POPs. The aim of this study is to explore the key genes involved in HCC when exposure to TCDD and α-endosulfan by weighted gene co-expression network analysis (WGCNA). First, we performed co-expressed analysis on HCC and normal condition, based on WGCNA. In results, seven co-expressed modules were identified from 56 human liver samples, and the brown module correlated with five stages of HCC. Subsequently, we predicted that human five liver diseases were associated with exposure to TCDD and/or α-endosulfan by Nextbio analysis. Functional enrichment analysis showed that the brown module enriched in oxidation-reduction process, DNA replication, oxidoreductase activity and aging, which were the same as the results when exposure to the mixture of TCDD and α-endosulfan. Lastly, based on the protein-protein interaction network, we identified three novel genes including HK2, EXO1 and PFKP as key genes in HCC associated with exposure to TCDD and α-endosulfan mixture. In addition, survival analysis of key genes in Kaplan-Meier plotter demonstrated that aberrant expression levels of all the three key genes were associated with poor prognosis of HCC. Finally, Western blot analysis confirmed that protein expression levels of PFKP and HK2 in the three exposed groups were significantly elevated, while EXO1 were significantly upregulated when exposure to TCDD and α-endosulfan mixture in HepaRG cells. This study provides a new perspective to the understanding of the genetic mechanism of HCC when exposure to POPs.