CTNNBIP1和ICMT基因共表达在膀胱尿路上皮癌中作用和功能分析

为了研究CTNNBIP1和ICMT基因共表达在膀胱尿路上皮癌中的作用及对患者预后的诊断价值,本研究选取癌症基因组图谱(TCGA)数据库中412例膀胱尿路上皮癌数据资料,利用cBioPortal数据库对膀胱尿路上皮癌CTNNBIP1基因、ICMT基因及其共表达基因作生存分析,将Pearson和Spearman相关系数均大于0.3的基因定义为中等程度以上共表达相关的基因。通过String数据库获取这些基因共表达的互作关系组。采用DAVID数据库、GO数据库、KEGG数据库分别对其进行信号通路和功能分析,生物学过程聚类分析,信号通路聚类分析。结果显示:TCGA数据库中的膀胱尿路上皮癌患者CTNNBIP1和ICMT存在共表达(Pearson相关系数=0.46和Spearman相关系数=0.47)。String数据库显示基因共表达有38组具有互作关系。KEGG数据库显示CTNNBIP1-ICMT基因共表达所富集的信号通路集中在细胞周期(p<0.05)。DAVID数据库分析其功能主要为调节细胞生长和有丝分裂、负性调控Wnt信号通路、蛋白激酶结合等。生存分析证实CTNNBIPI1和ICMT基因共表达与患者总生存率显著相关(p=0.002 72),CTNNBIPI1和ICMT共表达阳性患者的预后最差。由此得出结论:CTNNBIP1和ICMT在膀胱尿路上皮癌中存在共表达,对CTNNBIP1-ICMT基因共表达网络和生物信息学分析,可找到其在膀胱尿路上皮癌中的作用及信号通路,为深入研究膀胱尿路上皮癌的发病机制提供了理论基础,可提高膀胱尿路上皮癌患者的预后。     

三棱内酯B调控人冠状动脉内皮细胞基因表达谱的生物信息学分析

目的:分析三棱内酯B在人冠状动脉内皮细胞中的表达谱数据集,寻找三棱内酯B调控血管内皮功能的关键作用靶点。方法:基于GEO公共数据库,下载原始表达谱数据集(GSE44598),经过差异基因筛选,功能注释,通路富集,信号通路网络以及基因互作网络分析,找出三棱内酯B对人冠状动脉内皮细胞基因表达谱产生影响的关键基因和信号通路。结果:同对照组相比,三棱内酯B给药组共有5224个基因有显著性差异,包括2628个上调基因和2596个下调基因。基因功能注释和信号通路富集分析表明,差异基因主要参与了细胞周期过程。网络分析显示,MAPK信号通路、细胞周期通路以及PLCG2,PRKACA和ADCY4等为关键信号通路和基因。结论:三棱内酯B通过影响PLCG2,PRKACA,ADCY4等基因的表达,参与MAPK和细胞周期等信号通路,从而调节人冠状内皮细胞的功能。这些关键基因和信号通路是三棱内酯B在心血管疾病治疗应用中潜在的作用靶点。     

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

HBV感染者肝细胞样本的关键基因和信号通路的生物信息学挖掘

为探究乙型肝炎病毒(Hepatitis B virus,HBV)感染后肝细胞基因表达和信号通路的改变情况,从Gene Expression Omnibus(GEO)数据库中下载了HBV感染者肝细胞样本制作的基因表达谱数据集GSE83148,进行质量检测、数据标准化后筛选出差异表达基因,进一步做GO和KEGG富集分析以及基因网络相互作用分析,筛选关键基因和信号通路。从HBV感染样本中筛选出fold change≥2,p-value0. 05的上调差异表达基因44个,GO分析获得关键基因BAK1和TP63,差异表达基因网络互作获得5个位于枢纽位置的基因:NDUFS1、NDUFS2、COX7B、ATP5B、OPA1。KEGG分析获得关键信号通路有:乙型肝炎信号通路、病毒癌变信号通路、Fox O信号通路、PI3K-Akt信号通路。本研究筛选出的多数基因与线粒体和氧化呼吸链有关,造成这一现象的具体机制还需进一步探究。     

