禽巴氏杆菌四环素耐药株转录组测序与分析

为鉴定禽源多杀性巴氏杆菌(avian Pasteurella multocida,Pm)四环素(Tetracycline,Tet)敏感株与耐药株(MIC=32μg/mL)在转录组水平上的差异,本研究采用RNA-seq技术对Tet敏感株与耐药株进行双向测序,通过GO、KEGG和ARDB数据库对表达差异基因进行注释和富集性分析。结果显示,共筛选出442个差异表达基因,其中包括103个上调基因和339个下调基因。通过GO功能富集分析显示差异表达基因多数富集于核糖体结构成分、rRNA结合、跨膜运输和离子转运等生物过程;对差异表达基因进行信号通路富集性分析,大多数差异表达基因分布在ABC转运蛋白与核糖体代谢途径中;对差异表达基因进行耐药基因注释(ARDB),发现其中的耐药基因主要为核糖体保护蛋白和外排系统两类。本研究表明禽源Pm对四环素耐药的耐药机制可能与细胞膜上ABC转运蛋白的过量表达以及核糖体保护蛋白有关。     

禽巴氏杆菌链霉素耐药株转录组测序与分析

为鉴定禽源多杀性巴氏杆菌(Pm)链霉素(Str)敏感株与耐药株(MIC=1024μg/mL)在转录组水平上的差异,本研究采用RNA-seq技术对Str敏感株与耐药株进行双向测序,筛选得到差异表达基因,通过GO、KEGG和ARDB数据库对差异表达基因进行分组和富集性分析,并利用荧光定量PCR对筛选出的差异表达基因进行验证。结果显示,共筛选获得625个差异表达基因,其中包括581个上调基因和44个下调基因。通过GO功能富集分析显示差异表达基因多数发挥细胞膜转运功能;对差异表达基因进行信号通路富集性分析显示,大多数差异表达基因分布在ABC转运蛋白与双组分系统代谢途径中;采用ARDB对差异表达基因进行耐药基因分析,显示其中的耐药基因主要为氨基糖苷转移酶类和外排系统两类。本研究结果表明,禽源Pm对Str的耐药机制可能与细胞膜上多药外排泵的过量表达和氨基糖苷转移酶类灭活抗菌药物有关。同时,也为进一步研究禽源Pm对Str的耐药机制奠定了基础。     

捻转血矛线虫阿苯达唑敏感株和耐药株比较转录组学分析

为探明捻转血矛线虫阿苯达唑敏感株和耐药株在转录组水平的差异,挖掘耐药相关基因,本研究采用Illumina HiSeq 4000对敏感虫株和耐药虫株进行转录组测序,构建cDNA文库,通过质控后,对其进行生物信息学分析。结果显示:敏感和耐药虫株相比较共获得218个显著(P≤0.05)差异表达基因,对显著差异表达基因进行GO功能富集,有121、82和102个基因分别注释到生物学过程、细胞组分和分子功能三大类的206个GO terms上。KEGG数据库分析显示,有42个显著差异表达基因参与31条KEGG信号通路,显著富集在药物代谢-细胞色素P450、药物代谢-其他酶等通路。在敏感株和耐药株中,筛选出GST、CYPs和UGT等9个耐药相关基因,这些基因在捻转血矛线虫耐药中发挥重要的作用。本研究利用RNA-Seq预测出敏感虫株和耐药虫株差异表达基因在功能分类和代谢通路等方面的生物学特征,挖掘出耐药相关基因,进一步为捻转血矛线虫药物靶点预测及耐药性早期鉴别诊断的建立提供重要基础。     

美国短毛黑水貂快速生长过程中转录组差异表达分析

