基于生物信息学技术筛选影响胶质母细胞瘤化疗敏感性相关基因的研究

目的：应用生物信息学技术筛选影响胶质母细胞瘤（GBM）化疗敏感性的相关基因。方法：对2批胶质瘤患者BIOSTAR基因芯片进行分析。通过随访完善临床资料,筛选芯片中胶质母细胞瘤患者生存期长、短两组间的差异基因,明确差异基因参与的功能和通路,并构建与烷化剂相关基因的信号传导网络,结合芯片数据、患者预后和信号传导网络,筛选GBM化疗敏感性的相关基因。结果：两组芯片中间差异基因有503条。2批芯片的差异基因主要参与62项基因功能,主要参与31条信号传导通路。通过对差异基因功能、通路,烷化剂信号转导网络的分析,得到影响胶质母细胞瘤化疗敏感性的核心的差异基因IFNGR2、IL8、ITGA5、TNFRSF1B。结论：通过严谨的实验设计和科学的统计学判别,结合患者完整的生存资料,本研究成功地应用生物信息学技术对基因芯片的大量数据进行挖掘和分析,并筛选出了可能影响GBM患者预后和化疗药物敏感性的基因,为进一步功能实验和患者个体化治疗奠定了基础。     

基于生物信息学技术筛选影响胶质母细胞瘤化疗敏感性 相关基因的研究

目的:应用生物信息学技术筛选影响胶质母细胞瘤(GBM)化疗敏感性的相关基因。方法:对2批胶质瘤患者BIOSTAR基因芯片进行分析。通过随访完善临床资料,筛选芯片中胶质母细胞瘤患者生存期长、短两组间的差异基因,明确差异基因参与的功能和通路,并构建与烷化剂相关基因的信号传导网络,结合芯片数据、患者预后和信号传导网络,筛选GBM化疗敏感性的相关基因。结果:两组芯片中间差异基因有503条。2批芯片的差异基因主要参与62项基因功能,主要参与31条信号传导通路。通过对差异基因功能、通路,烷化剂信号转导网络的分析,得到影响胶质母细胞瘤化疗敏感性的核心的差异基因IFNGR2、IL8、ITGA5、TNFRSF1B。结论:通过严谨的实验设计和科学的统计学判别,结合患者完整的生存资料,本研究成功地应用生物信息学技术对基因芯片的大量数据进行挖掘和分析,并筛选出了可能影响GBM患者预后和化疗药物敏感性的基因,为进一步功能实验和患者个体化治疗奠定了基础。     

胶质母细胞瘤预后相关基因的数据挖掘分析

胶质母细胞瘤（glioblastoma, GBM）是恶性程度最高的颅内恶性肿瘤，目前临床上缺乏有效治疗药物，复发率高且预后差，开发新的抗GBM药物是目前临床上亟待解决的问题。为了筛选与GBM预后密切相关的基因，为寻找新的药物靶点提供线索，采用GEO2R工具从GEO数据库中的269个肿瘤组织和61个正常组织中初步筛选出差异表达基因，然后利用Cluster Profiler数据库进行基因功能富集分析，STRING及Cytoscape进一步筛选出37个差异表达基因，采用GEPIA交互分析对这37个基因在GBM肿瘤组织中的表达进行验证。为了进一步探索这些差异表达基因与患者预后的关系，研究中利用GEPIA工具对TCGA数据库中与患者预后相关的数据进行深入挖掘，最终发现PTTG1、RRM2、E2F7与患者中位生存期呈显著性负相关。研究筛选出的与患者预后密切相关的基因不仅可以为评估患者预后提供参考，同时也为开发新的抗GBM药物提供了潜在的靶点。     

家蝇Prohibitin基因的生物信息学分析

