Identification Of Putative Gene Signatures Associated With Diagnosis And Prognosis Of Breast Cancer

Purpose: This study aimed to identify potential diagnostic and therapeutic biomakers for the development ofbreast cancer (BC). Methods: GSE86374 dataset containing 159 samples was acquired from the Gene Expression Omnibus (GEO) database followed by differentially expressed genes (DEGs) identification and cluster analysis. Corresponding functional enrichment and protein-protein interaction (PPI) network analyses were performed to identify hub genes. Prognostic evaluation using clinical information obtained from TCGA database and hub genes was conducted to screen for crucial indicators for BC progression. The risk model was established and validated. Results: In total, 186 DEGs were identified and grouped into four clusters: 96 in cluster 1; 69 in cluster 2; 16 in cluster 3; and 5 in cluster 4. Functional enrichment analysis showed that DEGs, including ADH1B in cluster 1, were dramatically enriched in the tyrosine and drug metabolism pathways, while genes in cluster 2, including SPP1 and RRM2, played crucial roles in PI3K-Akt and p53 signalling pathway. SPP1 and RRM2 served as hub genes in the PPI network, resulting in an support vector machine classifier with good accuracy and specificity.Ad ditionally, the results of prognostic analysis suggest that age, metastasis stage, SPP1 and ADH1B were correlated with risk of BC, which was validated by using the established risk model analysis. Conclusion: SPP1, RRM2 and ADH1B appear to play vital roles in the development of BC. Age and TNM stage were also preferentially associated with risk of developing BC. Evaluation of the risk model based on larger sample size and further experimental validation are required.

Download Full-text

Identifcation of The Ferroptosis-Related Gene Signature In Placenta of Patients With Early-Onset Preeclampsia

10.21203/rs.3.rs-618453/v1 ◽

2021 ◽

Author(s):

Nana Yang ◽

Qianghua Wang ◽

Biao Ding ◽

Yinging Gong ◽

Yue Wu ◽

...

Keyword(s):

Iron Homeostasis ◽

Molecular Mechanisms ◽

Late Onset ◽

Expression Profiles ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Differentially Expressed ◽

Ppi Network ◽

Network Analyses

Abstract Background: The accumulation of ROS resulting from upregulated levels of oxidative stress is commonly implicated in preeclampsia (PE). Ferroptosis is a novel form of iron-dependent cell death instigated by lipid peroxidation likely plays important role in PE pathogenesis. This study aims to investigate expression profiles and functions of the ferroptosis-related genes (FRGs) in early- and late-onset preeclampsia.Methods: The gene expression data and clinical information were downloaded from GEO database. The “limma” R package was used for screening differentially expressed genes. GO(Gene Ontology), Kyoto Encyclopedia of Genes and Genomes(KEGG) and protein protein interaction (PPI) network analyses were conducted to investigate the bioinformatics functions and molecular interactions of significantly different FRGs. Quantitative real-time reverse transcriptase PCR was used to verify the expression of hub FRGs in PE.Results: A total number of 4,215 DEGs were identified between EOPE and preterm cases and 3,356 DEGs were found between EOPE and LOPE subtypes. 20 significantly different FRGs were identified in EOPE, while only 3 in LOPE. Functional enrichment analysis revealed that the differentially expressed FRGs was mainly involved in EOPE and enriched in hypoxia- and iron-related pathways, such as response to hypoxia, iron homeostasis and iron ion binding process. The PPI network analysis and verification by RT-qPCR resulted in the identification of the following six interesting FRGs: FTH1, HIF1A, FTL, IREB2, MAPK8 and PLIN2. Conclusions: EOPE and LOPE owned distinct underlying molecular mechanisms and ferroptosis may be mainly implicated in pathogenesis of EOPE. Further studies are necessary for deeper inquiry into placental ferroptosis and its role in the pathogenesis of EOPE.

Download Full-text

Xiaoyu Xiezhuo Drink Protects against Ischemia-Reperfusion Acute Kidney Injury in Aged Mice through Inhibiting the TGF-β1/Smad3 and HIF1 Signaling Pathways

BioMed Research International ◽

10.1155/2021/9963732 ◽

2021 ◽

Vol 2021 ◽

pp. 1-19

Author(s):

