Identification of Hub Genes of Keloid Fibroblasts by Coexpression Network Analysis and Degree Algorithm

Background. Keloid is a benign dermal tumor characterized by abnormal proliferation and invasion of fibroblasts. The establishment of biomarkers is essential for the diagnosis and treatment of keloids. Methods. We systematically identified coexpression modules using the weighted gene coexpression network analysis method (WGCNA). Differential expressed genes (DEGs) in GSE145725 and genes in significant modules were integrated to identify overlapping key genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were then performed, as well as protein-protein interaction (PPI) network construction for hub gene screening. Results. Using the R package of WGCNA, 22 coexpression modules consisting of different genes were identified from the top 5,000 genes with maximum mean absolute deviation in 19 human fibroblast samples. Blue-green and yellow modules were identified as the most important modules, where genes overlapping with DEGs were identified as key genes. We identified the most critical functions and pathways as extracellular structure organization, vascular smooth muscle contraction, and the cGMP-PKG signaling pathway. Hub genes from key genes as BMP4, MSX1, HAND2, TBX2, SIX1, IRX1, EDN1, DLX5, MEF2C, and DLX2 were identified. Conclusion. The blue-green and yellow modules may play an important role in the pathogenesis of keloid. 10 hub genes were identified as potential biomarkers and therapeutic targets for keloid.

ceRNA Network and Functional Enrichment Analysis of Preeclampsia by Weighted Gene Coexpression Network Analysis

Background. Preeclampsia (PE) is a multisystemic syndrome which has short- and long-term risk to mothers and children and has pluralistic etiology. Objective. This study is aimed at constructing a competitive endogenous RNA (ceRNA) network for pathways most related to PE using a data mining strategy based on weighted gene coexpression network analysis (WGCNA). Methods. We focused on pathways involving hypoxia, angiogenesis, and epithelial mesenchymal transition according to the gene set variation analysis (GSVA) scores. The gene sets of these three pathways were enriched by gene set enrichment analysis (GSEA). WGCNA was used to study the underlying molecular mechanisms of the three pathways in the pathogenesis of PE by analyzing the relationship among pathways and genes. The soft threshold power (β) and topological overlap matrix allowed us to obtain 15 modules, among which the red module was chosen for the downstream analysis. We chose 10 hub genes that satisfied ∣ log 2 Fold Change ∣ > 2 and had a higher degree of connectivity within the module. These candidate genes were subsequently confirmed to have higher gene significance and module membership in the red module. Coexpression networks were established for the hub genes to unfold the connection between the genes in the red module and PE. Finally, ceRNA networks were constructed to further clarify the underlying molecular mechanism involved in the occurrence of PE. 56 circRNAs, 17 lncRNAs, and 20 miRNAs participated in the regulation of the hub genes. Coagulation factor II thrombin receptor (F2R) and lumican (LUM) were considered the most relevant genes, and ceRNA networks of them were constructed. Conclusion. The microarray data mining process based on bioinformatics methods constructed lncRNA and miRNA networks for ten hub genes that were closely related to PE and focused on ceRNAs of F2R and LUM finally. The results of our study may provide insight into the mechanisms underlying PE occurrence.

Integrated weighted gene coexpression network analysis identifies Frizzled 2 (FZD2) as a key gene in invasive malignant pleomorphic adenoma

Background Invasive malignant pleomorphic adenoma (IMPA) is a highly malignant neoplasm of the oral salivary glands with a poor prognosis and a considerable risk of recurrence. Many disease-causing genes of IMPA have been identified in recent decades (e.g., P53, PCNA and HMGA2), but many of these genes remain to be explored. Weighted gene coexpression network analysis (WGCNA) is a newly emerged algorithm that can cluster genes and form modules based on similar gene expression patterns. This study constructed a gene coexpression network of IMPA via WGCNA and then carried out multifaceted analysis to identify novel disease-causing genes. Methods RNA sequencing (RNA-seq) was performed for 10 pairs of IMPA and normal tissues to acquire the gene expression profiles. Differentially expressed genes (DEGs) were screened out with the cutoff criteria of |log2 Fold change (FC)|> 1 and adjusted p value < 0.05. Then, WGCNA was applied to systematically identify the hidden diagnostic hub genes of IMPA. Results In this research, a total of 1970 DEGs were screened out in IMPA tissues, including 1056 upregulated DEGs and 914 downregulated DEGs. Functional enrichment analysis was performed for identified DEGs and revealed an enrichment of tumor-associated GO terms and KEGG pathways. We used WGCNA to identify gene module most relevant with the histological grade of IMPA. The gene FZD2 was then recognized as the hub gene of the selected module with the highest module membership (MM) value and intramodule connectivity in protein–protein interaction (PPI) network. According to immunohistochemistry (IHC) staining, the expression level of FZD2 was higher in low-grade IMPA than in high-grade IMPA. Conclusion FZD2 shows an expression dynamic that is negatively correlated with the clinical malignancy of IMPA and it plays a central role in the transcription network of IMPA. Thus, FZD2 serves as a promising histological indicator for the precise prediction of IMPA histological stages.