旨在对美国短毛黑水貂快速生长过程中转录组差异表达进行分析。本试验将美国短毛黑水貂45、90日龄各3只健康公貂的胸肌组织作为试验材料,利用Illumina HiSeq~(TM)2500高通量测序平台建立转录组文库,对两个生长阶段差异显著性的基因进行GO功能富集和KEGG Pathway分析,寻找调控水貂快速生长期的相关候选基因和代谢通路,并利用qRT-PCR对转录组测序结果进行验证。结果表明,以P0.05,|log_2FC|1为条件筛选到279个显著差异基因,对这些差异基因进行GO和KEGG分析,筛选到与肌肉生长相关显著富集GO条目有66条、相关差异表达基因44个,其中上调20个,下调24个;与肌肉生长相关显著富集KEGG通路12条,如p53信号通路、PPAR信号通路、FOXO信号通路等。这些差异基因可能在水貂的快速生长过程中起到调控作用,可以作为后续研究水貂生长发育的候选基因。本研究结果为探究水貂生长发育调控机制及培育大体型水貂品种奠定基础。     

鼠伤寒沙门菌伴侣蛋白Hfq调控基因的筛选

为筛选鼠伤寒沙门菌伴侣蛋白Hfq所调控的基因,本实验基于本研究室前期对指数生长期和稳定期的鼠伤寒沙门菌hfq基因缺失株转录组测序结果,从中随机挑选6个在转录水平差异表达的基因,并通过荧光定量PCR技术验证转录组数据的可靠性,然后在转录组数据可靠的前提下将两个转录组中所有显著差异表达(log_2FoldChang1)的基因GO富集细菌致病机制生物学过程,Pathway富集细菌ABC转运通路和群体感应通路,在富集得到的基因中筛选指数生长期和稳定期均表现为转录水平差异表达且趋势相同的基因。荧光定量PCR检测结果显示,6个基因转录水平差异上调或下调的趋势与转录组分析结果中一致,表明转录组数据可靠,可进行后续的研究。GO与Pathway富集结果显示,在鼠伤寒沙门氏菌hfq基因缺失株指数期和稳定期中均表现为转录水平差异上调或下调且差异倍数均大于2的基因共26个。其中,与细菌致病机制生物学过程相关的13个基因均上调;与细菌ABC转运通路相关的基因11个,其中9个基因转录水平均上调2个基因转录水平均下调;与细菌群体感应相关的2个基因转录水平均上调。26个基因中仅rpoE经验证受Hfq直接调控。以上结果表明,Hfq主要参与细菌致病机制生物学过程以及细菌ABC转运通路的负调控。本研究为后续Hfq靶基因的鉴定及其作用机理的研究提供了理论基础,同时丰富了Hfq调控基因的数目。     

初生及成年牦牛大脑皮质的转录组分析

旨在分析初生和成年牦牛大脑皮质转录组表达谱,筛选影响大脑皮质发育的候选基因。本研究选取初生(1～7日龄(NYB),n=3)和成年(3～4岁(AYB),n=3)牦牛,放血处死,采集大脑皮质,利用Illumina HiSeq平台对转录组进行测序与分析,筛选得到差异表达基因,在此基础上对差异基因进行GO功能注释、KEGG富集分析和实时荧光定量PCR(qRT-PCR)方法验证。测序结果共分析得到了4 790个显著差异表达的基因(P<0.05),其中包括2 670个上调表达的差异基因,2 120个下调表达的差异基因。随机选取了9个差异基因进行qRT-PCR验证,表达趋势与转录组测序结果一致,表明测序结果可靠。GO功能注释发现其中占比较多的条目为与蛋白合成及ATP利用相关的转录调控、RNA聚合酶II正调控转录、蛋白结合及ATP结合;KEGG富集分析发现占比较多与大脑发育相关的通路为调控轴突发育及重塑的轴突导向通路,与神经胶质细胞增殖相关的MAPK信号通路,与突触传递相关的神经活性配体-受体相互作用、 PI3K-Akt、钙信号及磷脂酶D信号通路,与免疫防御相关的Toll样受体通路,与低氧适应调...     

低海拔与高海拔斑头雁肺脏的比较转录组学分析

