irGSEA: the integration of single-cell rank-based gene set enrichment analysis.

irGSEA is an R package designed to assess the outcomes of various gene set scoring methods when applied to single-cell RNA sequencing data. This package incorporates six distinct scoring methods that rely on the expression ranks of genes, emphasizing relative expression levels over absolute values. The implemented methods include AUCell, UCell, singscore, ssGSEA, JASMINE and Viper. Previous studies have demonstrated the robustness of these methods to variations in dataset size and composition, generating enrichment scores based solely on the relative gene expression of individual cells. By employing the robust rank aggregation algorithm, irGSEA amalgamates results from all six methods to ascertain the statistical significance of target gene sets across diverse scoring methods. The package prioritizes user-friendliness, allowing direct input of expression matrices or seamless interaction with Seurat objects. Furthermore, it facilitates a comprehensive visualization of results. The irGSEA package and its accompanying documentation are accessible on GitHub (https://github.com/chuiqin/irGSEA).

Identification of Signature Genes of Dilated Cardiomyopathy Using Integrated Bioinformatics Analysis.

Dilated cardiomyopathy (DCM) is characterized by left ventricular or biventricular enlargement with systolic dysfunction. To date, the underlying molecular mechanisms of dilated cardiomyopathy pathogenesis have not been fully elucidated, although some insights have been presented. In this study, we combined public database resources and a doxorubicin-induced DCM mouse model to explore the significant genes of DCM in full depth. We first retrieved six DCM-related microarray datasets from the GEO database using several keywords. Then we used the "LIMMA" (linear model for microarray data) R package to filter each microarray for differentially expressed genes (DEGs). Robust rank aggregation (RRA), an extremely robust rank aggregation method based on sequential statistics, was then used to integrate the results of the six microarray datasets to filter out the reliable differential genes. To further improve the reliability of our results, we established a doxorubicin-induced DCM model in C57BL/6N mice, using the "DESeq2" software package to identify DEGs in the sequencing data. We cross-validated the results of RRA analysis with those of animal experiments by taking intersections and identified three key differential genes (including BEX1, RGCC and VSIG4) associated with DCM as well as many important biological processes (extracellular matrix organisation, extracellular structural organisation, sulphur compound binding, and extracellular matrix structural components) and a signalling pathway (HIF-1 signalling pathway). In addition, we confirmed the significant effect of these three genes in DCM using binary logistic regression analysis. These findings will help us to better understand the pathogenesis of DCM and may be key targets for future clinical management.

Integrated analysis of transcriptomics to identify hub genes in primary Sjögren's syndrome.

OBJECTIVE: The treatment of patients with primary Sjögren's syndrome is a clinical challenge. Gene expression profile analysis and comprehensive network methods for complex diseases can provide insight into molecular characteristics in the clinical context. MATERIALS AND METHODS: We downloaded gene expression datasets from the Gene Expression Omnibus (GEO) database. We screened differentially expressed genes (DEG) between the pSS patients and the controls by the robust rank aggregation (RRA) method. We explored DEGs' potential function using gene function annotation and PPI network analysis. RESULTS: GSE23117, GSE40611, GSE80805, and GSE127952 were included, including 38 patients and 30 controls. The RRA integrated analysis determined 294 significant DEGs (241 upregulated and 53 downregulated), and the most significant gene aberrantly expressed in SS was CXCL9 (p = 6.39E-15), followed by CXCL13 (p = 1.53E-13). Immune response (GO:0006955; p = 4.29E-32) was the most significantly enriched biological process in GO (gene ontology) analysis. KEGG pathway enrichment analysis showed that cytokine-cytokine receptor interaction (hsa04060; p = 6.46E-10) and chemokine signaling pathway (hsa04062; p = 9.54E-09) were significantly enriched. We defined PTPRC, CD86, and LCP2 as the hub genes based on the PPI results. CONCLUSION: Our integrated analysis identified gene signatures and helped understand molecular changes in pSS.

Identification of biomarkers, immune infiltration landscape, and treatment targets of ischemia-reperfusion acute kidney injury at an early stage by bioinformatics methods.

BACKGROUND: Mechanisms underlying ischemia/reperfusion injury-acute kidney injury (IRI-AKI) are not fully elucidated. We conducted an integrative analysis of IRI-AKI by bioinformatics methods. METHODS: We screened gene expression profiles of the IRI-AKI at early phase from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were identified and enrichment pathways were conducted based on gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) database, and Gene set enrichment analysis (GSEA). Immune cell infiltration analysis was performed to reveal the change of the microenvironment cell types. We constructed protein-protein interaction (PPI), and Cytoscape with plug-ins to find hub genes and modules. We performed robust rank aggregation (RRA) to combine DEGs and analyzed the target genes for miRNA/transcription factor (TF) and drug-gene interaction networks. RESULTS: A total of 239 and 384 DEGs were identified in GSE87024 and GSE34351 separately, with the 73 common DEGs. Enrichment analysis revealed that the significant pathways involve mitogen-activated protein kinase (MAPK) signaling, interleukin-17, and tumor necrosis factor (TNF) signaling pathway, etc. RRA analysis detected a total of 27 common DEGs. Immune cell infiltration analysis showed the plasma cells reduced and T cells increased in IRI-AKI. We identified JUN, ATF3, FOS, EGR1, HMOX1, DDIT3, JUNB, NFKBIZ, PPP1R15A, CXCL1, ATF4, and HSPA1B as hub genes. The target genes interacted with 23 miRNAs and 116 drugs or molecular compounds such as curcumin, staurosporine, and deferoxamine. CONCLUSION: Our study first focused on the early IRI-AKI adopting RRA analysis to combine DEGs in different datasets. We identified significant biomarkers and crucial pathways involved in IRI-AKI and first construct the immune landscape and detected the potential therapeutic targets of the IRI-AKI by drug-gene network.