Weighted Gene Coexpression Network Analysis in Mouse Livers following Ischemia-Reperfusion and Extensive Hepatectomy

In mouse models, the recovery of liver volume is mainly mediated by the proliferation of hepatocytes after partial hepatectomy that is commonly accompanied with ischemia-reperfusion. The identification of differently expressed genes in liver following partial hepatectomy benefits the better understanding of the molecular mechanisms during liver regeneration (LR) with appliable clinical significance. Briefly, studying different gene expression patterns in liver tissues collected from the mice group that survived through extensive hepatectomy will be of huge critical importance in LR than those collected from the mice group that survived through appropriate hepatectomy. In this study, we performed the weighted gene coexpression network analysis (WGCNA) to address the central candidate genes and to construct the free-scale gene coexpression networks using the identified dynamic different expressive genes in liver specimens from the mice with 85% hepatectomy (20% for seven-day survial rate) and 50% hepatectomy (100% for seven-day survial rate under ischemia-reperfusion condition compared with the sham group control mice). The WGCNA combined with Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analyses pinpointed out the apparent distinguished importance of three gene expression modules: the blue module for apoptotic process, the turquoise module for lipid metabolism, and the green module for fatty acid metabolic process in LR following extensive hepatectomy. WGCNA analysis and protein-protein interaction (PPI) network construction highlighted FAM175B, OGT, and PDE3B were the potential three hub genes in the previously mentioned three modules. This work may help to provide new clues to the future fundamental study and treatment strategy for LR following liver injury and hepatectomy.

Identification of Differentially Expressed and Prognostic lncRNAs for the Construction of ceRNA Networks in Lung Adenocarcinoma

Background. Long noncoding RNAs (lncRNAs) could function as competitive endogenous RNAs (ceRNAs) to competitively adsorb microRNAs (miRNAs), thereby regulating the expression of their target protein-coding mRNAs. In this study, we aim to identify more effective diagnostic and prognostic markers for lung adenocarcinoma (LUAD). Methods. We obtained differentially expressed lncRNAs (DElncRNAs), miRNAs (DEmiRNAs), and mRNAs (DEmRNAs) for LUAD by using The Cancer Genomes Atlas (TCGA) portal. Weighted gene coexpression network analysis (WGCNA) was performed to unveil core gene modules associated with LUAD. The Cox proportional hazards model was performed to determine the prognostic significance of DElncRNAs. The diagnostic and prognostic significance of DElncRNAs was further verified based on the receiver operating characteristic curve (ROC). Cytoscape was used to construct the ceRNA networks comprising the lncRNAs-miRNAs-mRNAs axis based on the correlation obtained from the miRcode, miRDB, and TargetScan. Results. Compared with normal lung tissues, 2355 DElncRNAs, 820 DEmiRNAs, and 17289 DEmRNAs were identified in LUAD tissues. We generated 8 WGCNA core modules in the lncRNAs coexpression network, 5 modules in the miRNAs, and 12 modules in the mRNAs coexpression network, respectively. One lncRNA module (blue) consisting of 441 lncRNAs, two miRNA modules (blue and turquoise) containing 563 miRNAs, and one mRNA module (turquoise), which consisted of 15162 mRNAs, were mostly significantly related to LUAD status. Furthermore, 67 DEmRNAs were found to be tumor-associated as well as the target genes of the DElncRNAs-DEmiRNAs axis. Survival analyses showed that 6 lncRNAs (LINC01447, WWC2-AS2, OGFRP1, LINC00942, LINC01168, and AC005863.1) were significantly correlated with the prognosis of LUAD patients. Ultimately, the potential ceRNA networks including 6 DElncRNAs, 4 DEmiRNAs, and 22 DEmRNAs were constructed. Conclusion. Our study indicated that 6 DElncRNAs had the possibilities as diagnostic and prognostic biomarkers for LUAD. The lncRNA-mediated ceRNA networks might provide novel insights into the molecular mechanisms of LUAD progression.