西伯利亚蝗卵发育过程的比较转录组分析及滞育关联基因筛选

【目的】通过筛选西伯利亚蝗Gomphocerus sibiricus卵滞育基因及代谢通路,初步明确西伯利亚蝗卵滞育分子机理。【方法】利用Illumina NovaSeq 6000测序平台对西伯利亚蝗早期发育(ES)、滞育(DS)和滞育后发育阶段(PS)的卵进行转录组测序;采用GO和KEGG富集分析并结合文献,预测滞育相关通路并聚类热图分析,筛选蝗卵滞育关联基因,并选取其中6个重要基因进行qRT-PCR验证。【结果】DSvs ES和PS vs DS两比较组分别富集到12 419和4 789个差异表达基因,且多上调表达;两比较组共表达差异基因为2 206个,主要与糖代谢、环境胁迫和生长发育相关。DS vs ES组富集最显著的GO条目为蛋白结合,PS vs DS组GO条目主要包含酶活性、细胞骨架构建及蛋白结合等过程。滞育关联基因主要集中在Wnt信号通路、胰岛素信号通路、细胞周期通路以及昆虫激素生物合成等通路。6个重要的滞育关联基因的表达量变化趋势与转录组数据结果基本一致。【结论】本研究初步明确了西伯利亚蝗卵滞育调控的重要代谢途径,并筛选出20个滞育关联基因,为后续深入研究该物种的适应机理奠定了基础。     

肾纤维化关键基因的生信筛选及分析

[目的]利用生物信息学技术筛选与肾纤维化相关基因，并预测与基因相关通路，疾病、细胞类型和有关药物。[方法]在公共基因芯片数据库(GEO)中获取与肾纤维化相关三个基因数据集，利用Venn图鉴定出共表达的差异基因，Metascape工具对差异基因进行功能(GO)和通路富集(KEGG)分析，还利用STRING构建蛋白相互作用网络(PPI)网络以及可视化工具Cytoscape筛选关键基因，最后用及Enrichr和Comparative Toxicogenomics Database(CTD)数据分析有关疾病、细胞类型和药物。[结果]与肾纤维化相关的数据集GSE148420、GSE38117、GSE54441经Venn图共得到83个DEGs。经Cytoscape对STRING工具得到的PPI网络可视化后，筛选出10个与肾纤维化相关的关键基因，分别是F13B、ALDH8A1、A1CF、PAH、KMO、ALDH6A1、SPP2、ACAT1、ABAT、CAT。GO和KEGG富集显示这些基因与氧化还原酶活性，血小板致密颗粒，色氨酸、结氨酸、亮氨酸和异亮氨酸等氨基酸代谢通路有关。利用Enrichr和CTD...     

