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

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

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

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

Identification of core genes and outcomes in hepatocellular carcinoma by bioinformatics analysis

Hepatocellular carcinoma (HCC) is the most common malignant liver disease in the world. However, the mechanistic relationships among various genes and signaling pathways are still largely unclear. In this study, we aimed to elucidate potential core candidate genes and pathways in HCC. The expression profiles GSE14520, GSE25097, GSE29721, and GSE62232, which cover 606 tumor and 550 nontumour samples, were downloaded from the Gene Expression Omnibus (GEO) database. Furthermore, HCC RNA-seq datasets were also downloaded from the Cancer Genome Atlas (TCGA) database. The differentially expressed genes (DEGs) were filtered using R software, and we performed gene ontology (GO) and Kyoto Encyclopedia of Gene and Genome (KEGG) pathway analysis using the online databases DAVID 6.8 and KOBAS 3.0. Furthermore, the protein-protein interaction (PPI) network complex of these DEGs was constructed by Cytoscape software, the molecular complex detection (MCODE) plug-in and the online database STRING. First, a total of 173 DEGs (41 upregulated and 132 downregulated) were identified that were aberrantly expressed in both the GEO and TCGA datasets. Second, GO analysis revealed that most of the DEGs were significantly enriched in extracellular exosomes, cytosol, extracellular region, and extracellular space. Signaling pathway analysis indicated that the DEGs had common pathways in metabolism-related pathways, cell cycle, and biological oxidations. Third, 146 nodes were identified from the DEG PPI network complex, and two important modules with a high degree were detected using the MCODE plug-in. In addition, 10 core genes were identified, TOP2A, NDC80, FOXM1, HMMR, KNTC1, PTTG1, FEN1, RFC4, SMC4, and PRC1. Finally, Kaplan-Meier analysis of overall survival and correlation analysis were applied to these genes. The abovementioned findings indicate that the identified core genes and pathways in this bioinformatics analysis could significantly enrich our understanding of the development and recurrence of HCC; furthermore, these candidate genes and pathways could be therapeutic targets for HCC treatment.     

Identifying potential biomarkers in hepatitis B virus infection and its response to the antiviral therapy by integrated bioinformatic analysis

The antiviral treatment efficacy varies among chronic hepatitis B (CHB) patients and the underlying mechanism is unclear. An integrated bioinformatics analysis was performed to investigate the host factors that affect the therapeutic responsiveness in CHB patients. Four GEO data sets (GSE54747, GSE27555, GSE66698 and GSE66699) were downloaded from the Gene Expression Omnibus (GEO) database and analysed to identify differentially expressed genes(DEGs). Enrichment analyses of the DEGs were conducted using the DAVID database. Immune cell infiltration characteristics were analysed by CIBERSORT. Upstream miRNAs and lncRNAs of hub DEGs were identified by miRWalk 3.0 and miRNet in combination with the MNDR platform. As a result, seventy-seven overlapping DEGs and 15 hub genes were identified including CCL5, CXCL9, MYH2, CXCR4, CD74, CCL4, HLA-DRB1, ACTA1, CD69, CXCL10, HLA-DRB5, HLA-DQB1, CXCL13, STAT1 and CKM. The enrichment analyses revealed that the DEGs were mainly enriched in immune response and chemokine signalling pathways. Investigation of immune cell infiltration in liver samples suggested significantly different infiltration between responders and non-responders, mainly characterized by higher proportions of CD8+ T cells and activated NK cells in non-responders. The prediction of upstream miRNAs and lncRNAs led to the identification of a potential mRNA-miRNA-lncRNA regulatory network composed of 2 lncRNAs (H19 and GAS5) and 5 miRNAs (hsa-mir-106b-5p, hsa-mir-17-5p, hsa-mir-20a-5p, hsa-mir-6720-5p and hsa-mir-93-5p) targeting CCL5 mRNA. In conclusion, our study suggested that host genetic factors could affect therapeutic responsiveness in CHB patients. The antiviral process might be associated with the chemokine-mediated immune response and immune cell infiltration in the liver microenvironment.     

Identification of a six-gene prognostic signature for oral squamous cell carcinoma