Identification of four hub genes in venous thromboembolism via weighted gene coexpression network analysis

Background The pathogenic mechanisms of venous thromboembolism (VT) remain to be defined. This study aimed to identify differentially expressed genes (DEGs) that could serve as potential therapeutic targets for VT. Methods Two human datasets (GSE19151 and GSE48000) were analyzed by the robust rank aggregation method. Gene ontology and Kyoto encyclopedia of genes and genomes pathway enrichment analyses were conducted for the DEGs. To explore potential correlations between gene sets and clinical features and to identify hub genes, we utilized weighted gene coexpression network analysis (WGCNA) to build gene coexpression networks incorporating the DEGs. Then, the levels of the hub genes were analyzed in the GSE datasets. Based on the expression of the hub genes, the possible pathways were explored by gene set enrichment analysis and gene set variation analysis. Finally, the diagnostic value of the hub genes was assessed by receiver operating characteristic (ROC) analysis in the GEO database. Results In this study, we identified 54 upregulated and 10 downregulated genes that overlapped between normal and VT samples. After performing WGCNA, the magenta module was the module with the strongest negative correlation with the clinical characteristics. From the key module, FECH, GYPA, RPIA and XK were chosen for further validation. We found that these genes were upregulated in VT samples, and high expression levels were related to recurrent VT. Additionally, the four hub genes might be highly correlated with ribosomal and metabolic pathways. The ROC curves suggested a diagnostic value of the four genes for VT. Conclusions These results indicated that FECH, GYPA, RPIA and XK could be used as promising biomarkers for the prognosis and prediction of VT.

Identification of crucial genes involved in pathogenesis of regional weakening of the aortic wall

Background The diameter of the abdominal aortic aneurysm (AAA) is the most commonly used parameter for the prediction of occurrence of AAA rupture. However, the most vulnerable region of the aortic wall may be different from the most dilated region of AAA under pressure. The present study is the first to use weighted gene coexpression network analysis (WGCNA) to detect the coexpressed genes that result in regional weakening of the aortic wall. Methods The GSE165470 raw microarray dataset was used in the present study. Differentially expressed genes (DEGs) were filtered using the “limma” R package. DEGs were assessed by Gene Ontology biological process (GO-BP) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. WGCNA was used to construct the coexpression networks in the samples with regional weakening of the AAA wall and in the control group to detect the gene modules. The hub genes were defined in the significant functional modules, and a hub differentially expressed gene (hDEG) coexpression network was constructed with the highest confidence based on protein–protein interactions (PPIs). Molecular compound detection (MCODE) was used to identify crucial genes in the hDEG coexpression network. Crucial genes in the hDEG coexpression network were validated using the GSE7084 and GSE57691 microarray gene expression datasets. Result A total of 350 DEGs were identified, including 62 upregulated and 288 downregulated DEGs. The pathways were involved in immune responses, vascular smooth muscle contraction and cell–matrix adhesion of DEGs in the samples with regional weakening in AAA. Antiquewhite3 was the most significant module and was used to identify downregulated hDEGs based on the result of the most significant modules negatively related to the trait of weakened aneurysm walls. Seven crucial genes were identified and validated: ACTG2, CALD1, LMOD1, MYH11, MYL9, MYLK, and TPM2. These crucial genes were associated with the mechanisms of AAA progression. Conclusion We identified crucial genes that may play a significant role in weakening of the AAA wall and may be potential targets for medical therapies and diagnostic biomarkers. Further studies are required to more comprehensively elucidate the functions of crucial genes in the pathogenesis of regional weakening in AAA.

Global Transcriptome and Coexpression Network Analyses Reveal New Insights Into Somatic Embryogenesis in Hybrid Sweetgum (Liquidambar styraciflua × Liquidambar formosana)

