Weighted Gene Co-expression Network Analysis for RNA-Sequencing Data of the Varicose Veins Transcriptome

Abstract

Objective: Varicose veins are a common problem worldwide and can cause significant impairments in health-related quality of life, but the etiology and pathogenesis remain not well defined. This study aims to elucidate transcriptomic regulations of varicose veins by detecting differentially expressed genes, pathways and regulator genes. Methods: We harvested great saphenous veins (GSV) from patients who underwent coronary artery bypass grafting (CABG) and varicose veins from conventional stripping surgery. RNA-Sequencing (RNA-Seq) technique was used to obtain the complete transcriptomic data of both GSVs from CABG patients and varicose veins. Weighted Gene Co-expression network analysis (WGCNA) and further analyses were then carried out with the aim to elucidate transcriptomic regulations of varicose veins by detecting differentially expressed genes, pathways and regulator genes. Results: From January 2015 to December 2016, 7 GSVs from CABG patients and 13 varicose veins were obtained. WGCNA identified 4 modules. In the brown module, gene ontology (GO) analysis showed that the biological processes were focused on response to stimulus, immune response and inflammatory response, etc. Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis showed that the biological processes were focused on cytokine-cytokine receptor interaction and TNF signaling pathway, etc. In the gray module, GO analysis showed that the biological processes were skeletal myofibril assembly related. The immunohistochemistry staining showed that the expression of ASC, Caspase-1 and NLRP3 were increased in GSVs from CABG patients compared with varicose veins. Histopathological analysis showed that in the varicose veins group, the thickness of vascular wall, tunica intima, tunica media and collagen/smooth muscle ratio were significantly increased, and that the elastic fiber/internal elastic lamina ratio was decreased. Conclusion: This study shows that there are clear differences in transcriptomic information between varicose veins and GSVs from CABG patients. Some inflammatory RNAs are down-regulated in varicose veins compared with GSVs from CABG patients. Skeletal myofibril assembly pathway may play a crucial role in the pathogenesis of varicose veins. Characterization of these RNAs may provide new targets for understanding varicose veins diagnosis, progression, and treatment.