金针菇单、双核菌丝差异表达基因分析

对金针菇Flammulina velutipes单核菌丝W23的菌丝体以及与L11质配后的双核菌丝H1123菌丝体进行了转录组测序,以本实验室已获得的W23基因组为参考基因组研究两样本间差异基因,并对这些差异基因进行了GO功能和Pathway显著性富集分析。差异基因分析显示,两个样本中共有显著性差异表达的基因3 504个,其中在双核菌丝中上调、下调的基因数分别为2 151和1 353个。研究发现差异表达基因含有很多的转录因子基因、蛋白激酶以及WD40 repeat-like蛋白。Gene Ontology（GO）功能分析结果表明,extracellular region和membrane-enclosed lumen条目下的差异基因全部为上调表达,而envelope下的差异基因全部为下调表达,以利于双核菌丝分裂时锁状联合的形成而便于核的迁移。Pathway 功能富集分析结果表明,脂肪酸、氨基酸以及大部分糖类合成相关基因具有比较活跃的上调表达。说明双核菌丝主要进行营养物质的富集,为下一步在合适条件下分化成原基,进入生殖生长阶段储备物质基础。     

胃癌组织中肿瘤相关成纤维细胞与正常成纤维细胞之间转录组学研究及验证

胃癌组织中肿瘤相关成纤维细胞(carcinoma associated fibroblasts, CAFs)是胃癌微环境的重要成分,主要来源于正常成纤维细胞(normal fibroblasts, NFs)的活化,对胃癌的发生发展有重要作用,但是两者之间的基因表达差异并不完全清楚。本研究选取从人胃癌组织中分离获得的CAFs及NFs 各3组,进行转录组学研究,筛选出3组细胞中交集且差异倍数较大的基因12个,用Omicsbean在线工具对差异基因进行Gene Ontology (GO)功能及KEGG通路富集,构建蛋白质相互作用调控网络;最后用RT-qPCR验证CAFs和NFs中差异基因的表达。结果显示,筛选出的12个差异表达基因主要参与NF-κB信号、炎症、细胞黏附、细胞表面受体和细胞因子等功能,上述功能均与肿瘤的发生发展密切相关。RT-qPCR检测发现,与NFs相比,CAFs中BCL2A1、NKX3-2、CXCL12、TNFAIP3、FOS、CDH4及CLDN1表达上调;ATF3、CYFIP2、CCL11、KLF2及GDF15基因表达下调,差异均具有统计学意义(P＜0.05)。结果提示,胃癌CAFs与NFs中存在肿瘤相关的差异表达基因,这些差异基因可能在胃癌微环境中发挥重要作用。     

京海黄鸡柔嫩艾美耳球虫感染后盲肠转录组分析

为了探究京海黄鸡感染柔嫩艾美尔球虫(E. tenella)后盲肠组织差异表达基因,以及球虫感染分子应答过程和免疫应答机制,试验采用RNA-seq技术对E. tenella感染和非感染组第7天的盲肠组织进行转录组测序,筛选差异表达基因,并进行差异基因的功能、通路富集分析。结果表明,在感染和非感染组中有显著差异的表达基因2 830个(P0.05),其中1 419个基因上调,1 411个基因下调。随机选取10个差异基因进行qRT-PCR验证,结果显示差异基因的表达倍数与RNA-seq检测结果显著相关(r=0.988,P0.000),决定系数达0.975。GO分析表明,有2 356个差异基因获得GO功能注释,显著富集的前30个GO terms主要涉及细胞交流、信号转导、血管生成、氧化还原酶活性等。KEGG分析发现差异基因显著富集的信号通路有黏着斑、细胞外基质-受体相互作用、过氧化物酶体增殖物激活受体等。这些通路中的差异基因有ANGPTL4、ACSL5、VEGFC、CD44和MAKP10等,提示这些基因在宿主柔嫩艾美耳球虫感染过程中发挥重要作用。     

不同海拔饲育的麦洼牦牛肺脏组织转录组学分析