Qingqing Ye ◽

Hongbo Chen ◽

Hongzhen Ma ◽

Xiaojun Xiang ◽

Shouci Hu ◽

...

Keyword(s):

Drug Targets ◽

Ischemia Reperfusion ◽

Kidney Injury ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Ppi Network ◽

Hub Genes ◽

Aged Mice ◽

Tgf Β1

Acute kidney injury (AKI) is responsible for significant mortality among hospitalized patients that is especially troubling aged people. An effective self-made Chinese medicine formula, Xiaoyu Xiezhuo Drink (XXD), displayed therapeutic effects on AKI. However, the compositions and underlying mechanisms of XXD remain to be elucidated. In this study, we used the ultra-high-performance liquid chromatography method coupled with hybrid triple quadrupole time-of-flight mass spectrometry (UHPLC-Q-TOF-MS) to investigate the chemical components in XXD. Then, the absorbable components of XXD were identified based on the five principles and inputted into the SwissTargetPrediction and STITCH databases to identify the drug targets. AKI-related targets were collected from the GenCLiP 3, GeneCards, and DisGeNET databases. The crossover genes of XXD and AKI were identified for functional enrichment analysis. The protein-protein interaction (PPI) network of crossover genes was constructed, followed by the identification of hub genes. Subsequently, the effects and potential mechanisms of XXD on AKI predicted by the network pharmacology and bioinformatics analyses were experimentally validated in ischemia-reperfusion (I/R) injury-induced AKI aged mouse models. A total of 122 components in XXD were obtained; among them, 58 components were found that could be absorbed in the blood. There were 800 potential drug targets predicted from the 58 absorbable components in AKI which shared 36 crossover genes with AKI-related targets. The results of functional enrichment analysis indicated that crossover genes mostly associated with the response to oxidative stress and the HIF1 signaling pathway. In the PPI network analysis, 12 hub genes were identified, including ALB, IL-6, TNF, TP53, VEGFA, PTGS2, TLR4, NOS3, EGFR, PPARG, HIF1A, and HMOX1. In AKI aged mice, XXD prominently alleviated I/R injury-induced renal dysfunction, abnormal renal pathological changes, and cellular senescence, inflammation, and oxidative damage with a reduction in the expression level of the inflammatory mediator, α-SMA, collagen-1, F4/80, TP53, VEGFA, PTGS2, TLR4, NOS3, EGFR, PPARG, HIF1A, ICAM-1, TGF-β1, Smad3, and p-Smad3 and an increase of nephridial tissue p-H3, Ki67, HMOX1, MMP-9, and Smad7 levels. In summary, our findings suggest that XXD has renoprotective effects against AKI in aged mice via inhibiting the TGF-β1/Smad3 and HIF1 signaling pathways.

Download Full-text

The Important Role of Perituberal Tissue in Epileptic Patients with Tuberous Sclerosis Complex by the Transcriptome Analysis

BioMed Research International ◽

10.1155/2020/4980609 ◽

2020 ◽

Vol 2020 ◽

pp. 1-11

Author(s):

Shuqiang Li ◽

Huijie Shao ◽

Liansheng Chang

Keyword(s):

Tuberous Sclerosis ◽

Tuberous Sclerosis Complex ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Receptor Interaction ◽

Ppi Network ◽

Hub Genes ◽

R Language ◽

Protein Protein Interaction

Epilepsy is most common in patients with tuberous sclerosis complex (TSC). However, in addition to the challenging treatment, the pathogenesis of epilepsy is still controversial. To determine the transcriptome characteristics of perituberal tissue (PT) and clarify its role in the pathogenesis of epilepsy, GSE16969 was downloaded from the GEO database for further study by comprehensive bioinformatics analysis. Identification of differentially expressed genes (DEGs), functional enrichment analysis, construction of protein-protein interaction (PPI) network, and selection of Hub genes were performed using R language, Metascape, STRING, and Cytoscape, respectively. Comparing with cortical tuber (CT), 220 DEGs, including 95 upregulated and 125 downregulated genes, were identified in PT and mainly enriched in collagen-containing extracellular matrix and positive regulation of receptor-mediated endocytosis, as well as the pathways of ECM-receptor interaction and neuroactive ligand-receptor interaction. As for normal cortex (NC), 1549 DEGs, including 30 upregulated and 1519 downregulated genes, were identified and mainly enriched in presynapse, dendrite and axon, and also the pathways of dopaminergic synapse and oxytocin signaling pathway. In the PPI network, 4 hub modules were found between PT and CT, and top 5 hub modules were selected between PT and NC. C3, APLNR, ANXA2, CD44, CLU, CP, MCHR2, HTR1E, CTSG, APP, and GNG2 were identified as Hub genes, of which, C3, CD44, ANXA2, HTR1E, and APP were identified as Hub-BottleNeck genes. In conclusion, PT has the unique characteristics different from CT and NC in transcriptome and makes us further understand its importance in the TSC-associated epilepsy.

