肺腺癌诊断标志物筛选及免疫细胞浸润分析

用生物信息学方法筛选肺腺癌(Lung adenocarcinoma,LUAD)的诊断生物标志物,并分析肺腺癌中免疫细胞浸润情况。从GEO和TCGA数据库下载肺腺癌的表达数据集,利用R软件筛选肺腺癌与正常肺组织间的差异表达基因(DEGs),使用DAVID网站对DEGs进行GO及KEGG富集分析,使用STRING及Cytoscape等工具对DEGs构建蛋白相互作用网络并筛选hub基因;利用Kaplan-Meier法对DEGs进行生存分析,并对hub基因进行ROC分析筛选诊断生物标志物,利用GSEA预测有预后价值的基因参与的信号通路;并用Cibersort软件反卷积算法分析肺腺癌中免疫细胞浸润情况。共得到肺腺癌的234个DEGs,这些基因主要参与信号转导、物质代谢、免疫反应等相关信号通路;构建PPI网络筛选出的20个hub基因中8个存在预后价值(CCNA2、DLGAP5、HMMR、MMP1、MMP9、MMP13、SPP1、TOP2A),ROC分析中DLGAP5、SPP1值分别是0.703、0.706;DLGAP5、SPP1基因表达水平与肺腺癌组织浆细胞、未活化的CD4⁺记忆细胞、调节T细胞、巨噬细胞M0、M1、M2及中性粒细胞浸润密切相关(P＜0.05)。肺腺癌中DLGAP5、SPP1具有较高诊断价值且参与肺腺癌组织免疫细胞浸润;DLGAP5、SPP1基因可作为肺腺癌诊断的生物标志物,可为肺腺癌的靶向治疗研究提供新思路。     

强直性脊柱炎差异表达基因的生物信息学分析

目的探讨强直性脊柱炎(AS)患者差异表达基因,并基于差异基因探讨强直性脊柱炎发病相关的可能生物学过程和信号通路。方法检索基因表达谱数据库(GEO)并筛选AS相关基因表达谱数据集。应用GEO在线分析功能GEO2R分析AS组和正常对照组的差异表达基因,用Cytoscape软件clueGO插件进行基因本体论和京都基因与基因组百科全书分析,采用String蛋白-蛋白相互作用(PPI)数据库分析差异表达基因编码蛋白间的相互作用;应用Cytoscape绘制蛋白相互作用网络图,并软件筛选信号通路关键基因分析。结果选取AS患者全血表达数据集GSE25101为研究对象,分析获得差异表达基因72个。72个差异表达基因分子功能主要为参与高迁移率族盒染色体蛋白1(HMGB1)转导机制;生物学过程主要富集于巨噬细胞迁移、骨髓细胞凋亡过程、线粒体呼吸链复合体装配、ATP合成偶联电子传输、线粒体ATP合成耦合电子输运等;细胞成分主要富集于呼吸链复合体、线粒体呼吸体等。信号通路富集于氧化磷酸化信号通路和帕金森综合征相关信号通路。PPI网络经过cytohubba插件筛选,ATP5J、NDUFS4、UQCRB、UQCRH、NDUFB3、COX7B、LSM3、ATP5EP2、ENY2、PSMA4被筛选为网络中的核心基因。结论通过生物信息学方法进行预测了AS的潜在机制,并筛选出10个潜在的与AS相关的重要分子,其中氧化磷酸化可能在AS发病机制中发挥了重要的作用。     

基于SNP共表达网络肝癌分子分型及预后分析

