Bioinformatic analysis of key pathways and genes involved in pediatric atopic dermatitis

The initiation of atopic dermatitis (AD) typically happens very early in life, but most of our understanding of AD is derived from studies on AD patients in adult. The aim of the present study was to identify gene signature speficic to pediatric AD comapred with adult AD. The gene expression profiles of four datasets ({"type":"entrez-geo","attrs":{"text":"GSE32924","term_id":"32924"}}GSE32924, {"type":"entrez-geo","attrs":{"text":"GSE36842","term_id":"36842"}}GSE36842, {"type":"entrez-geo","attrs":{"text":"GSE58558","term_id":"58558"}}GSE58558, and {"type":"entrez-geo","attrs":{"text":"GSE107361","term_id":"107361"}}GSE107361) were downloaded from the GEO database. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway (KEGG) enrichment analyses were performed, and protein–protein interaction (PPI) network was constructed by Cytoscape software. Total 654 differentially expressed genes (DEGs) (394 up-regulated and 260 down-regulated) were identified in pediatric AD samples with adult AD samples as control. The up-regulated DEGs were significantly enriched in the migration and chemotaxis of granulocyte and neutrophil, while down-regulated DEGs were significantly enriched in biological adhesion. KEGG pathway analysis showed that up-regulated DEGs participated in chemokine signaling pathway while down-regulated DEGs participated in adherens junction, focal adhesion, and regulation of actin cytoskeleton. The top 10 hub genes GAPDH, EGFR, ACTB, ESR1, CDK1, CXCL8, CD44, KRAS, PTGS2, and SMC3 were involved in chemokine signaling pathway, cytokine–cytokine receptor interaction, interleukin-17 signaling pathway, and regulation of actin cytoskeleton. In conclusion, we identified DEGs and hub genes involved in pediatric AD, which might be used as therapeutic targets and diagnostic biomarkers for pediatric AD.     

Bioinformatics analysis identifies COL1A1, THBS2 and SPP1 as potential predictors of patient prognosis and immunotherapy response in gastric cancer

Background: The present study aimed to use bioinformatics tools to explore pivotal genes associated with the occurrence of gastric cancer (GC) and assess their prognostic significance, and link with clinicopathological parameters. We also investigated the predictive role of COL1A1, THBS2, and SPP1 in immunotherapy.Materials and methods: We identified differential genes (DEGs) that were up- and down-regulated in the three datasets ({"type":"entrez-geo","attrs":{"text":"GSE26942","term_id":"26942"}}GSE26942, {"type":"entrez-geo","attrs":{"text":"GSE13911","term_id":"13911"}}GSE13911, and {"type":"entrez-geo","attrs":{"text":"GSE118916","term_id":"118916"}}GSE118916) and created protein–protein interaction (PPI) networks from the overlapping DEGs. We then investigated the potential functions of the hub genes in cancer prognosis using PPI networks, and explored the influence of such genes in the immune environment.Results: Overall, 268 overlapping DEGs were identified, of which 230 were up-regulated and 38 were down-regulated. CytoHubba selected the top ten hub genes, which included SPP1, TIMP1, SERPINE1, MMP3, COL1A1, BGN, THBS2, CDH2, CXCL8, and THY1. With the exception of SPP1, survival analysis using the Kaplan–Meier database showed that the levels of expression of these genes were associated with overall survival. Genes in the most dominant module explored by MCODE, COL1A1, THBS2, and SPP1, were primarily enriched for two KEGG pathways. Further analysis showed that all three genes could influence clinicopathological parameters and immune microenvironment, and there was a significant correlation between COL1A1, THBS2, SPP1, and PD-L1 expression, thus indicating a potential predictive role for GC response to immunotherapy.Conclusion: ECM–receptor interactions and focal adhesion pathways are of great significance in the progression of GC. COL1A1, THBS2, and SPP1 may help predict immunotherapy response in GC patients.     

Network-based approach to identify prognosis-related genes in tamoxifen-treated patients with estrogen receptor-positive breast cancer