Download Full-text

Identification of the Hub Genes in Alzheimer’s Disease

Computational and Mathematical Methods in Medicine ◽

10.1155/2021/6329041 ◽

2021 ◽

Vol 2021 ◽

pp. 1-8

Author(s):

Huiwen Gui ◽

Qi Gong ◽

Jun Jiang ◽

Mei Liu ◽

Huanyin Li

Keyword(s):

Alzheimer’S Disease ◽

Alzheimer's Disease ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

The Elderly ◽

Functional Enrichment ◽

Ppi Network ◽

Hub Genes ◽

Protein Protein Interaction ◽

Geo Database

Purpose. Alzheimer’s disease (AD) is considered to be the most common neurodegenerative disease and also one of the major fatal diseases affecting the elderly, thus bringing a huge burden to society. Therefore, identifying AD-related hub genes is extremely important for developing novel strategies against AD. Materials and Methods. Here, we extracted the gene expression profile GSE63061 from the National Center for Biotechnology Information (NCBI) GEO database. Once the unverified gene chip was removed, we standardized the microarray data after quality control. We utilized the Limma software package to screen the differentially expressed genes (DEGs). We conducted Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of DEGs. Subsequently, we constructed a protein-protein interaction (PPI) network using the STRING database. Result. We screened 2169 DEGs, comprising 1313 DEGs with upregulation and 856 DEGs with downregulation. Functional enrichment analysis showed that the response of immune, the degranulation of neutrophils, lysosome, and the differentiation of osteoclast were greatly enriched in DEGs with upregulation; peptide biosynthetic process, translation, ribosome, and oxidative phosphorylation were dramatically enriched in DEGs with downregulation. 379 nodes and 1149 PPI edges were demonstrated in the PPI network constructed by upregulated DEGs; 202 nodes and 1963 PPI edges were shown in the PPI network constructed by downregulated DEGs. Four hub genes, including GAPDH, RHOA, RPS29, and RPS27A, were identified to be the newly produced candidates involved in AD pathology. Conclusion. GAPDH, RHOA, RPS29, and RPS27A are expected to be key candidates for AD progression. The results of this study can provide comprehensive insight into understanding AD’s pathogenesis and potential new therapeutic targets.

Download Full-text

Single-cell RNA-sequencing Unravels Heterogeneity of Cardiomyocytes and Signaling Pathways in Pressure Overload

10.21203/rs.3.rs-37034/v1 ◽

2020 ◽

Author(s):

Rongjun Zou ◽

Wanting Shi ◽

Minghui Zou ◽

Weidan Chen ◽

Wenlei Li ◽

...

Keyword(s):

Cardiac Hypertrophy ◽

Mouse Models ◽

Expression Patterns ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Ppi Network ◽

Hub Genes ◽

Cell Clusters ◽

Single Cell Rna Sequencing

