Search

Background: Osteoarthritis (OA) is a degenerative disease requiring additional research. This study compared gene expression and immune infiltration between lesioned and preserved subchondral bone. The results were validated using multiple tissue datasets and experiments. Methods: Differentially expressed genes (DEGs) between the lesioned and preserved tibial plateaus of OA patients were identified in the GSE51588 dataset. Moreover, functional annotation and protein–protein interaction (PPI) network analyses were performed on the lesioned and preserved sides to explore potential therapeutic targets in OA subchondral bones. In addition, multiple tissues were used to screen coexpressed genes, and the expression levels of identified candidate DEGs in OA were measured by quantitative real-time polymerase chain reaction. Finally, an immune infiltration analysis was conducted. Results: A total of 1,010 DEGs were identified, 423 upregulated and 587 downregulated. The biological process (BP) terms enriched in the upregulated genes included "skeletal system development", "sister chromatid cohesion", and "ossification". Pathways were enriched in "Wnt signaling pathway" and "proteoglycans in cancer". The BP terms enriched in the downregulated genes included "inflammatory response", "xenobiotic metabolic process", and "positive regulation of inflammatory response". The enriched pathways included "neuroactive ligand–receptor interaction" and "AMP-activated protein kinase signaling". JUN, tumor necrosis factor α, and interleukin-1β were the hub genes in the PPI network. Collagen XI A1 and leucine-rich repeat-containing 15 were screened from multiple datasets and experimentally validated. Immune infiltration analyses showed fewer infiltrating adipocytes and endothelial cells in the lesioned versus preserved samples. Conclusion: Our findings provide valuable information for future studies on the pathogenic mechanism of OA and potential therapeutic and diagnostic targets. [ABSTRACT FROM AUTHOR]

Background Osteoarthritis (OA) is a degenerative disease requiring additional research. This study compared gene expression and immune infiltration between lesioned and preserved subchondral bone. The results were validated using multiple tissue datasets and experiments. Methods Differentially expressed genes (DEGs) between the lesioned and preserved tibial plateaus of OA patients were identified in the GSE51588 dataset. Moreover, functional annotation and protein–protein interaction (PPI) network analyses were performed on the lesioned and preserved sides to explore potential therapeutic targets in OA subchondral bones. In addition, multiple tissues were used to screen coexpressed genes, and the expression levels of identified candidate DEGs in OA were measured by quantitative real-time polymerase chain reaction. Finally, an immune infiltration analysis was conducted. Results A total of 1,010 DEGs were identified, 423 upregulated and 587 downregulated. The biological process (BP) terms enriched in the upregulated genes included “skeletal system development”, “sister chromatid cohesion”, and “ossification”. Pathways were enriched in “Wnt signaling pathway” and “proteoglycans in cancer”. The BP terms enriched in the downregulated genes included “inflammatory response”, “xenobiotic metabolic process”, and “positive regulation of inflammatory response”. The enriched pathways included “neuroactive ligand–receptor interaction” and “AMP-activated protein kinase signaling”. JUN, tumor necrosis factor α, and interleukin-1β were the hub genes in the PPI network. Collagen XI A1 and leucine-rich repeat-containing 15 were screened from multiple datasets and experimentally validated. Immune infiltration analyses showed fewer infiltrating adipocytes and endothelial cells in the lesioned versus preserved samples. Conclusion Our findings provide valuable information for future studies on the pathogenic mechanism of OA and potential therapeutic and diagnostic targets.

Both oxidative stress and the immune response are associated with heart failure (HF). In this study, our aim was to identify the hub genes associated with oxidative stress andimmune infiltration of HF by bioinformatics analysis and experimental verification. The expression profile of GSE36074 was obtained from the Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) were screened by GEO2R. The genes related to oxidative stress were extracted from GeneCards websites. Then, the functional enrichment analysis of oxidative stress-related DEGs (OSRDEGs) was performed using DAVID. In addition, we constructed a protein-protein interaction (PPI) network using the STRING database and screened for hub genes with Cytoscape software. We also used CIBERSORTx to analyze immune infiltration in mice heart tissues between the TAC and Sham groups and explored the correlation between immune cells and hub genes. Finally, the hub genes were carried out using reverse transcription quantitative PCR (RT-qPCR), immunohistochemistry (IHC) and western blot. A total of 136 OSRDEGs were found in GSE36074. Enrichment analysis revealed that these OSRDEGs were enriched in the mitochondrion, HIF-1, FoxO, MAPK and TNF signaling pathway. The five hub genes (Mapk14, Hif1a, Myc, Hsp90ab1, and Hsp90aa1) were screened by the cytoHubba plugin. The correlation analysis between immune cells and hub genes showed that Mapk14 was positively correlated with Th2 Cells, while Hif1a and Hsp90ab1exhibited a negative correlation with Th2 Cells; Myc exhibited a negative correlation with Monocytes; whereas, Hsp90aa1 was negatively correlated with NK Resting. Finally, five hub genes were validated by RT-qPCR, IHC and western blot. Mapk14, Hif1a, Myc, Hsp90ab1, and Hsp90aa1 are hub genes of HF and may play a critical role in the oxidative stress of HF. This study may provide new targets for the treatment of HF, and the potential immunotherapies are worthy of further study. [ABSTRACT FROM AUTHOR]

Background: The mechanism of immunoglobulin A nephropathy (IgAN) is still unknown. A bioinformatics analysis is a powerful method to identify the biomarkers and possible therapeutic targets of a certain disease from related datasets. Methods: The GSE93973 dataset, obtained from the Gene Expression Omnibus (GEO) database, was used to construct a weighted gene co-expression network (WGCNA) and filter differentially expressed genes (DEGs). The biological process (BP) enrichment among all the genes in the key modules was analyzed through a Gene Ontology (GO) enrichment analysis. We selected the overlap of hub genes in the WGCNA and Protein-Protein Interaction (PPI) network as the final hub genes in IgAN. We verified the final hub genes in two other datasets and in clinical kidney tissue specimens. A receiver operating characteristic (ROC) curve was used to evaluate the diagnostic efficacy of hub genes for IgAN. Results: The turquoise module, which contained 1,806 genes, was the module with the highest correlation coefficient with IgAN in the GSE93973 dataset. The GO enrichment analysis showed that these 1,806 genes were mainly enriched in inflammation and immune responses. There were five hub genes identified by WGCNA and 34 hub genes identified in a DEG analysis in the GSE93973 dataset. DUSP1 and FOSB were identified as the final hub genes in IgAN. The validation results of the final hub genes in two other databases and clinical kidney tissue specimens validated the result that, compared to the control group, FOSB and DUSP1 were expressed at lower levels in the glomerulus of IgAN patients. The ROC curve indicated that DUSP1 and FOSB were good diagnostic indicators for IgAN. Conclusions: Our analysis identified two hub genes that might be potential targets for the intervention and treatment of IgAN. [ABSTRACT FROM AUTHOR]

Background. Identification of accurate prognostic biomarkers is still particularly urgent for improving the poor survival of lung cancer patients. In this study, we aimed to identity the potential biomarkers in Chinese lung cancer population via bioinformatics analysis. Methods. In this study, the differentially expressed genes (DEGs) in lung cancer were identified using six datasets from Gene Expression Omnibus (GEO) database. Subsequently, enrichment analysis was conducted to evaluate the underlying molecular mechanisms involved in progression of lung cancer. Protein-protein interaction (PPI) and CytoHubba analysis were performed to determine the hub genes. The GEPIA, Human Protein Atlas (HPA), Kaplan-Meier plotter, and TIMER databases were used to explore the hub genes. The receiver operating characteristic (ROC) analysis was performed to evaluate the diagnostic value of hub genes. Reverse transcription quantitative PCR (qRT-PCR) was used to validate the expression levels of hub genes in 10 pairs of lung cancer paired tissues. Results. A total of 499 overlapping DEGs (160 upregulated and 339 downregulated genes) were identified in the microarray datasets. DEGs were mainly associated with pathways in cancer, focal adhesion, and protein digestion and absorption. There were nine hub genes (CDKN3, MKI67, CEP55, SPAG5, AURKA, TOP2A, UBE2C, CHEK1 and BIRC5) identified by PPI and module analysis. In GEPIA database, the expression levels of these genes in lung cancer tissues were significantly upregulated compared with normal lung tissues. The results of prognostic analysis showed that relatively higher expression of hub genes was associated with poor prognosis of lung cancer. In HPA database, most hub genes were highly expressed in lung cancer tissues. The hub genes have good diagnostic efficiency in lung cancer and normal tissues. The expression of any hub gene was associated with the infiltration of at least two immune cells. qRT- PCR confirmed that the expression level of CDKN3, MKI67, CEP55, SPAG5, AURKA, TOP2A were highly expressed in lung cancer tissues. Conclusions. The hub genes and functional pathways identified in this study may contribute to understand the molecular mechanisms of lung cancer. Our findings may provide new therapeutic targets for lung cancer patients. [ABSTRACT FROM AUTHOR]