Tamoxifen is an estrogen receptor (ER) antagonist that is most commonly used for the treatment of ER-positive breast cancer. However, tamoxifen resistance remains a major cause of cancer recurrence and progression. Here, we aimed to identify hub genes implicated in the progression and prognosis of ER-positive breast cancer following tamoxifen treatment. Microarray data ({"type":"entrez-geo","attrs":{"text":"GSE9893","term_id":"9893"}}GSE9893) for 155 tamoxifen-treated primary ER-positive breast cancer samples were obtained from the Gene Expression Omnibus database. In total, 1706 differentially expressed genes (DEGs), including 859 up-regulated and 847 down-regulated genes, were identified between relapse and relapse-free samples. Weighted correlation network analysis clustered genes into 13 modules, among which the tan and blue modules were the most significantly related to prognosis. From these two modules, we further identified and validated two prognosis-related hub genes (G-rich RNA sequence binding factor 1 (GRSF1) and microtubule-associated protein τ (MAPT)) via survival analysis based on several publicly available datasets. High expression of GRSF1 predicted poor prognosis, whereas MAPT indicated favorable outcomes in ER-positive breast cancer. Using breast cancer cell lines and tissue samples, we confirmed that GRSF1 was significantly up-regulated and MAPT was down-regulated in the tamoxifen-resistant group compared with the tamoxifen-sensitive group. The prognostic value of GRSF1 and MAPT was also verified in 48 tamoxifen-treated ER-positive breast cancer patients in our hospital. Gene set enrichment analysis (GSEA) suggested that GRSF1 was potentially involved in RNA degradation and cell cycle pathways, while MAPT was strongly linked to immune-related signaling pathways. Taken together, our findings established novel prognostic biomarkers to predict tamoxifen sensitivity, which may facilitate individualized management of breast cancer.     

Epigenetic alternations of microRNAs and DNA methylation contribute to gestational diabetes mellitus

This study aimed to identify epigenetic alternations of microRNAs and DNA methylation for gestational diabetes mellitus (GDM) diagnosis and treatment using in silico approach. Data of mRNA and miRNA expression microarray ({"type":"entrez-geo","attrs":{"text":"GSE103552","term_id":"103552"}}GSE103552 and {"type":"entrez-geo","attrs":{"text":"GSE104297","term_id":"104297"}}GSE104297) and DNA methylation data set ({"type":"entrez-geo","attrs":{"text":"GSE106099","term_id":"106099"}}GSE106099) were obtained from the GEO database. Differentially expressed genes (DEGs), differentially expressed miRNAs (DEMs) and differentially methylated genes (DMGs) were obtained by limma package. Functional and enrichment analyses were performed with the DAVID database. The protein‐protein interaction (PPI) network was constructed by STRING and visualized in Cytoscape. Simultaneously, a connectivity map (CMap) analysis was performed to screen potential therapeutic agents for GDM. In GDM, 184 low miRNA‐targeting up‐regulated genes and 234 high miRNA‐targeting down‐regulated genes as well as 364 hypomethylation–high‐expressed genes and 541 hypermethylation–low‐expressed genes were obtained. They were mainly enriched in terms of axon guidance, purine metabolism, focal adhesion and proteasome, respectively. In addition, 115 genes (67 up‐regulated and 48 down‐regulated) were regulated by both aberrant alternations of miRNAs and DNA methylation. Ten chemicals were identified as putative therapeutic agents for GDM and four hub genes (IGF1R, ATG7, DICER1 and RANBP2) were found in PPI and may be associated with GDM. Overall, this study identified a series of differentially expressed genes that are associated with epigenetic alternations of miRNA and DNA methylation in GDM. Ten chemicals and four hub genes may be further explored as potential drugs and targets for GDM diagnosis and treatment, respectively.     

Screening and identification of LMNB1 and DLGAP5, two key biomarkers in gliomas