为对家蝇Prohibitin蛋白序列进行生物信息学分析,从而为该基因功能研究奠定基础。利用在线分析程序和相关工具软件分析Prohibitin蛋白的理化性质、结构域、并预测其空间结构和功能。结果表明家蝇Prohibitin蛋白由277个氨基酸组成,分子量为30.54 k Da,理论等电点为5.26,为稳定蛋白,有跨膜区,但不含信号肽,该蛋白属于PHB保守结构域家族,亚细胞定位于细胞质,二级结构以α-螺旋为主。蛋白同源性比对结果显示,昆虫中的Prohibitin蛋白具有较高的同源性。这些分析结果可为今后深入研究该蛋白的结构特征和功能提供参考。     

自噬基因CTSL表达与胶质母细胞瘤患者预后相关

本研究旨在探讨自噬基因CTSL对胶质母细胞瘤(GBM)患者的预后影响。利用癌症基因组图谱(TCGA)、人类自噬数据库(HADB)、中国脑胶质瘤基因组图谱(CGGA)数据库、基因表达谱分析(GEPIA)获取数据信息,通过筛选差异表达基因及单因素和多因素COX分析确定GBM的独立预后危险因素,同时通过基因本体论(GO)、基因组百科全书途径(KEGG)、临床病理相关性、基因集富集分析(GSEA)、自噬基因网络分析CTSL的相关作用机制。结果显示:(1)富集分析显示胶质母细胞瘤中差异自噬基因(ARG)与自噬体的形成、细胞凋亡、血管生成、细胞化疗等相关;(2)GBM中CTSL的mRNA水平明显高于正常组织样本;(3)多因素COX回归分析显示自噬基因CTSL的高表达为GBM预后的独立危险因素,STUPP治疗(术后替莫唑胺[Tmz]同步放化疗+Tmz辅助化疗)为独立保护因素;(4)自噬基因CTSL在非GCIMP(CpG岛甲基化)型、间质型、IDH野生型、1p/19q无缺失型胶质母细胞瘤及化疗后表达量更高。综上所述,本研究分析了自噬基因在GBM中的作用,并表明自噬基因CTSL的过表达预示胶质母细胞瘤患者不良预后,显示自噬基因CTSL有作为有效靶标的潜质。     

甲基转移样酶1在胶质母细胞瘤中高表达并与促胶质母细胞瘤生长有关

目的 探讨甲基转移样酶1(methyltransferase-like 1,METTL1)在人胶质母细胞瘤中的表达及其临床意义.方法 通过免疫荧光双标方法观察正常人脑组织中METTL1蛋白的表达定位.采用TCGA数据库分析比较METTL1 mRNA在胶质母细胞瘤及正常人脑组织中的表达差异.选取60例人胶质母细胞瘤及10...     

基于单细胞测序构建胶质母细胞瘤预后模型

[目的]基于单细胞测序筛选胶质母细胞瘤特征基因并构建预后模型。[方法]分析GEO数据库单细胞RNA测序数据集GSE84465,筛选出GBM细胞分化相关的差异基因。下载TCGA数据库GBM的基因表达谱和临床数据,采用Lasso回归、Cox回归分析筛选出特征基因构建预后模型,根据独立预后因素构建列线图,GSE83300作为外部验证集。基于风险评分中位数将患者分组,比较两组生存差异。[结果]通过scRNA-seq得到492个分化差异基因,经过回归分析得到基于6个基因(PLAUR、RARRES2、G0S2、MDK、SERPINE2、CD81)的预后模型。其1、3、5年ROC曲线下面积均大于0.7;KM分析显示高低风险组预后存在差异(P<0.001),GSE83300验证结果与TCGA一致。多因素Cox回归分析表明年龄和风险评分可以作为独立影响因素(P<0.01);C-Index(0.679)、校准图显示列线图预测模型有良好的拟合度。GSEA分析示高低风险组差异基因集参与细胞因子受体相互作用、抗原处理与提呈等通路。[结论]由PLAUR、RARRES2、G0S2、MDK、SERPINE...     

基于生物信息学分析的肝癌核心基因的筛选及验证