Identification of four novel hub genes as monitoring biomarkers for colorectal cancer.

BACKGROUND: It must be admitted that the incidence of colorectal cancer (CRC) was on the rise all over the world, but the related treatment had not caught up. Further research on the underlying pathogenesis of CRC was conducive to improving the survival status of current CRC patients. METHODS: Differentially expressed genes (DEGs) screening were conducted based on "limma" and "RobustRankAggreg" package of R software. Weighted gene co-expression network analysis (WGCNA) was performed in the integrated DEGs that from The Cancer Genome Atlas (TCGA), and all samples of validation were from Gene Expression Omnlbus (GEO) dataset. RESULTS: The terms obtained in the functional annotation for primary DEGs indicated that they were associated with CRC. The MEyellow stand out whereby showed the significant correlation with clinical feature (disease), and 4 hub genes, including ABCC13, AMPD1, SCNN1B and TMIGD1, were identified in yellow module. Nine datasets from Gene Expression Omnibus database confirmed these four genes were significantly down-regulated and the survival estimates for the low-expression group of these genes were lower than for the high-expression group in Kaplan-Meier survival analysis section. MEXPRESS suggested that down-regulation of some top hub genes may be caused by hypermethylation. Receiver operating characteristic curves indicated that these genes had certain diagnostic efficacy. Moreover, tumor-infiltrating immune cells and gene set enrichment analysis for hub genes suggested that there were some associations between these genes and the pathogenesis of CRC. CONCLUSION: This study identified modules that were significantly associated with CRC, four novel hub genes, and further analysis of these genes. This may provide a little new insights and directions into the potential pathogenesis of CRC.

Identification of Robust and Key Differentially Expressed Genes during C2C12 Cell Myogenesis Based on Multiomics Data.