Glioma is the most common primary cancer in the central nervous system. Despite advances in surgery, radiotherapy and chemotherapy over the past decades, the prognosis of glioblastoma patients remains poor. We aim to identify robust gene signatures to better understand the complex molecular mechanisms and to discover potential novel molecular biomarkers for glioma. By exploring {"type":"entrez-geo","attrs":{"text":"GSE16011","term_id":"16011"}}GSE16011, {"type":"entrez-geo","attrs":{"text":"GSE4290","term_id":"4290"}}GSE4290 and {"type":"entrez-geo","attrs":{"text":"GSE50161","term_id":"50161"}}GSE50161 data in Gene Expression Omnibus (GEO) database, we screened out 380 differentially expressed genes between non-tumor and glioma tissues, and further selected 30 hub genes through the Molecular Complex Detection (MCODE) plug-in in Cytoscape. In addition, LMNB1 and DLGAP5 were selected for further analyses due to their high expression in gliomas and were verified by using our cohort. Our study confirmed that LMNB1 and DLGAP5 were up-regulated in gliomas, and patients with high expression of LMNB1 or DLGAP5 had poor survival rate. Furthermore, silence of LMNB1 and DLGAP5 inhibited the proliferation of glioma cells. Together, LMNB1 and DLGAP5 were two potentially novel molecular biomarkers for diagnosis and prognosis of glioma.     

Mining the key genes for ventilator-induced lung injury using co-expression network analysis

Mechanical ventilation is extensively adopted in general anesthesia and respiratory failure management, but it can also induce ventilator-induced lung injury (VILI). Therefore, it is of great urgency to explore the mechanisms involved in the VILI pathogenesis, which might contribute to its future prevention and treatment. Four microarray datasets from the GEO database were selected in our investigation, and were subjected to the Weighted Gene Co-Expression Network Analysis (WGCNA) to identify the VILI-correlated gene modules. The limma package in R software was used to identify the differentially expressed genes (DEGs) between the VILI and control groups. WGCNA was constructed by merging the {"type":"entrez-geo","attrs":{"text":"GSE9314","term_id":"9314"}}GSE9314, {"type":"entrez-geo","attrs":{"text":"GSE9368","term_id":"9368"}}GSE9368, {"type":"entrez-geo","attrs":{"text":"GSE11434","term_id":"11434"}}GSE11434 and {"type":"entrez-geo","attrs":{"text":"GSE11662","term_id":"11662"}}GSE11662 datasets. A total of 49 co-expression network modules were determined as associated with VILI. The intersected genes between hub genes screened from DEGs for VILI and those identified using WGCNA were as follows: Tlr2, Hmox1, Serpine1, Mmp9, Il6, Il1b, Ptgs2, Fos and Atf3, which were determined to be key genes for VILI. Those key genes were validated by {"type":"entrez-geo","attrs":{"text":"GSE86229","term_id":"86229"}}GSE86229 and quantitative PCR (qPCR) experiment to have significantly statistical difference in their expression between the VILI and control groups. In a nutshell, nine key genes with expression differences in VILI were screened by WGCNA by integrating multiple datasets.     

Classifying Integrated Signature Molecules in Macrophages of Rheumatoid Arthritis,Osteoarthritis, and Periodontal Disease: An Omics-Based Study

Rheumatoid arthritis (RA), osteoarthritis (OA), and periodontal disease (PD) are chronic inflammatory diseases that are globally prevalent, and pose a public health concern. The search for a potential mechanism linking PD to RA and OA continues, as it could play a significant role in disease prevention and treatment. Recent studies have linked RA, OA, and PD to Porphyromonas gingivalis (PG), a periodontal bacterium, through a similar dysregulation in an inflammatory mechanism. This study aimed to identify potential gene signatures that could assist in early diagnosis as well as gain insight into the molecular mechanisms of these diseases. The expression data sets with the series IDs {"type":"entrez-geo","attrs":{"text":"GSE97779","term_id":"97779"}}GSE97779, {"type":"entrez-geo","attrs":{"text":"GSE123492","term_id":"123492"}}GSE123492, and {"type":"entrez-geo","attrs":{"text":"GSE24897","term_id":"24897"}}GSE24897 for macrophages of RA, OA synovium, and PG stimulated macrophages (PG-SM), respectively, were retrieved and screened for differentially expressed genes (DEGs). The 72 common DEGs among RA, OA, and PG-SM were further subjected to gene–gene correlation analysis. A GeneMANIA interaction network of the 47 highly correlated DEGs comprises 53 nodes and 271 edges. Network centrality analysis identified 15 hub genes, 6 of which are DEGs (API5, ATE1, CCNG1, EHD1, RIN2, and STK39). Additionally, two significantly up-regulated non-hub genes (IER3 and RGS16) showed interactions with hub genes. Functional enrichment analysis of the genes showed that “apoptotic regulation” and “inflammasomes” were among the major pathways. These eight genes can serve as important signatures/targets, and provide new insights into the molecular mechanism of PG-induced RA, OA, and PD.     

