基于高通量测序技术分析家兔感染猪瘟兔化弱毒株后脾脏的转录组学变化

本研究基于高通量测序技术对感染猪瘟兔化弱毒株(C株)后0h,12h,24h,30h,36h,48h的家兔脾脏组织进行测序,以家兔基因组作为参考,筛选感染后各时间点和对照组(0h)之间的差异基因。利用Uniprot和NCBI数据库查找各组中变化水平显著的前10个差异基因的生物学功能。结果发现,感染后12h,24h,30h的前10个差异基因只有少数与炎症、免疫、凋亡等相关;感染后36h,48h的前10个差异基因完全相同,其中B2M、RLA-DR-ALPHA、CD74和IGJ参与抗病毒过程。GO功能注释结果显示:感染后各时间点的差异基因都和免疫、代谢、调控途径紧密相关。KEGG通路富集分析发现,感染后24h的差异基因富集到RIG-I样受体信号通路,感染后30h的差异基因富集到黏着斑和细胞外基质-受体相互作用通路,感染后36h和48h的差异基因共同富集到的通路有20个,其中蛋白酶体,溶酶体,核糖体,趋化因子信号通路,B细胞受体信号通路,抗原加工提呈通路和FcγR介导的吞噬通路与抗病毒相关。本研究建立的家兔感染C株后脾脏组织的差异表达基因数据库,将为进一步阐述家兔适应C株的分子机制奠定基础。     

红鳍东方鲀不同生长时期脑和肌肉的转录组比较分析

红鳍东方鲀(Takifugu rubripes)是中国北方重要的海水养殖鱼类.本研究通过Illumina测序技术对13月龄和15月龄红鳍东方鲀的脑组织和肌肉组织共12个样本进行了转录组测序与分析.结果表明:在脑组织样本中平均共得到5 510万条原始reads,在肌肉组织样本中平均共得到5 482万条原始reads.在15月龄相对13月龄脑(BB15 vs BB13)中筛选到554个差异基因,其中234个上调表达基因,320个下调表达基因.15月龄相对13月龄肌肉(MB15 vs MB13)中筛选到380差异基因,其中253个上调表达基因,127个下调表达基因.GO富集结果将差异基因分为生物过程、细胞组成、分子功能三个层面.根据KEGG代谢通路分析,在BB15 vs BB13中,注释到163个差异基因,共富集到82条通路上,其中MAPK信号通路和神经信号传递通路的富集程度最高.在MB15 vs MB13中注释到138个差异基因,富集到97条通路上,其中也是MAPK信号通路富集程度最高.共筛选出3个可能与生长相关基因MSTN2、IER2、C-FOS.研究结果为今后开展红鳍东方鲀的基因研究提供了有价值的参考信息,也为改善养殖红鳍东方鲀的生长、抗逆和抗病能力的研究提供科学依据.     

基于转录组的八角挥发油合成相关基因挖掘与分析

为更好地挖掘八角(Illicium verum)挥发油合成相关基因，该文对挥发油性状差异显著的优良无性系桂角69号及普通品种砧01号叶片进行了转录组测序及组装注释，并对差异表达基因进行了GO分类和KEGG通路分析。结果表明：(1)转录本经组装后获得84 182条序列，使用NR、NT、Swiss-Prot、KEGG、KOG、GO和Pfam数据库进行序列比对，共注释了59 161条序列，筛选出30 572个差异表达基因。与砧01号相比，桂角69号叶片中上调基因有15 025个，下调基因有15 547个。(2)GO分类结果显示共有20 287个差异基因被注释。KEGG分析结果表明，有21 600个差异基因被注释到133条KEGG通路上，其中挥发油合成相关的单萜生物合成通路、萜类骨架生物合成通路、苯丙素合成通路中的芳樟醇合酶、月桂烯合酶、香叶基香叶基焦磷酸合酶、肉桂酰辅酶A还原酶、咖啡酸3-O-甲基转移酶、肉桂醇脱氢酶等关键酶基因呈差异表达。(3)转录因子分析发现差异表达基因分布于31个转录因子家族，其中MYB家族序列数量最多。该文利用转录组测序技术分析八角优良无性系与普通品种叶片的差异基因及...     

基于高通量测序的文冠果转录组分析