Background. Hepatocellular carcinoma (HCC) is one of the most common malignancies around the world. Among the risk factors involved in liver carcinogenesis, hepatitis B virus (HBV) X protein (HBx) is considered to be a key regulator in hepatocarcinogenesis. Whether HBx promotes or protects against HCC remains controversial, therefore exploring new HBx-associated genes is still important. Methods. HBx was overexpressed in HepG2, HepG2.2.15 and SMMC-7721 cell lines, primary mouse hepatocytes and livers of C57BL/6N mice. High-throughput RNA sequencing profiling of HepG2 cells with HBx overexpression and related differentially-expressed genes (DEGs), pathway enrichment analysis, protein-protein interaction networks (PPIs), overlapping analysis were conducted. In addition, Gene Expression Omnibus (GEO) and proteomic datasets of HBV-positive HCC datasets were used to verify the expression and prognosis of selected DEGs. Finally, we also evaluated the known oncogenic role of HBx by oncogenic array analysis. Results. A total of 523 DEGs were obtained from HBx-overexpressing HepG2 cells. Twelve DEGs were identified and validated in cells transiently transfected with HBx and three datasets of HBV-positive HCC transcription profiles. In addition, using the Kaplan-Meier plotter database, the expression levels of the twelve different genes were further analyzed to predict patient outcomes. Conclusion. Among the 12 identified HBx-associated hub genes, HBV-positive HCC patients expressing ARG1 and TAT showed a good overall survival (OS) and relapse-free survival (RFS). Thus, ARG1 and TAT expression could be potential prognostic markers. [ABSTRACT FROM AUTHOR]

Background: Diffuse large B-cell lymphoma (DLBCL) is a highly heterogeneous malignancy with varied outcomes. However, the fundamental mechanisms remain to be fully defined. Aim: We aimed to identify core differentially co-expressed hub genes and perturbed pathways relevant to the pathogenesis and prognosis of DLBCL. Methods: We retrieved the raw gene expression profile and clinical information of GSE12453 from the Gene Expression Omnibus (GEO) database. We used integrated bioinformatics analysis to identify differentially co-expressed genes. The CIBERSORT analysis was also applied to predict tumor-infiltrating immune cells (TIICs) in the GSE12453 dataset. We performed survival and ssGSEA (single-sample Gene Set Enrichment Analysis) (for TIICs) analyses and validated the hub genes using GEPIA2 and an independent GSE31312 dataset. Results: We identified 46 differentially co-expressed hub genes in the GSE12453 dataset. Gene expression levels and survival analysis found 15 differentially co-expressed core hub genes. The core genes prognostic values and expression levels were further validated in the GEPIA2 database and GSE31312 dataset to be reliable (p < 0.01). The core genes’ main KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichments were Ribosome and Coronavirus disease-COVID-19. High expressions of the 15 core hub genes had prognostic value in DLBCL. The core genes showed significant predictive accuracy in distinguishing DLBCL cases from non-tumor controls, with the area under the curve (AUC) ranging from 0.992 to 1.00. Finally, CIBERSORT analysis on GSE12453 revealed immune cells, including activated memory CD4+ T cells and M0, M1, and M2-macrophages as the infiltrates in the DLBCL microenvironment. Conclusion: Our study found differentially co-expressed core hub genes and relevant pathways involved in ribosome and COVID-19 disease that may be potential targets for prognosis and novel therapeutic intervention in DLBCL. [ABSTRACT FROM AUTHOR]

Background: Chemotherapy resistance, especially platinum resistance, is the main cause of poor prognosis of ovarian cancer. It is of great urgency to find molecular markers and mechanism related to platinum resistance in ovarian cancer. Methods: One mRNA dataset (GSE28739) and one miRNA dataset (GSE25202) were acquired from Gene Expression Omnibus (GEO) database. The GEO2R tool was used to screen out differentially expressed genes (DEGs) and differentially expressed miRNAs (DE-miRNAs) between platinum-resistant and platinumsensitive ovarian cancer patients. Gene Ontology (GO) function and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis for DEGs were performed using the DAVID to present the most visibly enriched pathways. Protein–protein interaction (PPI) of these DEGs was constructed based on the information of the STRING database. Hub genes related to platinum resistance were visualized by Cytoscape software. Then, we chose seven interested hub genes to further validate using qRT-PCR in A2780 ovarian cancer cell lines. And, at last, the TF-miRNA-target genes regulatory network was predicted and constructed using miRNet software. Results: A total of 63 upregulated DEGs, 124 downregulated DEGs, four upregulated miRNAs and six downregulated miRNAs were identified. From the PPI network, the top 10 hub genes were identified, which were associated with platinum resistance. Our further qRT-PCR showed that seven hub genes (BUB1, KIF2C, NUP43, NDC80, NUF2, CCNB2 and CENPN) were differentially expressed in platinum-resistant ovarian cancer cells. Furthermore, the upstream transcription factors (TF) for upregulated DE-miRNAs were SMAD4, NFKB1, SMAD3, TP53 and HNF4A. Three overlapping downstream target genes (KIF2C, STAT3 and BUB1) were identified by miRNet, which was regulated by hsa-miR-494. Conclusions: The TF-miRNA–mRNA regulatory pairs, that is TF (SMAD4, NFKB1 and SMAD3)-miR-494-target genes (KIF2C, STAT3 and BUB1), were established. In conclusion, the present study is of great significance to find the key genes of platinum resistance in ovarian cancer. Further study is needed to identify the mechanism of these genes in ovarian cancer. [ABSTRACT FROM AUTHOR]

Background. Renal ischemia-reperfusion injury (IRI) is a disease with high incidence rate in kidney related surgery. Micro RNA (miRNA) and transcription factors (TFs) are widely involved in the process of renal IRI through regulation of their target genes. However, the regulatory relationships and functional roles of TFs, miRNAs and mRNAs in the progression of renal IRI are insufficiently understood. The present study aimed to clarify the underlying mechanism of regulatory relationships in renal IRI. Methods. Six gene expression profiles were downloaded from Gene Expression Omnibus (GEO). Differently expressed genes (DEGs) and differently expressed miRNAs (DEMs) were identified through RRA integrated analysis ofmRNAdatasets (GSE39548, GSE87025, GSE52004, GSE71647, and GSE131288) and miRNA datasets (GSE29495). miRDB and TransmiR v2.0 database were applied to predict target genes of miRNA and TFs, respectively. DEGs were applied for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, followed with construction of protein-protein interaction (PPI) network. Then, the TF-miRNA-mRNA network was constructed. Correlation coefficient and ROC analysis were used to verify regulatory relationship between genes and their diagnostic value in GSE52004. Furthermore, in independent mouse RNA-seq datasets GSE98622, human RNA-seq GSE134386 and in vitro, the expression of hub genes and genes from the network were observed and correlation coefficient and ROC analysis were validated. Results. A total of 21 DEMs and 187 DEGs were identified in renal IRI group compared to control group. The results of PPI analysis showed 15 hub genes. The TF-miRNAmRNA regulatory network was constructed and several important pathways were identified and further verified, including Junb-miR-223-Ranbp3l, Cebpb-miR-223- Ranbp3l, Cebpb-miR-21-Ranbp3l and Cebpb-miR-181b-Bsnd. Four regulatory loops were identified, including Fosl2-miR-155, Fosl2-miR-146a, Cebpb-miR-155 and MafkmiR- 25. The hub genes and genes in the network showed good diagnostic value in mice and human. [ABSTRACT FROM AUTHOR]