FRZB as a key molecule in abdominal aortic aneurysm progression affecting vascular integrity

Abdominal aortic aneurysm (AAA), when ruptured, results in high mortality. The identification of molecular pathways involved in AAA progression is required to improve AAA prognosis. The aim of the present study was to assess the key genes for the progression of AAA and their functional role. Genomic and clinical data of three independent cohorts were downloaded from the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) ({"type":"entrez-geo","attrs":{"text":"GSE57691","term_id":"57691"}}GSE57691, {"type":"entrez-geo","attrs":{"text":"GSE7084","term_id":"7084"}}GSE7084, and {"type":"entrez-geo","attrs":{"text":"GSE98278","term_id":"98278"}}GSE98278). To develop AAA diagnosis and progression-related differentially expressed genes (DEGs), we used a significance analysis of microarray (SAM). Spearman correlation test and gene set analysis were performed to identify potential enriched pathways for DEGs. Only the Frizzled-related protein (FRZB) gene and chromosome 1 open reading frame 24 (C1orf24) exhibited significant down-regulation in all analyses. With FRZB, the pathways were associated with RHO GTPase and elastin fiber formation. With C1orf24, the pathways were elastic fiber formation, extracellular matrix organization, and cell–cell communication. Since only FRZB was evolutionally conserved in the vertebrates, function of FRZB was validated using zebrafish embryos. Knockdown of frzb remarkably reduced vascular integrity in zebrafish embryos. We believe that FRZB is a key gene involved in AAA initiation and progression affecting vascular integrity.     

Identification of ribosomal protein family in triple-negative breast cancer by bioinformatics analysis

Triple-negative breast cancer (TNBC) accounts for ∼20% of all breast cancer (BC) cases. The management of TNBC represents a challenge due to its worse prognosis, heterogeneity and lack of targeted therapy. Moreover, its mechanisms are not fully clear. The aim of the study is to identify crucial genes between TNBC and non-TNBC for underlying targets for diagnostic and therapeutic methods of TNBC. The differentially expressed genes (DEGs) between TNBC and non-TNBC were selected from the Gene Expression Omnibus (GEO) database after the integrated analysis of two datasets ({"type":"entrez-geo","attrs":{"text":"GSE65194","term_id":"65194"}}GSE65194 and {"type":"entrez-geo","attrs":{"text":"GSE76124","term_id":"76124"}}GSE76124). Then Gene ontology (GO) and KEGG analysis were performed by DAVID database, protein–protein interaction (PPI) of DEGs was constructed by Search Tool for the Retrieval of Reciprocity Genes (STRING) database. Furthermore, centrality analysis and module analysis were carried out by Cytoscape to analyze the TNBC-related PPI. Subsequently, overall survival (OS) analysis was performed by GEPIA. Finally, the expressions of these key genes in TNBC and non-TNBC tissues were tested by qRT-PCR. The results showed that 955 DEGs were obtained, which were mainly enriched in ribosome, ribosomal subunit, and so on. Moreover, 19 candidate genes were focused on by centrality analysis and module analysis. Furthermore, we found the low expressions of ribosomal protein S9 (RPS9), ribosomal protein S14 (RPS14), ribosomal protein S27 (RPS27), ribosomal protein L11 (RPL11) and ribosomal protein L14 (RPL14) were related to a poor OS in BC patients. Additionally, qRT-PCR results suggested that these five genes were notably down-regulated in TNBC tissues. In summary, the present study suggests that ribosomal proteins are related to TNBC, and they may play an important role in the diagnosis, treatment and prognosis of TNBC.     

Identification of novel biomarkers and candidate small molecule drugs in rheumatoid arthritis and osteoarthritis based on bioinformatics analysis of high‐throughput data