Oral squamous cell carcinoma (OSCC) is one of the most common types of malignancies worldwide, and its morbidity and mortality have increased in the near term. Consequently, the purpose of the present study was to identify the notable differentially expressed genes (DEGs) involved in their pathogenesis to obtain new biomarkers or potential therapeutic targets for OSCC. The gene expression profiles of the microarray datasets GSE85195, GSE23558, and GSE10121 were obtained from the Gene Expression Omnibus (GEO) database. After screening the DEGs in each GEO dataset, 249 DEGs in OSCC tissues were obtained. Kyoto Encyclopedia of Genes and Genomes and Gene Ontology pathway enrichment analysis was employed to explore the biological functions and pathways of the above DEGs. A protein–protein interaction network was constructed to obtain a central gene. The corresponding total survival information was analyzed in patients with oral cancer from The Cancer Genome Atlas (TCGA). A total of six candidate genes (CXCL10, OAS2, IFIT1, CCL5, LRRK2, and PLAUR) closely related to the survival rate of patients with oral cancer were identified, and expression verification and overall survival analysis of six genes were performed based on TCGA database. Time-dependent receiver operating characteristic curve analysis yields predictive accuracy of the patient's overall survival. At the same time, the six genes were further verified by quantitative real-time polymerase chain reaction using samples obtained from the patients recruited to the present study. In conclusion, the present study identified the prognostic signature of six genes in OSCC for the first time via comprehensive bioinformatics analysis, which could become potential prognostic markers for OCSS and may provide potential therapeutic targets for tumors.     

应用GSEA和WGCNA方法分析细菌性败血症宿主关键差异基因及意义

【目的】采用生物信息学方法分析公共数据库来源的细菌性败血症患者全血转录组学表达谱,探讨细菌败血症相关的宿主关键差异基因及意义。【方法】基于GEO数据库中GSE80496和GSE72829全血转录组基因数据集,采用GEO2R、基因集富集分析(GSEA)联用加权基因共表达网络分析(WGCNA)筛选细菌性败血症患者相比健康人群显著改变的差异基因,通过R软件对交集基因进行GO功能分析和KEGG富集分析。同时,通过String 11.0和Cytoscape分析枢纽基因,验证枢纽基因在数据集GSE72809(Health组52例,Definedsepsis组52例)全血标本中的表达情况,并探讨婴儿性别、月(胎)龄、出生体重、是否接触抗生素等因素与靶基因表达谱间的关系。【结果】分析GSE80496和GSE72829数据集分别筛选得到932个基因和319个基因,联合WGCNA枢纽模块交集得到与细菌性败血症发病相关的10个枢纽基因(MMP9、ITGAM、CSTD、GAPDH、PGLYRP1、FOLR3、OSCAR、TLR5、IL1RN和TIMP1);GSEA分析获得关键通路(氨基酸糖类-核糖代谢、PPAR信号通路、聚糖生物合成通路、自噬调控通路、补体、凝血因子级联反应、尼古丁和烟酰胺代谢、不饱和脂肪酸生物合成和阿尔兹海默症通路)及生物学过程(类固醇激素分泌、腺苷酸环化酶的激活、细胞外基质降解和金属离子运输)。【结论】本项研究通过GEO2R、GSEA联用WGCNA分析,筛选出与细菌性败血症发病相关的2个枢纽模块、10个枢纽基因以及一些关键信号通路和生物学过程,可为后续深入研究细菌性败血症致病机制奠定理论依据。     

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.     

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 key candidate genes and pathways in multiple myeloma by integrated bioinformatics analysis

Multiple myeloma (MM) is a common hematologic malignancy for which the underlying molecular mechanisms remain largely unclear. This study aimed to elucidate key candidate genes and pathways in MM by integrated bioinformatics analysis. Expression profiles GSE6477 and GSE47552 were obtained from the Gene Expression Omnibus database, and differentially expressed genes (DEGs) with p < .05 and [logFC] > 1 were identified. Functional enrichment, protein–protein interaction network construction and survival analyses were then performed. First, 51 upregulated and 78 downregulated DEGs shared between the two GSE datasets were identified. Second, functional enrichment analysis showed that these DEGs are mainly involved in the B cell receptor signaling pathway, hematopoietic cell lineage, and NF-kappa B pathway. Moreover, interrelation analysis of immune system processes showed enrichment of the downregulated DEGs mainly in B cell differentiation, positive regulation of monocyte chemotaxis and positive regulation of T cell proliferation. Finally, the correlation between DEG expression and survival in MM was evaluated using the PrognoScan database. In conclusion, we identified key candidate genes that affect the outcomes of patients with MM, and these genes might serve as potential therapeutic targets.     

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

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

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.     