Background Diffuse large B-cell lymphoma (DLBCL) is a highly heterogeneous malignancy with varied outcomes. However, the fundamental mechanisms remain to be fully defined. Aim We aimed to identify core differentially co-expressed hub genes and perturbed pathways relevant to the pathogenesis and prognosis of DLBCL. Methods We retrieved the raw gene expression profile and clinical information of GSE12453 from the Gene Expression Omnibus (GEO) database. We used integrated bioinformatics analysis to identify differentially co-expressed genes. The CIBERSORT analysis was also applied to predict tumor-infiltrating immune cells (TIICs) in the GSE12453 dataset. We performed survival and ssGSEA (single-sample Gene Set Enrichment Analysis) (for TIICs) analyses and validated the hub genes using GEPIA2 and an independent GSE31312 dataset. Results We identified 46 differentially co-expressed hub genes in the GSE12453 dataset. Gene expression levels and survival analysis found 15 differentially co-expressed core hub genes. The core genes prognostic values and expression levels were further validated in the GEPIA2 database and GSE31312 dataset to be reliable (p < 0.01). The core genes’ main KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichments were Ribosome and Coronavirus disease-COVID-19. High expressions of the 15 core hub genes had prognostic value in DLBCL. The core genes showed significant predictive accuracy in distinguishing DLBCL cases from non-tumor controls, with the area under the curve (AUC) ranging from 0.992 to 1.00. Finally, CIBERSORT analysis on GSE12453 revealed immune cells, including activated memory CD4+ T cells and M0, M1, and M2-macrophages as the infiltrates in the DLBCL microenvironment. Conclusion Our study found differentially co-expressed core hub genes and relevant pathways involved in ribosome and COVID-19 disease that may be potential targets for prognosis and novel therapeutic intervention in DLBCL.

Background Alzheimer’s disease (AD) is an age-related neurodegenerative disease caused by central nervous system disorders. Late-onset Alzheimer disease (LOAD) is the most common neurodegenerative disorder worldwide. Differences at the expression level of certain genes, resulting from either genetic variations or environmental interactions, might be one of the mechanisms underlying differential risks for developing AD. Peripheral blood genome transcriptional profiling may provide a powerful and minimally invasive tool for the identification of novel targets beyond Aβ and tau for AD research. Methods This preliminary study explores molecular pathogenesis of LOAD-related inflammation through next generation sequencing, to assess RNA expression profiles in peripheral blood from five patients with LOAD and 10 healthy controls. Results The analysis of RNA expression profiles revealed 94 genes up-regulated and 147 down-regulated. Gene function analysis, including Gene Ontology (GO) and KOBAS-Kyoto Encyclopedia of DEGs and Genomes (KEGG) pathways indicated upregulation of interferon family (INF) signaling, while the down-regulated genes were mainly associated with the cell cycle process. KEGG metabolic pathways mapping showed gene expression alterations in the signaling pathways of JAK/STAT, chemokines, MAP kinases and Alzheimer disease. The results of this preliminary study provided not only a comprehensive picture of gene expression, but also the key processes associated with pathology for the regulation of neuroinflammation, to improve the current mechanisms to treat LOAD.

Background. Alzheimer’s disease (AD) is an age-related neurodegenerative disease caused by central nervous system disorders. Late-onset Alzheimer disease (LOAD) is the most common neurodegenerative disorder worldwide. Differences at the expression level of certain genes, resulting from either genetic variations or environmental interactions, might be one of the mechanisms underlying differential risks for developing AD. Peripheral blood genome transcriptional profiling may provide a powerful and minimally invasive tool for the identification of novel targets beyond Aβ and tau for AD research. Methods. This preliminary study explores molecular pathogenesis of LOAD-related inflammation through next generation sequencing, to assess RNA expression profiles in peripheral blood from five patients with LOAD and 10 healthy controls. Results. The analysis of RNA expression profiles revealed 94 genes up-regulated and 147 down-regulated. Gene function analysis, including Gene Ontology (GO) and KOBASKyoto Encyclopedia of DEGs and Genomes (KEGG) pathways indicated upregulation of interferon family (INF) signaling, while the down-regulated genes were mainly associated with the cell cycle process. KEGG metabolic pathways mapping showed gene expression alterations in the signaling pathways of JAK/STAT, chemokines, MAP kinases and Alzheimer disease. The results of this preliminary study provided not only a comprehensive picture of gene expression, but also the key processes associated with pathology for the regulation of neuroinflammation, to improve the current mechanisms to treat LOAD. [ABSTRACT FROM AUTHOR]

Chilling stress (CS) is an important limiting factor for the growth and development of passion fruit (Passiflora edulis) in winter in South China. However, little is known about how the passion fruit responds and adapts to CS. In this study, we performed transcriptome sequencing of cold-susceptible cultivar Huangjinguo (HJG) and cold-tolerant cultivar Tainong 1 (TN1) under normal temperature (NT) and CS conditions, and a total of 47,353 unigenes were obtained by seven databases. Using differentially expressed unigenes (DEGs) analysis, 3,248 and 4,340 DEGs were identified at two stages, respectively. The Gene Ontology (GO) enrichment analysis showed that the DEGs were mainly related to phosphorylation, membrane protein, and catalytic activity. In Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, the unigenes of plant-pathogen interaction, plant hormone signal transduction and fatty acid metabolism were enriched. Then, the 12,471 filtered unigenes were clustered into eight co-expression modules, and two modules were correlated with CS. In this two modules, 32 hub unigenes were obtained. Furthermore, the unigenes related to CS were validated using quantitative real-time PCR (RT-qPCR). This work showed that the expression levels of CS-related unigenes were very different in two passion fruit cultivars. The results provide information for the development of passion fruit with increased chilling tolerance.

Background Ischemia-reperfusion injury and drug-induced nephrotoxicity are the two most common reasons for acute kidney injury (AKI). However, little attention has been paid to early activation of fibroblasts in the progression of AKI to chronic kidney disease (CKD). The present study aimed to identify related genes and pathways on fibroblast activation in two mouse models of AKI: ischemia-reperfusion injury (IRI) model and folic acid (FA)-induced injury model. Methods The microarray expression profiles of GSE62732 and GSE121190 were downloaded from the GEO database, and the differentially expressed genes (DEGs) was analyzed using the Limma package of R software. Principal component analysis (PCA) was also performed using R. The functional information of gene products was annotated by Gene Ontology (GO) and DAVID online database, and the pathway analysis was carried out by using the Kyoto Encyclopedia of Genes and Genomes pathway (KEGG) database. Protein-protein interactions (PPI) network was constructed by STRING and Cytoscape. Furthermore, in the Hypoxia/Reoxygenation (H/R) model, the morphological changes of cells were observed under microscope and the expression of the hub genes in NRK-49F cells were validated by qRT-PCR assays. Results A total of 457 DEGs were identified. Among these, 215 DEGs were upregulated and 242 DEGs were downregulated in the acute injured samples compared with uninjured samples. The GO enrichment analysis indicated that these DEGs were mainly involved in transport, the oxidation-reduction process, the metabolic process, metal ion binding, hydrolase activity, and oxidoreductase activity. The KEGG analysis revealed that these DEGs were significantly enriched in the PI3K-Akt signaling pathway, protein digestion and absorption pathway, and focal adhesion pathway. The hub genes including Hnf4α, Pck1 and Timp1 were validated by the qRT-PCR assay in NRK-49F cells in the H/R model. Conclusions Hnf4α, Pck1 and Timp-1 may play a pivotal role in the early activation of fibroblasts, providing novel therapeutic strategies for early prediction and treatment of renal fibrosis.