Abstract Background This study aimed to unravel the heterogeneity of cardiomyocytes and probed out hub genes and hub pathways for cardiac hypertrophy based on transverse aortic constriction (TAC) mouse models using single-cell RNA sequencing (scRNA-seq). Methods scRNA-seq data of TAC mouse models were retrieved from the GSE95140 dataset. After filtering, cell clusters were detected using scRNA-seq data, followed by identification of differentially expressed genes (DEGs). Then, functional enrichment analysis of DEGs was presented. GSVA scores of hub pathways were calculated. After that, hub genes were detected by protein-protein interaction (PPI) network and expression association analysis. Cell subtypes were clustered using UMAP and the expression patterns of hub genes across different cell subtypes and different stages of cardiac hypertrophy were visualized. Finally, hub genes and hub pathways were verified using the GSE76 and GSE36074 datasets. Results Following data filtering and normalization, 3408 DEGs were identified between TAC and sham operation. As shown functional enrichment analysis, hub pathways were identified including cardiac hypertrophy, ion transport, myocardial remodeling, apoptosis, HIF pathway and metabolise. Eight hub genes (Vldlr, Ugp2, Tgm2, Pygm, Flnc, Ctsd, Clu and Atp1b1) with the highest degree in the PPI network and the strongest correlation with GSVA calculated score of hub pathways were identified for cardiac hypertrophy. Six cell subtypes were clustered, composed of fibroblast, CM-A, CM-V, trabecular CM and endothelial cell. There was a distinct heterogeneity in the expression patterns of hub genes and the GSVA scores of hub pathways across different cell clusters and different stages of cardiac hypertrophy. The hub genes and hub pathways were externally verified by the two independent datasets. Conclusion Our findings identified hub genes and hub pathways for cardiac hypertrophy, which had a distinct heterogeneity across different cell clusters and different stages of cardiac hypertrophy.

Download Full-text

Bioinformatics analysis of microarray data to identify the candidate biomarkers of lung adenocarcinoma

PeerJ ◽

10.7717/peerj.7313 ◽

2019 ◽

Vol 7 ◽

pp. e7313 ◽

Cited By ~ 6

Author(s):

Tingting Guo ◽

Hongtao Ma ◽

Yubai Zhou

Keyword(s):

Survival Analysis ◽

Lung Adenocarcinoma ◽

Molecular Mechanisms ◽

Prion Diseases ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Gene Expression Omnibus ◽

Functional Enrichment ◽

Ppi Network ◽

Hub Genes

Background Lung adenocarcinoma (LUAD) is the major subtype of lung cancer and the most lethal malignant disease worldwide. However, the molecular mechanisms underlying LUAD are not fully understood. Methods Four datasets (GSE118370, GSE85841, GSE43458 and GSE32863) were obtained from the gene expression omnibus (GEO). Identification of differentially expressed genes (DEGs) and functional enrichment analysis were performed using the limma and clusterProfiler packages, respectively. A protein–protein interaction (PPI) network was constructed via Search Tool for the Retrieval of Interacting Genes (STRING) database, and the module analysis was performed by Cytoscape. Then, overall survival analysis was performed using the Kaplan–Meier curve, and prognostic candidate biomarkers were further analyzed using the Oncomine database. Results Totally, 349 DEGs were identified, including 275 downregulated and 74 upregulated genes which were significantly enriched in the biological process of extracellular structure organization, leukocyte migration and response to peptide. The mainly enriched pathways were complement and coagulation cascades, malaria and prion diseases. By extracting key modules from the PPI network, 11 hub genes were screened out. Survival analysis showed that except VSIG4, other hub genes may be involved in the development of LUAD, in which MYH10, METTL7A, FCER1G and TMOD1 have not been reported previously to correlated with LUAD. Briefly, novel hub genes identified in this study will help to deepen our understanding of the molecular mechanisms of LUAD carcinogenesis and progression, and to discover candidate targets for early detection and treatment of LUAD.

Download Full-text

Identification of key gene in colorectal cancer using bioinformatics analysis.

10.21203/rs.3.rs-63457/v1 ◽

2020 ◽

Author(s):

Rongrong Xiao ◽

Ping Wang ◽

Tian Xia ◽

Chun-Yi Li ◽

Ting Jiang ◽

...

Keyword(s):

Colorectal Cancer ◽

Colon Cancer ◽

Tumor Microenvironment ◽

Differentially Expressed Genes ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Differentially Expressed ◽

Ppi Network ◽

Hub Genes

