基于转录组的肝豆状核变性调控网络的构建和分析

本研究旨在利用生物信息学方法构建经铜诱导的ATP7B基因敲除HepG2细胞系的转录调控网络。探讨关键转录因子在肝豆状核变性发生、发展中的潜在作用机制。收集公共基因表达数据库(gene expression omnibus, GEO)中包含野生型、ATP7B基因敲除型、铜诱导的野生型和铜诱导的ATP7B基因敲除型HepG2细胞系数据。筛选由铜诱导产生的差异表达基因(differentially expressed genes,DEGs)后进行基因本体论(gene ontology,GO)、京都基因和基因组百科全书(Kyoto encyclopedia of genes and genomes, KEGG)富集分析。基于蛋白相互作用网络,识别疾病关键基因和功能模块,并对关键功能模块中的基因进行富集分析。最后,构建转录调控网络,筛选核心转录因子。共筛选出1 034个差异表达基因,其中上调525个,下调509个。上、下调关键功能模块分别包括了3785个和3931个基因。关键功能模块中的基因主要定位于细胞-基质连接、染色体、剪接复合体、核糖体等区域,共同参与了mRNA加工、组蛋白修饰、RNA剪切...     

Identification and validation of key genes associated with non-small-cell lung cancer

Non-small-cell lung cancer (NSCLC) is one of the main causes of death induced by cancer globally. However, the molecular aberrations in NSCLC patients remain unclearly. In the present study, four messenger RNA microarray datasets (GSE18842, GSE40275, GSE43458, and GSE102287) were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) between NSCLC tissues and adjacent lung tissues were obtained from GEO2R and the overlapping DEGs were identified. Moreover, functional and pathway enrichment were performed by Funrich, while the protein–protein interaction (PPI) network construction were obtained from STRING and hub genes were visualized and identified by Cytoscape software. Furthermore, validation, overall survival (OS) and tumor staging analysis of selected hub genes were performed by GEPIA. A total of 367 DEGs (95 upregulated and 272 downregulated) were obtained through gene integration analysis. The PPI network consisted of 94 nodes and 1036 edges in the upregulated DEGs and 272 nodes and 464 edges in the downregulated DEGs, respectively. The PPI network identified 46 upregulated and 27 downregulated hub genes among the DEGs, and six (such as CENPE, NCAPH, MYH11, LRRK2, HSD17B6, and A2M) of that have not been identified to be associated with NSCLC so far. Moreover, the expression differences of the mentioned hub genes were consistent with that in lung adenocarcinoma and lung squamous cell carcinoma in the TCGA database. Further analysis showed that all the six hub genes were associated with tumor staging except MYH11, while only the upregulated DEG CENPE was associated with the worse OS of patients with NSCLC. In conclusion, the current study showed that CENPE, NCAPH, MYH11, LRRK2, HSD17B6, and A2M might be the key genes contributed to tumorigenesis or tumor progression in NSCLC, further functional study is needed to explore the involved mechanisms.     

Identification of lncRNA–miRNA–mRNA regulatory network associated with epithelial ovarian cancer cisplatin-resistant

To construct a long noncoding RNA (lncRNA)–microRNA (miRNA)–messenger RNA (mRNA) regulatory network related to epithelial ovarian cancer (EOC) cisplatin-resistant, differentially expressed genes (DEGs), differentially expressed lncRNAs (DELs), and differentially expressed miRNAs (DEMs) between MDAH and TOV-112D cells lines were identified. Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were conducted to analyze the biological functions of DEGs. Downstream mRNAs or upstream lncRNAs for miRNAs were analyzed at miRTarBase 7.0 or DIANA-LncBase V2, respectively. A total of 485 significant DEGs, 85 DELs, and 5 DEMs were identified. Protein–protein interaction (PPI) network of DEGs contrains 81 nodes and 141 edges was constructed, and 25 hub genes related to EOC cisplatin-resistant were identified. Subsequently, a lncRNA–miRNA–mRNA regulatory network contains 4 lncRNAs, 4 miRNAs, and 35 mRNAs was established. Taken together, our study provided evidence concerning the alteration genes involved in EOC cisplatin-resistant, which will help to unravel the mechanisms underlying drug resistant.     