Myogenesis is a central step in prenatal myofiber formation, postnatal myofiber hypertrophy, and muscle damage repair in adulthood. RNA-Seq technology has greatly helped reveal the molecular mechanism of myogenesis, but batch effects in different experiments inevitably lead to misinterpretation of differentially expressed genes (DEGs). We previously applied the robust rank aggregation (RRA) method to effectively circumvent batch effects across multiple RNA-Seq datasets from 3T3-L1 cells. Here, we also used the RRA method to integrate nine RNA-Seq datasets from C2C12 cells and obtained 3140 robust DEGs between myoblasts and myotubes, which were then validated with array expression profiles and H3K27ac signals. The upregulated robust DEGs were highly enriched in gene ontology (GO) terms related to muscle cell differentiation and development. Considering that the cooperative binding of transcription factors (TFs) to enhancers to regulate downstream gene expression is a classical epigenetic mechanism, differentially expressed TFs (DETFs) were screened, and potential novel myogenic factors (MAF, BCL6, and ESR1) with high connection degree in protein-protein interaction (PPI) network were presented. Moreover, KLF5 cooperatively binds with the three key myogenic factors (MYOD, MYOG, and MEF2D) in C2C12 cells. Motif analysis speculates that the binding of MYOD and MYOG is KLF5-independent, while MEF2D is KLF5-dependent. It was revealed that KLF5-binding sites could be exploited to filter redundant MYOD-, MYOG-, and MEF2D-binding sites to focus on key enhancers for myogenesis. Further functional annotation of KLF5-binding sites suggested that KLF5 may regulate myogenesis through the PI3K-AKt signaling pathway, Rap1 signaling pathway, and the Hippo signaling pathway. In general, our study provides a wealth of untapped candidate targets for myogenesis and contributes new insights into the core regulatory mechanisms of myogenesis relying on KLF5-binding signal.

Integrated analysis of multiple microarray studies to identify potential pathogenic gene modules in preeclampsia.

BACKGROUND: Preeclampsia is a life-threatening hypertensive disorder during pregnancy, while underlying pathogenesis and its diagnosis are incomplete. METHODS: In this study, we utilized the Robust Rank Aggregation method to integrate 6 eligible preeclampsia microarray datasets from Gene Expression Omnibus database. We used linear regression to assess the associations between significant differentially expressed genes (DEGs) and blood pressure. Functional annotation, protein-protein interaction, Gene Set Enrichment Analysis (GSEA) and single sample GSEA were employed for investigating underlying pathogenesis in preeclampsia. RESULTS: We filtered 52 DEGs and further screened for 5 hub genes (leptin, pappalysin 2, endoglin, fms related receptor tyrosine kinase 1, tripartite motif containing 24) that were positively correlated with both systolic blood pressure and diastolic blood pressure. Receiver operating characteristic indicated that hub genes were potential biomarkers for diagnosis and prognosis in preeclampsia. GSEA for single hub gene revealed that they were all closely related to angiogenesis and estrogen response in preeclampsia. Moreover, single sample GSEA showed that the expression levels of 5 hub genes were correlated with those of immune cells in immunologic microenvironment at maternal-fetal interface. CONCLUSIONS: These findings provide new insights into underlying pathogenesis in preeclampsia; 5 hub genes were identified as biomarkers for diagnosis and prognosis in preeclampsia.

Identification of hub genes in papillary thyroid carcinoma: robust rank aggregation and weighted gene co-expression network analysis.

BACKGROUND: Papillary thyroid carcinoma (PTC), which is the most common endocrine malignancy, has been steadily increasing worldwide in incidence over the years, while mechanisms underlying the pathogenesis and diagnostic for PTC are incomplete. The purpose of this study is to identify potential biomarkers for diagnosis of PTC, and provide new insights into pathogenesis of PTC. METHODS: Based on weighted gene co-expression network analysis, Robust Rank Aggregation, functional annotation, GSEA and DNA methylation, were employed for investigating potential biomarkers for diagnosis of PTC. RESULTS: Black and turquoise modules were identified in the gene co-expression network constructed by 1807 DEGs that from 6 eligible gene expression profiles of Gene Expression Omnibus database based on Robust Rank Aggregation and weighted gene co-expression network analysis. Hub genes were significantly down-regulated and the expression levels of the hub genes were different in different stages in hub gene verification. ROC curves indicated all hub genes had good diagnostic value for PTC (except for ABCA6 AUC = 89.5%, the 15 genes with AUC > 90%). Methylation analysis showed that hub gene verification ABCA6, ACACB, RMDN1 and TFPI were identified as differentially methylated genes, and the decreased expression level of these genes may relate to abnormal DNA methylation. Moreover, the expression levels of 8 top hub genes were correlated with tumor purity and tumor-infiltrating immune cells. These findings, including functional annotations and GSEA provide new insights into pathogenesis of PTC. CONCLUSIONS: The hub genes and methylation of hub genes may as potential biomarkers provide new insights for diagnosis of PTC, and all these findings may be the direction to study the mechanisms underlying of PTC in the future.