为了探讨牦牛适应高海拔低氧环境的基因表达特征与规律,对在高海拔（3 560 m）和低海拔地区（478 m）饲育4个月的2.5~3岁健康雄性麦洼牦牛肺组织进行转录组测序。转录组测序采用Illumina高通量测序平台(HiSeqTM2500/4000)进行,并以qRT-PCR验证差异表达基因的表达量。结果显示,高海拔组牦牛肺脏转录组平均每个测序样本得到约5.76亿条Clean Reads,低海拔组牦牛中得到约6.10亿条Clean Reads,比对到参考基因组上的Reads数分别占91.74%和91.28%以上,共发现了2 047个新转录本。低海拔组与高海拔组牦牛肺脏组织之间共有199个差异表达基因,其中含89个差异上调表达基因和110个差异下调表达基因。所得差异表达基因富集在297个GO条目和146个KEGG通路中,包含62个低氧适应相关的GO条目和35个低氧适应相关代谢通路。其中低氧适应相关GO条目在生物过程、细胞组成和分子功能三种类别中占比最多的分别为细胞粘附、蛋白复合物和钙离子结合。低氧适应相关KEGG通路中占比最多的为肿瘤坏死因子(TNF)信号通路,其次为低氧诱导因子1(HIF-1)信号通路。qRT-PCR验证结果显示,Ⅱ类人类白细胞抗原α链(HLA-DOA、HLA-DRA)、补体因子 (C2)和甘露糖结合凝集素相关丝氨酸蛋白酶1（MASP1）基因的表达量变化与转录组测序结果相符。本研究为全局和深入理解牦牛肺组织转录本表达对高海拔低氧的响应提供了有价值的切入点。     

基于转录组测序的羽衣甘蓝叶色相关基因分析

叶色是羽衣甘蓝重要的观赏性状之一。本研究采用Illumina Hi-Seq2500高通量测序技术,对基因型纯合的紫叶和白叶羽衣甘蓝叶片进行转录组测序,筛选差异基因并与GO和KEGG数据库比对进行注释分析,分析羽衣甘蓝叶色形成相关基因。结果显示,获得高质量短读序共104 608 770条,筛选出紫叶相对白叶的差异表达基因1 993个,其中上调表达基因1 094个,下调表达基因899个。根据GO功能分类可分为生物过程、细胞组分和分子功能3大类64功能组。根据KEGG代谢通路分析可以分为171类,在叶色相关的类黄酮生物合成途径中黄酮醇合成酶(FLS)、二氢黄酮醇4-还原酶(DFR)上调以及类胡萝卜素生物合成途径中的类胡萝卜素β-环化酶上调与紫叶形成关系密切。本研究丰富了羽衣甘蓝的转录组信息,获得了一些差异表达基因,为进一步研究羽衣甘蓝叶色形成的遗传机制提供了有价值的信息。     

MDCK细胞和MDCK-G1细胞感染H1N1流感病毒后病毒滴度和细胞转录组表达差异分析

目的通过分析细胞表达基因及其相关信号通路的差异等信息,探索MDCK和MDCK-G1细胞株在接种H1N1流感病毒后出现病毒滴度和细胞生长趋势差异的内在生物学特性等原因。方法通过Illumina测序平台对2株细胞进行转录组测序(RNA sequencing, RNA-seq)和测序结果生物信息学分析后,选择了17个差异基因进行实时定量聚合酶链反应(quantitative real-time polymerase chain reaction, qRT-PCR)验证。结果 MDCK和MDCK-G1细胞差异基因(differentially expressed genes, DEGs)|log_2(FoldChange)|0和padj≤0.05总数为2 786个。相比于MDCK细胞,MDCK-G1细胞有967个基因的表达显著上升,1 819个基因的表达显著下降。差异基因通过基因本体(Gene Ontology, GO)数据库富集功能注释和京都基因和基因组数据库(Kyoto Encyclopedia of Genes and Genomes, KEGG)进行通路分析发现,2株细胞在免疫调控和细胞生长周期调控上均存在差异,qRT-PCR验证的17个基因中有16个基因的表达趋势与转录组测序结果一致。结论 2株细胞的转录谱存在显著差异。     