为探究高低海拔斑头雁肺脏差异表达基因和差异代谢通路的变化,基于Illumina NovaSeq 6000测序平台对斑头雁肺脏组织进行转录组测序,分析差异表达基因,利用荧光定量PCR对部分差异表达基因进行验证,进一步对差异表达基因进行GO功能富集分析和KEGG信号通路富集分析。结果:与高海拔组斑头雁肺脏相比,筛选出差异表达基因84个,其中上调差异表达基因51个,下调差异表达基因33个。经GO功能注释后富集到21个显著性GO功能,"生物过程"显著富集到11个条目,涉及代谢过程以及内源性刺激的应答等。"细胞组成"显著富集到1个条目,涉及到细胞外区域。"分子功能"显著富集到9个条目,涉及到酶活性、激素、转录因子、钠离子跨膜转运等。注释到KEGG的显著性富集通路有3条,涉及到糖胺聚糖降解通路、丝裂原活化蛋白激酶信号通路以及丙氨酸、天冬氨酸和谷氨酸代谢通路。本研究结果丰富了斑头雁的基因资源,为该物种海拔适应相关基因的功能研究奠定了基础。     

禽呼肠孤病毒强毒株和弱毒疫苗株感染鸡肝癌细胞的转录组学分析

本研究旨在研究禽呼肠孤病毒(avian reovirus,ARV)强毒株S1133和弱毒疫苗株aS1133分别感染鸡肝癌(LMH)细胞后,LMH细胞基因组的转录水平变化。将1 MOI的ARV强毒株S1133和弱毒疫苗株aS1133分别感染LMH细胞,感染后6、12和24h,收集感染细胞和阴性对照细胞,制备总RNA样品,使用Illumina HiSeq~(TM)2000测序仪进行转录组学测序。以感染组之间样品中上下调转录差异≥2且P≤0.001的基因为差异基因,进行生物信息学分析并进行荧光定量PCR验证。分析结果表明,与ARV aS1133感染组相比,S1133感染组共有4 013个差异表达基因,其中6hpi时,58个上调表达差异基因,5个下调基因;12hpi时,1 429个上调基因,354个下调基因;24hpi时,1 680个上调基因,481个下调基因。K-平均值聚类分析结果说明这些差异基因的表达模式主要呈上调表达模式,并且有10种共表达趋势。GO功能注释和分析结果表明差异基因主要具有GTP酶调节活性、信号转导活性、蛋白激酶活性等生物功能;主要参与机体应激、细胞信号转导、细胞凋亡等生物过程。Pathway富集分析结果表明差异表达基因主要参与Toll样受体信号通路、TGF-β信号通路、病原体感染应答等细胞信号通路。随机选择部分差异表达基因进行荧光定量RT-PCR验证,其结果表明与转录组学测序结果基本一致。以上数据比较分析了ARV强毒株和弱毒疫苗株感染后LMH细胞的基因表达的变化差异,探讨了宿主细胞对ARV不同毒力毒株的可能不同应答机制,有利于进一步阐明ARV致病机制。     

捻转血矛线虫阿苯达唑耐药株给药前后比较转录组学分析

为探明捻转血矛线虫阿苯达唑耐药株给药前后在转录组水平上的差异,本研究采用Illumina Hiseq4000对给药前后的捻转血矛线虫进行转录组测序,筛选得到差异表达基因并通过GO和KEGG数据库对其进行功能注释及富集性分析,并利用荧光定量PCR验证部分差异表达基因。结果显示,共筛选获得851个显著差异表达基因,其中包括584个上调基因和267个下调基因。通过GO功能富集分析显示,有458、418和367个基因分别注释到生物学过程、分子功能和细胞组分三大类;对差异表达基因进行KEGG富集分析显示,有173个差异表达基因参与到75条KEGG通路中,显著富集在核糖体合成、细胞凋亡、过氧化物酶增殖体激活受体(PPAR)等信号通路。本试验初步筛选了阿苯达唑耐药株给药前后的差异表达基因,为深入探索捻转血矛线虫耐药性分子机制、筛选对耐药性检测的分子标记及耐药性早期鉴别诊断方法的建立提供重要数据。     

