Identification of molecular mechanisms causing skin lesions of cutaneous leishmaniasis using weighted gene coexpression network analysis (WGCNA)

Leishmaniasis is an infectious disease, caused by a protozoan parasite. Its most common form is cutaneous leishmaniasis, which leaves scars on exposed body parts from bites by infected female phlebotomine sandflies. Approximately 50% of cases of cutaneous leishmaniasis fail to respond to standard treatments, creating slow-healing wounds which cause permanent scars on the skin. We performed a joint bioinformatics analysis to identify differentially expressed genes (DEGs) in healthy skin biopsies and Leishmania cutaneous wounds. DEGs and WGCNA modules were analyzed based on the Gene Ontology function, and the Cytoscape software. Among almost 16,600 genes that had significant expression changes on the skin surrounding Leishmania wounds, WGCNA determined that one of the modules, with 456 genes, has the strongest correlation with the size of the wounds. Functional enrichment analysis indicated that this module includes three gene groups with significant expression changes. These produce tissue-damaging cytokines or disrupt the production and activation of collagen, fibrin proteins, and the extracellular matrix, causing skin wounds or preventing them from healing. The hub genes of these groups are OAS1, SERPINH1, and FBLN1 respectively. This information can provide new ways to deal with unwanted and harmful effects of cutaneous leishmaniasis.