Background: Rheumatoid arthritis (RA) and osteoarthritis (OA) are two major types of joint diseases. The present study aimed to identify hub genes involved in the pathogenesis and further explore the potential treatment targets of RA and OA.Methods: The gene expression profile of {"type":"entrez-geo","attrs":{"text":"GSE12021","term_id":"12021"}}GSE12021 was downloaded from Gene Expression Omnibus (GEO). Total 31 samples (12 RA, 10 OA and 9 NC samples) were used. The differentially expressed genes (DEGs) in RA versus NC, OA versus NC and RA versus OA groups were screened using limma package. We also verified the DEGs in {"type":"entrez-geo","attrs":{"text":"GSE55235","term_id":"55235"}}GSE55235 and {"type":"entrez-geo","attrs":{"text":"GSE100786","term_id":"100786"}}GSE100786. Functional annotation and protein–protein interaction (PPI) network construction of OA‐ and RA‐specific DEGs were performed. Finally, the candidate small molecules as potential drugs to treat RA and OA were predicted in CMap database.Results: 165 up-regulated and 163 down-regulated DEGs between RA and NC samples, 73 up-regulated and 293 down-regulated DEGs between OA and NC samples, 92 up-regulated and 98 down-regulated DEGs between RA and OA samples were identified. Immune response and TNF signaling pathway were significantly enriched pathways for RA‐ and OA‐specific DEGs, respectively. The hub genes were mainly associated with ‘Primary immunodeficiency’ (RA vs. NC group), ‘Ribosome’ (OA vs. NC group), and ‘Chemokine signaling pathway’ (RA vs. OA group). Arecoline and Cefamandole were the most promising small molecule to reverse the RA and OA gene expression.Conclusion: Our findings suggest new insights into the underlying pathogenesis of RA and OA, which may improve the diagnosis and treatment of these intractable chronic diseases.     

High expression of COL5A2, a member of COL5 family,indicates the poor survival and facilitates cell migration in gastric cancer

Background: Gastric cancer (GC) metastasis determines the prognosis of patients, and exploring the molecular mechanism of GC metastasis is expected to provide a theoretical basis for clinical treatment. Recent studies have shown that extracellular matrix protein is closely related to GC metastasis. The present study aimed to explore the expression profile and role of COL5A2, as an extracellular matrix protein, in GC.Methods: The expression, overall survival, and progression-free survival data of COL5 family members were extracted from The Cancer Genome Atlas (TCGA) database, respectively. Weighted gene co-expression network analysis of the {"type":"entrez-geo","attrs":{"text":"GSE62229","term_id":"62229"}}GSE62229 database was performed out to identify modules and associated genes.Results: COL5A2 was selected as our research target in the TCGA database, and was also verified in the {"type":"entrez-geo","attrs":{"text":"GSE62229","term_id":"62229"}}GSE62229 and {"type":"entrez-geo","attrs":{"text":"GSE15459","term_id":"15459"}}GSE15459 datasets. COL5A2 was up-regulated in GC tissues by paraffin immunohistochemistry and RT-qPCR. The prognosis of patients with low COL5A2 expression was better than that of patients with high COL5A2 expression. Scratch and migration experiments showed that knockdown of COL5A2 decreased the migration ability of gastric cancer cells compared with the control group. In vivo, mice with tail vein injection COL5A2 knockdown had fewer and smaller metastatic nodules in liver. GSEA results showed that the TCGA and {"type":"entrez-geo","attrs":{"text":"GSE62229","term_id":"62229"}}GSE62229 samples were significantly enriched in several well-known cancer-related pathways, such as the TGF-β, MAPK, and JAK2 signaling pathways.Conclusion: COL5A2 was most closely related to advanced GC among COL5 family members. High COL5A2 expression is associated with a poor prognosis, and may be a novel therapeutic target for GC.     

Early B-cell factors involve in the tumorigenesis and predict the overall survival of gastric cancer