基于SNP突变数据与mRNA表达谱关联分析,构建一种肝癌分子分型方法并对比不同分型预后的差异,并对不同分型肝癌的发生发展机制进一步研究。首先通过TCGA数据库收集359例肝细胞癌患者的SNP突变数据和mRNA表达数据,采用Wilcoxon秩和检验,筛选突变后差异表达基因,并通过生物信息学工具String和Cytoscape 构建差异表达基因的蛋白互作网络,筛选连接度最高的10个Hub基因。利用Consensus Cluster Plus软件包,基于Hub基因mRNA表达水平构建NMF分子分型模型,再结合生存数据评估各分型患者的预后。最后利用加权基因共表达网络分析(WGCNA),识别与肝癌分子分型相关的模块,并针对关键模块的基因进行通路富集,从而对不同分型肝癌的基因表达谱进行比较。结果:NMF模型将肝癌分为高危、低危2个分型,其中CDKN2A和FOXO1基因对分型贡献度高。生存分析显示低危组患者的生存情况显著优于高危组,高危组富集多个与肿瘤细胞侵蚀、转移、复发过程相关的信号通路,低危组则与细胞周期和胰液分泌相关。本研究在无先验性信息的前提下,基于突变后显著差异表达的Hub基因表达水平构建的肝癌分子分型对肝癌患者预后评估具有一定的指导意义,其中CDKN2A和FOXO1突变是肝癌患者的不良预后因素,针对二者的靶向药研发,可能为肝癌患者提供新的治疗策略。     

胰腺癌中CDK1的表达与预后的生物信息学分析

为探讨胰腺癌的发病机制并为胰腺癌的防治提供生物信息学依据,用GEO2R在线工具分析GSE16515中胰腺癌患者肿瘤组织和相应正常组织的差异表达基因(differentially expressed genes, DEGs),通过DAVID数据库对DEGs进行GO分析和KEGG通路富集分析,然后通过STRING数据库构建蛋白质相互作用(protein-protein interaction, PPI)网络,用Cytoscape软件进行关键基因(hub基因)筛选和功能模块分析,并在GEPIA数据库对hub基因进行验证,用CCLE数据库检测靶基因在胰腺癌组织及细胞系中的表达水平。分析结果显示胰腺癌中筛选出的376个DEGs主要涉及细胞周期、p53信号通路、蛋白质消化吸收、ECM-受体相互作用、PI3K-Akt信号通路、血小板激活信号通路。GEPIA数据库验证结果显示10个hub基因均在胰腺癌组织中高表达,其中8个hub基因与胰腺癌患者的不良预后有关。CCLE数据库检测结果显示周期蛋白依赖性激酶1 (cyclin-dependent kinase 1, CDK1)在胰腺癌组织和细胞中均有较高的表达水平。本研究结果表明CDK1可能与胰腺癌的发生发展最为相关,为进一步探究胰腺癌的发病机制提供了生物信息学依据。     

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的诊断和鉴别提供简便、可靠的依据补充,甚至提供治疗的新靶点。     

肾透明细胞癌Caki-1细胞系差异表达基因的生物信息学分析

目的比较肾透明细胞癌Caki-1细胞系与正常肾上皮细胞系ASE-5063中的差异表达基因（DEGs）,寻找潜在的肾透明细胞癌特异性分子标志物。 方法利用GEO数据库自带的GEO2R在线分析工具分析基因芯片GSE78179,将筛选出的DEGs分别导入Metascape、STRING以及Cytoscape进行综合分析并筛选出核心基因。最后使用FunRich等软件对筛选出的核心基因进行GO和KEGG富集分析。 结果共筛选出562个DEGs,其中上调基因345个,下调基因217个。进一步使用MCODE筛选出36个关键基因,GO功能分析发现这些基因与细胞粘附分子活性、趋化因子活性、细胞通讯和信号转导等密切相关;KEGG通路富集结果则表明差异基因主要集中在趋化因子信号通路、TNF信号通路以及NF-κB信号通路等多种与肿瘤相关的通路上。 结论运用生物信息学方法筛选出肾透明细胞癌Caki-1细胞系中DEGs,其中数个核心基因广泛参与多种肿瘤的病理进程,但尚未在肾透明细胞癌有相关研究报道,提示其可能是治疗肾透明细胞癌的潜在靶点。     

GWGS芯片整合分析识别外周血细胞中与骨质疏松症相关的枢纽基因