Background Hepatocellular carcinoma (HCC) is one of the most commonly diagnosed cancers and the fourth leading cause of cancer-related deaths in the world. Although the treatment of HCC has made great progress in recent years, the therapeutic effects on HCC are still unsatisfactory due to difficulty in early diagnosis, chemoresistance and high recurrence rate post-surgery. Methods In this study, we identified differentially expressed genes (DEGs) based on four Gene Expression Omnibus (GEO) datasets (GSE45267, GSE98383, GSE101685 and GSE112790) between HCC and normal hepatic tissues. A protein–protein interaction (PPI) network was established to identify the central nodes associated with HCC. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of the central nodes were conducted to find the hub genes. The expression levels of the hub genes were validated based on the ONCOMINE and Gene Expression Profiling Interactive Analysis (GEPIA) databases. Additionally, the genetic alterations of the hub genes were evaluated by cBioPortal. The role of the hub genes on the overall survival (OS) and relapse survival (RFS) of HCC patients was evaluated by Kaplan-Meier plotter. At last, the mechanistic role of the hub genes was illustrated by in vitro experiments. Results We found the following seven hub genes: BUB1B, CCNB1, CCNB2, CDC20, CDK1, MAD2L1 and RRM2 using integrated bioinformatics analysis. All of the hub genes were significantly upregulated in HCC tissues. And the seven hub genes were associated with the OS and RFS of HCC patients. Finally, in vitro experiments indicated that BUB1B played roles in HCC cell proliferation, migration, invasion, apoptosis and cell cycle by partially affecting mitochondrial functions. Conclusions In summary, we identified seven hub genes that were associated with the expression and prognosis of HCC. The mechanistic oncogenic role of BUB1B in HCC was first illustrated. BUB1B might play an important role in HCC and could be potential therapeutic targets for HCC.

Background Idiopathic pulmonary fibrosis (IPF) is a chronic and progressive interstitial lung disease, characterized by a decline in lung function. To date, the pathophysiologic mechanisms associated with lung dysfunction remain unclear, and no effective therapy has been identified to improve lung function. Methods In the present study, we used weighted gene co-expression network analysis (WGCNA) to identify key modules and hub genes associated with lung function in IPF. Three datasets, containing clinical information, were downloaded from Gene Expression Omnibus. WGCNA was performed on the GSE32537 dataset. Differentially expressed gene s (DEGs) between IPF patients and healthy controls were also identified to filter hub genes. The relationship between hub genes and lung function was then validated using the GSE47460 and GSE24206 datasets. Results The red module, containing 267 genes, was positively correlated with the St. George’s Respiratory Questionnaire score (r = 0.37, p

Background Colorectal cancer (CRC) is the third most common cancer in the world. The present study is aimed at identifying hub genes associated with the progression of CRC. Method The data of the patients with CRC were obtained from the Gene Expression Omnibus (GEO) database and assessed by weighted gene co-expression network analysis (WGCNA), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses performed in R by WGCNA, several hub genes that regulate the mechanism of tumorigenesis in CRC were identified. Differentially expressed genes in the data sets GSE28000 and GSE42284 were used to construct a co-expression network for WGCNA. The yellow, black and blue modules associated with CRC level were filtered. Combining the co-expression network and the PPI network, 15 candidate hub genes were screened. Results After validation using the TCGA-COAD dataset, a total of 10 hub genes (MT1X, MT1G, MT2A, CXCL8, IL1B, CXCL5, CXCL11, IL10RA, GZMB, KIT) closely related to the progression of CRC were identified. The expressions of MT1G, CXCL8, IL1B, CXCL5, CXCL11 and GZMB in CRC tissues were higher than normal tissues (p-value < 0.05). The expressions of MT1X, MT2A, IL10RA and KIT in CRC tissues were lower than normal tissues (p-value < 0.05). Conclusions By combinating with a series of methods including GO enrichment analysis, KEGG pathway analysis, PPI network analysis and gene co-expression network analysis, we identified 10 hub genes that were associated with the progression of CRC.

Background Pancreatic cancer is one of the most common malignant cancers worldwide. Currently, the pathogenesis of pancreatic cancer remains unclear; thus, it is necessary to explore its precise molecular mechanisms. Methods To identify candidate genes involved in the tumorigenesis and proliferation of pancreatic cancer, the microarray datasets GSE32676, GSE15471 and GSE71989 were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) between Pancreatic ductal adenocarcinoma (PDAC) and nonmalignant samples were screened by GEO2R. The Database for Annotation Visualization and Integrated Discovery (DAVID) online tool was used to obtain a synthetic set of functional annotation information for the DEGs. A PPI network of the DEGs was established using the Search Tool for the Retrieval of Interacting Genes (STRING) database, and a combination of more than 0.4 was considered statistically significant for the PPI. Subsequently, we visualized the PPI network using Cytoscape. Functional module analysis was then performed using Molecular Complex Detection (MCODE). Genes with a degree ≥10 were chosen as hub genes, and pathways of the hub genes were visualized using ClueGO and CluePedia. Additionally, GenCLiP 2.0 was used to explore interactions of hub genes. The Literature Mining Gene Networks module was applied to explore the cocitation of hub genes. The Cytoscape plugin iRegulon was employed to analyze transcription factors regulating the hub genes. Furthermore, the expression levels of the 13 hub genes in pancreatic cancer tissues and normal samples were validated using the Gene Expression Profiling Interactive Analysis (GEPIA) platform. Moreover, overall survival and disease-free survival analyses according to the expression of hub genes were performed using Kaplan-Meier curve analysis in the cBioPortal online platform. The relationship between expression level and tumor grade was analyzed using the online database Oncomine. Lastly, the eight snap-frozen tumorous and adjacent noncancerous adjacent tissues of pancreatic cancer patients used to detect the CDK1 and CEP55 protein levels by western blot. Conclusions Altogether, the DEGs and hub genes identified in this work can help uncover the molecular mechanisms underlying the tumorigenesis of pancreatic cancer and provide potential targets for the diagnosis and treatment of this disease.

Background Lung squamous cell carcinoma (LUSC) is a major subtype of lung cancer with limited therapeutic options and poor clinical prognosis. Methods Three datasets (GSE19188, GSE33532 and GSE33479) were obtained from the gene expression omnibus (GEO) database. Differentially expressed genes (DEGs) between LUSC and normal tissues were identified by GEO2R, and functional analysis was employed using the Database for Annotation, Visualization and Integrated Discovery (DAVID) online tool. Protein–protein interaction (PPI) and hub genes were identified via the Search Tool for the Retrieval of Interacting Genes (STRING) and Cytoscape software. Hub genes were further validated in The Cancer Genome Atlas (TCGA) database. Subsequently, survival analysis was performed using the Kapla–Meier curve and Cox progression analysis. Based on univariate and multivariate Cox progression analysis, a gene signature was established to predict overall survival. Receiver operating characteristic curve was used to evaluate the prognostic value of the model. Results A total of 116 up-regulated genes and 84 down-regulated genes were identified. These DEGs were mainly enriched in the two pathways: cell cycle and p53 signaling way. According to the degree of protein nodes in the PPI network, 10 hub genes were identified. The mRNA expression levels of the 10 hub genes in LUSC were also significantly up-regulated in the TCGA database. Furthermore, a novel seven-gene signature (FLRT3, PPP2R2C, MMP3, MMP12, CAPN8, FILIP1 and SPP1) from the DEGs was constructed and acted as a significant and independent prognostic signature for LUSC. Conclusions The 10 hub genes might be tightly correlated with LUSC progression. The seven-gene signature might be an independent biomarker with a significant predictive value in LUSC overall survival.

Purpose Pterygium results from a variety of biological pathways that are involved in the formation of ocular surface diseases. However, the exact pathogenesis of pterygium is still unclear. Our study focused on gene expression profiles to better understand the potential mechanisms of pterygium. Methods RNA sequencing experiments were performed on clinical pterygium tissues and normal conjunctival tissues. To identify the hub genes for the development of pterygium, we further conducted weighted gene co-expression network analysis (WGCNA). qRT-PCR was utilized to validate the dysregulation of the most significant differentially expressed genes (DEGs) and key hub genes in the independent subjects. Results A total of 339 DEGs (P-adjusted < 0.05 and log2 fold change [log2FC] ≥ 1.0) were obtained that reached statistical significance with p-values < 0.05. Among them, 200 DEGs were upregulated; these genes were mainly associated with the extracellular matrix and with cell adhesion or migration. In contrast, the 139 downregulated genes were enriched for endocrine and inflammation pathways. With regard to WGCNA, five modules were assigned based on the DEG profiles, and the biological functions of each module were verified with previously published GO terms. The functions included ECM-receptor interactions, the PI3K-Akt signalling pathway and an endoplasmic reticulum (ER)-related pathway. The five hub genes with the highest connectivity in each module and the five most significant DEGs showed dysregulated expression in the independent cohort samples. Conclusions RNA sequencing and WGCNA provided novel insights into the potential regulatory mechanisms of pterygium. The identified DEGs and hub genes, which were classified into two groups according to different functions or signalings, may provide important references for further research on the molecular biology of pterygium.

