基于GEO数据库对儿童哮喘急性发作靶基因的生物信息学分析

为了通过生物信息学技术筛选并分析关键基因以探究儿童哮喘急性发作的病理机制,从GEO数据库下载儿童哮喘急性发作相关基因芯片数据集(GSE103166),用R软件包工具筛选出哮喘急性发作患儿与健康儿童的差异表达基因(differentially expressed genes, DEGs),用DAVID (version 6.8)在线工具获取DEGs列表的GO和KEGG注释结果,同时通过STRING (version 11.0)获取其蛋白质互作数据,应用Cytoscape及其插件MCODE构建蛋白质互作网络、鉴定关键基因。结果显示,共筛选出78个DEGs,其中上调基因共49个,下调基因共29个。这些DEGs主要富集于核质转运调控、免疫系统过程调节等生物过程;核小体、DNA包装复合物等细胞组分; Rho GTP酶结合、离子型谷氨酸受体结合等分子功能。KEGG通路分析结果提示其主要富集在系统性红斑狼疮和哮喘信号通路。在这些DEGs中,筛选得到一个基因模块, GO分析显示该模块主要富集在γ干扰素介导的信号通路、通过MHCⅡ类分子进行抗原加工和外源性抗原肽的呈递、通过MHCⅡ类分子进行抗原加工和抗原肽的呈递等生物过程; MHCⅡ类蛋白复合物、内质网膜腔侧、内质网膜腔侧组成部分等细胞组分;MHCⅡ类受体活性、抗原肽结合、跨膜信号受体活性等分子功能。KEGG通路分析提示该模块主要富集在哮喘、移植物抗宿主病、同种异体排斥反应等通路。该基因模块涵盖5个关键基因,分别为HLA-DPB1、HLADQB1、HLA-DQB2、MT2A和KIF11。上述分析结果表明,文中筛选的DEGs和关键基因有助于加深对儿童哮喘急性发作病理机制的理解,其中HLA-DPB1、HLA-DQB1、HLA-DQB2等关键基因已得到相关研究的验证,未经证实的MT2A、KIF11关键基因,可能是儿童哮喘急性发作研究的新靶点。     

慢性鼻-鼻窦炎伴鼻息肉基因表达谱的生物信息学分析

本研究从美国国立生物信息中心(NCBI)的子数据库基因表达数据库(GEO)中选择基因表达谱GSE36830数据集,采用GEO2R筛选正常钩突和鼻息肉组织之间的差异表达基因(DEGs),对关键通路和差异表达基因进行数据库挖掘和分析,经筛选后的差异表达基因采用戴维在线工具对其进行基因本体富集分析(GO)、京都基因和基因组百科全书(KEGG)分析,然后将DEGs导入String数据库进行蛋白质互作网络分析,绘制差异表达基因互作网络图,将其数据导入Cytoscape软件中,筛选网络中心节点和关键基因,分析关键子网络。共筛选出699个DEGs,其中475个基因为上调表达基因,224个基因为下调表达基因。在GO分析中,针对生物过程,上调的DEGs包括:炎症反应、免疫反应、细胞趋化性、炎症反应的正向调节和细胞的粘附等;下调的DEGs主要参与:唾液分泌、生物矿物组织发展、细胞氨基酸生物合成过程、视网膜内稳态及离子跨膜转运等。在KEGG分析中,上调的DEGs主要在参与造血细胞系、细胞因子-细胞因子受体相互作用、破骨细胞分化、趋化因子信号通路、癌症中的转录失调、哮喘、金黄色葡萄球菌感染等信号通路中富集,而下调的DEGs在唾液腺分泌及胆汁分泌信号通路中富集。差异表达基因互作网络图筛选出前10个关键基因:ITGAM、IL10、CD86、TLR8、ITGAX、CCL2、CCR7、SRC、EGF及ITGB2。本研究得到了一组鼻息肉差异表达基因的生物信息学分析结果,但仍需进一步用基础试验来验证。本文分析的结论为慢性鼻-鼻窦炎、鼻息肉的研究提供了新的研究方向,也为鼻息肉发病机制研究的思路提供了一定的建设性作用。     

胰腺癌中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可能与胰腺癌的发生发展最为相关,为进一步探究胰腺癌的发病机制提供了生物信息学依据。     

COPD 差异表达基因的生物信息学分析及在LUSC样本中的表达