应用Illumina Hi-seqTM2000高通量测序技术对文冠果花芽进行转录组分析。共获得N50为1 180bp、平均长度为686bp的unigene 58 311条。与公共数据库Nr和Swiss-Prot同源性比较后发现37 047条unigene获得基因注释,另有21 264条unigene未被注释。利用COG数据库将unigene分成25类。通过GO分类和KEGG Pathway富集性分析,将unigene分别归类于55个GO term和128个代谢途径。此外,在9 794条unigene中共搜索到12 213个SSR位点,单核苷酸重复基元出现频率最高(34.95%),其次分别为二核苷酸(32.74%)和三核苷酸(28.64%)。在获得的unigene中发掘出涉及4个开花调控途径(光周期途径、春化途径、GA途径和自主途径)多个基因的同源序列。研究结果可在一定程度上解析文冠果花芽形态分化的分子调控模式与机制。     

草地贪夜蛾雄性成虫和5龄幼虫的转录组比较分析

【目的】草地贪夜蛾Spodoptera frugiperda是一种新近入侵我国的重要害虫。本研究旨在对草地贪夜蛾雄性成虫和5龄幼虫两个不同发育阶段的转录组进行比较分析。【方法】利用高通量测序技术对草地贪夜蛾雄性成虫和5龄幼虫进行转录组测序和数据组装,并对转录组数据进行功能注释和比较分析。【结果】经de novo组装共获得209 002条转录本,平均长度为687.55 bp,N50为982 bp。共有46 198条(57.43%) unigene在至少一个数据库中获得功能注释,其中1 713条(2.13%) unigene在所有数据库中均能获得注释。在GO数据库中获得205 269条unigene的注释,主要包括68个功能分类;在KEGG数据库中共有3 408条unigene得到注释,涉及277个代谢通路。共鉴定到424个嗅觉相关的基因,并且在雄性成虫和5龄幼虫之间的表达存在差异。通过比较转录组分析,在雄性成虫中鉴定到9 162个上调和6 399个下调差异表达基因(DEGs);功能富集分析发现在上调DEGs中涉及信息素以及信号转导的代谢通路显著富集,而下调DEGs中涉及解毒相关的通路显著富集。【结论】这些转录组数据为探究草地贪夜蛾的生长发育、嗅觉相关功能基因以及候选分子靶标提供资源信息。     

不同大小规格墨西哥湾扇贝(Argopecten irradians concentricus)转录组分析及生长相关基因筛选

墨西哥湾扇贝(Argopecten irradians concentricus)是中国双壳贝类养殖的重要经济品种,具有较高的研究价值,然而,转录组学及基因组学信息的匮乏,严重制约了该贝的遗传改良进程。本研究基于Illumina Hi SeqTM2000平台,对5月龄大规格(平均体质量16.05 g)和小规格(平均体质量5.57 g)墨西哥湾扇贝进行转录组测序,经de nove组装,大、小规格墨西哥湾扇贝分别获得120 563和118 794条高质量unigenes,unigene的平均长度分别为944 bp和969 bp;按照FDR≤0.001且|log2Ratio|≥1的原则,共筛选出9 901个差异表达基因,其中,大规格组相对于小规格组,上调基因4 976个(50.26%),下调基因4 925个(49.74%);六大公共数据库Blast比对结果显示,差异表达基因在Nr数据库的注释率最高,其次是GO数据库;从已获功能注释的差异表达基因中筛选得到47个生长相关候选基因,涉及多个生长因子及其受体;GO功能富集结果显示,膜、L-氨基酸运输和糖基键水解酶活性分别是细胞组分、生物过程和分子功能富集unigene最多的条目;Pathway富集分析发现,3 699条差异表达基因被成功注释到252条代谢通路中,其中33条为显著性富集通路(p<0.05),富集度最高的是焦点粘附;大规格组中共鉴定出10 870个SSRs和242 551个SNPs,小规格组中共鉴定出10 857个SSRs和228 121个SNPs。研究结果为墨西哥湾扇贝分子标记批量开发、功能相关基因挖掘以及分子遗传育种提供了理论数据支撑。     

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