Osteoporosis is a major public health problem that is associated with high morbidity and mortality, and its prevalence is increasing as the world’s population ages. Therefore, understanding the molecular basis of the disease is becoming a high priority. In this regard, studies have shown that an imbalance in adipogenic and osteogenic differentiation of bone marrow mesenchymal stem cells (MSCs) is associated with osteoporosis. In this study, we conducted a Weighted Gene Co-Expression Network Analysis to identify gene modules associated with the differentiation of bone marrow MSCs. Gene Ontology and Kyoto Encyclopedia of Genes and Genome enrichment analysis showed that the most significant module, the brown module, was enriched with genes involved in cell cycle regulation, which is in line with the initial results published using these data. In addition, the Cytoscape platform was used to identify important hub genes and lncRNAs correlated with the gene modules. Furthermore, differential gene expression analysis identified 157 and 40 genes that were upregulated and downregulated, respectively, after 3 h of MSCs differentiation. Interestingly, regulatory network analysis, and comparison of the differentially expressed genes with those in the brown module identified potential novel biomarker genes, including two transcription factors (ZNF740, FOS) and two hub genes (FOXQ1, SGK1), which were further validated for differential expression in another data set of differentiation of MSCs. Finally, Gene Set Enrichment Analysis suggested that the two most important candidate hub genes are involved in regulatory pathways, such as the JAK-STAT and RAS signaling pathways. In summary, we have revealed new molecular mechanisms of MSCs differentiation and identified novel genes that could be used as potential therapeutic targets for the treatment of osteoporosis.

Purpose Mantle cell lymphoma (MCL) is a rare and aggressive subtype of non-Hodgkin lymphoma that is incurable with standard therapies. The use of gene expression analysis has been of interest, recently, to detect biomarkers for cancer. There is a great need for systemic coexpression network analysis of MCL and this study aims to establish a gene coexpression network to forecast key genes related to the pathogenesis and prognosis of MCL. Methods The microarray dataset GSE93291 was downloaded from the Gene Expression Omnibus database. We systematically identified coexpression modules using the weighted gene coexpression network analysis method (WGCNA). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis were performed on the modules deemed important. The protein–protein interaction networks were constructed and visualized using Cytoscape software on the basis of the STRING website; the hub genes in the top weighted network were identified. Survival data were analyzed using the Kaplan–Meier method and were compared using the log-rank test. Results Seven coexpression modules consisting of different genes were applied to 5,000 genes in the 121 human MCL samples using WGCNA software. GO and KEGG enrichment analysis identified the blue module as one of the most important modules; the most critical pathways identified were the ribosome, oxidative phosphorylation and proteasome pathways. The hub genes in the top weighted network were regarded as real hub genes (IL2RB, CD3D, RPL26L1, POLR2K, KIF11, CDC20, CCNB1, CCNA2, PUF60, SNRNP70, AKT1 and PRPF40A). Survival analysis revealed that seven genes (KIF11, CDC20, CCNB1, CCNA2, PRPF40A, CD3D and PUF60) were associated with overall survival time (p < 0.05). Conclusions The blue module may play a vital role in the pathogenesis of MCL. Five real hub genes (KIF11, CDC20, CCNB1, CCNA2 and PUF60) were identified as potential prognostic biomarkers as well as therapeutic targets with clinical utility for MCL.

Background Juvenile dermatomyositis (JDM) is an immune-mediated disease characterized by chronic organ inflammation. The pathogenic mechanisms remain ill-defined. Methods Raw microarray data of JDM were obtained from the gene expression omnibus (GEO) database. Based on the GSE3307 dataset with 39 samples, weighted correlation network analysis (WGCNA) was performed to identify key modules associated with pathological state. Functional enrichment analyses were conducted to identify potential mechanisms. Based on the criteria of high connectivity and module membership, candidate hub genes were selected. A protein-protein interaction network was constructed to identify hub genes. Another dataset (GSE11971) was used for the validation of real hub genes. Finally, the real hub genes were used to screen out small-molecule compounds via the Connectivity map database. Results Three modules were considered as key modules for the pathological state of JDM. Functional enrichment analysis indicated that responses to interferon and metabolism were dysregulated. A total of 45 candidate hub genes were selected according to the pre-established criteria, and 20 genes could differentiate JDM from normal controls by validation of another external dataset (GSE11971). These real hub genes suggested the pivotal role of mitochondrial dysfunction and interferon signature in JDM. Furthermore, drug repositioning highlighted the importance of acacetin, helveticoside, lanatoside C, deferoxamine, LY-294002, tanespimycin and L01AD from downregulated genes with the potential to perturb the development of JDM, while betonicine, felodipine, valproic acid, trichostatin A and sirolimus from upregulated genes provided potentially therapeutic goals for JDM. Conclusions There are 20 real hub genes associated with the pathological state of JDM, suggesting the pivotal role of mitochondrial dysfunction and interferon signature in JDM. This analysis predicted several kinds of small-molecule compounds to treat JDM.

Background Adrenocortical carcinoma (ACC) is a rare and aggressive malignant cancer in the adrenal cortex with poor prognosis. Though previous research has attempted to elucidate the progression of ACC, its molecular mechanism remains poorly understood. Methods Gene transcripts per million (TPM) data were downloaded from the UCSC Xena database, which included ACC (The Cancer Genome Atlas, n = 77) and normal samples (Genotype Tissue Expression, n = 128). We used weighted gene co-expression network analysis to identify gene connections. Overall survival (OS) was determined using the univariate Cox model. A protein–protein interaction (PPI) network was constructed by the search tool for the retrieval of interacting genes. Results To determine the critical genes involved in ACC progression, we obtained 2,953 significantly differentially expressed genes and nine modules. Among them, the blue module demonstrated significant correlation with the “Stage” of ACC. Enrichment analysis revealed that genes in the blue module were mainly enriched in cell division, cell cycle, and DNA replication. Combined with the PPI and co-expression networks, we identified four hub genes (i.e., TOP2A, TTK, CHEK1, and CENPA) that were highly expressed in ACC and negatively correlated with OS. Thus, these identified genes may play important roles in the progression of ACC and serve as potential biomarkers for future diagnosis.

Background Understanding hub genes involved in gastric cancer (GC) metastasis could lead to effective approaches to diagnose and treat cancer. In this study, we aim to identify the hub genes and investigate the underlying molecular mechanisms of GC. Methods To explore potential therapeutic targets for GC,three expression profiles (GSE54129, GSE33651, GSE81948) of the genes were extracted from the Gene Expression Omnibus (GEO) database. The GEO2R online tool was applied to screen out differentially expressed genes (DEGs) between GC and normal gastric samples. Database for Annotation, Visualization and Integrated Discovery was applied to perform Gene Ontology (GO) and KEGG pathway enrichment analysis. The protein-protein interaction (PPI) network of these DEGs was constructed using a STRING online software. The hub genes were identified by the CytoHubba plugin of Cytoscape software. Then, the prognostic value of these identified genes was verified by gastric cancer database derived from Kaplan-Meier plotter platform. Results A total of 85 overlapped upregulated genes and 44 downregulated genes were identified. The majority of the DEGs were enriched in extracellular matrix organization, endodermal cell differentiation, and endoderm formation. Moreover, five KEGG pathways were significantly enriched, including ECM-receptor interaction, amoebiasis, AGE-RAGE signaling pathway in diabetic complications, focal adhesion, protein digestion and absorption. By combining the results of PPI network and CytoHubba, a total of nine hub genes including COL1A1, THBS1, MMP2, CXCL8, FN1, TIMP1, SPARC, COL4A1, and ITGA5 were selected. The Kaplan-Meier plotter database confirmed that overexpression levels of these genes were associated with reduced overall survival, except for THBS1 and CXCL8. Conclusions Our study suggests that COL1A1, MMP2, FN1, TIMP1, SPARC, COL4A1, and ITGA5 may be potential biomarkers and therapeutic targets for GC. Further study is needed to assess the effect of THBS1 and CXCL8 on GC.