Protein–protein interaction and SNP analysis in intraductal papillary mucinous neoplasm

Intraductal papillary mucinous neoplasm (IPMN) is a type of tumor that grows within the pancreatic ducts. It is a progress from hyperplasia to intraductal adenoma (IPMA), to noninvasive carcinoma, and ultimately to invasive carcinoma (IPMC). The objective of this study was to explore the molecular mechanism of the progression from IPMA to IPMC. By using the GSE19650 affymetrix microarray data accessible from Gene Expression Omnibus (GEO) database, we first identified the differentially expressed genes (DEGs) between IPMA and IPMC, followed by the protein–protein interaction and single-nucleotide polymorphism (SNP) analysis of the DEGs. Our study identified thousands of DEGs which involved regulation of cell cycle and apoptosis in this progression from IPMA to IPMC. Protein–protein interaction network construction found that MYC, IL6ST, NR3C1, CREBBP, GATA1 and LRP1 might play an important role in the progression. Furthermore, the SNP analysis confirmed the association between BRAC1 and pancreas cancer. In conclusion, our data provide a comprehensive bioinformatics analysis of genes and pathways which may be involved in the progression of IPMN from IPMA to IPMC.     

Identification of Biomarkers for Early Diagnosis of Acute Myocardial Infarction

Acute myocardial infarction (AMI) is a common disease with serious consequences in mortality and cost. Here we aim to screen the differentially expressed genes (DEGs) as biomarkers for early diagnosis of AMI. The microarray data of AMI was downloaded from Gene Expression Omnibus (GEO), including two independent types of research GSE66360 and GSE62646. The DEGs between control and processed samples were screened out by using limma package. Meanwhile, we performed functional analysis of GO terms and/or KEGG pathways to observe the enriched pathways of the DEGs. Finally, regression analysis of raw data was performed by using packet affyPLM in R language. Dataset GSE62646 contained 53 DEGs (FC log2>1 and P value <0.05) between first‐day samples from 28 STEMI patients and control samples from 14 patients without myocardial infarction history. There were 12 up‐regulated and 41 down‐regulated genes, 35 DEGs annotated. In GSE66360, a total of 3034 DEGs between 32 AMI patients and 38 healthy persons were obtained, including 1861 up‐regulated and 1173 down‐regulated DEGs. The comparison of the DEGs in two datasets revealed four common up‐regulated genes (EGR1, STAB1, SOCS3, TMEM176A). In enrichment analysis, STAB1, SOCS3, EGR1 involved in metabolic regulation and signaling pathways related to coronary artery disease with a characteristic change (P < 0.05). The DEGs, especially the four up‐regulated common genes, could serve as biomarkers for early diagnosis of AMI. Additionally, the relative biological pathways these DEGs enriched in might provide a good reference to explore the molecular expression mechanism of AMI. J. Cell. Biochem. 119: 650–658, 2018. © 2017 Wiley Periodicals, Inc.     

IgA 肾病患者与膜性肾病患者外周血单核细胞mRNA分析

为了寻找诊断、鉴别IgA肾病(IgAN)和膜性肾病(MN)的血液特异性标记物,利用公共数据库中的IgAN和MN患者的外周血单核细胞(PBMCs)的转录组表达谱数据集识别特异性生物标记物,为诊断和鉴别提供简便、可靠的依据补充。从公共基因表达数据库(GEO)下载IgAN患者组(n=15)和MN患者组(n=8)芯片数据集,筛选前250个差异表达基因(DEGs)。通过分析筛选关键基因和途径,进行基因本体(GO)富集分析、京都基因与基因组百科全书(KEGG)通路分析和蛋白质与蛋白质相互作用关系(PPI)分析等进一步了解DEGs。通过分析共发现75个显著DEGs,其中73个上调基因,2个下调基因。GO富集分析的生物学过程(BP)主要包括蛋白质转运、内溶酶体到溶酶体转运、趋化因子介导的信号通路作用等。显著富集差异表达基因KEGG通路分析包括Endocytosis和Hepatitis B的相关信号通路。PPI筛选出EPS15、STAT4、CCL2、SUN2、SEC24C、SEC31A、GOLGB1、F2R,RAB12和PTK2B等关键基因。成功筛选出核心差异表达基因,为IgAN和MN的诊断和鉴别提供简便、可靠的依据补充,甚至提供治疗的新靶点。     