Abstract Background Tumor microenvironment plays important roles in the development of cancer. The aim of our study was to examine the expression of genes in colorectal cancer and also to evaluate the association value between expression level of these genes and clinical features. Methods We combined The Cancer Genome Atlas (TCGA) datasets to identify differentially expressed genes in colon cancer. Using these differentially expressed genes, we constructed protein-protein interaction network and conducted functional enrichment analysis. Genes with degree beyond 10 in the PPI network were regarded as hub genes. Then, we verified of the expression of molecules in Oncomine datasets and conducted Kaplan-Meier curve and log-rank test and functional enrichment analysis on these hub genes. Finally, we analyzed the relationship clinicopathological features analysis with the key gene. Results There were 719 differentially expressed genes identified to be associated with colon cancer microenvironment. We screened out 10 hub genes by construction of PPI network. The functions of these hub genes were enriched in cytokine-cytokine receptor interaction, alcoholism and systemic lupus erythematosus, which provided further insight into the roles of these genes in the tumor microenvironment. GNG4, with the highest degrees in the PPI network, were highly exprepressed in metastasis(P = 9.5-05) ,N1(P = 0.0025) and N2(,0.037).It was a relationship with stage. It was significantly different between with stage I and IV, II and III, II and IV,III and IV (P = 0.0015,0.029,3.9-05,0.00074,0.01,respectively) Conclusions We identified GNG4 can be regarded as a prognostic biomarker in colon cancer.

Download Full-text

Identification of miRNA-target gene regulatory networks in liver fibrosis based on bioinformatics analysis

PeerJ ◽

10.7717/peerj.11910 ◽

2021 ◽

Vol 9 ◽

pp. e11910

Author(s):

Yang Tai ◽

Chong Zhao ◽

Jinhang Gao ◽

Tian Lan ◽

Huan Tong

Keyword(s):

Immune Response ◽

Liver Cirrhosis ◽

Bioinformatics Analysis ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Ppi Network ◽

Hub Genes ◽

Differentially Expressed Mirnas ◽

And Control

Background Liver cirrhosis is one of the leading causes of death worldwide. MicroRNAs (miRNAs) can regulate liver fibrosis, but the underlying mechanisms are not fully understood, and the interactions between miRNAs and mRNAs are not clearly elucidated. Methods miRNA and mRNA expression arrays of cirrhotic samples and control samples were acquired from the Gene Expression Omnibus database. miRNA-mRNA integrated analysis, functional enrichment analysis and protein-protein interaction (PPI) network construction were performed to identify differentially expressed miRNAs (DEMs) and mRNAs (DEGs), miRNA-mRNA interaction networks, enriched pathways and hub genes. Finally, the results were validated with in vitro cell models. Results By bioinformatics analysis, we identified 13 DEMs between cirrhotic samples and control samples. Among these DEMs, six upregulated (hsa-miR-146b-5p, hsa-miR-150-5p, hsa-miR-224-3p, hsa-miR-3135b, hsa-miR-3195, and hsa-miR-4725-3p) and seven downregulated (hsa-miR-1234-3p, hsa-miR-30b-3p, hsa-miR-3162-3p, hsa-miR-548aj-3p, hsa-miR-548x-3p, hsa-miR-548z, and hsa-miR-890) miRNAs were further validated in activated LX2 cells. miRNA-mRNA interaction networks revealed a total of 361 miRNA-mRNA pairs between 13 miRNAs and 245 corresponding target genes. Moreover, PPI network analysis revealed the top 20 hub genes, including COL1A1, FBN1 and TIMP3, which were involved in extracellular matrix (ECM) organization; CCL5, CXCL9, CXCL12, LCK and CD24, which participated in the immune response; and CDH1, PECAM1, SELL and CAV1, which regulated cell adhesion. Functional enrichment analysis of all DEGs as well as hub genes showed similar results, as ECM-associated pathways, cell surface interaction and adhesion, and immune response were significantly enriched in both analyses. Conclusions We identified 13 differentially expressed miRNAs as potential biomarkers of liver cirrhosis. Moreover, we identified 361 regulatory pairs of miRNA-mRNA and 20 hub genes in liver cirrhosis, most of which were involved in collagen and ECM components, immune response, and cell adhesion. These results would provide novel mechanistic insights into the pathogenesis of liver cirrhosis and identify candidate targets for its treatment.

Download Full-text

Construction of Gene Modules and Analysis of Prognostic Biomarkers for Cervical Cancer by Weighted Gene Co-Expression Network Analysis

Frontiers in Oncology ◽

10.3389/fonc.2021.542063 ◽

2021 ◽

Vol 11 ◽