为确定慢性阻塞性肺病(COPD)的分子标记物及COPD与肺鳞状细胞癌(LUSC)共存的差异表达基因,探寻COPD合并肺癌的预测因子,发现新的治疗靶点。本研究采用生物信息学方法,从GEO数据库中筛选3套基因芯片数据集,挖掘COPD患者小气道上皮细胞(SAEC)的差异表达基因(DEG)以及潜在的生物标记物,并通过基因本体(GO)、京都基因与基因组百科全书(KEGG)富集分析预测DEGs的功能及参与的代谢途径。继而对DEGs构建PPI网络,使用Cytoscape软件筛选子模块和Hub基因,并将Hub基因通过TCGA数据库分析其在LUSC中的差异表达情况及差异基因间的相关性。结果共获得52个上调基因和24个下调基因,代谢通路主要集中在细胞色素P450对外源物质的代谢、化学致癌、花生四烯酸代谢及甲状腺激素合成四条途径上,通过Cytoscape软件从PPI网络中筛选得到2个功能模块和10个Hub基因,进一步验证发现其中5个基因在TCGA数据库中的LUSC样本中同样差异表达。由此推测SPP1、ALDH3A1、SPRR3、KRT6A和SPRR1B 可能为COPD 分子标记物及COPD与LUSC共存的DEGs,从而为研究COPD和LUSC的发病机制及二者潜在关系奠定良好的基础。     

肾透明细胞癌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,其中数个核心基因广泛参与多种肿瘤的病理进程,但尚未在肾透明细胞癌有相关研究报道,提示其可能是治疗肾透明细胞癌的潜在靶点。     

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

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