Gastric cancer (GC) is a heavy health burden around the world, which is the fifth most frequent tumor and leads to the third most common cancer-related deaths. It is urgent to identify prognostic markers as the guideline for personalized treatment and follow-up. We accessed the prognostic value of Early B-cell factors (EBFs) in GC. A total of 415 GC tissues and 34 normal tissues from The Cancer Genome Atlas Stomach Adenocarcinoma (TCGA-STAD) cohort, 616 external patients from {"type":"entrez-geo","attrs":{"text":"GSE15459","term_id":"15459"}}GSE15459, {"type":"entrez-geo","attrs":{"text":"GSE22377","term_id":"22377"}}GSE22377, {"type":"entrez-geo","attrs":{"text":"GSE51105","term_id":"51105"}}GSE51105, {"type":"entrez-geo","attrs":{"text":"GSE62245","term_id":"62245"}}GSE62245 were enrolled for analysis. Univariate and multivariate Cox regression analyses were employed to evaluate the sole and integrative prognostic value of EBFs, respectively. Genetic alterations, DNA methylation of EBFs were also evaluated, as well as the involved signaling pathways. We revealed that increased EBFs associated with the poor prognosis of GC patients, the prognostic model was established in TCGA-STAD cohort, and validated in Gene Expression Omnibus (GEO) cohorts, with effectiveness in both HER2 positive and negative patients. DNA methylation was involved in the impact on prognosis. Cell cycle, immune-associated, and MAPK pathways were influenced by EBFs. Anti-CTLA4 immunotherapy is more suitable for EBFs determining high-risk groups, but not anti-PD-1/PD-L1 therapy. 5-Fluorouracil, methotrexate, vorinostat are suitable to inhibit the function of EBFs. Our new findings provide novel insight into the prediction of prognosis and clinical treatment of GC patients based on EBFs.     

Tumor microenvironment characterization in stage IV gastric cancer

