Interactions between multiple genes are involved in the development of complex diseases. However, there are few analyses of gene interactions associated with papillary thyroid cancer (PTC). Weighted gene co-expression network analysis (WGCNA) is a novel and powerful method that detects gene interactions according to their co-expression similarities. In the present study, WGCNA was performed in order to identify functional genes associated with PTC using R package. First, differential gene expression analysis was conducted in order to identify the differentially expressed genes (DEGs) between PTC and normal samples. Subsequently, co-expression networks of the DEGs were constructed for the two sample groups, respectively. The two networks were compared in order to identify a poorly preserved module. Concentrating on the significant module, validation analysis was performed to confirm the identified genes and combined functional enrichment analysis was conducted in order to identify more functional associations of these genes with PTC. As a result, 1062 DEGs were identified for network construction. A brown module containing 118 highly related genes was selected as it exhibited the lowest module preservation. After validation analysis, 61 genes in the module were confirmed to be associated with PTC. Following the enrichment analysis, two PTC-related pathways were identified: Wnt signal pathway and transcriptional misregulation in cancer. LRP4, KLK7, PRICKLE1, ETV4 and ETV5 were predicted to be candidate genes regulating the pathogenesis of PTC. These results provide novel insights into the etiology of PTC and the identification of potential functional genes.
Keywords: ETV4; ETV5; KLK7; LRP4; PRICKLE1; co-expression; gene interactions; papillary thyroid cancer.