Identification of key biomarkers in diabetic nephropathy via bioinformatic analysis

Diabetic nephropathy (DN) is a major cause of end-stage renal disease. Although intense efforts have been made to elucidate the pathogenesis, the molecular mechanisms of DN remain to be clarified. To identify the candidate genes in the progression of DN, microarray datasets GSE30122, GSE30528, and GSE47183 were downloaded from the Gene Expression Omnibus database. The differentially expressed genes (DEGs) were identified, and function enrichment analyses were performed. The protein-protein interaction network was constructed and the module analysis was performed using the Search Tool for the Retrieval of Interacting Genes and Cytoscape. A total of 61 DEGs were identified. The enriched functions and pathways of the DEGs included glomerulus development, extracellular exosome, collagen binding, and the PI3K-Akt signaling pathway. Fifteen hub genes were identified and biological process analysis revealed that these genes were mainly enriched in acute inflammatory response, inflammatory response, and blood vessel development. Correlation analysis between unexplored hub genes and clinical features of DN suggested that COL6A3, MS4A6A,PLCE1, TNNC1, TNNI1, TNN2, and VSIG4 may involve in the progression of DN. In conclusion, DEGs and hub genes identified in this study may deepen our understanding of molecular mechanisms underlying the progression of DN, and provide candidate targets for diagnosis and treatment of DN.     

Analysis of the potential pathways and target genes in spinal cord injury using bioinformatics methods

Spinal cord injury (SCI) remains to be the most devastating type of trauma for patients because of long lasting disability and limited response to the acute drug administration and efforts at rehabilitation. With the purpose to identify potential targets for SCI treatment and to gain more insights into the mechanisms of SCI, the microarray data of GSE2270, including 119 raphe magnus (RM) samples and 125 sensorimotor cortex (SMTC) samples, was downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were screened in RM group and SMTC group compared with their corresponding controls, respectively. A protein–protein interaction (PPI) network was constructed based on the common DEGs identified in both RM group and SMTC group. Gene ontology (GO) and pathway enrichment analyses of the overlapping DEGs were performed. Furthermore, the common DEGs enriched in each pathway were analyzed to identify significant regulatory elements. Totally, 173 overlapping DEGs (130 up-regulated and 43 down-regulated) were identified in both RM and SMTC samples. These overlapping DEGs were enriched in different GO terms. Pathway enrichment analysis revealed that DEGs were mainly related to inflammation and immunity. CD68 molecule (CD68) was a hub protein in the PPI network. Moreover, the regulatory network showed that ras-related C3 botulinum toxin substrate 2 (RAC2), CD44 molecule (CD44), and actin related protein 2/3 complex (ARPC1B) were hub genes. RAC2, CD44, and ARPC1B may be significantly involved in the pathogenesis of SCI by participating significant pathways such as extracellular matrix-receptor signaling pathway and Toll-like receptor signaling pathway.     

Identification of key genes and pathways in benign prostatic hyperplasia