Immunotherapy is remarkably affected by the immune environment of the principal tumor. Nonetheless, the immune environment’s clinical relevance in stage IV gastric cancer (GC) is largely unknown. The gene expression profiles of 403 stage IV GC patients in the three cohorts: GEO (Gene Expression Omnibus, {"type":"entrez-geo","attrs":{"text":"GSE84437","term_id":"84437"}}GSE84437 (n=292) and {"type":"entrez-geo","attrs":{"text":"GSE62254","term_id":"62254"}}GSE62254 (n=77), and TCGA (The Cancer Genome Atlas, n=34) were used in the present study. Using four publicly available stage IV GC expression datasets, 29 immune signatures were expression profiled, and on this basis, we classified stage IV GC. The classification was conducted using the hierarchical clustering method. Three stage IV GC subtypes L, M, and H were identified representing low, medium, and high immunity, respectively. Immune correlation analysis of these three types revealed that Immune H exhibited a better prognostic outcome as well as a higher immune score compared with other subtypes. There was a noticeable difference in the three subgroups of HLA genes. Further, on comparing with other subtypes, CD86, CD80, CD274, CTLA4, PDCD1, and PDCD1LG2 had higher expression in the Immunity H subtype. In stage IV GC, potentially positive associations between immune and pathway activities were displayed, due to the enrichment of pathways including TNF signaling, Th-17 cell differentiation, and JAK-STAT signaling pathways in Immunity H vs Immunity L subtypes. External cohorts from TCGA cohort ratified these results. The identification of stage IV GC subtypes has potential clinical implications in stage IV GC treatment.     

An autophagy-related prognostic signature associated with immune microenvironment features of uveal melanoma

Autophagy is involved in cancer initiation and progression but its role in uveal melanoma (UM) was rarely investigated. Herein, we built an autophagy-related gene (ARG) risk model of UM patients by univariate Cox regression and least absolute shrinkage and selection operator (Lasso) regression model and filtrated out nine prognostic ARGs in The Cancer Genome Atlas (TCGA) cohort. Survival and Receiver Operating Characteristic (ROC) Curve analysis in the TCGA and other four independent UM cohorts ({"type":"entrez-geo","attrs":{"text":"GSE22138","term_id":"22138"}}GSE22138, {"type":"entrez-geo","attrs":{"text":"GSE27831","term_id":"27831"}}GSE27831, {"type":"entrez-geo","attrs":{"text":"GSE44295","term_id":"44295"}}GSE44295 and {"type":"entrez-geo","attrs":{"text":"GSE84976","term_id":"84976"}}GSE84976) proved that the ARG-signature possessed robust and steady prognosis predictive ability. We calculated risk scores for patients included in our study and patients with higher risk scores showed worse clinical outcomes. We found the expressions of the nine ARGs were significantly associated with clinical and molecular features (including risk score) and overall survival (OS) of UM patients. Furthermore, we utilized univariate and multivariate Cox regression analyses to determine the independent prognostic ability of the ARG-signature. Functional enrichment analysis showed the ARG-signature was correlated with several immune-related processes and pathways like T-cell activation and T-cell receptor signaling pathway. Gene set enrichment analysis (GSEA) found tumor hallmarks including angiogenesis, IL6-JAK-STAT3-signaling, reactive oxygen species pathway and oxidative phosphorylation were enriched in high-risk UM patients. Finally, infiltrations of several immune cells and immune-related scores were found significantly associated with the ARG-signature. In conclusion, the ARG-signature might be a strong predictor for evaluating the prognosis and immune infiltration of UM patients.     

A panel of two miRNAs correlated to systolic blood pressure is a good diagnostic indicator for stroke

Background: We aimed to develop a diagnostic indicator of stroke based on serum miRNAs correlated to systolic blood pressure.Methods: Using miRNA expression profiles in {"type":"entrez-geo","attrs":{"text":"GSE117604","term_id":"117604"}}GSE117604 from the Gene Expression Omnibus (GEO), we utilized the WGCNA to identify hub miRNAs correlated to systolic blood pressure (SBP). Differential analysis was applied to highlight hub differentially expressed miRNAs (DE-miRNAs), whereby we built a miRNA-based diagnostic indicator for stroke using bootstrap ranking Least Absolute Shrinkage and Selection Operator (LASSO) regression with 10-fold cross-validation. The classification value of the indicator was validated with receiver operating characteristic (ROC) analysis in both the training set and test set, as well as quantitative real-time PCR (qRT-PCR) for the feature miRNAs. Further, target genes of hub miRNAs and hub DE-miRNAs were retrieved for functional enrichment.Results: A total of 447 hub miRNAs in the blue modules were significantly correlated with systolic blood pressure (r = 0.32, false discovery rate = 10⁻⁶). Target genes predicted with the hub miRNAs were mostly implicated in the Kyoto Encyclopedia of Genes and Genomes (KEGG) terms including mitogen-activated protein kinase (MAPK) pathway, senescence, and TGF-β signaling pathway. The diagnostic indicator with miR-4420 and miR-6793-5p showed remarkable performance in the training set (area under curve [AUC]= 0.953), as well as in the test set (AUC = 0.894). Results of qRT-PCR validated the diagnostic value of the two miRNAs embedded in the proposed indicator.Conclusions: We developed a panel of two miRNAs, which is a good diagnostic indicator for stroke. These results require further investigation.     

Identification of potential biomarkers in dengue via integrated bioinformatic analysis

Dengue fever virus (DENV) is a global health threat that is becoming increasingly critical. However, the pathogenesis of dengue has not yet been fully elucidated. In this study, we employed bioinformatics analysis to identify potential biomarkers related to dengue fever and clarify their underlying mechanisms. The results showed that there were 668, 1901, and 8283 differentially expressed genes between the dengue-infected samples and normal samples in the {"type":"entrez-geo","attrs":{"text":"GSE28405","term_id":"28405"}}GSE28405, {"type":"entrez-geo","attrs":{"text":"GSE38246","term_id":"38246"}}GSE38246, and {"type":"entrez-geo","attrs":{"text":"GSE51808","term_id":"51808"}}GSE51808 datasets, respectively. Through overlapping, a total of 69 differentially expressed genes (DEGs) were identified, of which 51 were upregulated and 18 were downregulated. We identified twelve hub genes, including MX1, IFI44L, IFI44, IFI27, ISG15, STAT1, IFI35, OAS3, OAS2, OAS1, IFI6, and USP18. Except for IFI44 and STAT1, the others were statistically significant after validation. We predicted the related microRNAs (miRNAs) of these 12 target genes through the database miRTarBase, and finally obtained one important miRNA: has-mir-146a-5p. In addition, gene ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment were carried out, and a protein–protein interaction (PPI) network was constructed to gain insight into the actions of DEGs. In conclusion, our study displayed the effectiveness of bioinformatics analysis methods in screening potential pathogenic genes in dengue fever and their underlying mechanisms. Further, we successfully predicted IFI44L and IFI6, as potential biomarkers with DENV infection, providing promising targets for the treatment of dengue fever to a certain extent.     

Bioinformatics Analysis of Key micro-RNAs and mRNAs under the Hand,Foot, and Mouth Disease Virus Infection

To clarify crucial key micro-RNAs and mRNAs associated with hand, foot, and mouth disease (HFMD) virus infection, we conducted this bioinformatics analysis from four GEO datasets. The following datasets were used for the analysis: {"type":"entrez-geo","attrs":{"text":"GSE85829","term_id":"85829"}}GSE85829, {"type":"entrez-geo","attrs":{"text":"GSE94551","term_id":"94551"}}GSE94551, {"type":"entrez-geo","attrs":{"text":"GSE52780","term_id":"52780"}}GSE52780, and {"type":"entrez-geo","attrs":{"text":"GSE45589","term_id":"45589"}}GSE45589. Differentially expressed genes (DEGs) were acquired, and the analysis of functional and pathway enrichment and the relative regulatory network were conducted. After screening common differentially expressed miRNAs (DE-miRNAs), five key miRNAs were acquired: miR-100-3p, miR-125a-3p, miR-1273g-3p, miR-5585-3p, and miR-671-5p. There were three common enriched GO terms between miRNA-derived prediction and mRNA-derived analysis: biosynthetic process, cytosol, and nucleoplasm. There was one common KEGG pathway, i.e., cell cycle shared between miRNA-based and mRNA-based enrichment. Using TarBase V8 in DIANA tools, we acquired 1,520 potential targets (mRNA) from the five key DE-miRNAs, among which the159 DE-mRNAs also included 11 DEGs. These common DEGs showed a PPI network mainly connected by SMC1A, SMARCC1, SF3B3, LIG1, and BRMS1L. Together, changes in five key miRNAs and 11 key mRNAs may play crucial roles in HFMD progression. A combination of these roles may benefit the early diagnosis and treatment of HFMD.Key words: HFMD, micro-RNA, protein-protein interaction, microarray, regulatory network     

Exploring the Novel Computational Drug Target and Associated Key Pathways of Oral Cancer

Oral cancer (OC) is a serious health concern that has a high fatality rate. The oral cavity has seven kinds of OC, including the lip, tongue, and floor of the mouth, as well as the buccal, hard palate, alveolar, retromolar trigone, and soft palate. The goal of this study is to look into new biomarkers and important pathways that might be used as diagnostic biomarkers and therapeutic candidates in OC. The publicly available repository the Gene Expression Omnibus (GEO) was to the source for the collection of OC-related datasets. {"type":"entrez-geo","attrs":{"text":"GSE74530","term_id":"74530"}}GSE74530, {"type":"entrez-geo","attrs":{"text":"GSE23558","term_id":"23558"}}GSE23558, and {"type":"entrez-geo","attrs":{"text":"GSE3524","term_id":"3524"}}GSE3524 microarray datasets were collected for analysis. Minimum cut-off criteria of |log fold-change (FC)| > 1 and adjusted p < 0.05 were applied to calculate the upregulated and downregulated differential expression genes (DEGs) from the three datasets. After that only common DEGs in all three datasets were collected to apply further analysis. Gene ontology (GO) and pathway analysis were implemented to explore the functional behaviors of DEGs. Then protein–protein interaction (PPI) networks were built to identify the most active genes, and a clustering algorithm was also implemented to identify complex parts of PPI. TF-miRNA networks were also constructed to study OC-associated DEGs in-depth. Finally, top gene performers from PPI networks were used to apply drug signature analysis. After applying filtration and cut-off criteria, 2508, 3377, and 670 DEGs were found for {"type":"entrez-geo","attrs":{"text":"GSE74530","term_id":"74530"}}GSE74530, {"type":"entrez-geo","attrs":{"text":"GSE23558","term_id":"23558"}}GSE23558, and {"type":"entrez-geo","attrs":{"text":"GSE3524","term_id":"3524"}}GSE3524 respectively, and 166 common DEGs were found in every dataset. The GO annotation remarks that most of the DEGs were associated with the terms of type I interferon signaling pathway. The pathways of KEGG reported that the common DEGs are related to the cell cycle and influenza A. The PPI network holds 88 nodes and 492 edges, and CDC6 had the highest number of connections. Four clusters were identified from the PPI. Drug signatures doxorubicin and resveratrol showed high significance according to the hub genes. We anticipate that our bioinformatics research will aid in the definition of OC pathophysiology and the development of new therapies for OC.