Bioinformatics analysis of the biological changes involved in the osteogenic differentiation of human mesenchymal stem cells

The mechanisms underlying the osteogenic differentiation of human bone marrow mesenchymal stem cells (hBMSCs) remain unclear. In the present study, we aimed to identify the key biological processes during osteogenic differentiation. To this end, we downloaded three microarray data sets from the Gene Expression Omnibus (GEO) database: GSE12266, GSE18043 and GSE37558. Differentially expressed genes (DEGs) were screened using the limma package, and enrichment analysis was performed. Protein‐protein interaction network (PPI) analysis and visualization analysis were performed with STRING and Cytoscape. A total of 240 DEGs were identified, including 147 up‐regulated genes and 93 down‐regulated genes. Functional enrichment and pathways of the present DEGs include extracellular matrix organization, ossification, cell division, spindle and microtubule. Functional enrichment analysis of 10 hub genes showed that these genes are mainly enriched in microtubule‐related biological changes, that is sister chromatid segregation, microtubule cytoskeleton organization involved in mitosis, and spindle microtubule. Moreover, immunofluorescence and Western blotting revealed dramatic quantitative and morphological changes in the microtubules during the osteogenic differentiation of human adipose‐derived stem cells. In summary, the present results provide novel insights into the microtubule‐ and cytoskeleton‐related biological process changes, identifying candidates for the further study of osteogenic differentiation of the mesenchymal stem cells.     

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.     

Novel key genes in triple-negative breast cancer identified by weighted gene co-expression network analysis

Triple-negative breast cancer (TNBC) is a special subtype of breast cancer (BC) with poor prognosis. Although some molecular mechanisms of TNBC have been elucidated, the efficacy of current treatments is limited. Therefore, it is urgently demanded to screen for novel biomarkers and drug targets for TNBC. In this study, we obtained four independent data sets (GSE76250, GSE31448, GSE43358, and METABRIC) from the Gene Expression Omnibus (GEO) database and the cBioPortal website. In the GSE76250 data set, 890 differentially expressed genes were identified and weighted gene co-expression network analysis was performed based on them. Then, two preserved modules associated with the KI67 score were detected. Gene ontology and pathway enrichment analyses showed genes in the modules participated in some cancer-related biological processes or pathways. Non-SMC condensin I complex subunit G (NCAPG) and ATP-binding cassette subfamily A member 9 (ABCA9) were identified as hub genes of the modules, and the significance of hub genes was validated in the GSE43358 data set. Finally, their prognostic value was assessed by survival analysis. These findings suggested that NCAPG and ABCA9 may be the key genes of TNBC. Moreover, ABCA9 was first reported in TNBC. They deserved further studies.     

Immune-related chemotactic factors were found in acute coronary syndromes by bioinformatics

DNA microarray data for thrombus-related leukocyte from patients with acute coronary syndrome (ACS) was analyzed to acquire key genes associated with ACS. Microarray data set GSE19339, including four ACS patients’ samples and four normal samples, were downloaded from Gene Expression Omnibus database. Raw data was pre-processed and differentially expressed genes (DEGs) were identified by Affy packages of R. The interaction network was established with STRING. DrugBank was retrieved to obtain relevant small molecules. A total of 487 differentially expressed genes were identified as DEGs between normal and disease samples. Among which, ten up-regulated genes belonging to chemokine family (CCL2, CCR1, CXCL3, CXCL2, CCL8, CXCL11, CCL7, IL10, CCL22 and CCL20) were related to inflammatory response. In addition, two inhibitors of CCL2 (L-Mimosine) were retrieved from the DrugBank database. Considering the roles of inflammatory response in the progression of ACS and the functions of the ten up-regulated genes, we speculated that these genes might be related to ACS. Moreover, the inhibitors could provide guidelines for future drug design acting on these genes.