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

Abstract

ObjectiveVaricose 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.MethodsWe 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.ResultsFrom 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.ConclusionThis 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.

Investigating Transcriptional Dynamics Changes and Time-Dependent Marker Gene Expression in the Early Period After Skeletal Muscle Injury in Rats

“…We then used dynamic tree cutting to divide the gene co-expression modules according to the common gene expression file. The clearest gene coexpression module distribution obtained at early post-SMI period using the Dynamic Merge tool included 24 modules comprising 23 meaningful modules and one undefined gray module (Figure 1D) (Zhang et al, 2019b).…”

Deep Transcriptomic Analysis Reveals the Dynamic Developmental Progression during Early Development of Channel Catfish (Ictalurus punctatus)

“…WGCNA [17] is a powerful method to identify co-expressed groups of genes from large RNA expression datasets [21], and is widely used to explore the correlation among transcriptomic datasets, identify hub genes and discover new pathways in both model and non-model species [22][23][24]. WGCNA has proven its superiority over partial correlation methods and provides a powerful tool for identifying higher-order correlation in complex traits of interest, by presenting a simplified network on the integrated function of gene modules [25,26].…”

Identification of circulating hub long noncoding RNAs associated with hypertrophic cardiomyopathy using weighted correlation network analysis

“…after relating modules to clinical traits, biologically relevant modules can be identified in any given disease. compared with standard comparative analysis, WGcna can identify critical genes with key roles in the phenotype and development of a disease from interesting modules associated with important clinical traits (11). despite the strengths and popularity of network analysis, WGcna has not been employed to analyse lncrna microarray data in HcM.…”