Somatic embryogenesis (SE) is a process of somatic cells that dedifferentiate to totipotent embryonic stem cells and generate embryos in vitro. Despite recent scientific headway in deciphering the difficulties of somatic embryogenesis, the overall picture of key genes, pathways, and co-expression networks regulating SE is still fragmented. Therefore, deciphering the molecular basis of somatic embryogenesis of hybrid sweetgum remains pertinent. In the present study, we analyzed the transcriptome profiles and gene expression regulation changes via RNA sequencing from three distinct developmental stages of hybrid sweetgum: non-embryogenic callus (NEC), embryogenic callus (EC), and redifferentiation. Comparative transcriptome analysis showed that 19,957 genes were differentially expressed in ten pairwise comparisons of SE. Among these, plant hormone signaling-related genes, especially the auxin and cytokinin signaling components, were significantly enriched in NEC and EC early. The K-means method was used to identify multiple transcription factors, including HB-WOX, B3-ARF, AP2/ERF, and GRFs (growth regulating factors). These transcription factors showed distinct stage- or tissue-specific expression patterns mirroring each of the 12 superclusters to which they belonged. For example, the WOX transcription factor family was expressed only at NEC and EC stages, ARF transcription factor was expressed in EC early, and GRFs was expressed in late SE. It was noteworthy that the AP2/ERF transcription factor family was expressed during the whole SE process, but almost not in roots, stems and leaves. A weighted gene co-expression network analysis (WGCNA) was used in conjunction with the gene expression profiles to recognize the genes and modules that may associate with specific tissues and stages. We constructed co-expression networks and revealed 22 gene modules. Four of these modules with properties relating to embryonic potential, early somatic embryogenesis, and somatic embryo development, as well as some hub genes, were identified for further functional studied. Through a combination analysis of WGCNA and K-means, SE-related genes including AUX22, ABI3, ARF3, ARF5, AIL1, AIL5, AGL15, WOX11, WOX9, IAA29, BBM1, MYB36, LEA6, SMR4 and others were obtained, indicating that these genes play an important role in the processes underlying the progression from EC to somatic embryos (SEs) morphogenesis. The transcriptome information provided here will form the foundation for future research on genetic transformation and epigenetic control of plant embryogenesis at a molecular level. In follow-up studies, these data could be used to construct a regulatory network for SE; Key genes obtained from coexpression network analysis at each critical stage of somatic embryo can be considered as potential candidate genes to verify these networks.

Identification of Hub Gene GRIN1 Correlated with Histological Grade and Prognosis of Glioma by Weighted Gene Coexpression Network Analysis

The function of glutamate ionotropic receptor NMDA type subunit 1 (GRIN1) in neurodegenerative diseases has been widely reported; however, its role in the occurrence of glioma remains less explored. We obtained clinical data and transcriptome data from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA). Hub gene’s expression differential analysis and survival analysis were conducted by browsing the Gene Expression Profiling Interactive Analysis (GEPIA) database, Human Protein Atlas database, and LOGpc database. We conducted a variation analysis of datasets obtained from GEO and TCGA and performed a weighted gene coexpression network analysis (WGCNA) using the R programming language (3.6.3). Kaplan-Meier (KM) analysis was used to calculate the prognostic value of GRIN1. Finally, we conducted Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Using STRING, we constructed a protein–protein interaction (PPI) network. Cytoscape software, a prerequisite of visualizing core genes, was installed, and CytoHubba detected the 10 most tumor-related core genes. We identified 185 differentially expressed genes (DEGs). GO and KEGG enrichment analyses illustrated that the identified DEGs are imperative in different biological functions and ascertained the potential pathways in which the DEGs may be enriched. The overall survival calculated by KM analysis showed that patients with lower expression of GRIN1 had worse prognoses than patients with higher expression of GRIN1 ( p = 0.004 ). The GEPIA and LOGpc databases were used to verify the expression difference of GRIN1 among GBM, LGG, and normal brain tissues. Ultimately, immunohistochemical assay results showed that GRIN1 was detected in normal tissue and not in the tumor specimens. Our results highlight a potential target for glioma treatment and will further our understanding of the molecular mechanisms underlying the treatment of glioma.