Benign prostatic hyperplasia (BPH) is one of the most common causes of lower urinary tract symptoms (LUTS) in elderly man. However, the underlying molecular mechanisms of BPH have not been completely elucidated. We identified the key genes and pathways by using analysis of Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were identified using edgeR. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed for the DEGs by Database for Annotation, Visualization and Integrated Discovery (DAVID) database and ConsensusPathDB, respectively. Then, protein–protein interaction (PPI) networks were established by the Search Tool for the Retrieval of Interacting Genes (STRING) database and visualized by Cytoscape software. Finally, we identified 660 DEGs ultimately including 268 upregulated genes and 392 downregulated genes. GO analysis revealed that DEGs were mainly enriched in extracellular exosome, identical protein binding, mitochondrial adenosine triphosphate (ATP) synthesis coupled proton transport, extracelluar matrix, focal adhesion, cytosol, Golgi apparatus, cytoplasm, protein binding, and Golgi membrane. Focal adhesion pathway, FoxO signaling pathway, and autophagy pathway were selected. Ubiquitin-conjugating enzyme E2 C (UBE2C), serine/threonine kinase (AKT1), mitogen-activated protein kinase 1 (MAPK1), cyclin B1 (CCNB1), polo-like kinase 1 (PLK1) were filtrated as the hub genes according to the degree of connectivity from the PPI network. The five hub genes including UBE2C, AKT1, MAPK1, CCNB1, PLK1 may play key roles in the pathogenesis of benign prostatic hyperplasia (BPH). Focal adhesion pathway, FoxO signaling pathway, and autophagy pathway may be crucial for the progression of BPH.     

Analyzing time-series microarray data reveals key genes in spinal cord injury

Although many scholars have utilized high-throughput microarrays to delineate gene expression patterns after spinal cord injury (SCI), no study has evaluated gene changes in raphe magnus (RM) and somatomotor cortex (SMTC), two areas in brain primarily affected by SCI. In present study, we aimed to analyze the differentially expressed genes (DEGs) of RM and SMTC between SCI model and sham injured control at 4, 24 h, 7, 14, 28 days, and 3 months using microarray dataset GSE2270 downloaded from gene expression omnibus and unpaired significance analysis of microarray method. Protein–protein interaction (PPI) network was constructed for DEGs at crucial time points and significant biological functions were enriched using DAVID. The results indicated that more DEGs were identified at 14 days in RM and at 4 h/3 months in SMTC after SCI. In the PPI network for DEGs at 14 days in RM, interleukin 6, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), FBJ murine osteosarcoma viral oncogene homolog (FOS), tumor necrosis factor, and nuclear receptor subfamily 3, group C, member 1 (glucocorticoid receptor) were the top 5 hub genes; In the PPI network for DEGs at 3 months in SMTC, the top 5 hub genes were ubiquitin B, Ras‐related C3 botulinum toxin substrate 1 (rho family, small GTP binding protein Rac1), FOS, Janus kinase 2 and vascular endothelial growth factor A. Hedgehog and Wnt signaling pathways were the top 2 significant pathways in RM. These hub DEGs and pathways may be underlying therapeutic targets for SCI.     

Identification of hub genes and pathways in adrenocortical carcinoma by integrated bioinformatic analysis

Adrenocortical carcinoma (ACC), a rare malignant neoplasm originating from adrenal cortical cells, has high malignancy and few treatments. Therefore, it is necessary to explore the molecular mechanism of tumorigenesis, screen and verify potential biomarkers, which will provide new clues for the treatment and diagnosis of ACC. In this paper, three gene expression profiles (GSE10927, GSE12368 and GSE90713) were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were obtained using the Limma package. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were enriched by DAVID. Protein‐protein interaction (PPI) network was evaluated by STRING database, and PPI network was constructed by Cytoscape. Finally, GEPIA was used to validate hub genes’ expression. Compared with normal adrenal tissues, 74 up‐regulated DEGs and 126 down‐regulated DEGs were found in ACC samples; GO analysis showed that up‐regulated DEGs were enriched in organelle fission, nuclear division, spindle, et al, while down‐regulated DEGs were enriched in angiogenesis, proteinaceous extracellular matrix and growth factor activity; KEGG pathway analysis showed that up‐regulated DEGs were significantly enriched in cell cycle, cellular senescence and progesterone‐mediated oocyte maturation; Nine hub genes (CCNB1, CDK1, TOP2A, CCNA2, CDKN3, MAD2L1, RACGAP1, BUB1 and CCNB2) were identified by PPI network; ACC patients with high expression of 9 hub genes were all associated with worse overall survival (OS). These hub genes and pathways might be involved in the tumorigenesis, which will offer the opportunities to develop the new therapeutic targets of ACC.     