脂多糖调控对大鼠心脏微血管内皮细胞的转录组分析

目的探索脂多糖(LPS)对大鼠心脏微血管内皮细胞(rCMECs)转录组的调节作用。方法对照组(正常培养rCMECs),LPS组(100 ng/mL LPS处理6 h的rCMECs),每组进行3个生物学重复转录组测序。得到差异基因后,使用实时定量PCR对部分差异基因mRNA的表达进行验证。分别对上调和下调基因进行GO和KEGG富集,并对差异基因进行共表达网络分析。采用独立t检验进行统计学分析。结果LPS处理后,265个基因表达上调,118个基因的表达下调。前10个最显著上调基因为:Mt2a、Cyp7b1、Sod2、Icam1、Ccl2、AC128848.1、Mt1、Cebpd、Serpinb2和Tnfrsf11b。前10个最显著下调基因为:Cavin2、Ankrd1、Edn1、Prss35、Lmod1、Dhrs3、Ttc22、Sema6a、Map2k3和Sema7a。定量PCR的结果表明Mt2a、Sod、Ccl2、Cxcl1、Icam1和Vcaml基因的表达得到了上调(P均<0.01);而Cavin2、Ankrd1、Edn1和Prss35基因表达下调(P均<0.05)。GO和KEGG富集的结果表明,上调基因与内皮细胞对炎性免疫细胞的趋化作用和黏附作用密切相关;而下调基因则是与钙离子信号和G蛋白相关通路以及内皮通透性增加有关。此外,差异基因进行共表达网络分析发现Sod2处于核心位置,提示其可能与LPS诱导的rCMECs的各种变化密切相关。结论LPS调控了rCMECs中大量与炎性免疫细胞进入心肌组织相关基因的表达。     

急性氨氮胁迫下大口黑鲈的肝脏转录组特征分析