Author(s):

Jiamei Liu ◽

Shengye Liu ◽

Xianghong Yang

Keyword(s):

Gene Expression ◽

Cervical Cancer ◽

Network Analysis ◽

Cancer Progression ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Functional Enrichment ◽

Ppi Network ◽

Hub Genes ◽

Gene Modules

BackgroundDespite advances in the understanding of neoplasm, patients with cervical cancer still have a poor prognosis. Identifying prognostic markers of cervical cancer may enable early detection of recurrence and more effective treatment.MethodsGene expression profiling data were acquired from the Gene Expression Omnibus database. After data normalization, genes with large variation were screened out. Next, we built co-expression modules by using weighted gene co-expression network analysis to investigate the relationship between the modules and clinical traits related to cervical cancer progression. Functional enrichment analysis was also applied on these co-expressed genes. We integrated the genes into a human protein-protein interaction (PPI) network to expand seed genes and build a co-expression network. For further analysis of the dataset, the Cancer Genome Atlas (TCGA) database was used to identify seed genes and their correlation to cervical cancer prognosis. Verification was further conducted by qPCR and the Human Protein Atlas (HPA) database to measure the expression of hub genes.ResultsUsing WGCNA, we identified 25 co-expression modules from 10,016 genes in 128 human cervical cancer samples. After functional enrichment analysis, the magenta, brown, and darkred modules were selected as the three most correlated modules for cancer progression. Additionally, seed genes in the three modules were combined with a PPI network to identify 31 tumor-specific genes. Hierarchical clustering and Gepia results indicated that the expression quantity of hub genes NDC80, TIPIN, MCM3, MCM6, POLA1, and PRC1 may determine the prognosis of cervical cancer. Finally, TIPIN and POLA1 were further filtered by a LASSO model. In addition, their expression was identified by immunohistochemistry in HPA database as well as a biological experiment.ConclusionOur research provides a co-expression network of gene modules and identifies TIPIN and POLA1 as stable potential prognostic biomarkers for cervical cancer.

Download Full-text

Underlying Genes and Molecular Mechanism of Keloids Investigated by Integrated Bioinformatics Analysis

10.21203/rs.3.rs-148691/v1 ◽

2021 ◽

Author(s):

chuying LI ◽

Mei-Tong Jin ◽

Yin-Li Luo ◽

Zhe-Hu Jin ◽

Long-Quan PI

Keyword(s):

Expression Profiles ◽

Normal Skin ◽

Enrichment Analysis ◽

Functional Enrichment Analysis ◽

Underlying Mechanism ◽

Cholinergic Receptor ◽

Gene Expression Omnibus ◽

Functional Enrichment ◽

Ppi Network ◽

Hub Genes

Abstract Background: We aimed to identify the overlapping differentially expressed genes (DEGs) of keloids distinguished from normal scar and normal skin and relevant underlying mechanism using integrated bioinformatics methods.Methods: The expression profiles of 18 keloid samples, 7 normal skin and 5 normal scar, were obtained from the GSE7890, GSE44270, GSE92566, and GSE3189 datasets in the Gene Expression Omnibus database. DEGs were identified using the LIMMA package in R. Gene ontology (GO) functional enrichment analysis was performed using the R software. A DEG-associated protein–protein interaction (PPI) network was constructed using STRING and MCODE was used for module analysis of the PPI network. Moreover, the hub genes were verified by qRT-PCR. The predicted DEGs, their regulatory miRNA and TF regulation network was analyzed using miRnet. Results: A total of 978 common DEGs were identified in the keloid samples. Genes with more than 45 interaction degrees, including neuropeptide Y (NPY), opioid receptor mu 1 (OPRM1), cholinergic receptor muscarinic 2 (CHRM2), and proopiomelanocortin (POMC), were found in the PPI network. Hsa-miR-335 and Sp1 as upstream-regulators regulated CHRM2, NPY, and POMC. Functional enrichment analysis revealed that hub genes were commonly enriched in the “G protein-coupled receptor signaling pathway” GO_BP termConclusion: Taken together, CHRM2, NPY, POMC, and OPRM1 potentially have crucial roles in keloid disease. Furthermore, miR-335 and Sp1 are potential targets for preventing keloid formation.

Download Full-text