本研究基于GEO数据库,选取由慢性乙型肝炎诱导的肝细胞癌芯片数据GSE121248为研究对象,利用GEO2R软件分析数据,筛选出差异表达基因,利用DAVID数据库进行GO分析和KEGG pathway富集分析.利用STRING数据库构建PPI网络,分析筛选核心基因.利用GEPIA对核心基因的表达进行验证,Kaplan ...     

多形胶质母细胞瘤细胞系端粒酶活性的测定

端粒缩短见于星形细胞瘤发展过程中,但其长度在胶质母细胞瘤／细胞系相对稳定,提示胶质瘤细胞内存在端粒修复机制的可能性．为证实此点,利用端粒重复片段扩增技术（ＴＲＡＰ）,对８株人／大鼠多形胶质母细胞系的蛋白提取液中端粒酶活性加以测定．结果显示：８例胶质瘤样本的反应液均可见端粒ＰＣＲ扩增片段;用无ＤＮａｓｅ的ＲＮａｓｅ事先处理蛋白提取液,可明显降低或消除ＰＣＲ产物的出现,说明ＴＲＡＰ反应中的ＰＣＲ扩增是在端粒酶的介导下进行而非ＤＮＡ污染或其它端粒修复因子所致．从而不但建立起检测人癌细胞内端粒酶活性的可靠方法,也为针对端粒酶的胶质母细胞瘤生物／药物治疗提供了实验依据．     

单纯疱疹病毒胸苷激酶基因对人胶质母细胞瘤细胞的杀伤效应研究

目的　研究单纯疱疹病毒胸苷激酶 (HSV1 TK)基因转染并联合抗病毒药更昔洛韦 (Ganciclovir ,GCV)对人胶质母细胞瘤细胞的杀伤效应。方法　采用基因工程技术构建带HSV1-TK基因的逆转录病毒重组体pLX SN TK ,采用脂质体介导入PA317包装细胞 ,建立重组逆转录病毒载体分泌细胞株PA317 TK ;用该细胞上清液转导人胶质母细胞瘤细胞 ,用不同浓度GCV作用人胶质母细胞瘤细胞SW0 38 C2 TK和野生型SW0 38 C2 ,采用噻唑蓝(MTT)比色法检测 72h后细胞存活率 ,并求出半杀伤浓度。结果　重组逆转录病毒载体pLXSN TK能有效地将TK基因导入SW0 38-C2细胞内 ,并使其获得对GCV的敏感性 ;体外细胞毒实验 :在 2 0 μmol LGCV存在的情况下 ,野生型细胞无明显改变 ,而实验组细胞明显死亡 ,半杀伤浓度IC50 为 0 8μmol L ,与SW0 38 C2相比 ,对GCV的敏感性提高了 6 0 0倍 ,旁杀效果也较明显。结论　SW0 38 C2转染TK基因后 ,能有效地被GCV杀灭。显示了其潜在的临床应用价值     

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

Identifying breast cancer subtypes associated modules and biomarkers by integrated bioinformatics analysis

Breast cancer is the most common form of cancer afflicting women worldwide. Patients with breast cancer of different molecular classifications need varied treatments. Since it is known that the development of breast cancer involves multiple genes and functions, identification of functional gene modules (clusters of the functionally related genes) is indispensable as opposed to isolated genes, in order to investigate their relationship derived from the gene co-expression analysis. In total, 6315 differentially expressed genes (DEGs) were recognized and subjected to the co-expression analysis. Seven modules were screened out. The blue and turquoise modules have been selected from the module trait association analysis since the genes in these two modules are significantly correlated with the breast cancer subtypes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment show that the blue module genes engaged in cell cycle, DNA replication, p53 signaling pathway, and pathway in cancer. According to the connectivity analysis and survival analysis, 8 out of 96 hub genes were filtered and have shown the highest expression in basal-like breast cancer. Furthermore, the hub genes were validated by the external datasets and quantitative real-time PCR (qRT-PCR). In summary, hub genes of Cyclin E1 (CCNE1), Centromere Protein N (CENPN), Checkpoint kinase 1 (CHEK1), Polo-like kinase 1 (PLK1), DNA replication and sister chromatid cohesion 1 (DSCC1), Family with sequence similarity 64, member A (FAM64A), Ubiquitin Conjugating Enzyme E2 C (UBE2C) and Ubiquitin Conjugating Enzyme E2 T (UBE2T) may serve as the prognostic markers for different subtypes of breast cancer.     