Both oxidative stress and the immune response are associated with heart failure (HF). In this study, our aim was to identify the hub genes associated with oxidative stress andimmune infiltration of HF by bioinformatics analysis and experimental verification. The expression profile of GSE36074 was obtained from the Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) were screened by GEO2R. The genes related to oxidative stress were extracted from GeneCards websites. Then, the functional enrichment analysis of oxidative stress-related DEGs (OSRDEGs) was performed using DAVID. In addition, we constructed a protein-protein interaction (PPI) network using the STRING database and screened for hub genes with Cytoscape software. We also used CIBERSORTx to analyze immune infiltration in mice heart tissues between the TAC and Sham groups and explored the correlation between immune cells and hub genes. Finally, the hub genes were carried out using reverse transcription quantitative PCR (RT-qPCR), immunohistochemistry (IHC) and western blot. A total of 136 OSRDEGs were found in GSE36074. Enrichment analysis revealed that these OSRDEGs were enriched in the mitochondrion, HIF-1, FoxO, MAPK and TNF signaling pathway. The five hub genes (Mapk14, Hif1a, Myc, Hsp90ab1, and Hsp90aa1) were screened by the cytoHubba plugin. The correlation analysis between immune cells and hub genes showed that Mapk14 was positively correlated with Th2 Cells, while Hif1a and Hsp90ab1exhibited a negative correlation with Th2 Cells; Myc exhibited a negative correlation with Monocytes; whereas, Hsp90aa1 was negatively correlated with NK Resting. Finally, five hub genes were validated by RT-qPCR, IHC and western blot. Mapk14, Hif1a, Myc, Hsp90ab1, and Hsp90aa1 are hub genes of HF and may play a critical role in the oxidative stress of HF. This study may provide new targets for the treatment of HF, and the potential immunotherapies are worthy of further study., Competing Interests: The authors declare there are no competing interests., (©2023 Gu et al.)

Background: The mechanism of immunoglobulin A nephropathy (IgAN) is still unknown. A bioinformatics analysis is a powerful method to identify the biomarkers and possible therapeutic targets of a certain disease from related datasets., Methods: The GSE93973 dataset, obtained from the Gene Expression Omnibus (GEO) database, was used to construct a weighted gene co-expression network (WGCNA) and filter differentially expressed genes (DEGs). The biological process (BP) enrichment among all the genes in the key modules was analyzed through a Gene Ontology (GO) enrichment analysis. We selected the overlap of hub genes in the WGCNA and Protein-Protein Interaction (PPI) network as the final hub genes in IgAN. We verified the final hub genes in two other datasets and in clinical kidney tissue specimens. A receiver operating characteristic (ROC) curve was used to evaluate the diagnostic efficacy of hub genes for IgAN., Results: The turquoise module, which contained 1,806 genes, was the module with the highest correlation coefficient with IgAN in the GSE93973 dataset. The GO enrichment analysis showed that these 1,806 genes were mainly enriched in inflammation and immune responses. There were five hub genes identified by WGCNA and 34 hub genes identified in a DEG analysis in the GSE93973 dataset. DUSP1 and FOSB were identified as the final hub genes in IgAN. The validation results of the final hub genes in two other databases and clinical kidney tissue specimens validated the result that, compared to the control group, FOSB and DUSP1 were expressed at lower levels in the glomerulus of IgAN patients. The ROC curve indicated that DUSP1 and FOSB were good diagnostic indicators for IgAN., Conclusions: Our analysis identified two hub genes that might be potential targets for the intervention and treatment of IgAN., Competing Interests: The authors declare that they have no competing interests., (© 2022 Jianping et al.)

Background: Identification of accurate prognostic biomarkers is still particularly urgent for improving the poor survival of lung cancer patients. In this study, we aimed to identity the potential biomarkers in Chinese lung cancer population via bioinformatics analysis., Methods: In this study, the differentially expressed genes (DEGs) in lung cancer were identified using six datasets from Gene Expression Omnibus (GEO) database. Subsequently, enrichment analysis was conducted to evaluate the underlying molecular mechanisms involved in progression of lung cancer. Protein-protein interaction (PPI) and CytoHubba analysis were performed to determine the hub genes. The GEPIA, Human Protein Atlas (HPA), Kaplan-Meier plotter, and TIMER databases were used to explore the hub genes. The receiver operating characteristic (ROC) analysis was performed to evaluate the diagnostic value of hub genes. Reverse transcription quantitative PCR (qRT-PCR) was used to validate the expression levels of hub genes in 10 pairs of lung cancer paired tissues., Results: A total of 499 overlapping DEGs (160 upregulated and 339 downregulated genes) were identified in the microarray datasets. DEGs were mainly associated with pathways in cancer, focal adhesion, and protein digestion and absorption. There were nine hub genes (CDKN3, MKI67, CEP55, SPAG5, AURKA, TOP2A, UBE2C, CHEK1 and BIRC5) identified by PPI and module analysis. In GEPIA database, the expression levels of these genes in lung cancer tissues were significantly upregulated compared with normal lung tissues. The results of prognostic analysis showed that relatively higher expression of hub genes was associated with poor prognosis of lung cancer. In HPA database, most hub genes were highly expressed in lung cancer tissues. The hub genes have good diagnostic efficiency in lung cancer and normal tissues. The expression of any hub gene was associated with the infiltration of at least two immune cells. qRT-PCR confirmed that the expression level of CDKN3, MKI67, CEP55, SPAG5, AURKA, TOP2A were highly expressed in lung cancer tissues., Conclusions: The hub genes and functional pathways identified in this study may contribute to understand the molecular mechanisms of lung cancer. Our findings may provide new therapeutic targets for lung cancer patients., Competing Interests: Mengzhen Li and Zhouyu Wang are employed by MyGene Diagnostics Co., Ltd., (©2022 Liu et al.)

Background: Hepatocellular carcinoma (HCC) is one of the most common malignancies around the world. Among the risk factors involved in liver carcinogenesis, hepatitis B virus (HBV) X protein (HBx) is considered to be a key regulator in hepatocarcinogenesis. Whether HBx promotes or protects against HCC remains controversial, therefore exploring new HBx-associated genes is still important., Methods: HBx was overexpressed in HepG2, HepG2.2.15 and SMMC-7721 cell lines, primary mouse hepatocytes and livers of C57BL/6N mice. High-throughput RNA sequencing profiling of HepG2 cells with HBx overexpression and related differentially-expressed genes (DEGs), pathway enrichment analysis, protein-protein interaction networks (PPIs), overlapping analysis were conducted. In addition, Gene Expression Omnibus (GEO) and proteomic datasets of HBV-positive HCC datasets were used to verify the expression and prognosis of selected DEGs. Finally, we also evaluated the known oncogenic role of HBx by oncogenic array analysis., Results: A total of 523 DEGs were obtained from HBx-overexpressing HepG2 cells. Twelve DEGs were identified and validated in cells transiently transfected with HBx and three datasets of HBV-positive HCC transcription profiles. In addition, using the Kaplan-Meier plotter database, the expression levels of the twelve different genes were further analyzed to predict patient outcomes., Conclusion: Among the 12 identified HBx-associated hub genes, HBV-positive HCC patients expressing ARG1 and TAT showed a good overall survival (OS) and relapse-free survival (RFS). Thus, ARG1 and TAT expression could be potential prognostic markers., Competing Interests: The authors declare there are no competing interests., (©2021 Ni et al.)