Identification of potential key genes associated with osteosarcoma based on integrated bioinformatics analyses

Due to high rates of metastasis and poor clinical outcomes for patients, it is important to study the pathomechanisms of osteosarcoma. However, due to the fact that osteosarcoma shows significant interindividual variation and high heterogeneity, the identification of differentially expressed genes (DEGs) at the population level cannot answer many important questions related to osteosarcoma tumorigenesis. Therefore, a new strategy to identify dysregulated genes in osteosarcoma samples is required. The aim of this study was to improve our understanding of osteosarcoma pathogenesis by identifying genes with universal aberrant expression in osteosarcoma samples. Because the relative expression ordering of genes is stable in normal bone tissues but is disrupted in osteosarcoma tissues, we used the RankComp algorithm to identify DEGs in normal and osteosarcoma tissue samples. We then calculated the dysregulation frequency for each gene. Genes with deregulation frequencies above 80% were deemed to be universal DEGs. Next, coexpression, pathway enrichment, and protein-protein interaction network analyses were performed to characterize the functions of these genes. From 188 samples of osteosarcoma obtained from four datasets measured on different platforms, 51 universal DEGs were identified, including 4 universally upregulated genes and 47 universally downregulated genes. Genes that were differentially coexpressed with these universal DEGs were found to be enriched in 46 cancer-related pathways. In addition, functional and network analyses showed that genes with high dysregulation frequencies were involved in cancer-related functions. Thus, the commonly aberrant genes identified in osteosarcoma tissues may be important targets for osteosarcoma diagnosis and therapy.     

Identification of EPHX2 and RMI2 as two novel key genes in cervical squamous cell carcinoma by an integrated bioinformatic analysis

Cervical cancer is the fourth most common malignancy in women worldwide and cervical squamous cell carcinoma (CESC) is the most common histological type of cervical cancer. The dysregulation of genes plays a significant role in cancer. In the present study, we screened out differentially expressed genes (DEGs) of CESC in the GSE63514 data set from the Gene Expression Omnibus database. An integrated bioinformatics analysis was used to select hub genes, as well as to investigate their related prognostic signature, functional annotation, methylation mechanism, and candidate molecular drugs. As a result, a total of 1907 DEGs were identified (944 were upregulated and 963 were downregulated). In the protein–protein interaction network, three hub modules and 30 hub genes were identified. And two hub modules and 116 hub genes were screened out from four CESC-related modules by the weighted gene coexpression network analysis. The gene ontology term enrichment analysis and Kyoto encyclopedia of genes and genomes pathway analysis were performed to better understand functions and pathways. Genes with a significant prognostic value were found by prognostic signature analysis. And there were five genes (EPHX2, CHAF1B, KIAA1524, CDC45, and RMI2) identified as significant CESC-associated genes after expression validation and survival analysis. Among them, EPHX2 and RMI2 were noted as two novel key genes for the CESC-associated methylation and expression. In addition, four candidate small molecule drugs for CESC (camptothecin, resveratrol, vorinostat, and trichostatin A) were defined. Further studies are required to explore these significant CESC-associated genes for their potentiality in diagnosis, prognosis, and targeted therapy.     

Fifteen hub genes associated with progression and prognosis of clear cell renal cell carcinoma identified by coexpression analysis