齐兴肉兔转录组测序及肌肉生长发育相关基因的筛选

为了研究不同发育阶段齐兴肉兔肌肉生长发育相关的差异表达基因,试验选取胚胎期18日龄、26日龄的齐兴肉兔骨骼肌进行转录组测序,与兔基因组数据库进行比对,筛选差异表达基因并进行GO分析,进一步筛选与肌肉生长发育相关的差异表达基因,并通过qRT-PCR对筛选基因进行验证。结果表明:胚胎期26日龄共获得5 320个差异表达基因,其中上调基因2 607个、下调基因2 713个;KEGG功能注释发现差异基因显著富集到癌症通路、细胞骨架调控通路、pi3k-akt信号通路和RAP1信号通路;筛选获得的6个基因与肌肉生长发育关系较为密切,其中上调基因分别为TPM1、CASQ1、MYOZ1、CKM、MYL3基因,下调基因为IGF2BP1基因。     

基于RNA-Seq分析食欲素A对绵羊卵巢黄体化颗粒细胞基因表达的影响

食欲素A(orexin A)作为一种重要的神经肽类激素,通过与G蛋白耦联受体结合影响生殖功能。为了研究orexin A对绵羊黄体化颗粒细胞基因表达的影响及其信号网络,培养绵羊卵巢黄体化颗粒细胞进行鉴定,提取黄体化颗粒细胞组和添加Orexin A的黄体化颗粒细胞组的RNA,采用lllumina测序技术进行测序,利用生物信息学方法对转录组数据进行分析,在测序数据质量控制的基础上,对差异表达的基因进行筛选、功能注释和富集分析。结果表明,黄体化颗粒细胞组的孕酮分泌量显著高于未黄体化颗粒细胞组;转录组测序共获得了50个差异基因,其中上调表达20个,下调30个;发现多个与细胞周期调控、繁殖及代谢相关的高表达基因,如PRRT2、ID1、RRM2、ATP6、ATP8、KIRREL、SOX4、TBX3等基因,结合GO分析和KEGG通路分析表明,差异基因主要参与代谢途径、氧化磷酸化、癌症信号通路、GnRH信号通路、cGMP-PKG信号通路、Rap1信号通路、Wnt信号通路以及Ca离子信号通路等。本试验为进一步阐明orexin A影响绵羊的生殖调控提供了依据。     

猪源大肠杆菌O157∶H7毒力差异株的转录组测序分析

为了解大肠杆菌O157∶H7毒力差异株转录组差异,丰富O157∶H7转录组数据信息,本研究采用Illumina HiSeq^TM 2000平台对两株大肠杆菌O157∶H7毒力差异株进行高通量测序,测序数据采用测序评估、基因功能注释等生物信息学方法进行分析。结果发现,经过测序,两个菌株分别获得3 113 118和2 944 912条reads,比对到参考基因组上的reads分别占总reads的83.76%和78.97%。以中等毒力株为参考,在高毒力株中共获得941个差异表达基因,其中上调基因637个,下调基因304个。GO功能注释分析表明,差异表达基因主要与催化活性功能、黏附、转运活性、受体活性、酶调节活性、定位、生化调节、运动等诸多生理生化过程相关;KEGG富集分析发现共有425个基因注释到160个代谢通路中,其中新陈代谢、核糖体、鞭毛合成、嘧啶代谢、糖类代谢、细菌趋化等通路显著富集。此次通过大肠杆菌O157∶H7毒力差异株转录组研究对差异表达基因涉及的信号调控及可能的功能基因进行了探索,丰富了转录组信息,为进一步开展大肠杆菌O157∶H7毒力相关基因的研究及分子调控机制奠定了基础。     

基于转录组学技术分析绿原酸对鸡源大肠杆菌抑菌及耐药消除的作用机制