Background: Chemotherapy resistance, especially platinum resistance, is the main cause of poor prognosis of ovarian cancer. It is of great urgency to find molecular markers and mechanism related to platinum resistance in ovarian cancer., Methods: One mRNA dataset (GSE28739) and one miRNA dataset (GSE25202) were acquired from Gene Expression Omnibus (GEO) database. The GEO2R tool was used to screen out differentially expressed genes (DEGs) and differentially expressed miRNAs (DE-miRNAs) between platinum-resistant and platinum-sensitive ovarian cancer patients. Gene Ontology (GO) function and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis for DEGs were performed using the DAVID to present the most visibly enriched pathways. Protein-protein interaction (PPI) of these DEGs was constructed based on the information of the STRING database. Hub genes related to platinum resistance were visualized by Cytoscape software. Then, we chose seven interested hub genes to further validate using qRT-PCR in A2780 ovarian cancer cell lines. And, at last, the TF-miRNA-target genes regulatory network was predicted and constructed using miRNet software., Results: A total of 63 upregulated DEGs, 124 downregulated DEGs, four upregulated miRNAs and six downregulated miRNAs were identified. From the PPI network, the top 10 hub genes were identified, which were associated with platinum resistance. Our further qRT-PCR showed that seven hub genes (BUB1, KIF2C, NUP43, NDC80, NUF2, CCNB2 and CENPN) were differentially expressed in platinum-resistant ovarian cancer cells. Furthermore, the upstream transcription factors (TF) for upregulated DE-miRNAs were SMAD4, NFKB1, SMAD3, TP53 and HNF4A. Three overlapping downstream target genes (KIF2C, STAT3 and BUB1) were identified by miRNet, which was regulated by hsa-miR-494., Conclusions: The TF-miRNA-mRNA regulatory pairs, that is TF (SMAD4, NFKB1 and SMAD3)-miR-494-target genes (KIF2C, STAT3 and BUB1), were established. In conclusion, the present study is of great significance to find the key genes of platinum resistance in ovarian cancer. Further study is needed to identify the mechanism of these genes in ovarian cancer., Competing Interests: The authors declare that they have no competing interests., (© 2021 Wang et al.)

Background: Renal ischemia-reperfusion injury (IRI) is a disease with high incidence rate in kidney related surgery. Micro RNA (miRNA) and transcription factors (TFs) are widely involved in the process of renal IRI through regulation of their target genes. However, the regulatory relationships and functional roles of TFs, miRNAs and mRNAs in the progression of renal IRI are insufficiently understood. The present study aimed to clarify the underlying mechanism of regulatory relationships in renal IRI., Methods: Six gene expression profiles were downloaded from Gene Expression Omnibus (GEO). Differently expressed genes (DEGs) and differently expressed miRNAs (DEMs) were identified through RRA integrated analysis of mRNA datasets (GSE39548, GSE87025, GSE52004, GSE71647, and GSE131288) and miRNA datasets (GSE29495). miRDB and TransmiR v2.0 database were applied to predict target genes of miRNA and TFs, respectively. DEGs were applied for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, followed with construction of protein-protein interaction (PPI) network. Then, the TF-miRNA-mRNA network was constructed. Correlation coefficient and ROC analysis were used to verify regulatory relationship between genes and their diagnostic value in GSE52004. Furthermore, in independent mouse RNA-seq datasets GSE98622, human RNA-seq GSE134386 and in vitro , the expression of hub genes and genes from the network were observed and correlation coefficient and ROC analysis were validated., Results: A total of 21 DEMs and 187 DEGs were identified in renal IRI group compared to control group. The results of PPI analysis showed 15 hub genes. The TF-miRNA-mRNA regulatory network was constructed and several important pathways were identified and further verified, including Junb-miR-223-Ranbp3l, Cebpb-miR-223-Ranbp3l, Cebpb-miR-21-Ranbp3l and Cebpb-miR-181b-Bsnd. Four regulatory loops were identified, including Fosl2-miR-155, Fosl2-miR-146a, Cebpb-miR-155 and Mafk-miR-25. The hub genes and genes in the network showed good diagnostic value in mice and human., Conclusions: In this study, we found 15 hub genes and several TF-miRNA-mRNA pathways, which are helpful for understanding the molecular and regulatory mechanisms in renal IRI. Junb-miR-223-Ranbp3l, Cebpb-miR-223-Ranbp3l, Cebpb-miR-21-Ranbp3l and Cebpb-miR-181b-Bsnd were the most important pathways, while Spp1, Fos, Timp1, Tnc, Fosl2 and Junb were the most important hub genes. Fosl2-miR-155, Fosl2-miR-146a, Cebpb-miR-155 and Mafk-miR-25 might be the negative feedback loops in renal IRI., Competing Interests: The authors declare there are no competing interests., (©2021 Ke et al.)

Background: Ischemia-reperfusion injury and drug-induced nephrotoxicity are the two most common reasons for acute kidney injury (AKI). However, little attention has been paid to early activation of fibroblasts in the progression of AKI to chronic kidney disease (CKD). The present study aimed to identify related genes and pathways on fibroblast activation in two mouse models of AKI: ischemia-reperfusion injury (IRI) model and folic acid (FA)-induced injury model., Methods: The microarray expression profiles of GSE62732 and GSE121190 were downloaded from the GEO database, and the differentially expressed genes (DEGs) was analyzed using the Limma package of R software. Principal component analysis (PCA) was also performed using R. The functional information of gene products was annotated by Gene Ontology (GO) and DAVID online database, and the pathway analysis was carried out by using the Kyoto Encyclopedia of Genes and Genomes pathway (KEGG) database. Protein-protein interactions (PPI) network was constructed by STRING and Cytoscape. Furthermore, in the Hypoxia/Reoxygenation (H/R) model, the morphological changes of cells were observed under microscope and the expression of the hub genes in NRK-49F cells were validated by qRT-PCR assays., Results: A total of 457 DEGs were identified. Among these, 215 DEGs were upregulated and 242 DEGs were downregulated in the acute injured samples compared with uninjured samples. The GO enrichment analysis indicated that these DEGs were mainly involved in transport, the oxidation-reduction process, the metabolic process, metal ion binding, hydrolase activity, and oxidoreductase activity. The KEGG analysis revealed that these DEGs were significantly enriched in the PI3K-Akt signaling pathway, protein digestion and absorption pathway, and focal adhesion pathway. The hub genes including Hnf4α, Pck1 and Timp1 were validated by the qRT-PCR assay in NRK-49F cells in the H/R model., Conclusions: Hnf4α, Pck1 and Timp-1 may play a pivotal role in the early activation of fibroblasts, providing novel therapeutic strategies for early prediction and treatment of renal fibrosis., Competing Interests: The authors declare there are no competing interests., (©2021 Deng et al.)

Chilling stress (CS) is an important limiting factor for the growth and development of passion fruit ( Passiflora edulis ) in winter in South China. However, little is known about how the passion fruit responds and adapts to CS. In this study, we performed transcriptome sequencing of cold-susceptible cultivar Huangjinguo (HJG) and cold-tolerant cultivar Tainong 1 (TN1) under normal temperature (NT) and CS conditions, and a total of 47,353 unigenes were obtained by seven databases. Using differentially expressed unigenes (DEGs) analysis, 3,248 and 4,340 DEGs were identified at two stages, respectively. The Gene Ontology (GO) enrichment analysis showed that the DEGs were mainly related to phosphorylation, membrane protein, and catalytic activity. In Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, the unigenes of plant-pathogen interaction, plant hormone signal transduction and fatty acid metabolism were enriched. Then, the 12,471 filtered unigenes were clustered into eight co-expression modules, and two modules were correlated with CS. In this two modules, 32 hub unigenes were obtained. Furthermore, the unigenes related to CS were validated using quantitative real-time PCR (RT-qPCR). This work showed that the expression levels of CS-related unigenes were very different in two passion fruit cultivars. The results provide information for the development of passion fruit with increased chilling tolerance., Competing Interests: The authors declare there are no competing interests., (©2021 Wu et al.)

Background: Hepatocellular carcinoma (HCC) is one of the most commonly diagnosed cancers and the fourth leading cause of cancer-related deaths in the world. Although the treatment of HCC has made great progress in recent years, the therapeutic effects on HCC are still unsatisfactory due to difficulty in early diagnosis, chemoresistance and high recurrence rate post-surgery., Methods: In this study, we identified differentially expressed genes (DEGs) based on four Gene Expression Omnibus (GEO) datasets (GSE45267, GSE98383, GSE101685 and GSE112790) between HCC and normal hepatic tissues. A protein-protein interaction (PPI) network was established to identify the central nodes associated with HCC. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of the central nodes were conducted to find the hub genes. The expression levels of the hub genes were validated based on the ONCOMINE and Gene Expression Profiling Interactive Analysis (GEPIA) databases. Additionally, the genetic alterations of the hub genes were evaluated by cBioPortal. The role of the hub genes on the overall survival (OS) and relapse survival (RFS) of HCC patients was evaluated by Kaplan-Meier plotter. At last, the mechanistic role of the hub genes was illustrated by in vitro experiments., Results: We found the following seven hub genes: BUB1B, CCNB1, CCNB2, CDC20, CDK1, MAD2L1 and RRM2 using integrated bioinformatics analysis. All of the hub genes were significantly upregulated in HCC tissues. And the seven hub genes were associated with the OS and RFS of HCC patients. Finally, in vitro experiments indicated that BUB1B played roles in HCC cell proliferation, migration, invasion, apoptosis and cell cycle by partially affecting mitochondrial functions., Conclusions: In summary, we identified seven hub genes that were associated with the expression and prognosis of HCC. The mechanistic oncogenic role of BUB1B in HCC was first illustrated. BUB1B might play an important role in HCC and could be potential therapeutic targets for HCC., Competing Interests: The authors declare there are no competing interests., (©2021 Fu et al.)