Dysregulated microRNAs in laryngeal cancer: a comprehensive meta-analysis using a robust rank aggregation approach.

We need a reasonable method of compiling data from different studies regarding the expression of microRNA (miRNA) in laryngeal squamous cell carcinoma (LSCC). The robust rank aggregation method was used to integrate the rank lists of miRNAs from 11 studies. The enrichment analysis was performed on target genes of meta-signature miRNAs. The Cancer Genome Atlas database was used to confirm the results of meta-analysis. Three meta-signature miRNAs (miR-21-5p, miR-196a-5p and miR-145-5p) were obtained. All three miRNAs could be prognostic for LSCC. The enrichment analysis showed that these miRNAs were associated significantly with multiple cancer-related signaling pathways. The robust rank aggregation approach is an effective way to identify important miRNAs from different studies. All identified miRNAs could be candidates for LSCC diagnostic and prognostic biomarkers.

Identification of key genes and pathways of thyroid cancer by integrated bioinformatics analysis.

Thyroid cancer is a common endocrine malignancy with a rapidly increasing incidence worldwide. Although its mortality is steady or declining because of earlier diagnoses, its survival rate varies because of different tumour types. Thus, the aim of this study was to identify key biomarkers and novel therapeutic targets in thyroid cancer. The expression profiles of GSE3467, GSE5364, GSE29265 and GSE53157 were downloaded from the Gene Expression Omnibus database, which included a total of 97 thyroid cancer and 48 normal samples. After screening significant differentially expressed genes (DEGs) in each data set, we used the robust rank aggregation method to identify 358 robust DEGs, including 135 upregulated and 224 downregulated genes, in four datasets. Gene Ontology and Kyoto Encyclopaedia of Genes and Genomes pathway enrichment analyses of DEGs were performed by DAVID and the KOBAS online database, respectively. The results showed that these DEGs were significantly enriched in various cancer-related functions and pathways. Then, the STRING database was used to construct the protein-protein interaction network, and modules analysis was performed. Finally, we filtered out five hub genes, including LPAR5, NMU, FN1, NPY1R, and CXCL12, from the whole network. Expression validation and survival analysis of these hub genes based on the The Cancer Genome Atlas database suggested the robustness of the above results. In conclusion, these results provided novel and reliable biomarkers for thyroid cancer, which will be useful for further clinical applications in thyroid cancer diagnosis, prognosis and targeted therapy.

Identification of differentially expressed genes and biological characteristics of colorectal cancer by integrated bioinformatics analysis.

Colorectal cancer (CRC) ranks as one of the most common malignant tumors worldwide. Its mortality rate has remained high in recent years. Therefore, the aim of this study was to identify significant differentially expressed genes (DEGs) involved in its pathogenesis, which may be used as novel biomarkers or potential therapeutic targets for CRC. The gene expression profiles of GSE21510, GSE32323, GSE89076, and GSE113513 were downloaded from the Gene Expression Omnibus (GEO) database. After screening DEGs in each GEO data set, we further used the robust rank aggregation method to identify 494 significant DEGs including 212 upregulated and 282 downregulated genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed by DAVID and the KOBAS online database, respectively. These DEGs were shown to be significantly enriched in different cancer-related functions and pathways. Then, the STRING database was used to construct the protein-protein interaction network. The module analysis was performed by the MCODE plug-in of Cytoscape based on the whole network. We finally filtered out seven hub genes by the cytoHubba plug-in, including PPBP, CCL28, CXCL12, INSL5, CXCL3, CXCL10, and CXCL11. The expression validation and survival analysis of these hub genes were analyzed based on The Cancer Genome Atlas database. In conclusion, the robust DEGs associated with the carcinogenesis of CRC were screened through the GEO database, and integrated bioinformatics analysis was conducted. Our study provides reliable molecular biomarkers for screening and diagnosis, prognosis as well as novel therapeutic targets for CRC.

Identification of dysregulated miRNAs in triple negative breast cancer: A meta-analysis approach.