为了对骨质疏松症基因芯片数据集进行整合分析并识别出外周血细胞中与骨质疏松症相关的枢纽基因,通过检索GEO和ArrayExpress数据库获得骨质疏松症相关的表达谱芯片数据集;运用GWGS (genome-wide global significance)方法对纳入的数据集进行整合分析,筛选出差异表达基因(differentially expressed genes, DEGs);然后,运用GO (gene ontology)富集分析和KEGG (kyoto encyclopedia of genes and genomes)通路富集分析对差异表达基因进行功能注释,并建立蛋白质相互作用(protein-protein interaction, PPI)网络,筛选出骨质疏松症相关的枢纽基因。公共数据库检索得到3个符合纳入排除标准的研究集, GWGS整合分析筛选出排序前200的DEGs,这些基因主要富集的GO条目为脂多糖的细胞反应、凋亡过程和炎症反应,与骨质疏松症相关的KEGG富集通路为破骨细胞分化等。PPI分析进一步检测到与骨质疏松症相关的10个枢纽基因,其中9个基因已有研究报道和骨质疏松症的发生发展相关,而ELANE基因还未有研究报道与骨质疏松症有关。ELANE基因同时在人的骨髓组织、小鼠骨髓和骨组织中高表达,这个基因很可能与骨质疏松症有潜在的联系。本研究的结果将有助于进一步理解骨质疏松症的分子致病机理。     

基于生物信息学分析的宫颈鳞状细胞癌预后相关基因的筛选

为了分析宫颈鳞状细胞癌(cervical squamous cell carcinoma, CESC)与正常组织中的差异表达基因(differentially expressed genes, DEGs),鉴定与CESC预后相关的关键基因,从GEO和TCGA数据库下载CESC的基因表达谱数据,利用R软件筛选CESC组织与正常组织中的DEGs,并对这些DEGs开展功能和通路富集分析;然后构建蛋白质-蛋白质相互作用(protein-protein interaction, PPI)网络,筛选出关键(hub)基因;最后对hub基因进行LASSO COX回归及总体生存率(overall survival, OS)分析。研究共筛选出167个DEGs,这些基因主要涉及染色体分离、DNA复制等生物过程,介导染色质结合、G蛋白偶联受体结合等分子功能,富集于染色体区域、纺锤体和MCM复合体。GSEA分析结果显示,富集的通路主要涉及DNA复制和细胞周期信号通路。此外,从PPI网络中筛选出20个hub基因, LASSO COX回归结果显示MAD2L1、ZWINT、RRM2、TTK、CDC6、PBK、TOP2A、KIF11、KIF20A、NCAPG、NUSAP1、CCNB1及CDK1与CESC患者的预后相关; Kaplan-Meier曲线显示, ZWINT、DTL、CCNB1、CDC6、TOP2A、CDK1、PBK、RFC4及NUSAP1的m RNA表达水平与CESC患者生存预后相关。本研究结果表明, ZWINT、CDC6、PBK、TOP2A、NUSAP1、CCNB1和CDK1为CESC的预后关键基因,为阐明CESC的分子机制提供了理论依据。     

基于生物信息学分析的非小细胞肺癌诊断预后相关基因的筛选

为了从基因层面探讨非小细胞肺癌(non-small cell lung cancer, NSCLC)发生发展的内在机制,筛选与NSCLC诊断、预后相关的基因,为NSCLC分子机制的进一步研究提供生物信息学依据,利用生物信息学方法对GEO数据库和TCGA数据库的数据集进行合并分析,筛选NSCLC组织与正常肺组织之间的差异表达基因(differentially expressed genes, DEGs),并对所取交集的DEGs进行基因集富集分析(gene set enrichment analysis, GSEA)、基因本体论(gene ontology, GO)分析、KEGG (kyoto encyclopedia of genes and genomes)通路富集分析、蛋白质相互作用(protein-protein interaction, PPI)分析、ROC曲线诊断效能分析及LASSO生存分析。文中共筛选出240个DEGs,主要涉及核分裂、染色体分离等生物学过程。GSEA分析结果显示,富集的通路主要涉及DNA修复和细胞周期。从PPI网络中筛选出20个hub基因, ROC结果显示, UBE2C (AUC=0.939)、TOP2A(AUC=0.927)、RRM2 (AUC=0.927)、CCNB1 (AUC=0.928)、MKI67 (AUC=0.930)、AURKA (AUC=0.931)、MELK(AUC=0.950)相对具有较高的诊断价值, LASSO COX回归结果则显示IL6、KIAA0101、MKI67、TPX2、AURKA、CDKN3及CDCA5与NSCLC患者的预后强相关。本研究结果表明, ZWINT、KIF2C、MELK、CDCA5可能在NSCLC中发挥着重要的作用,为阐明NSCLC的分子机制提供了新思路。     