Identification of key genes and pathways of thyroid cancer by integrated bioinformatics analysis

Thyroid cancer is a common endocrine malignancy with a rapidly increasing incidence worldwide. Although its mortality is steady or declining because of earlier diagnoses, its survival rate varies because of different tumour types. Thus, the aim of this study was to identify key biomarkers and novel therapeutic targets in thyroid cancer. The expression profiles of GSE3467, GSE5364, GSE29265 and GSE53157 were downloaded from the Gene Expression Omnibus database, which included a total of 97 thyroid cancer and 48 normal samples. After screening significant differentially expressed genes (DEGs) in each data set, we used the robust rank aggregation method to identify 358 robust DEGs, including 135 upregulated and 224 downregulated genes, in four datasets. Gene Ontology and Kyoto Encyclopaedia of Genes and Genomes pathway enrichment analyses of DEGs were performed by DAVID and the KOBAS online database, respectively. The results showed that these DEGs were significantly enriched in various cancer-related functions and pathways. Then, the STRING database was used to construct the protein–protein interaction network, and modules analysis was performed. Finally, we filtered out five hub genes, including LPAR5, NMU, FN1, NPY1R, and CXCL12, from the whole network. Expression validation and survival analysis of these hub genes based on the The Cancer Genome Atlas database suggested the robustness of the above results. In conclusion, these results provided novel and reliable biomarkers for thyroid cancer, which will be useful for further clinical applications in thyroid cancer diagnosis, prognosis and targeted therapy.     

Screening and identification of hub genes in pancreatic cancer by integrated bioinformatics analysis

Pancreatic cancer (Pa) is a malignant tumor of the digestive tract with high degree of malignancy, this study aimed to obtain the hub genes in the tumorigenesis of Pa. Microarray datasets GSE15471, GSE16515, and GSE62452 were downloaded from Gene Expression Omnibus (GEO) database, GEO2R was conducted to screen the differentially expressed genes (DEGs), and functional enrichment analyses were carried out by Database for Annotation, Visualization and Integrated Discovery (DAVID). The protein-protein interaction (PPI) network was constructed with the Search Tool for the Retrieval of Interacting Genes (STRING), and the hub genes were identified by Cytoscape. Totally 205 DEGs were identified, consisting of 51 downregulated genes and 154 upregulated genes enriched in Gene Ontology terms including extracellular matrix (ECM) organization, collagen binding, cell adhesion, and pathways associated with ECM-receptor interaction, focal adhesion, and protein digestion. Two modules in the PPI were chosen and biological process analyses showed that the module genes were mainly enriched in ECM and cell adhesion. Twenty-four hub genes were confirmed, the survival analyses from the cBioPortal online platform revealed that topoisomerase (DNA) II α (TOP2A), periostin (POSTN), plasminogen activator, urokinase (PLAU), and versican (VCAN) may be involved in the carcinogenesis and progression of Pa, and the receiver-operating characteristic curves indicated their diagnostic value for Pa. Among them, TOP2A, POSTN, and PLAU have been previously reported as biomarkers for Pa, and far too little attention has been paid to VCAN. Analysis from R2 online platform showed that Pa patients with high VCAN expression were more sensitive to gemcitabine than those with low level, suggesting that VCAN may be an indicator to guide the use of the chemotherapeutic drug. In vitro experiments also showed that the sensitivity of the VCAN siRNA group to gemcitabine was lower than that of the control group. In conclusion, this study discerned hub genes and pathways related to the development of Pa, and VCAN was identified as a novel biomarker for the diagnose and therapy of Pa.     

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.