Triple negative breast cancer (TNBC) is an aggressive subtype of breast cancer with poor clinical outcomes and lack of approved targeted therapy. Dysregulated microRNAs (miRNAs) have been considered a promising biomarker, which plays an important role in the tumorigenesis of human cancer. Due to the increase in miRNA profiling datasets of TNBC, a proper analysis is required for studying. Therefore, this study used meta-analysis to amalgamate ten miRNA profiling studies of TNBC. By the robust rank aggregation method, metasignatures of six miRNAs (4 upregulated: hsa-miR-135b-5p, hsa-miR-18a-5p, hsa-miR-9-5p and hsa-miR-522-3p; 2 downregulated: hsa-miR-190b and hsa-miR-449a) were obtained. The gene ontology analysis revealed that target genes regulated by miRNAs were associated with processes like the regulation of transcription, DNA dependent, and signal transduction. Also, it is noted from the pathway analysis that signaling and cancer pathways were associated with the progression of TNBC. A Naïve Bayes-based classifier built with miRNA signatures discriminates TNBC and non-TNBC samples in test data set with high diagnostic sensitivity and specificity. From the analysis carried out by the study, it is suggested that the identified miRNAs are of great importance in improving the diagnostics and therapeutics for TNBC.

Prognostic value of gastric cancer-associated gene signatures: Evidence based on a meta-analysis using integrated bioinformatics methods.

Selecting differentially expressed genes (DEGs) based on integrated bioinformatics analyses has been used in previous studies to explore potential biomarkers in gastric cancer (GC) with microarray and RNA sequencing data. However, the genes obtained may be inaccurate because of noisy data and errors, as well as insufficient clinical sample sizes. Thus, we aimed to find robust and strong DEGs with prognostic value for GC, where the robust rank aggregation method was employed to select significant DEGs from eight Gene Expression Omnibus data sets with a total of 140 up-regulated and 206 down-regulated genes. Network data mining was then used to screen hub genes, and 11 genes were filtered using Fisher's exact test. Based on these results, we built a prognostic signature with seven genes (FBN1, MMP1, PLAU, SPARC, COL1A2, COL2A1 and ATP4A) using stepwise multivariate Cox proportional hazard regression. According to the risk score for each patient, we found that high-risk group patients had significantly worse survival results compared with those in the low-risk group (log-rank test P-value < 0.001). This seven-gene signature was then validated with an external data set. Thus, we established a signature based on seven DEGs with prognostic value for GC patients using multi-steps bioinformatics methods, which may provide novel insights and potential biomarkers for prognosis, as well as possibly serving as new therapeutic targets in clinical applications.

Identification of Potential Key Genes Associated with Adipogenesis through Integrated Analysis of Five Mouse Transcriptome Datasets.

Adipose tissue is the most important energy metabolism and secretion organ, and these functions are conferred during the adipogenesis process. However, the cause and the molecular events underlying adipogenesis are still unclear. In this study, we performed integrated bioinformatics analyses to identify vital genes involved in adipogenesis and reveal potential molecular mechanisms. Five mouse high-throughput expression profile datasets were downloaded from the Gene Expression Omnibus (GEO) database; these datasets contained 24 samples of 3T3-L1 cells during adipogenesis, including 12 undifferentiated samples and 12 differentiated samples. The five datasets were reanalyzed and integrated to select differentially expressed genes (DEGs) during adipogenesis via the robust rank aggregation (RRA) method. Functional annotation of these DEGs and mining of key genes were then performed. We also verified the expression levels of some potential key genes during adipogenesis. A total of 386 consistent DEGs were identified, with 230 upregulated genes and 156 downregulated genes. Gene Ontology (GO) analysis showed that the biological functions of the DEGs primarily included fat cell differentiation, lipid metabolic processes, and cell adhesion. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that these DEGs were mainly associated with metabolic pathways, the peroxisome proliferator-activated receptor (PPAR) signaling pathway, regulation of lipolysis in adipocytes, the tumor necrosis factor (TNF) signaling pathway, and the FoxO signaling pathway. The 30 most closely related genes among the DEGs were identified from the protein⁻protein interaction (PPI) network and verified by real-time quantification during 3T3-L1 preadipocyte differentiation. In conclusion, we obtained a list of consistent DEGs during adipogenesis through integrated analysis, which may offer potential targets for the regulation of adipogenesis and treatment of adipose dysfunction.

Predicting biomarkers related to idiopathic pulmonary fibrosis: Robust ranking aggregation analysis and animal experiment verification.

