WGCNA reveals key gene modules regulated by the combined treatment of colon cancer with PHY906 and CPT11

Irinotecan (CPT11) is one of the most effective drugs for treating colon cancer, but its severe side effects limit its application. Recently, a traditional Chinese herbal preparation, named PHY906, has been proved to be effective for improving therapeutic effect and reducing side effects of CPT11. The aim of the present study was to provide novel insight to understand the molecular mechanism underlying PHY906-CPT11 intervention of colon cancer. Based on the GSE25192 dataset, for different three treatments (PHY906, CPT11, and PHY906-CPT11), we screened out differentially expressed genes (DEGs) and constructed a co-expression network by weighted gene co-expression network analysis (WGCNA) to identify hub genes. The key genes of the three treatments were obtained by merging the DEGs and hub genes. For the PHY906-CPT11 treatment, a total of 18 key genes including Eif4e, Prr15, Anxa2, Ddx5, Tardbp, Skint5, Prss12 and Hnrnpa3, were identified. The results of functional enrichment analysis indicated that the key genes associated with PHY906-CPT11 treatment were mainly enriched in ‘superoxide anion generation’ and ‘complement and coagulation cascades’. Finally, we validated the key genes by Gene Expression Profiling Interactive Analysis (GEPIA) and RT-PCR analysis, the results indicated that EIF4E, PRR15, ANXA2, HNRNPA3, NCF1, C3AR1, PFDN2, RGS10, GNG11, and TMSB4X might play an important role in the treatment of colon cancer with PHY906-CPT11. In conclusion, a total of 18 key genes were identified in the present study. These genes showed strong correlation with PHY906-CPT11 treatment in colon cancer, which may help elucidate the underlying molecular mechanism of PHY906-CPT11 treatment in colon cancer.