Background: Idiopathic pulmonary fibrosis (IPF) is a chronic and progressive interstitial lung disease, characterized by a decline in lung function. To date, the pathophysiologic mechanisms associated with lung dysfunction remain unclear, and no effective therapy has been identified to improve lung function., Methods: In the present study, we used weighted gene co-expression network analysis (WGCNA) to identify key modules and hub genes associated with lung function in IPF. Three datasets, containing clinical information, were downloaded from Gene Expression Omnibus. WGCNA was performed on the GSE32537 dataset. Differentially expressed gene s (DEGs) between IPF patients and healthy controls were also identified to filter hub genes. The relationship between hub genes and lung function was then validated using the GSE47460 and GSE24206 datasets., Results: The red module, containing 267 genes, was positively correlated with the St. George's Respiratory Questionnaire score ( r  = 0.37, p  < 0.001) and negatively correlated with the percent predicted forced vital capacity (FVC% predicted) ( r  =  - 0.46, p  < 0.001) and the percent predicted diffusion capacity of the lung for carbon monoxide (Dlco% predicted) ( r  =  - 0.42, p  < 0.001). Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis suggested that the genes in the red module were primarily involved in inflammation and immune pathways. Based on Module Membership and Gene Significance, 32 candidate hub genes were selected in the red module to construct a protein-protein interaction network . Based on the identified DEGs and the degree of connectivity in the network, we identified three hub genes, including interleukin 6 ( IL6 ), suppressor of cytokine signaling-3 ( SOCS3 ), and serpin family E member 1 ( SERPINE1 ). In the GSE47460 dataset, Spearman correlation coefficients between Dlco% predicted and expression levels of IL6 , SERPINE1 , SOCS3 were -0.32, -0.41, and -0.46, respectively. Spearman correlation coefficients between FVC% predicted and expression levels of IL6 , SERPINE1 , SOCS3 were -0.29, -0.33, and -0.27, respectively. In the GSE24206 dataset, all three hub genes were upregulated in patients with advanced IPF., Conclusion: We identified three hub genes that negatively correlated with the lung function of IPF patients. Our results provide insights into the pathogenesis underlying the progressive disruption of lung function, and the identified hub genes may serve as biomarkers and potential therapeutictargets for the treatment of IPF patients., Competing Interests: The authors declare there are no competing interests., (©2020 Xia et al.)

Background: Colorectal cancer (CRC) is the third most common cancer in the world. The present study is aimed at identifying hub genes associated with the progression of CRC., Method: The data of the patients with CRC were obtained from the Gene Expression Omnibus (GEO) database and assessed by weighted gene co-expression network analysis (WGCNA), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses performed in R by WGCNA, several hub genes that regulate the mechanism of tumorigenesis in CRC were identified. Differentially expressed genes in the data sets GSE28000 and GSE42284 were used to construct a co-expression network for WGCNA. The yellow, black and blue modules associated with CRC level were filtered. Combining the co-expression network and the PPI network, 15 candidate hub genes were screened., Results: After validation using the TCGA-COAD dataset, a total of 10 hub genes (MT1X, MT1G, MT2A, CXCL8, IL1B, CXCL5, CXCL11, IL10RA, GZMB, KIT) closely related to the progression of CRC were identified. The expressions of MT1G, CXCL8, IL1B, CXCL5, CXCL11 and GZMB in CRC tissues were higher than normal tissues ( p -value < 0.05). The expressions of MT1X, MT2A, IL10RA and KIT in CRC tissues were lower than normal tissues ( p -value < 0.05)., Conclusions: By combinating with a series of methods including GO enrichment analysis, KEGG pathway analysis, PPI network analysis and gene co-expression network analysis, we identified 10 hub genes that were associated with the progression of CRC., Competing Interests: The authors declare that they have no competing interests., (© 2020 Meng et al.)

Background: Pancreatic cancer is one of the most common malignant cancers worldwide. Currently, the pathogenesis of pancreatic cancer remains unclear; thus, it is necessary to explore its precise molecular mechanisms., Methods: To identify candidate genes involved in the tumorigenesis and proliferation of pancreatic cancer, the microarray datasets GSE32676, GSE15471 and GSE71989 were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) between Pancreatic ductal adenocarcinoma (PDAC) and nonmalignant samples were screened by GEO2R. The Database for Annotation Visualization and Integrated Discovery (DAVID) online tool was used to obtain a synthetic set of functional annotation information for the DEGs. A PPI network of the DEGs was established using the Search Tool for the Retrieval of Interacting Genes (STRING) database, and a combination of more than 0.4 was considered statistically significant for the PPI. Subsequently, we visualized the PPI network using Cytoscape. Functional module analysis was then performed using Molecular Complex Detection (MCODE). Genes with a degree ≥10 were chosen as hub genes, and pathways of the hub genes were visualized using ClueGO and CluePedia. Additionally, GenCLiP 2.0 was used to explore interactions of hub genes. The Literature Mining Gene Networks module was applied to explore the cocitation of hub genes. The Cytoscape plugin iRegulon was employed to analyze transcription factors regulating the hub genes. Furthermore, the expression levels of the 13 hub genes in pancreatic cancer tissues and normal samples were validated using the Gene Expression Profiling Interactive Analysis (GEPIA) platform. Moreover, overall survival and disease-free survival analyses according to the expression of hub genes were performed using Kaplan-Meier curve analysis in the cBioPortal online platform. The relationship between expression level and tumor grade was analyzed using the online database Oncomine. Lastly, the eight snap-frozen tumorous and adjacent noncancerous adjacent tissues of pancreatic cancer patients used to detect the CDK1 and CEP55 protein levels by western blot., Conclusions: Altogether, the DEGs and hub genes identified in this work can help uncover the molecular mechanisms underlying the tumorigenesis of pancreatic cancer and provide potential targets for the diagnosis and treatment of this disease., Competing Interests: The authors declare there are no competing interests., (©2020 Jin et al.)

Background: Lung squamous cell carcinoma (LUSC) is a major subtype of lung cancer with limited therapeutic options and poor clinical prognosis., Methods: Three datasets (GSE19188, GSE33532 and GSE33479) were obtained from the gene expression omnibus (GEO) database. Differentially expressed genes (DEGs) between LUSC and normal tissues were identified by GEO2R, and functional analysis was employed using the Database for Annotation, Visualization and Integrated Discovery (DAVID) online tool. Protein-protein interaction (PPI) and hub genes were identified via the Search Tool for the Retrieval of Interacting Genes (STRING) and Cytoscape software. Hub genes were further validated in The Cancer Genome Atlas (TCGA) database. Subsequently, survival analysis was performed using the Kapla-Meier curve and Cox progression analysis. Based on univariate and multivariate Cox progression analysis, a gene signature was established to predict overall survival. Receiver operating characteristic curve was used to evaluate the prognostic value of the model., Results: A total of 116 up-regulated genes and 84 down-regulated genes were identified. These DEGs were mainly enriched in the two pathways: cell cycle and p53 signaling way. According to the degree of protein nodes in the PPI network, 10 hub genes were identified. The mRNA expression levels of the 10 hub genes in LUSC were also significantly up-regulated in the TCGA database. Furthermore, a novel seven-gene signature (FLRT3, PPP2R2C, MMP3, MMP12, CAPN8, FILIP1 and SPP1) from the DEGs was constructed and acted as a significant and independent prognostic signature for LUSC., Conclusions: The 10 hub genes might be tightly correlated with LUSC progression. The seven-gene signature might be an independent biomarker with a significant predictive value in LUSC overall survival., Competing Interests: The authors declare that they have no competing interests., (© 2020 Ma et al.)