Idiopathic pulmonary fibrosis (IPF) is a progressive and incurable lung disease characterized by unknown etiology. This study employs robust ranking aggregation to identify consistent differential genes across multiple datasets, aiming to enhance prognostic evaluation and facilitate the development of more effective immunotherapy strategies for IPF. Using the GSE10667, GSE110147, and GSE24206 datasets, the analysis identifies 92 robust differentially expressed genes (DEGs), including SPP1, IGF1, ASPN, and KLHL13, highlighted as potential biomarkers through machine learning and experimental validation. Additionally, significant differences in immune cell types between IPF samples and controls, such as Plasma cells, Macrophages M0, Mast cells resting, T cells CD8, and NK cells resting, inform the construction of diagnostic and survival prediction models, demonstrating good applicability. These findings provide insights into IPF pathophysiology and suggest potential therapeutic targets.

Identification of hub genes in calcific aortic valve disease.

Calcific aortic valve disease (CAVD) is a heart valve disorder characterized primarily by calcification of the aortic valve, resulting in stiffness and dysfunction of the valve. CAVD is prevalent among aging populations and is linked to factors such as hypertension, dyslipidemia, tobacco use, and genetic predisposition, and can result in becoming a growing economic and health burden. Once aortic valve calcification occurs, it will inevitably progress to aortic stenosis. At present, there are no medications available that have demonstrated effectiveness in managing or delaying the progression of the disease. In this study, we mined four publicly available microarray datasets (GSE12644 GSE51472, GSE77287, GSE233819) associated with CAVD from the GEO database with the aim of identifying hub genes associated with the occurrence of CAVD and searching for possible biological targets for the early prevention and diagnosis of CAVD. This study provides preliminary evidence for therapeutic and preventive targets for CAVD and may provide a solid foundation for subsequent biological studies.

Uncovering hub genes and immunological characteristics for heart failure utilizing RRA, WGCNA and Machine learning.

Background: Heart failure (HF) is a major public health issue with high mortality and morbidity. This study aimed to find potential diagnostic markers for HF by the combination of bioinformatics analysis and machine learning, as well as analyze the role of immune infiltration in the pathological process of HF. Methods: The gene expression profiles of 124 HF patients and 135 nonfailing donors (NFDs) were obtained from six datasets in the NCBI Gene Expression Omnibus (GEO) public database. We applied robust rank aggregation (RRA) and weighted gene co-expression network analysis (WGCNA) method to identify critical genes in HF. To discover novel diagnostic markers in HF, three machine learning methods were employed, including best subset regression, regularization technique, and support vector machine-recursive feature elimination (SVM-RFE). Besides, immune infiltration was investigated in HF by single-sample gene set enrichment analysis (ssGSEA). Results: Combining RRA with WGCNA method, we recognized 39 critical genes associated with HF. Through integrating three machine learning methods, FCN3 and SMOC2 were determined as novel diagnostic markers in HF. Differences in immune infiltration signature were also found between HF patients and NFDs. Moreover, we explored the potential associations between two diagnostic markers and immune response in the pathogenesis of HF. Conclusions: In summary, FCN3 and SMOC2 can be used as diagnostic markers of HF, and immune infiltration plays an important role in the initiation and progression of HF.

IGFBP2 functions as an endogenous protector against hepatic steatosis via suppression of the EGFR-STAT3 pathway.