为探究绿原酸对鸡源抗生素耐药大肠杆菌抑菌及耐药消除的作用机制，以影印法分离绿原酸浓度为1.25 mg/mL（亚抑菌浓度，1/2 MIC）的LB肉汤培养的第3代禽源大肠杆菌耐药逆转菌株，以微板法测定耐药逆转菌株的左氧氟沙星最小抑菌浓度（MIC），并进一步通过转录组测序方法分析绿原酸消除大肠杆菌耐药性的分子机制。结果显示：耐药逆转菌株对左氧氟沙星的MIC由16 μg/mL降低至4~8 μg/mL，说明1.25 mg/mL的绿原酸可在一定程度上消除耐药大肠杆菌的左氧氟沙星耐药性。为进一步解析绿原酸对鸡源大肠杆菌耐药消除作用的分子机制，通过转录组测序方法对耐药性消除前后鸡源大肠杆菌基因表达水平进行对比分析发现：绿原酸作用后共有12个基因的表达量发生显著下调，通过荧光定量PCR验证结果基本一致。经过GO功能富集分析和KEGG代谢通路富集分析发现，差异表达基因在功能上集中于膜组成成分和DNA重组，在代谢通路富集中差异基因均与代谢相关，差异基因中bhsA、cysP为Coli-HB染色质转录RNA（GenBank登录号：CP020933），oqxA、oqxR、emaA、pinR、xerD、folA、repA、repC、adhP、IS26为Coli-HB质粒1转录RNA（GenBank登录号：CP020934）。从差异表达基因的分布可以看出，绿原酸可以抑制细菌DNA重组、使细菌抗逆性减弱、抑制耐药基因活性等，在一定程度上消除大肠杆菌的左氧氟沙星耐药性，起到抑菌及耐药消除作用。     

基于转录组测序挖掘广西麻鸡生长发育相关基因

【目的】通过对1和60日龄广西麻鸡腿肌进行转录组测序(RNA-Seq),筛选出广西麻鸡生长发育关键基因。【方法】选择1和60日龄健康母鸡各3只,分别采集腿肌组织,利用Illumina Novaseq^TM 6000平台进行mRNA转录组测序,利用DESeq2软件进行差异表达基因的分析,并对差异表达基因进行GO功能和KEGG通路富集分析,进一步筛选出生长发育相关基因,并用实时荧光定量PCR验证转录组测序结果。【结果】转录组测序结果显示,6个样本分别获得34 296 198～41 647 642条clean reads。与1日龄腿肌相比,60日龄腿肌共获得2 304个差异表达基因,其中998个基因上调,1 306个基因下调。GO功能富集分析发现,共获得富集条目572条,其中生物过程381条,细胞组分70条,分子功能121条。KEGG通路富集分析发现,共获得34个显著富集的信号通路,其中心肌收缩、MAPK信号通路、紧密连接、肌动蛋白细胞骨架信号通路、心肌细胞中的肾上腺素能信号传导等与生长发育相关。通过GO功能和KEGG通路富集分析获得5个生长发育相关基因,分别为肌球蛋白重...     

MSTN基因编辑牛肌肉转录组测序及生物信息学分析

为探究肌生长抑制素（myostatin,MSTN）基因对牛肌肉发育的具体调控机制,本研究选取同一牛场健康鲁西黄牛10头,其中通过转基因技术得到的基因编辑牛MSTN^-/-和同种非转基因野生型牛各5头。分别采集两组牛腿臀肌肉样品,利用IIlumina HiSeq高通量测序技术进行转录组测序分析,通过生物信息学方法比较两组样本间的差异表达基因,并进行GO和KEGG富集分析,最后利用实时荧光定量PCR验证转录组测序数据。结果显示,基因编辑型牛和野生型牛之间共检测到18 071个基因。在log₂|FoldChange|≥ 1.48条件下,筛选出406个差异表达基因,其中347个显著上调,59个显著下调。GO功能富集分析显示,MSTN基因编辑后显著影响915个功能类别（P<0.05）,差异基因主要参与结合、生物系统调节、免疫系统等相关功能。KEGG通路富集分析结果共涉及211个通路,差异基因主要富集在细胞黏附分子、趋化因子信号通路、细胞因子互作等信号通路上,进一步从中筛选出可能参与细胞生长、肌肉发育的差异基因（CD14、KIT、CSF1R、FBP1、DUSP4、ULBP21、PRKCB、SPN、CHAD、SRC）。实时荧光定量PCR检测结果显示,所选差异基因表达水平与转录组表达水平一致,证明测序结果的可靠性。本研究结果表明,MSTN基因发挥作用后可以介导多个下游基因表达,从而影响相关信号通路及生物学过程;同时,所筛选出的差异表达基因可作为进一步研究骨骼肌调控机制的候选靶标。     