为了探究京海黄鸡感染柔嫩艾美尔球虫(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等,提示这些基因在宿主柔嫩艾美耳球虫感染过程中发挥重要作用。     

绿豆象触角转录组及嗅觉相关基因的分析

【目的】建立绿豆象Callosobruchus chinensis触角转录组数据库,挖掘绿豆象的基因数据。【方法】采用高通量测序平台(Illumina Hiseq)对绿豆象成虫触角进行转录组测序、序列组装及生物信息学分析。【结果】共获得51.10 Gb的clean data(NCBI SRA数据库登录号:SRP119884)及83 535条unigenes,长度在1 kb以上的unigenes有15 075条,unigenes的N50为1 492 bp。使用BLAST软件将绿豆象触角Unigene序列与公共数据库比对,共注释到22 148条unigenes,其中NR数据库注释的unigenes最多,为18 744条,且与赤拟谷盗Tribolium castaneum的相似基因最多,达39.57%。通过GO数据库注释,将unigenes功能分为细胞组分、分子功能与生物学过程三大类54个分支,其中参与代谢进程以及结合与催化活性的unigenes比例较大。KEGG代谢途径分析表明,9 084条unigenes形成289条代谢通路,其中核糖体代谢通路所含unigene最多,为516条。进一步基因注释分析筛选到140个嗅觉相关基因,包括12个气味结合蛋白(odorant binding protein,OBP)基因,4个化学感受蛋白(chemosensory protein,CSP)基因,116个气味受体(odorant receptor,Or)基因,1个离子型受体(ionotropic receptor,IR)基因和7个感觉神经元膜蛋白(sensory neuron membrane protein,SNMP)基因,并用FPKM值对基因表达量进行评估。【结论】本研究获得了绿豆象触角转录组数据,为进一步研究绿豆象的基因功能分析及嗅觉感受机制奠定分子基础。     

假单胞菌减毒活疫苗免疫后的大黄鱼脾脏转录组分析

大黄鱼是中国东部沿海养殖的重要经济鱼类品种,多年来受到假单胞菌引起的内脏白点病的危害,疫苗免疫是预防该病的理想手段,减毒活疫苗(live attenuated vaccine,LAV)或DNA疫苗等能有效刺激宿主细胞免疫应答的疫苗可能是该菌疫苗发展的重要方向。为了评估减毒活疫苗的有效性,解悉大黄鱼对疫苗的免疫应答机制,本研究以重要毒力基因缺失的减毒活菌株ΔpcrV NB2011免疫鱼体,在免疫后0 h和48 h分别采取实验鱼的脾脏组织进行转录组测序。通过美吉云平台在线分析,一共注释了30 438条Unigene和70 618条转录本。筛选出上调和下调表达基因分别为311个和313个,其中免疫相关的显著性差异基因28个,涉及抗原加工呈递、补体系统、细胞因子、模式识别受体和吞噬作用,尤其是MHCI型抗原加工和递呈通路显著上调。对上调表达基因进行KEGG富集分析,涉及的最主要通路包括谷胱甘肽代谢通路、精氨酸-脯氨酸代谢通路、抗原提呈与加工通路和溶菌酶通路等。选取5个显著性差异表达的免疫相关基因进行qRT-PCR验证,验证结果与测序结果变化趋势基本一致,提示了LAV免疫后有效激发鱼体的细胞免疫应答,有助于快速清除胞内病原。本研究部分揭示了大黄鱼对致病性假单胞菌LAV的免疫应答机制,可以为疫苗的研制提供科学依据。     

东亚特有珍稀蕨类植物岩穴蕨(碗蕨科)高通量转录组测序及分析

岩穴蕨(Monachosorum maximowiczii)隶属于碗蕨科稀子蕨属, 是东亚中高海拔地区所特有的珍稀濒危植物。为了在分子水平对该物种有进一步的认识, 本文首次利用二代高通量测序技术(RNA-seq)对岩穴蕨进行了转录组测序分析。通过Illumina Hiseq 2500测序平台, 共计获得4.95 Gb原始数据(raw data), 经过滤后得到4.83 Gb有效数据(clean reads), 并进行从头组装得到了101,448条unigene。其中, 54,106条unigene预测到完整的开放阅读框。我们利用目前已知的51个植物基因组数据, 对岩穴蕨的unigene进行了详尽的功能注释, 并通过GO、COG、KEGG注释进一步了解了这些编码基因的作用方式、特征以及所参与的代谢通路。同时, 转录因子分析结果也为岩穴蕨的环境适应机制研究提供了初步线索。     

多年宿根甘蔗的转录组分析及差异基因挖掘

甘蔗(Saccharum hybirds)宿根性直接关系到甘蔗生产成本和种植效益,品种、环境和栽培措施均会影响甘蔗宿根能力,但品种的种性是影响宿根性最关键因素。国内外从分子生物学层面解释甘蔗宿根性差异的文献鲜见报道,从转录水平分析不同宿根年限GR2号和ROC22号基因表达的差异。结果表明：转录组测序共得到100 558条转录本和25 582条Unigene。获得53 790条Unigene的注释结果,分别在NR、Swiss-Prot、KEGG、COG、KOG、GO和Pfam数据库进行比对。GO功能注释共分成三大类,51小类,6年宿根蔗注释到1 029个差异基因,3年宿根蔗注释到3 391个差异基因。主要KEGG代谢通路有8条,分别涉及植物激素信号转导,淀粉和蔗糖代谢,甘氨酸、丝氨酸和苏氨酸代谢,苯丙素生物合成,同源重组,DNA复制,错配修复及植物病原体相互作用。筛选出脱落酸(ABA)相关差异基因6个,分别为脱落酸不敏感蛋白2基因(ABI2)、bZIP转录因子超家族蛋白、碱性亮氨酸拉链型转录因子基因(ABI5)、aba响应元件结合因子1基因(ABF1)、G-box结合因子基因(GBFs)...     

鸭肠炎病毒感染对鸭肠道黏膜组织转录组变化分析

为了探究鸭肠炎病毒感染对鸭十二指肠黏膜组织相关免疫基因的变化影响,本试验采用鸭肠炎病毒接种35日龄的健康樱桃谷鸭,将其分为3组(2个试验组和1个对照组),在感染后24 h、48 h采集十二指肠,提取总RNA进行浓度、纯度及完整性检测后构建cDNA文库,使用BGISEQ-500平台对2个实验组和对照组黏膜组织测序,筛选差异表达基因,并用GO和KEGG数据库进行分析。结果显示,2个实验组感染鸭肠炎病毒后T24和T48分别筛选出共720个差异表达基因,其中T24表达上调58个,下调163个;T48表达上调77个,表达下调422个。GO功能注释结果显示,感染后差异基因都与代谢、免疫、炎症和调控途径等密切相关。KEGG通路富集分析发现,感染T24的差异表达基因富集到神经活性配体-受体相互作用、补体和凝血级联、氮代谢、cAMP信号通路和IL-17信号转导通路;感染T48的主要差异表达基因富集到神经活性配体-受体相互作用、细胞因子-细胞因子-受体相互作用、趋化因子信号通路、造血细胞谱系和补体和凝血级联。借助BGISEQ-500平台对鸭肠炎病毒感染后十二指肠黏膜组织相关因子的变化来探究其致病机制及其他生物学通路,弥补了非模式生物转录组研究中缺乏基因信息的不足。     

藏茵陈基源植物皱边喉毛花的全长转录组信息分析

皱边喉毛花为藏药藏茵陈基源植物之一,其包含丰富的药用成分。为进一步了解皱边喉毛花转录组,丰富其基因注释、代谢通路等遗传信息,该研究利用PacBio测序平台对皱边喉毛花叶片进行全长转录组测序。结果表明:(1)全长转录组测序共获得17 Gb的高质量数据,对795 698 个CCS序列进行聚类和去冗余,最终获得87 814 条高质量的全长转录本。(2)与7个数据库比对后,共有277 451 条转录本注释成功,其中注释到NR数据库的转录本最多,有39 214 条。26 396 条转录本成功注释到KOG数据库中,共有26 个子类。39 104 条转录本注释到KEGG数据库中,涉及6 个主要通路和40 个子通路。39 102 条转录本注释到GO数据库中,按分子功能、生物学过程和细胞成分3大类对注释成功的转录本进行分类。(3)SSR分析共鉴定到22 861 个SSR,其中单碱基重复最为丰富; 共检测到1 874 个转录因子和15 166 个长非编码RNA(LncRNA),而注释到转录本最多的转录因子家族是C3H。(4)筛选出55 条与单萜类及黄酮类化合物合成相关的转录本。该研究结果丰富了皱边喉毛花的转录组信息,为进一步筛选皱边喉毛花药用成分合成相关的关键基因提供了重要的遗传资源。     

火龙果转录组测序、基因表达与功能分析

利用高通量测序技术对火龙果(Hylocereus undulatus Britt)红肉品种‘大红二号’的花芽、果实和枝条不同发育阶段的基因表达进行研究。结果显示,转录组测序共获得468.68 Gb原始数据(Raw data),从头组装获得239 152条转录本和162 519条unigene,约53.74%的unigene得到注释。分别在43 506条和16 251条unigene中检测到600 283个SNP位点和56 147个SSR位点。基因表达分析结果表明,在火龙果不同组织Fl510、Fl513、Fl514、Fl518、F711、F715、S513、S419中分别有31、7、5、152、17、63、17、8个特异表达的unigene。通过GO和KEGG富集分析,发现了一些组织特异的GO条目和代谢通路,如在Fl510中富集的类萜骨架生物合成代谢通路等。本研究还对参与花发育的候选基因进行了鉴定和表达分析,他们包括COL基因、FT-like基因、分生组织决定基因和器官决定基因等。     

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株细胞的转录谱存在显著差异。     

禽致病性大肠杆菌pagP基因缺失株相关毒力基因的RNA-Seq分析

本研究为筛选禽致病性大肠杆菌pagP基因缺失后表达量变化的相关毒力基因,对禽致病性大肠杆菌AE17及pagP基因缺失株进行RNA-Seq测序,从中筛选出与毒力相关的差异基因并对其表达量的变化进行GO和KEGG分析。以|Fold-change|≥2且FDR≤0.05为条件从测序结果中筛选得到372个差异表达基因,其中339个基因表达上调,33个基因表达下调。GO功能分类结果显示差异基因主要集中在物质跨膜运输、转Class录调控以及生物膜的形成等功能注释;KEGG数据库显示,差异表达基因主要富集在鞭毛装配、双组分系统、脂多糖的生物合成以及ABC转运蛋白等通路。pagP基因缺失会引起相关毒力基因表达量的变化,根据RNA-Seq测序结果筛选出7个对APEC具有重要作用的毒力基因,为深入研究APEC的致病机理提供依据。     

大弹涂鱼皮肤转录组测序及抗菌肽基因分析

抗菌肽是鱼类用于抵御外界微生物入侵的天然防御多肽,也是开发新型药物的重要先导分子。大弹涂鱼(Boleophthalmus pectinirostris)是一种特殊的可以营两栖生活的鱼类,其皮肤表面具有丰富的粘液,其中所含抗菌肽对其免疫防御和适应两栖生活具有重要意义。为深入了解大弹涂鱼皮肤组织基因表达谱并从中筛选抗菌肽相关基因,采用新一代Illumina高通量测序平台对大弹涂鱼皮肤组织进行了转录组测序。利用Trinity软件从头组装,从测得的78 608 366条双端测序读长(paired-end read)共计6 GB的序列数据中获得119 848条高质量的蛋白质编码基因(unigene)。经公共数据库序列检索和比对,发现7个unigene编码的多肽与已知的鱼类5大家族的抗菌肽高度同源,即β-防御素(β-defensin)、hepcidin、NK-lysin、piscidin和肝表达抗菌肽-2(liver-expressed antimicrobial peptide-2,LEAP-2)。最后,对上述抗菌肽相关unigene开展了组织表达差异分析、序列比对及进化树分析。研究结果为进一步了解大弹涂鱼适应两栖生活的免疫防御机制和利用新鉴定的抗菌肽开发新型抗菌药物奠定了基础。     

青藏高原地区3种蒿属植物转录组比较分析

米蒿（Artemisia dalai-lamae）、冷蒿（Artemisia frigida）和臭蒿（Artemisia hedinii）作为西北干旱半干旱地区常见植物,具有防风固沙、阻止草原退化、杀菌消炎和治疗多种疾病的功效。通过Illumina测序技术分别对3种蒿属植物进行转录组测序和分析,分别获得54 268 322,46 434 864和43 971 646条Clean reads,并对测序结果进行了SNP和可变剪接分析。在米蒿和冷蒿对比组（ADL vs AF）中筛选到6 107个差异基因,在米蒿和臭蒿对比组（ADL vs AH）中筛选到4 822个差异基因,在冷蒿和臭蒿对比组（AF vs AH）中筛选到3 755个差异基因。GO富集将差异基因注释到生物过程、细胞组成和分子功能这3个类别共47个条目。KEGG代谢通路分析将3个对比组的差异基因分别富集到198、198和197条通路。对每个对比组的前10个高表达差异基因进行了分析,并筛选出25个生物碱、萜类和类黄酮的活性成分可能相关的差异基因。本研究为蒿属植物的物种鉴定、抗逆性研究和资源利用提供一定的科学依据。     

松树蜂毒腺转录组分析

【目的】构建入侵种松树蜂Sirex noctilio毒腺转录组数据库,筛选并分析松树蜂毒腺基因数据。【方法】采用新一代高通量测序平台Illumina HiSeq^TM 4000对松树蜂雌成虫毒腺进行转录组测序、数据组装和生物信息学分析。【结果】共获得12.7 Gb松树蜂雌成虫毒腺有效转录组数据,并组装到37 098条unigenes,平均长度968 bp,N50长度为2 364 bp。将所得的unigenes数据使用BlastX与各大数据库比对,共注释到13 515条unigenes,并且在NR数据库中注释的unigenes最多,共11 108条(占总数的29.94%),其中相似基因占比最高的物种为丽蝇蛹集金小蜂Nasonia vitripennis,达815条(占总数的7.29%)。在GO数据库中注释到5 726条unigenes,根据功能被分为生物学进程、细胞组分和分子功能3大类63个亚类。KEGG代谢通路分析表明,7 602条unigenes注释到357个代谢通路。根据基因注释信息进一步筛选到43条嗅觉相关基因,包括嗅觉受体(odorant receptor, Or)基因25条、化学感受蛋白(chemosensory protein, CSP)基因10条、离子型受体(ionotropic receptor, IR)基因5条和气味结合蛋白(odorant binding protein, OBP)基因3条。此外,还筛选出11条漆酶基因,包括漆酶1(laccase1, LAC1)基因5条、漆酶2(laccase2, LAC2)基因4条、漆酶4(laccase4, LAC4)基因1条和漆酶9(laccase9, LAC9)基因1条,且其中1条LAC2基因在所有被注释的基因中表达量最高(FPKM值=21 126)。【结论】本研究获得的松树蜂毒腺转录组数据为松树蜂毒液组分的鉴定和生物学功能的研究奠定了一定的理论基础。     

幽门螺杆菌感染与胃癌相关基因的筛选及预后分析

目的通过分析幽门螺杆菌感染胃黏膜组织和胃癌细胞系后的差异基因变化,并在癌症基因组图谱(The Cancer Genome Atlas,TCGA)数据库和肿瘤基因芯片(Oncomine)数据库进行验证,探究幽门螺杆菌导致胃癌发生、发展的分子机制。方法分析基因表达汇编(Gene Expression Omnibus,GEO)数据库幽门螺杆菌感染相关芯片集GSE5081与GSE70394,绘制维恩图查找幽门螺杆菌感染后共同上调的差异基因。对共同上调的差异基因进行功能富集分析。通过TCGA和Oncomine数据库验证差异基因在胃癌中的表达。利用Kaplan-Meier Plotter数据库和GEPIA数据库分析差异基因表达高低与胃癌患者预后是否存在相关性。结果通过差异基因筛选和维恩分析,两个芯片集共有21个共同上调差异基因。GO分析发现共同上调差异基因主要富集在对细菌来源分子的反应、趋化因子CXCR受体结合、中性粒细胞趋化作用等相关的基因功能上;KEGG通路主要富集在癌症通路、TNF信号通路、趋化因子信号通路等。STRING以及PPI数据库分析发现21个基因中PRDM1、IL10、NRP1、BIRC3、GNG13、CXCL1、CXCL2、CXCL3、CXCL8基因存在有网络关系,属于关键枢纽基因。通过TCGA和Oncomine数据库筛选及验证,发现在胃癌组织中NRP1、CXCL1、CXCL8基因明显上调,结果差异有统计学意义(TCGA数据库中,三者P值均小于0.05,Oncomine数据库中,NRP1:t=4.607,P0.01;CXCL1:t=5.854,P0.01;CXCL8:t=5.316,P0.01)。在Kaplan-Meier Plotter数据库(210615-at芯片:P0.01;210510-s-at芯片:P0.01;212298-at芯片:P0.01)以及GEPIA数据库(P0.01)两个数据库中,NRP1的高表达均与胃癌的预后负相关。结论不同的数据库均显示NRP1、CXCL1、CXCL8三个基因与幽门螺杆菌感染相关,同时在胃癌中高表达,并且NRP1的高表达与胃癌的不良预后相关,这些结果为进一步探究幽门螺杆菌导致胃癌发生、发展的分子机制提供了重要基础。