OBJECTIVE: Non-alcoholic fatty liver disease (NAFLD) is deemed as an emerging global epidemic, whereas the underlying pathogenic mechanism remains to be clarified. We aimed to systemically analyze all the NAFLD-related gene expression datasets from published human-based studies, by which exploring potential key factors and mechanisms accounting for the pathogenesis of NAFLD. METHODS: Robust rank aggregation (RRA) method was used to integrate NAFLD-related gene expression datasets. For fatty liver study, adeno-associated virus (AAV) delivery and genetic knockout mice were used to create IGFBP2 (Insulin-like growth factor binding protein 2) gain- or loss-of function models. Western blot, Co-immunoprecipitation (Co-IP), immunofluorescent (IF) staining, luciferase assay, molecular docking simulation were performed to reveal the IGFBP2-EGFR-STAT3 axis involved. Key axis protein levels in livers from healthy donors and patients with NAFLD were assessed via immunohistochemical staining. RESULTS: By using RRA method, the present study identified IGFBP2 being the most significantly down-regulated gene in all NAFLD subjects. The decreased IGFBP2 expression was further confirmed in the liver tissues from patients and animal models of NAFLD. IGFBP2 deficiency aggravated hepatic steatosis and NASH phenotypes and promoted lipogenic gene expression both in vivo and in vitro. Mechanistically, IGFBP2 directly binds to and regulates EGFR, whereas blockage of the IGFBP2-EGFR complex by knockdown of IGFBP2 resulted in the EGFR-STAT3 pathway activation, which in turn promoted the promoter activity of Srebf1. By using molecular docking simulation and protein-protein interaction analysis, the sequence of 233-257 amino acids in IGFBP2 was characterized as a key motif responding for its specific binding to EGFR and the protective effect against hepatic steatosis. CONCLUSIONS: The current study has, for the first time, identified IGFBP2 as a novel protector against hepatosteatosis. The protective effect is mediated by its specific interaction with EGFR and thereby suppressing the EGFR-STAT3 pathway. Therefore, pharmaceutically targeting the IGFBP2-EGFR-STAT3 axis may provide a theoretical basis for for the treatment of NAFLD/NASH and the associated diseases.

Identification of potential diagnostic biomarkers of atherosclerosis based on bioinformatics strategy.

BACKGROUND: Atherosclerosis is the main pathological change in atherosclerotic cardiovascular disease, and its underlying mechanisms are not well understood. The aim of this study was to explore the hub genes involved in atherosclerosis and their potential mechanisms through bioinformatics analysis. METHODS: Three microarray datasets from Gene Expression Omnibus (GEO) identified robust differentially expressed genes (DEGs) by robust rank aggregation (RRA). We performed connectivity map (CMap) analysis and functional enrichment analysis on robust DEGs and constructed a protein‒protein interaction (PPI) network using the STRING database to identify the hub gene using 12 algorithms of cytoHubba in Cytoscape. Receiver operating characteristic (ROC) analysis was used to assess the diagnostic potency of the hub genes.The CIBERSORT algorithm was used to perform immunocyte infiltration analysis and explore the association between the identified biomarkers and infiltrating immunocytes using Spearman's rank correlation analysis in R software. Finally, we evaluated the expression of the hub gene in foam cells. RESULTS: A total of 155 robust DEGs were screened by RRA and were revealed to be mainly associated with cytokines and chemokines by functional enrichment analysis. CD52 and IL1RN were identified as hub genes and were validated in the GSE40231 dataset. Immunocyte infiltration analysis showed that CD52 was positively correlated with gamma delta T cells, M1 macrophages and CD4 memory resting T cells, while IL1RN was positively correlated with monocytes and activated mast cells. RT-qPCR results indicate that CD52 and IL1RN were highly expressed in foam cells, in agreement with bioinformatics analysis. CONCLUSIONS: ​This study has established that CD52 and IL1RN may play a key role in the occurrence and development of atherosclerosis, which opens new lines of thought for further research on the pathogenesis of atherosclerosis.

Identification of hub genes and biological pathways in glioma via integrated bioinformatics analysis.

OBJECTIVE: Glioma is the most common intracranial primary malignancy, but its pathogenesis remains unclear. METHODS: We integrated four eligible glioma microarray datasets from the gene expression omnibus database using the robust rank aggregation method to identify a group of significantly differently expressed genes (DEGs) between glioma and normal samples. We used these DEGs to explore key genes closely associated with glioma survival through weighted gene co-expression network analysis. We then constructed validations of prognosis and survival analyses for the key genes via multiple databases. We also explored their potential biological functions using gene set enrichment analysis (GSEA) and gene set variation analysis (GSVA). RESULTS: We selected DLGAP5, CDCA8, NCAPH, and CCNB2, as four genes that were abnormally up-regulated in glioma samples, for verification. They showed high levels of isocitrate dehydrogenase gene mutation and tumor grades, as well as good prognostic and diagnostic value for glioma. Their methylation levels were generally lower in glioma samples. GSEA and GSVA analyses suggested the genes were closely involved with glioma proliferation. CONCLUSION: These findings provide new insights into the pathogenesis of glioma. The hub genes have the potential to be used as diagnostic and therapeutic markers.