为研究大口黑鲈(Micropterus salmoides)受到氨氮胁迫时基因表达的变化规律,文章采用Illumina平台的HiSeq测序策略分析了氨氮胁迫12h和24h后大口黑鲈肝脏的基因表达谱。氨氮胁迫后在大口黑鲈肝脏共获得111.66 Gb的有效数据,组装获得133486个单基因簇(Unigenes), N50为1000 bp。通过比较转录组分析在2个时间点共获得2072个差异表达基因,其中氨氮胁迫12h获得1516个差异表达基因, 907个上调, 609个下调;氨氮胁迫24h获得556个差异表达基因,其中330个上调, 226个下调。GO功能注释分析结果表明,在氨氮胁迫12h差异基因富集数目明显多于24h,但在“氧化还原酶活性”条目中富集的差异基因数目则是氨氮胁迫24h多于12h。此外, KEGG功能富集分析发现,氨氮胁迫对大口黑鲈肝脏氧化应激、细胞自噬和凋亡相关等途径产生影响。选取10个差异表达基因进行qPCR验证,其表达水平与转录组差异基因分析结果一致,证明了测序结果的可靠性。其中,与氧化应激相关基因缺氧诱导因子1α (HIF1α)、生长停滞DNA损伤可诱导蛋白β (Gad...     

眼优角蚱两性成体转录组差异分析

为了解蚱类昆虫两性间基因转录水平上的整体差异,本研究以眼优角蚱Eucriotettix oculatus为研究对象,通过高通量转录组(RNA-Seq)测序和分析技术对其雌、雄成体的转录组进行了测定和分析。结果显示,在雌、雄转录组中共获得13 432条转录本,两性间差异表达基因共3 597个,其中雌性相对于雄性显著上调基因1 295个,显著下调基因2 302个。对其中9个差异基因进行实时荧光定量PCR验证,其结果与RNA-Seq分析一致。差异基因显著富集的GO条目大多是与细胞的形成有关,结合、细胞组分和细胞过程等条目相关基因显著高表达。KEGG富集显示差异基因主要是与性腺发育、遗传信息加工、能量代谢和消化等通路有关。研究结果对从分子水平进一步研究眼优角蚱两性差异提供一定的理论基础。     

茶树种子发育过程的转录组分析

该研究以‘铁观音’茶树品种的种子为试验材料,采用转录组测序技术分析种子发育的3个时期(幼果期、膨大期、成熟期)的表达差异,探究茶树种子油脂代谢的分子机制。结果表明:(1)经转录组测序、组装后共获得30 940 581个clean reads,经数据合并拼接最终得到36 951条非冗余Unigene序列,其中28 476个Unigene可得到功能注释;在转录本中能够被注释到GO分类的Unigene有11 201条(30.3%),KEGG分析发现共有17 172个基因参与了127个代谢通路。(2)经KEGG通路筛选出14条与脂肪酸代谢相关的通路,且随着茶籽的发育,大部分脂肪酸调控途径相关基因呈下调趋势,其中上调基因数最多的有α-亚麻酸代谢途径和脂肪酸降解途径(有17个基因表达量上调),下调基因数最多的是甘油磷脂代谢途径(有58个基因表达量下调);在茶籽发育幼果期α-亚麻酸代谢途径中表达量上调的基因数超过表达量下调的基因数。(3)研究发现茶籽脂肪酸合成相关的基因涉及14个脂类调控途径,共409条差异基因;随着茶树种子发育到成熟期,上调的差异表达基因数量在减少,下调的差异表达基因数量增加,其中α-亚麻酸途径中的基因PLA2G16、DAD1、pldA、FabF、FabI表达量上调显著,随后表达量下调。(4)qRT-PCR检测结果表明,7个茶树FAD和1个ACP差异表达基因的水平与转录组测序结果基本一致;随着茶籽的发育,基因CsFAD7和Δ6-CsFAD从幼果期、果实膨大期至果实成熟期都为差异下调表达,CsFAD2、CsFAD6和Δ7-CsFAD为差异上调表达,CsFAD8、Δ8-CsFAD和CsACP在幼果期至果实膨大期差异上调表达,在果实膨大期至果实成熟期差异下调表达。     

基于高通量转录组测序的草鱼雌雄性腺差异表达基因分析

草鱼雌性个体生长优势明显高于雄性,挖掘精巢和卵巢差异表达的功能基因对草鱼生殖调控具有重要的价值及应用前景。采用Illumina HiSep 2500转录组测序技术分别对12月龄的3尾雌性草鱼和3尾雄性草鱼的性腺组织的mRNA进行测序分析,精巢和卵巢共得到8 363个差异表达基因。与卵巢相比,精巢上调基因6 446个(77%),下调基因1 917个(23%),差异明显。将上述差异基因在Nr、GO和KEGG数据库进行比对,获得48种Go功能注释分类和313条KEGG代谢通路。分析发现,Go功能分类中"绑定"功能富集差异基因的数量最多,占总数的15.8%;KEGG代谢通路中"信号转导"通路富集差异基因的数量最多,占总数的10.2%。结合转录本数据,挑选与草鱼性腺发育相关的关键差异基因Dmrt1、amh、foxl2、cyp19a1a进行分析,研究显示,与卵巢相比,Dmrt1和amh基因在草鱼精巢中极显著高表达(p0.01);与精巢相比,cyp19a1a和foxl2基因在草鱼卵巢中显著高表达(p0.05),与转录组数据分析结论一致,说明Dmrt1和amh基因与早期草鱼精巢发育调控相关,cyp19a1a和foxl2基因与早期草鱼卵巢发育调控相关。本实验为进一步了解草鱼性别基因的调控机理提供数据依据。     

应用CRISPR-Cas9技术构建NFIB敲除的膀胱癌细胞株及其转录组学分析

为了探究NFIB(nuclear factor I/B)基因在膀胱癌中的生物学功能,该研究应用CRISPR-Cas9技术构建NFIB敲除的膀胱癌细胞株,然后通过转录组学技术分析对照细胞NFIB-NC与NFIB敲除细胞株NFIB-KO-I和NFIB-KO-II之间的差异表达基因,同时对其进行KEGG富集分析,并通过Western blot和qRT-PCR技术对转录组学结果进行验证。结果显示,与对照细胞NFIB-NC相比,NFIB敲除细胞株NFIB-KO-I中有134个差异表达基因,其中上调基因62个,下调基因72个;NFIB-KO-II中有131个差异表达基因,表达上调和下调的基因分别为50个和81个。KEGG分析结果显示,差异表达基因与PI3K-AKT信号通路密切相关,Western blot结果证实敲除NFIB基因后,膀胱癌细胞中的p-AKT水平显著上调。qRT-PCR结果表明,敲除NFIB后,ITGA4基因表达水平上调,TNC和ANGPT4基因表达水平下调。该研究初步揭示了NFIB基因在膀胱癌细胞中调控的分子信号通路,也为后续研究NFIB基因介导膀胱癌发生发展的分子机制提供了依据。     

烟草脉斑驳病毒侵染本氏烟的转录组研究

为了研究烟草脉斑驳病毒（TVMV）侵染本氏烟的抗病分子机制,本文建立了TVMV转录组数据库,挖掘抗病相关基因以及代谢和信号通路,采用Illumina Novaseq 6000高通量测序平台对侵染TVMV的本氏烟进行转录组测序,并进行生物信息学分析;利用R package:edge R等软件分析了有关抗病的差异表达基因,对差异表达基因进行基因本体论（GO）及京都基因和基因组百科全书（KEGG）基于途径的通路富集分析。基于TVMV侵染本氏烟转录组的测序结果中共获得4 593个差异表达基因,其中3 564个基因上调表达,1 029个基因下调表达。共筛选出10个抗病相关的差异表达显著基因,6个上调,4个下调。GO数据库中注释到生物学过程、细胞组分及分子功能等三大类共50个功能组。其中,在第三大类分子功能中,蛋白质结合功能以及ATP结合功能所涉及的差异表达基因最显著,分别为501个和453个差异基因。KEGG共富集20条通路,注释到核糖体途径的差异基因富集程度最高,基因最显著,差异表达基因数为307个。蛋白质结合、ATP结合、核糖体、真核生物核糖体的生物反应、DNA复制（DNA replicat...     

镧对铜胁迫下水稻转录组模式的调控分析

探究稀土镧对铜胁迫下水稻转录组的影响,鉴定镧在调控水稻铜胁迫应答中的关键基因和功能。利用抗氧化酶活性筛选最适铜胁迫及镧处理浓度,进行转录组测序和差异表达分析,并用qRT-PCR验证差异表达水平。对测序数据进行差异表达分析发现3 222个基因显著上调,3 798个显著下调,qRT-PCR验证了4个细胞壁防御相关基因CHIT13、Laccase、Expansin和GH3.4。功能和通路分析获得95组GO条目和112条KEGG通路,富集于激酶和氧分子结合等分子功能、细胞壁及质膜等细胞组分、次级代谢和脂质代谢等生物学过程以及苯丙烷生物合成和植物激素信号转导通路。镧通过调节细胞壁形成和组分提高水稻铜胁迫耐受性。     

基于单细胞转录组的基因富集方法分析星形胶质细胞与阿尔茨海默病的关系*

目的:通过生物信息学方法分析阿尔茨海默病(Alzheimer disease, AD)中与星形胶质细胞相关的糖代谢通路,为揭示AD患者的星形胶质细胞在大脑中的糖代谢过程提供理论基础。方法:首先根据细胞特异性表达基因将AD患者和健康人脑组织单细胞转录组学测序结果进行降维分析,再根据星形胶质细胞不同亚型的基因表达特征进行细胞分群,对星形胶质细胞差异表达基因进行基因注释(Gene Ontology. GO)、信号通路分析(Kyoto Encyclopedia of Genes and Genomes, KEGG)以及基因集富集分析(Gene Set Enrichment Analysis, GSEA),采用转录调控网络分析与AD的星形胶质细胞相关的转录辅助因子。结果:所有细胞降维分析结果显示AD患者脑内星形胶质细胞和兴奋性神经元数量显著减少;星形胶质细胞降维分析结果显示其可以被进一步分为6个亚群,其中在AD患者中减少的星形胶质细胞主要为RASGEF1B⁺SLC26A3⁺亚群和NRGN⁺CALM1⁺亚群;GO分析结果显示AD患者与健康对照星形胶质细胞差异表达基因主要与轴突发生、神经元的迁移、胶质细胞分化、体内锌离子稳态、突触传递的正调控、血管运输有关。KEGG结果显示,上述差异基因主要与PI3K-Akt信号通路、AMPK信号通路、钙信号通路有关。GSEA分析结果显示,AD患者差异基因在糖酵解/糖异生通路中得到富集,其中丙酮酸激酶PKM、PFKL、ACSS1、乳酸脱氢酶LDHB在AD患者星形胶质细胞中下调。转录调控网络分析结果显示,星形胶质细胞中差异表达转录辅助因子有5个,其中PKM、SOX2、SOX9在AD患者星形胶质细胞中下调。SREBF1和BCL6在AD患者星形胶质细胞中上调。结论:AD患者脑内兴奋性神经元和星形胶质细胞数量降低,以及星形胶质细胞糖酵解相关基因下调。结合星形胶质细胞作为神经元的主要乳酸供应细胞,其数量减少和糖酵解能力减低提示星形胶质细胞供能不足可能是AD发生的机制之一。     

TYRO3~(-/-)AXL~(-/-)MER~(-/-)小鼠骨髓细胞mRNA的差异表达谱

TYRO3、AXL和MER受体是受体酪氨酸激酶亚家族之一,广泛地表达在哺乳动物神经、免疫、生殖、血液等系统多种细胞中,促进细胞存活、增殖和分化。本研究利用基因芯片技术筛选基因敲除的TYRO3~(~(-/-))AXL~(-/-)MER~(-/-)小鼠骨髓细胞中差异性表达的m RNA,并进行生物信息学分析。与正常小鼠骨髓细胞相比,TYRO3~(-/-)AXL~(-/-)MER~(-/-)小鼠骨髓细胞中差异表达基因探针共1 363条,其中表达上调的基因探针499条、下调的基因探针864条。GO功能富集分析发现上调的基因主要参与免疫反应,下调的基因主要参与造血。KEGG Pathway富集分析发现上调的基因主要参与细胞粘附分子和抗原呈递等信号通路。Real-time PCR验证5个差异基因,与芯片中表达变化趋势一致。因此TYRO3、AXL和MER受体共同参与调控机体的免疫反应和血细胞分化。     

小鼠正常黑色素细胞与黑色素瘤细胞(B16)mRNA表达差异分析

黑色素瘤是一种高侵袭性的恶性皮肤肿瘤,转移率高、预后差。研究黑素瘤细胞生物学特性对黑素瘤的治疗和控制具有重要的意义。本研究以C57BL/6J小鼠的正常黑色素细胞及B16黑色素瘤细胞为研究对象,采用二代测序技术分析两种细胞间的转录组表达差异,筛选差异基因,为后续黑色素瘤的形成机制研究提供理论依据。采用差异倍数及错误率分析测序数据,鉴定出1 436个新的mRNA和4 086个差异表达的已知mRNA。GO数据库和KEGG数据库分析显示,差异表达的mRNAs参与了149个调控途径,主要集中在疾病调控、细胞周期调节和环境信息调控方面。qRT-PCR及Western印迹检测发现,调节细胞增殖、迁移的Pdgf-B、Integrinβ1和Integrinβ5以及调节黑色素颗粒增加的Mitf、Tyr、Tyrp1和Tyrp2在B16细胞中的表达量显著高于在正常黑色素细胞中的表达。本研究获得的差异基因为后续黑色素瘤的研究提供了新的候选基因。     

ABCA1对滋养细胞基因表达谱的影响

该文从基因学角度探讨了ABCA1对人滋养细胞功能的影响。在人滋养细胞中调控ABCA1的表达,通过表达谱芯片检测ABCA1表达改变后滋养细胞中基因及相关信号通路的改变,Western blot及qRT-PCR验证芯片结果。结果显示,上调ABCA1表达时,有197个基因表达升高,有190个基因表达下降;而下调ABCA1表达时,有335个基因表达升高, 459个基因表达下降。GO和KEGG分析表明, ABCA1表达上升或下降后可导致滋养细胞内多个信号通路发生改变。qRT-PCR及Western blot检测发现,与阴性对照组相比, ABCA1上调后细胞中S1PR1的m RNA及蛋白水平明显升高,而CCL8、CXCL10与CXCL11的mRNA及蛋白水平明显降低;下调ABCA1的表达后,细胞中S1PR1的mRNA及蛋白水平明显降低,而CCL8、CXCL10与CXCL11的mRNA及蛋白水平明显升高,与芯片的结果完全一致。这些结果表明, ABCA1可通过调控多个信号通路而影响滋养细胞功能。     

小鼠正常黑素细胞与黑色素瘤细胞（B16）mRNA表达差异分析

黑色素瘤是一种高侵袭性的恶性皮肤肿瘤,转移率高、预后差。研究黑素瘤细胞生物学特性对黑素瘤的治疗和控制具有重要的意义。本研究以C57BL/6J小鼠的正常黑色素细胞及B16黑色素瘤细胞为研究对象,采用二代测序技术分析两种细胞间的转录组表达差异,筛选差异基因,为后续黑色素瘤的形成机制研究提供理论依据。采用差异倍数及错误率分析测序数据,鉴定出1 436个新的mRNA和4 086个差异表达的已知mRNA。GO数据库和KEGG数据库分析显示,差异表达的mRNAs参与了149个调控途径, 主要集中在疾病调控、细胞周期调节和环境信息调控方面。qRT-PCR及Western印迹检测发现,调节细胞增殖、迁移的Pdgf-B、Integrin-β1和Integrin-β5以及调节黑色素颗粒增加的Mitf、Tyr、Tyrp1和Tyrp2在B16细胞中的表达量显著高于在正常黑色素细胞中的表达。本研究获得的差异基因为后续黑色素瘤的研究提供了新的候选基因。     

阿尔泰蝠蛾(鳞翅目:蝙蝠蛾科)幼虫低温胁迫转录组分析

[目的]本研究旨在通过转录组测序技术分析低温胁迫引起的阿尔泰蝠蛾Hepialus altaicola Wang幼虫基因转录差异及上调表达基因的主要功能类群和参与的主要代谢通路.[方法]采用Illumina HiSeq TM 2500测序平台对阿尔泰蝠蛾幼虫进行转录组测序、组装,利用Blast软件进行数据库比对和基因功能注释,用DEGSeq R软件包分析4℃低温处理与室内适温饲养试虫的差异表达基因,并对上调表达基因进行GO和KEGG代谢途径富集分析.[结果]经序列拼接后共获得100300个unigenes,总长度81600309 bp,平均长度813 bp,N50长度1719 bp.与7大数据库同源比对,共获得34691(34.59％)条unigenes.低温胁迫转录组分析得到11569个差异表达基因(DEGs),7158条基因上调,4411条基因下调.富集到47个GO类群,217个KEGG途径.其中代谢过程、催化、结合活性类群占有重要比例,核糖体、碳代谢、剪接体等途径显著富集.另外,热激蛋白、昆虫表皮蛋白、海藻糖酶、超氧化物歧化酶等非生物胁迫相关基因显著上调表达.[结论]阿尔泰蝠蛾幼虫低温胁迫转录组分析揭示,代谢过程、细胞过程、生物调节、对刺激的反应等生物学过程相关基因和部分非生物胁迫响应基因显著上调表达,提示蝠蛾幼虫可能从抗氧化防御、分子伴侣、体温调节和维持细胞的渗透平衡等多方面应对低温胁迫.