Objective: Atopic eczema (AE) is a chronic relapsing inflammatory skin disease. This study aims to identify key genes related to the development of AE.
Materials and methods: The GSE6012 dataset was obtained from the Gene Expression Omnibus (GEO) database. The limma package was used to analyze differentially expressed genes (DEGs). Then, the weighted gene co-expression network analysis (WGCNA) package was utilized to generate weighted correlation networks of up-and downregulated genes. Additionally, the WGCNA package was used for enrichment analyses to explore the underlying functions of DEGs in modules (weighted correlation sub-networks) significantly associated with AE.
Results: A total of 515 DEGs were identified between lesional and non-lesional skin samples. For the upregulated genes, the blue module was found to have a significant positive correlation with AE. Importantly, small proline-rich protein 2C (SPRR2C) and defensin, beta 4A (DEFB4A) exhibited higher |log fold change (FC) values and were the key nodes of the network. Moreover, KEGG pathway analysis revealed that the upregulated genes in the blue module were primarily involved in cytokine-cytokine receptor interaction. Additionally, for the downregulated genes, the brown module was found to have a significant positive correlation with AE. Further, WNT inhibitory factor 1 (WIF1), cryptochrome 2 (CRY2), and keratin 19 (KRT19) had higher |log FC| values and were key nodes of the network.
Conclusions: SPRR2C, DEFB4A, WIF1, CRY2, KRT19 and cytokine-cytokine receptor interaction might be correlated with the development of AE.