目的探讨强直性脊柱炎(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发病机制中发挥了重要的作用。     

家族性双侧大结节性肾上腺增生症生物标志物的生信分析

为寻找与家族性双侧大结节性肾上腺皮质增生症发展有关的潜在治疗靶点和生物标志物。从GEO数据库中下载GSE171558数据集,筛选受家族影响的肾上腺结节与正常的肾上腺组织之间的差异表达基因(Differentially expressed genes, DEGs),并进行基因功能富集分析和蛋白质-蛋白质相互作用网络分析。通过Cytoscape v3.9.1软件中的插件cytoHubba筛选出关键基因,进一步经NetworkAnalyst分析TF-miRNA共调控网络和蛋白质-化合物相互作用。共鉴定出336个DEGs,这些基因主要富集在细胞粘附过程、细胞增殖的正调节过程和RNA加工过程等生物过程,并涉及钙信号通路、PI3K-Akt信号通路和cAMP信号通路等。通过cytoHubba插件获得5个hub基因,经验证分析,多功能蛋白聚糖(Versican,VCAN)、双糖链蛋白聚糖(Biglycan,BGN)被认为是家族性双侧大结节性肾上腺皮质增生症的潜在生物标志物。进一步的GSEA分析结果显示,VCAN主要与丁酸代谢、ECM-受体相互作用和类固醇生物合成等有关。BGN主要涉及剪接体、皮质醇的合...     

髓母细胞瘤基因表达谱芯片关键基因的生物信息学分析

筛选髓母细胞瘤(medulloblastoma, MB)发生发展的关键基因,可为MB分子机制的进一步研究提供生物信息学依据。本文通过下载GEO (Gene Expression Omnibus)数据库GSE50161原始数据,利用R语言对正常脑组织与髓母细胞瘤组织中差异表达的基因进行分析;通过生物信息学分析工具(DAVID、STRING和Cytoscape)对差异基因进行生物学功能和蛋白质相互作用(protein-protein interaction, PPI)分析,并通过PPI筛选互作调控的关键基因。结果显示,总共筛选出999个差异表达的基因,鉴定了CCNB1、AURKB、MAD2L1、CENPE、KIF2C、BUB1、BUB1B、NDC80、CENPF、CDC20十个关键基因。差异基因生物学功能主要富集于有丝分裂的核分裂、染色体分离、微管蛋白结合、RAGE受体结合等生物过程。KEGG信号通路分析结果显示差异基因主要富集于细胞周期、NF-κB、IL-17和T细胞受体等信号通路。10个关键基因的生物学功能和信号通路主要富集于细胞有丝分裂和细胞周期通路。因此,细胞周期通路对MB的增殖和分裂起着关键性的作用,相关的分子机制值得进一步深入研究。     

基于生物信息学方法识别非小细胞肺癌(NSCLC)潜在治疗靶基因

本研究运用生物信息学方法识别非吸烟女性非小细胞肺癌(NSCLC)潜在的靶基因,并从分子水平探索其潜在的发病机制。从GEO数据库下载非吸烟女性非小细胞肺癌相关基因芯片数据集,经癌症组和癌旁对照组差异表达基因识别,并利用R软件对差异基因进行层次聚类分析,DAVID进行基因本体(gene ontology)和KEGG通路富集分析,STRING和Cytoscape软件构建蛋白-蛋白交互(PPI)网络,以及运用PASTAA分析,识别NSCLC相关转录因子,构建转录因子-基因共表达网络。结果表明,185个基因在NSCLC中差异表达,其中40个上调,145个下调;通过PASTAA分析识别出5个NSCLC基因相关转录因子。差异基因与胶原分解代谢过程、炎症反应的正调控等生物过程密切相关,基因的产物主要参与蛋白质细胞外基质、胶原三聚体等细胞组分,且主要发挥调节金属内肽酶活性、肝素结合和调节受体活性等分子功能;KEGG通路富集分析表明差异基因显著富集到胞外基质-受体信号通路、粘着斑信号通路、PPAR信号通和PI3K-Akt信号通路等,与非小细胞肺癌的发生发展密切相关。通过生物信息学方法,最终筛选到4个NSCLC关键基因:IL6、MMP1、COL1A1、CD36,其可能是非吸烟女性NSCLC潜在的治疗靶点。     

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

本研究旨在利用生物信息学方法构建经铜诱导的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个。上、下调关键功能模块分别包括了3 785个和3 931个基因。关键功能模块中的基因主要定位于细胞-基质连接、染色体、剪接复合体、核糖体等区域,共同参与了mRNA加工、组蛋白修饰、RNA剪切、DNA代谢调节、蛋白磷酸化等生物学过程,且与转录共调控活性、DNA转录因子结合、泛素样蛋白连接酶结合等分子功能相关。KEGG分析表明功能模块中的基因显著富集的通路包括乙型肝炎、有丝分裂原活化蛋白激酶（mitogen-activated protein kinase,MAPK）信号通路、细胞衰老和凋亡、神经营养信号通路和神经变性途径等。肝豆状核变性转录调控网络包括11个差异表达转录因子和96个差异表达基因,其中U2AF1、NFRKB、FUS、MAX、SRSF1、CEBPA和RXRA为核心差异表达转录因子。该研究为肝豆状核变性转录调控相关分子的生物学功能研究提供了重要的参考依据。     

Screening of hub genes and prediction of putative drugs in arsenic-related bladder carcinoma: An in silico study

BackgroundEvidence showed that inorganic arsenic (iAs) can trigger malignant transformation in cells with complex mechanisms. Thus, we aimed to investigate the possible molecules, pathways and therapeutic drugs for iAs-induced bladder cancer (BC) by using bioinformatics approaches.MethodsMicroarray-based data were analyzed to screen the differentially expressed genes (DEGs) between iAs-related BC cells and controls. Then, the roles of DEGs were annotated and the hub genes were screened out by protein-protein interaction network. The key genes were further selected from the hub genes through an assessment of the prognostic values. Afterward, potential drugs were predicted by using CMAP analysis.ResultsAnalysis of a dataset (GSE90023) generated 21 upregulated and 47 downregulated DEGs, which were enriched in various signaling pathways. Among the DEGs, four hub genes including WNT7B, SFRP1, DNAJB2, and ATF3, were identified as the key genes because they might predict poor prognosis in BC patients. Lastly, Cantharidin was predicted to be a potential drug reversing iAs-induced malignant transformation in urinary epithelium cells.ConclusionThe present study found several hub genes involved in iAs-induced malignant transformation in urinary epithelium cells, and predicted several small agents for iAs toxicity prevention or therapy.     

应用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个枢纽基因以及一些关键信号通路和生物学过程,可为后续深入研究细菌性败血症致病机制奠定理论依据。     

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.

     

Uncovering the pathogenesis and identifying novel targets of pancreatic cancer using bioinformatics approach

Pancreatic cancer is a uniformly lethal disease that can be difficult to diagnose at its early stage. Thus, our present study aimed to explore the underlying mechanism and identify new targets for this disease. The data GSE16515, including 36 tumor and 16 normal samples were available from Gene Expression Omnibus. Differentially expressed genes (DEGs) were screened out using Robust Multichip Averaging and LIMMA package. Moreover, gene ontology and pathway enrichment analyses were performed to DEGs. Followed with protein–protein interaction (PPI) network construction by STRING and Cytoscape, module analysis was conducted using ClusterONE. Finally, based on PubMed, text mining about these DEGs was carried out. Total 274 up-regulated and 93 down-regulated genes were identified as the common DEGs and these genes were discovered significantly enriched in cell adhesion and extracellular region terms, as well as ECM-receptor interaction pathway. In addition, five modules were screened out from the up-regulated PPI network with none in down-regulated network. Finally, the up-regulated genes, including MIA, MET and CEACAMS, and down-regulated genes, such as FGF, INS and LAPP, had the most references in text mining analysis. Our findings demonstrate that the up- and down-regulated genes play important roles in pancreatic cancer development and might be new targets for the therapy.     

Identification of prognostic biomarkers for breast cancer brain metastases based on the bioinformatics analysis

PurposeThe prognosis of breast cancer (BC) patients who develop into brain metastases (BMs) is very poor. Thus, it is of great significance to explore the etiology of BMs in BC and identify the key genes involved in this process to improve the survival of BC patients with BMs.Patients and methodsThe gene expression data and the clinical information of BC patients were downloaded from TCGA and GEO database. Differentially expressed genes (DEGs) in TCGA-BRCA and GSE12276 were overlapped to find differentially expressed metastatic genes (DEMGs). The protein-protein interaction (PPI) network of DEMGs was constructed via STRING database. ClusterProfiler R package was applied to perform the gene ontology (GO) enrichment analysis of DEMGs. The univariate Cox regression analysis and the Kaplan-Meier (K-M) curves were plotted to screen DEMGs associated with the overall survival and the metastatic recurrence survival, which were identified as the key genes associated with the BMs in BC. The immune infiltration and the expressions of immune checkpoints for BC patients with brain relapses and BC patients with other relapses were analyzed respectively. The correlations among the expressions of key genes and the differently infiltrated immune cells or the differentially expressed immune checkpoints were calculated. The gene set enrichment analysis (GSEA) of each key gene was conducted to investigate the potential mechanisms of key genes involved in BC patients with BMs. Moreover, CTD database was used to predict the drug-gene interaction network of key genes.ResultsA total of 154 DEGs were identified in BC patients at M0 and M1 in TCGA database. A total of 667 DEGs were identified in BC patients with brain relapses and with other relapses. By overlapping these DEGs, 17 DEMGs were identified, which were enriched in the cell proliferation related biological processes and the immune related molecular functions. The univariate Cox regression analysis and the Kaplan-Meier curves revealed that CXCL9 and GPR171 were closely associated with the overall survival and the metastatic recurrence survival and were identified as key genes associated with BMs in BC. The analyses of immune infiltration and immune checkpoint expressions showed that there was a significant difference of the immune microenvironment between brain relapses and other relapses in BC. GSEA indicated that CXCL9 and GPR171 may regulate BMs in BC via the immune-related pathways.ConclusionOur study identified the key genes associated with BMs in BC patients and explore the underlying mechanisms involved in the etiology of BMs in BC. These findings may provide a promising approach for the treatments of BC patients with BMs.     

Bioinformatics Analysis of Microarray Data to Reveal Novel Genes Related to Cold-Resistance of Maize

Low temperature has become a major abiotic stress factor that can reduce maize yield and cause a number of economic loss. This study was designed to identify key genes and pathways associated with coldresistance of maize. The gene expression profile GSE46704, including 4 control temperature treated plants and 4 low temperature treated plants, was downloaded from the Gene Expression Omnibus database. Differentially-expressed genes (DEGs) were identified by limma package. Then, protein-protein interaction (PPI) network and module selection were constructed using Cytoscape. Moreover, the DEGs were re-matched based on the Zea mays L. gene ID and symbol data from PlantRegMap. Finally, the re-matched DEGs were performed functional and pathway enrichment analyses by the DAVID online tool. A total of 750 DEGs were screened (including 387 up-regulated and 363 down-regulated genes) In the PPI network, GRMZM2G070837_P01 and GRMZM2G114578_P01 had higher degrees. Besides, carbohydrate metabolic process, starch and sucrose metabolism and biosynthesis of secondary metabolites were significantly enriched in functional and pathway enrichment analysis. GRMZM2G070837_P01 and GRMZM2G114578_P01 might play a critical role in cold-resistance of maize. Meanwhile, carbohydrate metabolic process, starch and sucrose metabolism and biosynthesis of secondary metabolites might function in cold-resistance of maize.