Identification of KIF4A and its effect on the progression of lung adenocarcinoma based on the bioinformatics analysis

Background: Lung adenocarcinoma (LUAD) is the most frequent histological type of lung cancer, and its incidence has displayed an upward trend in recent years. Nevertheless, little is known regarding effective biomarkers for LUAD.Methods: The robust rank aggregation method was used to mine differentially expressed genes (DEGs) from the gene expression omnibus (GEO) datasets. The Search Tool for the Retrieval of Interacting Genes (STRING) database was used to extract hub genes from the protein–protein interaction (PPI) network. The expression of the hub genes was validated using expression profiles from TCGA and Oncomine databases and was verified by real-time quantitative PCR (qRT-PCR). The module and survival analyses of the hub genes were determined using Cytoscape and Kaplan–Meier curves. The function of KIF4A as a hub gene was investigated in LUAD cell lines.Results: The PPI analysis identified seven DEGs including BIRC5, DLGAP5, CENPF, KIF4A, TOP2A, AURKA, and CCNA2, which were significantly upregulated in Oncomine and TCGA LUAD datasets, and were verified by qRT-PCR in our clinical samples. We determined the overall and disease-free survival analysis of the seven hub genes using GEPIA. We further found that CENPF, DLGAP5, and KIF4A expressions were positively correlated with clinical stage. In LUAD cell lines, proliferation and migration were inhibited and apoptosis was promoted by knocking down KIF4A expression.Conclusion: We have identified new DEGs and functional pathways involved in LUAD. KIF4A, as a hub gene, promoted the progression of LUAD and might represent a potential therapeutic target for molecular cancer therapy.     

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 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.     

Restoration of miR-29b exerts anti-cancer effects on glioblastoma

Background

Methylation plays an important role in the etiology and pathogenesis of colorectal cancer (CRC). This study aimed to identify aberrantly methylated-differentially expressed genes (DEGs) and pathways in CRC by comprehensive bioinformatics analysis.

Methods

Data of gene expression microarrays (GSE68468, GSE44076) and gene methylation microarrays (GSE29490, GSE17648) were downloaded from GEO database. Aberrantly methylated-DEGs were obtained by GEO2R. Functional and enrichment analyses of selected genes were performed using DAVID database. Protein–protein interaction (PPI) network was constructed by STRING and visualized in Cytoscape. MCODE was used for module analysis of the PPI network.

Results

Totally 411 hypomethylation-high expression genes were identified, which were enriched in biological processes of response to wounding or inflammation, cell proliferation and adhesion. Pathway enrichment showed cytokine–cytokine receptor interaction, p53 signaling and cell cycle. The top 5 hub genes of PPI network were CAD, CCND1, ATM, RB1 and MET. Additionally, 239 hypermethylation-low expression genes were identified, which demonstrated enrichment in biological processes including cell–cell signaling, nerve impulse transmission, etc. Pathway analysis indicated enrichment in calcium signaling, maturity onset diabetes of the young, cell adhesion molecules, etc. The top 5 hub genes of PPI network were EGFR, ACTA1, SST, ESR1 and DNM2. After validation in TCGA database, most hub genes still remained significant.

Conclusion

In summary, our study indicated possible aberrantly methylated-differentially expressed genes and pathways in CRC by bioinformatics analysis, which may provide novel insights for unraveling pathogenesis of CRC. Hub genes including CAD, CCND1, ATM, RB1, MET, EGFR, ACTA1, SST, ESR1 and DNM2 might serve as aberrantly methylation-based biomarkers for precise diagnosis and treatment of CRC in the future.

     

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.     

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.     

Gene microarray analysis of expression profiles in Suberoyllanilide hyroxamic acid-treated Dendritic cells

Objective

The purpose of this study is to provide a further theoretical basis for the role of Suberoyllanilide hyroxamic acid (SAHA) affect on Dendritic cells (DCs).