干扰lnc23后牛骨骼肌卫星细胞转录组测序的生物信息学分析

试验旨在分析lnc23在调控牛骨骼肌卫星细胞分化过程中转录组水平的变化,筛选出可能参与调节骨骼肌卫星细胞生长的调控通路及相关mRNAs。采用第二代测序技术(next-generation sequencing,NGS)基于Illumia hiSeq测序平台对成功构建的lnc23抑制模型及相应对照组的牛骨骼肌卫星细胞进行转录组测序,将数据整理、过滤及评估后,通过生物信息学方法分析样品间基因转录差异表达,对这些差异基因进行GO和KEGG富集分析,并应用实时荧光定量PCR技术对转录组数据结果进行验证。结果显示,转录组测序共鉴定到19 358个基因,筛选到1 297个差异表达基因,其中上调856个,下调441个。差异基因的GO功能分析共包括细胞组分、分子功能和生物学过程3大类222个分支,其中有大量的差异基因与催化活性、分子功能调节、信号转导、生物黏附、新陈代谢相关。KEGG分析显示,差异基因共涉及190个通路,显著富集在PI3K-Akt、P53、TNF、RIG-Ⅰ样受体、神经活性配体受体互作、细胞因子互作等信号通路上,进一步从中筛选出可能参与细胞生长、肌肉发育过程的差异基因,如CCNA2、RRM2、IL-6、MYL2等。实时荧光定量PCR结果显示,所选的11个差异基因中有9个与转录组测序结果变化一致,表明了测序结果的可靠性。本试验完成了干扰lnc23及其对照后牛骨骼肌卫星细胞的转录组测序分析,获取差异基因的功能注释信息,初步揭示lnc23调控牛骨骼肌卫星细胞分化的潜在基因和通路,为深入探究lnc23调控牛骨骼肌卫星细胞的分化机制打下良好的基础。     

LPS诱导的奶牛蹄真皮细胞基因表达谱变化分析

通过转录组学测序技术,分析LPS诱导的奶牛蹄真皮细胞基因表达谱变化,筛选差异基因并进行生物信息学分析,为从细胞和分子水平上研究奶牛蹄叶炎提供必要的依据。结果显示,10 mg/L LPS作用24 h是诱导细胞模型的最佳条件。转录组测序结果显示,LPS作用前后,共有2 194个基因的表达水平发生变化,其中上调基因1 138个,下调基因1 056个。经过GO功能富集分析和KEGG代谢通路富集分析,发现差异基因主要富集于TNF、PI3K-Akt和MAPK信号通路等与免疫反应、炎症反应相关的通路中,推测LPS主要通过激活蹄真皮细胞相关炎性通路,诱导炎症介质的产生,从而导致奶牛蹄真皮细胞炎性损伤。     

基于转录组测序的油酸诱导延边牛骨骼肌卫星细胞成脂分化差异表达基因的分析