基于生物信息学方法识别非小细胞肺癌(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潜在的治疗靶点。     

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

肠道病毒71型感染相关调控枢纽基因筛选

肠道病毒71型(enterovirus 71)是引起婴幼儿手足口病的重要病原体,目前尚无特效抗病毒药物,并且宿主对病毒感染的应答研究有限。为深入研究病毒感染后宿主调控机制,分析了EV71感染RD细胞的基因表达芯片数据,筛选到与感染时间相关的6 642个差异基因;同时使用加权基因共表达网络分析方法(WGCNA)构建基因调控网络,得到和病毒感染时间呈强正负关联的pink模块和darkgreen模块。结果表明:GO分析发现目标模块分别富集于前剪切体的反式组装与翻译起始,KEGG Pathway富集分析发现pink模块没有显著的通路富集,darkgreen模块通路富集于核糖体相关通路;Visant展示目标模块的调控网络并筛选出pink模块枢纽基因RPS4X、HSPA13、CXCL9、CD55与darkgreen模块枢纽基因ABLIM1、MPDU1、SUPT7L、CTSS。q PCR验证的3个持续上调的基因TXNIP、EGR1、c-FOS和8个枢纽基因的转录水平与芯片数据一致,其中TXNIP、EGR1、c-FOS表达上调超过2.5倍,CXCL9下调4.5倍。这些关键基因的发现可能为EV71感染机理的研究和药物研发提供参考。     

基于RNA-seq研究胶质母细胞瘤中差异表达基因及蛋白相互作用网络

本研究基于RNA-seq数据分析了胶质母细胞瘤中的差异表达基因,并对差异表达基因进行了功能(GO term)和通路(KEGG)富集分析。进一步通过蛋白相互作用网络挖掘了胶质母细胞瘤的调控机理。结果表明,405个基因在肿瘤组织中差异表达(p-value≤0.05,|log FC|≥1.5),其中216个(53.3%)基因上调,189个(46.6%)基因下调。基因本体(gene ontology,GO)富集结果表明,这些差异表达基因参与了离子转运,神经冲动传递,细胞信号转导和细胞粘附等。此外,KEGG通路富集结果表明,差异基因参与了许多重要的生物学通路,包括ECM受体相互作用、黏着和钙信号等通路。进一步的蛋白相互作用网络分析鉴定了5个关键的hub基因,包括PTK2B、CDK1、FN1、CCND1和FOS。这5个关键基因对于胶质母细胞瘤的发生和发展可能起到了关键作用,可以作为潜在的调控位点和筛选的标志物。     

家兔前体脂肪细胞分化不同时期基因表达谱分析

脂肪的过度积累严重危害人类健康。前体脂肪细胞分化是脂肪发育的关键过程,研究前体脂肪细胞分化相关基因的表达有助于认识脂肪沉积的机理。尽管家兔是一种理想的研究脂肪发育的动物模型,但是针对其前体脂肪细胞分化不同时期基因表达谱的研究鲜见报道。本研究通过诱导家兔前体脂肪细胞分化,在分化第0 d、3 d和9 d收集脂肪细胞,利用转录组测序(RNA-seq),在分化第3 d样本与第0 d样本的比较中筛选出1352个差异表达基因(differentially expressed genes, DEGs),在分化第9 d样本与第3 d样本的比较中筛选出888个DEGs。GO (gene ontology)功能富集和KEGG (kyoto encyclopedia of genes and genomes)通路分析发现,0~3 d分化期上调的DEGs显著富集在PPAR信号通路和PI3K-Akt信号通路上,3~9d分化期上调的DEGs显著富集到与细胞周期调控有关的GO条目和KEGG信号通路,0~3d和3~9d阶段特异上调的DEGs可能分别作用于细胞质和细胞核。通过DEGs的蛋白-蛋白互作(protein-protein interaction, PPI)网络分析发现,筛选出的核心节点(hub node)基因可能通过调控细胞周期而影响家兔前体脂肪细胞分化。     

单细胞测序分析人类胚胎干细胞神经分化的分子机制

目的探讨人类胚胎干细胞 (ESCs)分化为神经细胞的关键性靶基因及分子机制,为临床靶向治疗神经康复患者提供分子理论依据。方法基于GEO数据平台芯片,采用单细胞测序方法 (scRNA-seq),利用R语言从多分子维度(单细胞差异基因、蛋白互作网络和基因通路等)分析人类ESCs分化过程中的关键Marker基因并利用质控和数据过滤、PCA、TSNE分析、细胞轨迹分析、GO富集分析、KEGG富集分析、KEGG通路分析等论证Marker基因调控人类ESCs分化作用机制。结果 GO功能富集分析结果为:Marker基因在胚层分化、细胞外基质和信号转导中作用显著;Marker基因互作网络及KEGG通路显示了特征性纤维连接蛋白1(FN1)、同源域蛋白 (NANOG)、生长因子受体结合蛋白2 (GRB2)为关键基因;KEGG通路分析显示:(1)FN1在调控细胞外基质通路中作用显著;(2)NANOG、GRB2在调控干细胞多能性信号通路作用显著。结论 FN1可能通过调控细胞外基质通路介导人类ESCs基质环境的变化;NANOG、GRB2上调在干细胞多能性信号通路中高表达介导人类ESCs定向分化为神经组织。     

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.     

miR-29a靶基因预测及其相关生物信息学分析

目的:通过对miR-29a进行靶基因预测及相关生物信息学分析,为miR-29a靶基因的实验验证提供数据支持,以期为深入研究miR-29a的生物学功能和调控机制提供理论指导。方法:利用PubMed检索miR-29a相关文章,通过miRBase在线工具分析miR-29a序列。应用TargetScan及miRNAda两种计算方法预测miR-29a靶基因并取其交集作为分析的基因集合,分别进行基因本体(gene ontology,GO)中的分子功能和生物学过程以及KEGG(Kyoto Encyclopedia of Genes and Genomes)生物通路富集分析。结果:(1)miR-29a序列在多物种间具有高度保守性。(2)两种方法预测miR-29a靶基因交集共191个。(3)miR-29a靶基因GO分子功能集中于转录因子活性、DNA结合和钙离子结合等(P0.05);miR-29a靶基因GO生物学过程集中于调控转录、细胞粘附、细胞增殖与凋亡等(P0.05);KEGG生物通路主要富集于PI3K-AKT信号通路、JAK-STAT信号通路、T细胞受体信号通路和胰岛素信号通路等信号转导通路,以及肺小细胞癌和子宫内膜癌等疾病通路(P0.05)。结论:miR-29a可能通过参与多个靶基因信号通路的调控,在机体的多种生理病理过程中发挥重要作用,是一个颇有研究价值的生物学靶标。     

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

本研究旨在利用生物信息学方法构建经铜诱导的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为核心差异表达转录因子。该研究为肝豆状核变性转录调控相关分子的生物学功能研究提供了重要的参考依据。     

Identification of Key Genes and Pathways Involved in the Heterogeneity of Intrinsic Growth Ability Between Neurons After Spinal Cord Injury in Adult Zebrafish

In the adult central nervous system (CNS), axon regeneration is a major hurdle for functional recovery after trauma. The intrinsic growth potential of an injured axon varies widely between neurons. The underlying molecular mechanisms of such heterogeneity are largely unclear. In the present study, the adult zebrafish dataset GSE56842 were downloaded. Differentially expressed genes (DEGs) were sorted and deeply analyzed by bioinformatics methods. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of DEGs were performed with the DAVID. A DEGs-associated protein–protein interaction network was constructed from the STRING database and visualized with Cytoscape software. In total, 621 DEGs were identified. GO analysis showed that the biological processes of DEGs focused mainly on the Notch signaling pathway, cell differentiation and positive regulation of neuron differentiation. The molecular functions mainly included calcium-transporting ATPase activity and calcium ion binding and structural constituents of the cytoskeleton. The cellular components included the plasma membrane, spectrin, and cytoplasmic and membrane-bound vesicles. KEGG pathway analysis showed that these DEGs were mainly involved in the metabolic pathway and Notch signaling pathway, and subnetworks revealed that genes within modules were involved in the metabolic pathway, Wnt signaling pathway, and calcium signaling pathway. This study identified DEG candidate genes and pathways involved in the heterogeneity of the intrinsic growth ability between neurons after spinal cord injury in adult zebrafish, which could facilitate our understanding of the molecular mechanisms underlying axon regeneration, and these candidate genes and pathways could be therapeutic targets for the treatment of CNS injury.

     

金针菇菌柄发育的转录组与蛋白组分析

菌柄是金针菇等食用菌的主要商品部位,但其生长机制仍不明确。本研究对金针菇伸长期和成熟期菌柄进行了转录组联合蛋白组分析,结果显示,两样本显著性差异表达基因和蛋白分别为721个和61个,均以上调表达为主。GO（gene ontology）功能聚类分析表明：有72.41%的差异表达基因富集在催化活性（catalytic activity）条目下。细胞组分（cell part）和绑定结合（binding）条目同时富集了较多的差异表达基因和蛋白。KEGG通路富集分析显示：碳水化合物代谢通路（carbohydrate metabolism）和氨基酸代谢通路（amino acid metabolism）富集的差异表达基因较多。差异表达蛋白富集较多的通路是单环菌素生物合成（monobactam biosynthesis,ko00261）、链霉素生物合成（streptomycin biosynthesis,ko00521）和有机含硒化合物代谢（selenocompound metabolism,ko00450）等。内质网蛋白质加工（protein processing in endoplasmic reticulum,ko04141）和MAPK信号通路（MAPK signaling pathway-yeast,ko04011）在转录组和蛋白组的KEGG富集分析中均为差异通路。本研究联合转录组和蛋白组数据筛选了40个金针菇菌柄发育中差异表达基因,为深入研究揭示食用菌菌柄发育过程提供候选基因。     

家族性高胆固醇血症HDL-miRNAs调节血单核细胞功能机制

本研究通过生物信息学方法分析家族性高胆固醇血症患者外周血单核细胞差异表达基因、HDL载体差异表达miRNA及其生物学功能,研究差异HDL-miRNA与单核细胞差异基因的相关性,探讨HDL-miRNA调控外周血单核细胞功能机制,寻找动脉粥样硬化防治新靶点。运用R语言分析GEO数据库共享平台家族性高胆固醇血症外周血单核细胞基因及HDL-miRNA探针芯片得到差异基因及差异miRNA,利用miRwalk2. 0预测miRNA靶基因,并利用STRING进行蛋白互作分析,构建差异miRNA与差异基因之间的调控网络。运用GO及KEGG分析研究基因功能。利用GEO数据(GSE6054)筛选出834个差异表达基因,利用GEO数据(GSE25108)筛选出HDL上差异miRNA28个。交叉匹配得到由19个差异miRNA和56个差异基因组配对的74对miRNA-靶基因。GO富集分析56个差异基因主要富集于肾上腺素受体信号等分子功能。KEGG分析56个差异基因主要富集于造血谱系通路上。家族性高胆固醇血症差异HDL-miRNA与外周血单核细胞差异mRNA具有相关性,HDL-miRNA有通过调控血单核细胞功能的可能性,可能参与高胆固醇血症导致动脉粥样硬化过程。     

低温胁迫下意大利蜜蜂预蛹的差异表达基因趋势分析

【目的】蜜蜂是典型的具有发育狭温性的全变态昆虫。本研究以对低温最敏感的意大利蜜蜂Apis mellifera ligustica预蛹为研究对象,通过低温胁迫不同时间的差异表达基因(differentially expressed genes, DEGs)趋势分析,探讨低温胁迫对蜜蜂发育影响的关键基因。【方法】对3日龄意大利蜜蜂封盖子预蛹进行20℃低温胁迫18 h(T18)和36 h(T36),以未经低温胁迫的预蛹为对照(CK),通过Illumina HiSeq~(TM)平台进行转录组学测定。利用Short Time-series Expression Miner(STEM)软件对2个处理组与对照组比较共有的差异表达基因进行趋势分析,再进一步对显著富集趋势模式中富集的差异表达基因进行GO分类和KEGG pathway分析。利用RT-qPCR对随机挑选的5个DEGs的表达模式进行验证。【结果】对检测到1 062个T18 vs CK和T36 vs CK共有的DEGs进行趋势分析,发现3个显著的基因表达模式,包括2个上调表达模式(Profile 6,有539个基因;Profile 7,有271个基因),1个下调表达模式(Profile 1,有183个基因)。对3个显著富集趋势模式DEGs分别进行GO富集分析和KEGG pathway分析,在Profile 6中找到同时富集在FoxO信号通路和寿命调节信号通路上的胰岛素样肽基因ILP和叉头蛋白O基因FoxO持续上调,说明低温胁迫会影响蜜蜂预蛹蜕皮激素信号传递;在Profile 7中CREB结合蛋白基因CBP持续上调,蜜蜂预蛹受到低温胁迫会影响细胞发育;在Profile 1中细胞色素P450基因CYP450 306a1 (phm)和WNT1显著下调,说明化蛹时受到低温胁迫会影响蜕皮激素的合成。通过RT-qPCR分析挑选的5个基因的表达结果和高通量测序结果一致。【结论】通过对蜜蜂响应低温胁迫的差异表达基因趋势分析发现,蜜蜂胰岛素和蜕皮激素共同调控的FoxO可能是低温胁迫抑制蜜蜂化蛹的一个关键基因。本研究为探索低温胁迫影响蜜蜂发育变态的分子调控机制奠定了基础。