Renal cell carcinoma (RCC) is the most common type of renal tumor, and the clear cell renal cell carcinoma (ccRCC) is the most frequent subtype. In this study, our aim is to identify potential biomarkers that could effectively predict the prognosis and progression of ccRCC. First, we used The Cancer Genome Atlas (TCGA) RNA-sequencing (RNA-seq) data of ccRCC to identify 2370 differentially expressed genes (DEGs). Second, the DEGs were used to construct a coexpression network by weighted gene coexpression network analysis (WGCNA). Moreover, we identified the yellow module, which was strongly related to the histologic grade and pathological stage of ccRCC. Then, the functional annotation of the yellow module and single-samples gene-set enrichment analysis of DEGs were performed and mainly enriched in cell cycle. Subsequently, 18 candidate hub genes were screened through WGCNA and protein–protein interaction (PPI) network analysis. After verification of TCGA’s ccRCC data set, Gene Expression Omnibus (GEO) data set (GSE73731) and tissue validation, we finally identified 15 hub genes that can actually predict the progression of ccRCC. In addition, by using survival analysis, we found that patients of ccRCC with high expression of each hub gene were more likely to have poor prognosis than those with low expression. The receiver operating characteristic curve showed that each hub gene could effectively distinguish between localized and advanced ccRCC. In summary, our study indicates that 15 hub genes have great predictive value for the prognosis and progression of ccRCC, and may contribute to the exploration of the pathogenesis of ccRCC.     

Gene regulation network analysis reveals core genes associated with survival in glioblastoma multiforme

Glioblastoma multiforme (GBM) is a very serious mortality of central nervous system cancer. The microarray data from GSE2223 , GSE4058 , GSE4290 , GSE13276 , GSE68848 and GSE70231 (389 GBM tumour and 67 normal tissues) and the RNA‐seq data from TCGA‐GBM dataset (169 GBM and five normal samples) were chosen to find differentially expressed genes (DEGs). RRA (Robust rank aggregation) method was used to integrate seven datasets and calculate 133 DEGs (82 up‐regulated and 51 down‐regulated genes). Subsequently, through the PPI (protein‐protein interaction) network and MCODE/ cytoHubba methods, we finally filtered out ten hub genes, including FOXM1, CDK4, TOP2A, RRM2, MYBL2, MCM2, CDC20, CCNB2, MYC and EZH2, from the whole network. Functional enrichment analyses of DEGs were conducted to show that these hub genes were enriched in various cancer‐related functions and pathways significantly. We also selected CCNB2, CDC20 and MYBL2 as core biomarkers, and further validated them in CGGA, HPA and CCLE database, suggesting that these three core hub genes may be involved in the origin of GBM. All these potential biomarkers for GBM might be helpful for illustrating the important role of molecular mechanisms of tumorigenesis in the diagnosis, prognosis and targeted therapy of GBM cancer.     

Hub genes identification and association of key pathways with hypoxia in cancer cells: A bioinformatics analysis

Three human cancer cell lines (A549, HCT116, and HeLa) were used to investigate the molecular mechanisms and potential prognostic biomarkers associated with hypoxia. We obtained gene expression data from Gene Expression Omnibus (GEO) datasets GSE11704, GSE147384, and GSE38061, which included 5 hypoxic and 8 control samples. Using the GEO2R tool and Venn diagram software, we identified common differentially expressed genes (cDEGs). The cDEGs were then subjected to Gene ontology (GO) and Kyoto Encyclopedia of Gene and Genome (KEGG) pathway analysis by employing DAVID. The hub genes were identified from critical PPI subnetworks through CytoHuba plugin and these genes' prognostic significance and expression were verified using Kaplan-Meier analysis and Gene Expression Profiling Interactive Analysis (GEPIA), respectively. The research showed 676 common DEGs (cDEGs), with 207 upregulated and 469 downregulated genes. The STRING analysis showed 673 nodes and 1446 edges in the PPI network. We identified 4 significant modules and 19 downregulated hub genes. GO analysis revealed all of them were majorly involved in ribosomal large subunit assembly and biogenesis, rRNA processing, ribosome biogenesis, translation, RNA & protein binding frequently at the sites of nucleolus and nucleoplasm while 11 were significantly associated with a better prognosis of hypoxic tumors. Our research sheds light on the molecular mechanisms that underpin hypoxia in human cancer cell lines and identifies potential prognostic biomarkers for hypoxic tumors.