本研究旨在研究油酸(OA)引起延边牛骨骼肌卫星细胞成脂分化途径的分子调控机制。对从12日龄延边牛中分离得到的骨骼肌卫星细胞进行体外培养,加入不同浓度OA的诱导分化培养液,进行分化处理,试验设1个空白对照组[CON组,5%马血清(HS)],3个不同浓度OA诱导组:OAL组(5%HS+50μmol/L OA)、OAM组(5%HS+100μmol/L OA)、OAH组(5%HS+200μmol/L OA),诱导96 h后,对延边牛骨骼肌卫星细胞进行了转录组测序及分析。结果表明:共得到13 951条单基因,对4个组差异表达基因进行比较发现,与CON组相比,OAL组有3 412个差异表达基因,OAM组有1 045个差异表达基因,OAH组有1 428个差异表达基因,各组之间共有278个差异表达基因。GO富集分析显示,差异基因参与了多种生物过程,包括代谢过程、细胞过程、生物调节过程等;在分子功能类别上,大多数的基因功能与结合活性、转录调节活性、催化活性等相关;在细胞组分范畴内,大部分的基因被富集到细胞、细胞器、膜等。KEGG富集分析表明,差异表达基因主要富集通路为单磷酸腺苷激活的蛋白激酶(AMPK)信号通路、过氧化物酶体增殖物激活受体(PPAR)信号通路、脂肪酸代谢通路、脂肪酸降解通路等。不同浓度OA差异基因荧光定量PCR结果显示,所选的5个差异基因与转录组测序结果基本一致,表明了测序结果的可靠性。本试验完成了OA诱导延边牛骨骼肌卫星细胞成脂分化的转录组测序分析,获取差异基因的功能注释信息,初步揭示了OA诱导牛骨骼肌卫星细胞成脂分化的潜在基因和通路。     

牦牛发情期卵巢比较转录组学研究

旨在进一步了解牦牛发情期卵巢的分子机制,解析牦牛繁殖的特殊性。本研究应用RNA-seq技术对牦牛和平原黄牛发情期卵巢进行转录组高通量测序及全基因组差异表达模式比对分析。通过比较分析牦牛和黄牛卵巢转录组数据,共筛选出1 307个差异表达基因,其中661个基因表达量上调和646个基因表达量下调。进一步功能分析表明,这些差异基因涉及多种GO分类及KEGG通路。其中GO分类注释显示,差异基因与细胞粘附、激素调控等生物学过程存在密切关联,同时钙离子结合、阳离子跨膜转运等分子事件表现活跃。KEGG通路分析显示,补体和凝血级联通路的富集水平最高,其次为细胞色素P450相关通路。昼夜节律等一些新型通路,尽管与生殖功能没有明显关联,但也表现出显著富集。本研究首次对比牦牛和黄牛发情期卵巢转录组数据,筛选并分析相关差异基因。该研究结果为进一步阐述牦牛卵巢的基本分子机理提供基础,同时也为全面理解牦牛繁殖特异性提供新的思路。     

sRNA伴侣蛋白Hfq敲除条件下沙门菌的转录组分析

为探明沙门菌(Salmonella)在sRNA伴侣蛋白Hfq敲除条件下转录组变化情况。本试验采用HiSeq测序平台对沙门菌LT2菌株及其Hfq敲除株进行高通量测序,通过DESeq2差异分析方法筛选敲除sRNA伴侣蛋白Hfq条件下的差异表达基因,对其进行生物信息学GO功能显著性富集分析和Pathway显著性富集分析。结果显示,对照组和试验组分别得到12 753 534条、8 254 308条Clean Reads。通过DESeq2差异分析获得差异基因1 055个,其中表达上调基因516个,表达下调基因539个。GO功能显著性富集分析表明显著差异表达基因主要富集在钴胺素代谢过程、钴胺素生物合成过程、碳水化合物代谢过程等生物过程中。Pathway显著性富集分析表明显著差异表达基因富集到细菌趋化性、丙酸代谢和碳代谢等11个信号通路中。结果表明,本试验应用RNA-seq技术丰富了sRNA伴侣蛋白Hfq调控基因数量,筛选出20个受伴侣蛋白Hfq调控的显著差异表达基因,其中4个基因功能及编码蛋白未知,注释了差异表达基因的生物学功能及信号通路,推断Hfq在细菌对数期主要通过影响营养提供和趋化性等途径,继而影响细菌的正常生理活动,为Hfq协同sRNA的调控机理研究及sRNA靶基因的筛选奠定了基础。