持续感染口蹄疫病毒VP1蛋白对细胞转录组的影响

【目的】通过转录组测序（RNA-Seq）方法分析持续感染口蹄疫病毒结构蛋白VP1诱导的PK-15细胞基因表达谱的变化,探索VP1蛋白对宿主细胞信号通路的影响。为后续开展VP1影响FMDV持续感染相关机制研究提供基础数据和理论支持。【方法】PK-15细胞分别转染pCAGGS-HA（对照载体）、pCAGGS-HA-Pi-VP1质粒,转染36 h后,Western blot检测VP1蛋白表达,同时提取两组细胞总RNA进行转录组测序（每组3个重复,2组6个样本）。将测序所得的序列进行过滤比对,获得差异表达基因后,对差异表达基因进行了生物信息学分析。采用qRT-PCR验证关键差异表达基因的表达。【结果】RNA-Seq共鉴定出3 376个差异表达基因（DEGs）|log2(Fold Change)|>0和padj<0.05,其中上调基因1 744个,下调基因1 632个。许多参与免疫反应和炎症的效应基因差异表达。GO和KEGG富集分析显示,差异基因GO条目绝大部分与核糖体有关,而KEGG富集中涉及了与天然免疫相关的信号通路,如TNF信号通路、Toll受体信号通路。【结论】pCAGGS-...     

水稻与稻瘟病菌不同小种互作的基因差异表达谱分析

【目的】分析水稻接种稻瘟病菌（Magnaporthe grisea）48 h后亲和互作与非亲和互作反应的基因表达谱,探索水稻对不同稻瘟病菌抗性差异的分子机理。【方法】运用Affymetrix表达谱芯片分析差异表达mRNA;通过分子注释系统平台（MAS 3.0）对差异表达基因进行了基因注释及GO分析;应用实时定量PCR 对部分差异表达基因进行验证。【结果】在49 824个转录本中共检测到大约24 000个转录本,Fold change大于2.5的基因共1 028个,其中非亲和互作上调基因460个,下调基因568个。所验证的4个基因的荧光定量PCR结果与芯片结果基本一致。经GO分析差异基因对应蛋白的功能主要涉及信号转导、酶的调节、转录及分子转运等。【结论】水稻与稻瘟病菌非亲和与亲和互作的基因表达谱存在较大差异,基因芯片筛选到的差异表达基因通过GO注释明确了差异基因的分子功能及信号通路,有利于进一步了解植物抗病机制,并可能为稻瘟病防治措施提供新的途径。     

盐胁迫下转 DcCIPK24 拟南芥的转录组比较分析

【 目 的】 探 究 铁 皮 石 斛（Dendrobium catenatum Lindl.）DcCIPK24 的 下 游 响 应 基 因， 为 研 究DcCIPK24 提高耐盐性的分子调控途径提供指导。【方法】利用 Illumina NovaSeq 6000 测序平台对 200 mmol/L NaCl 处理和对照组的转 DcCIPK24 拟南芥、野生型拟南芥进行测序，分析差异表达基因的富集通路，并对通路中的差异基因进行荧光定量 PCR 验证。【结果】盐胁迫下，从 2 个转基因拟南芥株系与野生型拟南芥的比较组合中共鉴定出 96 个差异表达基因；GO 注释分析发现 96 个差异基因主要被富集在分子功能相关通路；KEGG 富集分析发现差异基因主要被注释在氨基糖和核苷酸糖代谢、植物 MAPK 信号通路和甘油酯代谢通路中，选取上调表达的差异基因 AtGPAT5、AtSIRK、AtMEE25、AtUGE5 和下调表达基因 AtAMT1-2、AtATTI3、AtCYP71A25、AtLHCB4.3 进行实时荧光定量 PCR 验证，结果表明盐胁迫下 8 个差异基因表达趋势与转录组数据基本一致。【结论】盐胁迫下，DcCIPK24 主要通过氨基糖和核苷酸糖代谢途径、植物 MAPK 信号途径和甘油酯代谢通路响应耐盐。     

基于转录组测序对紫花苜蓿细胞质雄性不育系相关代谢通路的鉴定

【目的】研究紫花苜蓿(Medicago sative L.)细胞质雄性不育(CMS)的机理,为紫花苜蓿利用不育系的育种工作提供理论基础。【方法】本研究以紫花苜蓿细胞质雄性不育系MSGN1A及其保持系MSGN1B为材料,进行转录组测序(RNA-Seq)以及基因功能注释,筛选与不育系机理相关的差异表达基因及代谢通路。【结果】Illumina测序共得到紫花苜蓿95 679条功能基因,经过差异表达分析筛选出187个显著差异表达的基因(DEGs),其中包括18个上调基因,169个下调基因。采用GO、KEGG、COG注释库分别对187个DEGs进行功能注释,其主要涉及催化活性、代谢活动、信号传导机制、合成、膜功能、碳水化合物运输和代谢、能量代谢等相关通路。【结论】筛选出与紫花苜蓿细胞质雄性不育性相关的代谢通路,初步探讨了紫花苜蓿细胞质雄性不育系的分子机理。     

卵形鲳鲹卵巢不同发育时期的转录组测序分析

【目的】鉴定筛选出与卵形鲳鲹卵巢发育相关的候选基因及信号通路,为揭示其卵巢性成熟过程的分子机制打下基础。【方法】挑选卵巢发育处于?期和Ш期的雌性卵形鲳鲹,分别构建卵形鲳鲹卵巢?期和Ш期的cDNA文库,采用Illumina HiSeqTM 2500进行转录组测序,经过滤、质量控制及拼接组装后获得的Unigenes在七大数据库（Nr、Nt、Pfam、KOG/COG、Swiss-Prot、KEGG和GO）中进行比对;通过FPKM及DEGseq筛选出差异表达基因,以GOseq和KOBAS对差异表达基因分别进行功能注释及信号通路富集分析,并采用MISA和GATK3进行SSR鉴定及SNP分析。【结果】卵形鲳鲹卵巢组织转录组测序获得的325156432条Raw reads,经过滤筛选得到317206752条Clean reads,拼接组装后得到59554条Unigenes;69.65%的Unigenes在Nr、Nt、Pfam、KOG/COG、Swiss-Prot、KEGG和GO等七大数据库中注释成功,其中有24599条Unigenes被注释到GO数据库,15997条Unigenes被注释到KEGG数据库。在卵形鲳鲹卵巢组织的2个发育时期共鉴定获得56115个基因,经差异表达分析后获得17737个差异基因,其中8169个基因在卵巢Ш期上调表达、9568个基因在卵巢Ш期下调表达。GO功能注释分析发现,卵形鲳鲹卵巢差异表达基因主要注释在细胞过程、氮化合物代谢过程、初级代谢过程、核、核部分、离子结合及水解酶活性等条目上;而KEGG信号通路富集分析结果显示,17737个差异表达基因显著富集在318条代谢途径上,其中前20条KEGG信号通路包括2-氧代羧酸代谢、PI3K-Akt信号通路、甲状腺激素信号通路、磷脂酶D信号通路、Fc εRI信号通路和细胞周期等。卵形鲳鲹卵巢转录组（59554条Unigenes）中共存在30133个SSRs和82490个SNPs。【结论】GnRHR、FSHR、FSHβ、CYP11A、SIRT3和PEG3等差异表达基因及PI3K-Akt信号通路和VEGF信号通路等与卵形鲳鲹卵巢的发育密切相关,共同调节卵巢的发育与成熟,在卵巢性成熟过程中发挥重要作用。     

2种不同生长规格墨瑞鳕肌肉组织转录组测序分析

【目的】掌握不同生长规格墨瑞鳕个体间差异表达基因的表达特点,为其功能相关基因深度挖掘及分子遗传育种提供科学依据。【方法】挑选同一养殖条件下极大个体和极小个体的墨瑞鳕,构建肌肉组织cDNA文库后,采用Illumina HiSeqTM 4000测序平台对存在生长差异的墨瑞鳕肌肉组织进行转录组测序分析,获得的Unigenes在Nr、Nt、Pfam、KOG/COG、Swiss-Prot、KEGG和GO等数据库中进行比对;通过FPKM及DEGseq筛选出差异表达基因,以GOseq和KOBAS对差异表达基因分别进行GO功能注释及KEGG信号通路富集分析,并采用MISA进行SSR鉴定分析。【结果】从墨瑞鳕肌肉组织中共测序获得39749条Unigenes,其长度范围在301~55230 bp,平均长度为1705 bp。注释到Nt、Nr、Swiss-Prot、Pfam数据库的Unigenes分别有27046、20824、18268和17772条,在7个数据库中均得到注释的Unigenes共计6742条,占Unigenes总数的16.96%。根据差异表达基因筛选条件P<0.05且|log₂ Fold Change|>1,共筛选出722个差异表达基因,其中上调基因308个、下调基因414个。差异表达基因GO功能注释分析结果表明,注释基因数目较多的GO功能条目包括细胞过程、代谢过程、膜、细胞器及结合等;KEGG信号通路富集分析发现,差异表达基因被成功富集到234条信号通路上,主要涉及磷脂酰肌醇3激酶/蛋白激酶信号通路、MAPK信号通路、胰岛素信号通路及FoxO信号通路等。在39749条Unigenes中鉴定筛选出22120个SSRs,占Unigene总数的55.65%,SSR的平均间距为3063 bp。【结论】基于转录组测序分析获得的墨瑞鳕肌肉组织差异表达基因以发挥结合、细胞过程及代谢过程等功能为主,且主要富集在PI3K-Akt信号通路、核糖体信号通路、FoxO信号通路及细胞凋亡等能量代谢相关通路上,通过共同协调而对墨瑞鳕的生长发育起调控作用。     

基于高通量测序的山羊卵巢基质、卵泡转录组分析

【目的】比较山羊Capra hircus卵巢基质、大卵泡和小卵泡间的m RNA表达图谱,为探究卵泡发育机制提供一定的研究基础。【方法】采用高通量测序技术对发情期川中黑山羊的卵巢基质、大卵泡和小卵泡进行转录组测序,利用生物信息学方法检测mRNA表达谱,筛选差异基因,对其进行GO和KEGG通路分析,最后随机抽取5个差异基因进行荧光定量PCR验证。【结果】卵巢基质vs大卵泡、卵巢基质vs小卵泡和小卵泡vs大卵泡的差异表达基因分别为524、180和403个。筛选出INHA、TNFRSF19等15个与卵泡发育相关的基因,其主要富集在内质网蛋白质加工、类固醇生物合成、卵母细胞减数分裂等信号通路。通过q RT-PCR验证,定量结果与测序结果基本一致。【结论】川中黑山羊卵巢基质、大卵泡和小卵泡的基因表达模式不同,小卵泡与卵巢基质的基因表达模式更相近。     

西方蜜蜂工蜂不同虫态发育的转录组学分析

【目的】基于转录组学对西方蜜蜂工蜂不同虫态间的差异表达基因（DEGs）进行筛选和功能注释分析,揭示与工蜂生长发育相关的信号通路,为深入解析工蜂生长发育的分子调控机理提供基础数据。【方法】以西方蜜蜂工蜂的3日龄幼虫、1日龄蛹和1日龄羽化工蜂3个虫态为研究对象,利用llumina NovaSeq 6000平台进行转录组测序,采用DESeq2筛选不同虫态样品间的表达差异基因,然后分别进行GO功能注释分析及KEGG信号通路富集分析,并通过实时荧光定量PCR进行验证。【结果】经转录组测序,在西方蜜蜂工蜂3日龄幼虫与1日龄蛹间筛选出4823个差异表达基因（51.86%上调,48.14%下调）,在1日龄蛹与1日龄羽化工蜂间筛选出3295个差异表达基因（57.51%上调,42.49%下调）,在3日龄幼虫与1日龄羽化工蜂间筛选出5267个差异表达基因（52.95%上调,47.05%下调）。GO功能注释分析结果显示,3日龄幼虫与1日龄蛹间的差异表达基因注释到43个GO功能条目,1日龄蛹与1日龄羽化工蜂间的差异表达基因注释到45个GO功能条目,3日龄幼虫与1日龄羽化工蜂间的差异表达基因注释到44个GO功能条目,主要涉及细胞过程、细胞部分及结合等。KEGG信号通路富集分析发现,3日龄幼虫与1日龄蛹间有2905个差异表达基因富集到332条KEGG信号通路上,其中17条KEGG信号通路呈显著富集,涉及核糖体、氧化磷酸化和昆虫激素生物合成等;1日龄蛹与1日龄羽化工蜂间有1644个差异表达基因富集到331条KEGG信号通路上,其中45条KEGG信号通路呈显著富集,涉及氧化磷酸化、生热作用和胰岛素分泌等;3日龄幼虫与1日龄羽化工蜂间有2958个差异表达基因富集到337条KEGG信号通路上,其中14条KEGG信号通路呈显著富集,涉及核糖体、蛋白酶体和胰岛素分泌等。6个随机挑选差异表达基因的实时荧光定量PCR检测结果与转录组测序结果相符,即转录组测序结果可靠。【结论】昆虫激素生物合成通路相关差异表达基因调控与西方蜜蜂工蜂各虫态JH滴度变化规律一致,氧化磷酸化信号通路则与各虫态的营养摄入和活动行为相关,而胰岛素分泌通路涉及各虫态的营养调控、脂肪体合成及细胞凋亡。可见,昆虫激素生物合成、胰岛素分泌和氧化磷酸化3种信号通路在西方蜜蜂工蜂幼虫、蛹和成虫的发育调控中发挥着重要作用。     

FST基因的过表达对猪成纤维细胞内基因及生物学通路的影响

应用Illumina HiSeq测序平台,对过表达FST基因的细胞和正常细胞进行转录组测序,并对筛选的差异表达基因进行基因功能注释及信号通路分析。结果表明：过表达FST后共有286个基因的转录发生了变化,其中176个上调,110个下调。这些差异基因共得到44个GO功能注释,其中单一生物过程、代谢过程、细胞过程;细胞器、细胞区域、细胞和黏合terms下聚集的差异基因数量最多,均>100个。共发现9条显著富集的Pathway,包括与细胞增殖相关的细胞周期、DNA复制、碱基切除修复、p53信号通路等,与卵母细胞发育相关的卵母细胞减数分裂和孕激素介导的卵母细胞成熟通路。说明在成纤维细胞中过表达FST后,最显著的特征是与细胞增殖分化相关的基因和生物学通路被富集。     

半番鸭与番鸭精巢组织差异表达转录组测序分析

【目的】研究半番鸭与番鸭精巢组织转录组差异表达基因,为进一步阐明半番鸭不育的遗传机制奠定理论基础。【方法】利用转录组测序方法对半番鸭和番鸭的精巢组织进行研究,筛选其差异表达基因,并对其功能进行注释和荧光定量PCR(quantitative real-time PCR,QRT-PCR)验证。【结果】测序共获得43.84Gb Clean Data,组装后共获得193 535条Unigene。DESeq分析发现3 597个基因在两个鸭品种间差异表达,其中上调基因1 194个和下调基因2 403个,包括与生殖功能相关的基因,如成纤维细胞生长因子(fibroblast growth factor,FGF)、蛋白激酶A(protein kinase A,PKA)、丝裂原活化蛋白激酶7(mitogen-activated protein kinase 7-like,partial,BMK)、生长因子受体结合蛋白2(growth factor receptor-bound protein 2,GRB2)、c-Jun氨基末端激酶(c-Jun N-terminal kinase,JNK)和肿瘤坏死因子受体超家族成员6(tumor necrosis factor receptor superfamily member 6,FAS)等。GO(gene Ontology)分析发现382个差异基因获得功能注释,其中97个基因涉及发育繁殖生物学过程。KEGG(kyoto Encyclopedia of Genes and Genomes)通路分析表明差异表达基因共富集到50条信号通路中,其中17个通路显著富集,包括丝裂原活化蛋白激酶信号转导通路(mitogen-activated protein kinase signaling pathway,MAPK)、甘油酯代谢(glycerolipid metabolism)以及钙信号途径(calcium signaling pathway)信号通路等,与生殖功能密切相关的有促性腺激素释放激素信号通路(gonadotropin releasing hormone(Gn RH)signaling pathway)和MAPK信号转导通路。经QRT-PCR验证,差异基因表达变化模式与转录组测序结果一致,测序结果可靠。【结论】在转录组水平上筛选出半番鸭和番鸭精巢组织差异表达基因,揭示了Gn RH和MAPK信号通路在鸭的生殖活动中发挥了重要作用,为进一步探索半番鸭生殖系统的分化机理提供可靠的参考依据。     

光周期影响罗非鱼脑组织转录组基因表达分析

【目的】发掘尼罗罗非鱼响应光周期变化的重要功能基因,为研究光周期变化对罗非鱼产生影响的分子机制提供参考。【方法】采用新一代高通量测序RNA-seq技术分析24L∶0D、12L∶12D和0L∶24D 3个光周期条件下尼罗罗非鱼脑组织基因表达变化情况。【结果】转录组测序结果显示:3个组别分别产生55 524 504、61 410 946和58 855 762个Clean reads。12L∶12D光周期与0L∶24D光周期文库比较发现279个差异表达基因,12L∶12D光周期与24L∶0D光周期文库比较发现304个差异表达基因,24L∶0D光周期与0L∶24D光周期文库比较发现383个差异表达基因。GO分类分析表明,差异表达基因属于生物学过程、细胞定位和分子功能的42个类别。KEGG通路显著性富集分析差异表达基因共涉及143条代谢途径,包括单纯疱疹感染、ECM-受体相互作用、细胞因子-细胞因子与受体相互作用、肠道IgA生成免疫网络、细胞粘附分子(CAMs)、Wnt信号等通路。【结论】光周期变化影响尼罗罗非鱼脑组织基因表达水平,表达差异基因主要富集在单纯疱疹感染、细胞因子-细胞因子与受体相互作用、Wnt信号通路等信号通路。     

不同糖源对葡萄试管苗蛋白激酶相关基因表达的影响

【目的】探究不同外源糖对葡萄试管苗生长发育及蛋白激酶基因转录调控的影响,应用转录组测序挖掘蛋白质磷酸化过程中的基因,为葡萄蛋白激酶相关基因功能的验证奠定一定基础。【方法】在基本培养基中分别添加2%的蔗糖、葡萄糖和果糖,以无糖为对照,分别命名为S20、G20、F20和CK,经过37 d培养后,测定不同处理的地上和地下鲜重,并采用Illumina HiSeq ^TM 2000对各处理叶片进行转录组测序,通过综合生物信息学分析（参考基因组比对、差异基因（DEGs）筛选、COG（Cluster of Orthologous Groups of proteins）注释、GO（Gene Ontology）注释等）筛选出蛋白激酶相关基因,通过qRT-PCR分析该蛋白激酶相关基因的表达特性。【结果】F20、G20和S20处理下的葡萄（‘红地球’）试管苗与CK相比,地上鲜重具有明显差异,且F20最高,而G20地下鲜重最高。SNP统计发现,转换是主要的变异类型,颠换次之,且发生在基因间区的SNP数量最多,其次是下游;剪接位点供体和同义终止发生的基因数量最少且相等。4个样品中共获得了2 633个差异基因,3个处理与CK相比,共有差异基因180个且被聚类为3组,第一组中127个基因仅在CK中高表达,第二组19个基因仅在G20下高表达,而第三组34个基因在3个处理下表达模式不尽相同。这些共有的差异基因在COG中注释到了26个基因并分在11个功能类别中,且主要注释在一般的功能类别中。在GO分类中,共有的基因分别被注释在分子功能、生物学过程和细胞组分的14、22和13个功能类别中。共筛选出7种蛋白激酶,分别为葡萄糖激酶（Glucokinase,GK）、丝裂原活化蛋白激酶（Mitogen-activated protein kinases,MAPKs）、钙调蛋白激酶（Calcineurin protein kinase,CBL）、蛋白磷酸酶2（Protein phosphatase 2,PP2A）、己糖激酶（Hexokinase,HXK）、组氨酸蛋白激酶（Histidine protein kinase,HPK）和酪氨酸激酶（Tyrosine kinase,TK）,其不同激酶的基因在不同处理中具有各自的表达模式,经qRT-PCR验证,选择的20个差异基因中有17个基因表达与转录组测序结果相一致。【结论】在葡萄试管苗培养中,果糖较葡萄糖和蔗糖相比对生长较好。测序得出180个差异基因对3种不同糖均作出响应,这些基因在COG数据库中主要富集在膜酯转运和代谢、次级代谢物和碳水化合物的合成、转运和分解;GO中大多注释在蛋白激酶和氧化还原酶的活性中;筛选出了7种蛋白激酶,这些差异基因在数量、功能分类和代谢通路上对糖的响应各不相同。     

海马齿根系响应盐胁迫的转录组分析

【目的】对盐胁迫下海马齿根系进行转录组测序分析,挖掘海马齿根系耐盐相关基因,为揭示海马齿耐盐的分子机制提供参考。【方法】利用Illumina测序技术对0 mmol/L NaCl （对照组）和400 mmol/L NaCl胁迫处理（盐胁迫处理组）下的海马齿根系进行转录组测序分析,从中筛选出差异表达基因,选取13个基因进行实时荧光定量PCR （qRTPCR）检测,以验证转录组数据的可靠性。【结果】在海马齿根系转录组中共鉴定出305145个转录本,平均长度为622 bp,其中,对照组有146177个长度>300 bp的转录本,盐胁迫处理组有72173个长度>300 bp的转录本;共有65535条Unigenes在Nr、GO、Swiss-Prot、COG和KEGG五大数据库注释成功,占Unigenes总数的52.36%。对照组和盐胁迫处理组共有65535个差异Unigenes,其中,有182个热休克蛋白基因。对照组和盐胁迫处理组间共有24042个差异表达基因,从中选取13个基因进行qRT-PCR检测,结果显示,9个基因表达上调,其余4个基因表达下调,与转录组测序结果一致。24042个差异表达基因中,共有10106个显著差异基因富集到129条代谢通路,其中富集程度排名前10的代谢途径为核糖体、次级代谢生物合成、RNA转运、内吞作用、剪接体、甘油磷脂代谢、内质网加工、吞噬、醚脂类代谢和植物-病原体相互作用,参与盐胁迫相关的硫代谢、脯氨酸积累、活性氧（ROS）代谢、与盐胁迫相关的钙信号通路和过氧化氢代谢等途径的差异基因上调。【结论】在盐胁迫下海马齿差异表达基因如小分子量热激蛋白基因、抗氧化酶相关基因及与离子交换相关基因发挥了重要调控作用。     

小麦响应禾谷镰刀菌侵染的核心基因筛选与鉴定

【目的】筛选鉴定小麦品种连麦12响应禾谷镰刀菌侵染过程中的核心基因,为小麦赤霉病防御分子网络提供理论参考。【方法】对成功感染禾谷镰刀菌的连麦12小穗进行转录组测序,通过DESeq筛选差异表达基因（DEGs）,对DEGs进行GO功能注释和KEGG信号通路富集分析,构建DEGs的蛋白质互作网络,通过MCODE和cytoHubba筛选核心基因。【结果】对小麦品种连麦12接种禾谷镰刀菌（处理组）和接种无菌去离子水（对照,CK）72 h后的DEGs进行筛选,结果共筛选出14180个DEGs,其中有10966个DEGs上调表达,有3214个DEGs下调表达。经GO功能注释分成三大类,共51个条目,其中,细胞组分包含14个条目,主要富集在细胞、细胞部件和膜3个条目;分子功能包含11个条目,主要富集在催化活性和结合2个条目;生物学过程包含26个条目,主要富集在细胞过程、代谢过程和刺激反应3个条目。KEGG信号通路富集分析结果显示,富集到代谢、次级代谢物生物合成、苯丙烷类生物合成、植物信号激素转导、植物MAPK信号通路的DEGs数量最多;富集到光合作用—天线蛋白通路,苯丙氨酸、酪氨酸和色氨酸生物合成,谷...     

黔北麻羊毛色差异主要调控基因及其通路分析

【目的】明确调控黔北麻羊毛色差异的候选基因及其通路,为揭示黔北麻羊特征毛色的形成机制提供理论依据。【方法】以1岁黔北麻羊公羊为研究对象,分别采集其肩部黑色羊毛下的皮肤组织块和腹部白色羊毛下的皮肤组织块,以RNA为模板合成双链cDNA,经末端修复等处理后进行PCR富集构建cDNA文库,采用Illumina HiSeqTM4000测序平台完成转录组测序,经过滤及质控后以DESeq2筛选出差异表达基因（DEGs）,然后分别进行GO功能注释分析和KEGG信号通路富集分析。【结果】与白色羊毛对照组（White）相比,黑色羊毛试验组（Black）存在509个显著上调表达基因、195个显著下调表达基因;在黑色羊毛试验组中的509个显著上调基因中筛选获得5个与黑色素生成相关且表达差异显著（基因表达差异倍数log2Fold Change>1、校正后的P-adjust<0.05）的基因,分别是酪氨酸酶基因（TYR）、酪氨酸相关蛋白1基因（TYRP1）、黑素皮质素受体1基因（MC1R）、角蛋白19基因（KRT19）和角蛋白79基因（KRT79）。显著差异表达基因GO功能注释分析结果表明,除了色素...     

禽致病性大肠杆菌phoP/Q基因缺失对雏鸡脾脏microRNAs功能的影响

从microRNA(miRNA)水平探究禽致病性大肠杆菌(APEC)phoP/Q基因缺失对雏鸡脾脏功能表达的影响,为APEC的防制提供参考资料。用禽致病性大肠杆菌和其phoP/Q基因的缺失株分别攻毒14日龄雏鸡,采集其脾脏组织为材料进行高通量测序获得miRNA表达谱,选择差异倍数大于8的miRNA进行Real-time PCR,对差异表达的miRNA进行靶基因预测以及GO、KEGG富集分析。测序结果获得10个差异表达的miRNA,其中1个miRNA表达量下调,9个miRNA表达量上调,Real-time PCR结果表明,miRNA变化趋势与测序结果一致,GO分析结果表明,差异表达基因主要富集在能量代谢、免疫系统、细胞生长与凋亡、糖合成与代谢等生物学过程。KEGG分析结果显示预测靶基因参与MAPK信号通路、糖代谢通路、Wnt信号通路等重要通路。本试验分析禽致病性大肠杆菌phoP/Q基因缺失后感染雏鸡脾脏miRNA的表达差异,从miRNA水平探讨了phoP/Q基因缺失对APEC致病性的影响,为APEC的防治提供参考资料。     

拟轮枝镰孢与玉米籽粒互作的差异基因筛选及代谢通路分析

【目的】拟轮枝镰孢(Fusarium verticillioides)引起的玉米穗腐病是我国玉米产区发生严重的病害之一，论文旨在了解病原菌与玉米籽粒互作过程中的基因表达差异，为揭示病原菌的致病机制和玉米抗病机制提供依据。【方法】对拟轮枝镰孢侵染玉米籽粒0、4、12和72 h的样品进行转录组测序，之后采用生物信息学分析，分别以玉米和拟轮枝镰孢基因组为参考，以|log₂FC|≥1,P-adjust<0.05为阈值筛选互作过程中玉米和拟轮枝镰孢的差异表达基因，利用GO和KEGG对其进行功能注释及富集分析。Goatools软件分析植物-病原互作、MAPK途径和植物激素信号转导通路相关差异基因的表达变化，采用实时荧光定量PCR(qRT-PCR)方法对测序差异基因进行验证。【结果】在互作4、12和72 h后拟轮枝镰孢分别有140、400和1 945个基因上调表达，有9、302和1 784个基因下调表达；玉米分别有293、692和1 426个基因上调表达，320、482和153个基因下调表达。GO和KEGG富集分析显示，侵染早期拟轮枝镰孢在细胞间隙生长，差异基因主要富集在...     

山麻鸭开产期和产蛋高峰期卵巢组织转录组分析

【目的】探明山麻鸭开产期和产蛋高峰期卵巢组织的转录组差异,丰富蛋鸭转录组数据信息。【方法】选取6份山麻鸭卵巢组织样品（开产期和产蛋高峰期各3份）,利用Illumina HiSeq^TM 2000进行高通量测序,构建山麻鸭开产期和产蛋高峰期卵巢组织转录组文库,并用测序评估、基因功能注释等生物信息学方法进行分析。【结果】经过测序获得质量不低于20的碱基比例（Q20）均高于93%;开产期获得了43 489 798条reads,4 380 847 061 bp数据量;高峰期获得了42 782 676条reads,4 307 499 083 bp数据量;分别有70.92%和72.70%的reads能比对到鸭参考基因组序列上。对开产期与产蛋高峰期转录组进行比较,发现开产期有26 808个表达基因,产蛋高峰期有27 013个表达基因,与开产期为参考,在高峰期卵巢组织中共获得1 929个差异表达基因,其中上调基因989个,下调基因940个;进一步将这些差异基因在Nr数据库进行注释,发现获得注释的基因有1 423个,其中上调基因695个,下调基因728个;COG功能注释分析发现这些差异表达基因共获得663个功能注释,涉及25个功能分类;GO 功能注释分析表明,有1 122个基因获得GO 功能注释,可以分为61个功能分类,这些分类主要涉及到分子绑定、催化活性、细胞过程、生物调节等诸多生理生化过程,其中涉及到发育繁殖生物学过程的相关注释基因有406个;KEGG分析发现共有425个基因注释到160个代谢通路中,其中GnRH信号通路、GPI-anchor生物合成、TGF-β信号通路、核糖体生物合成等通路显著富集。GnRH信号通路和TGF-β信号通路在卵巢的生理生殖活动发挥了重要作用,而GPI-anchor生物合成和核糖体生物合成,这两类代谢通路在卵巢生理生殖活动的作用尚不明确。【结论】利用高通量测序技术对山麻鸭开产期和产蛋高峰期卵巢组织的转录组进行测序和分析,揭示了山麻鸭不同生理状态下卵巢组织差异表达基因的数量,获得了差异表达基因的功能、分类和代谢通路。为丰富蛋鸭卵巢组织转录组信息,为开展蛋鸭产蛋性状相关基因的研究及分子调控机制研究奠定基础。     

藏黄牛与宣汉黄牛心脏miRNA表达谱比较

【目的】miRNA作为一类非编码RNA,广泛参与机体多种生命活动,研究旨在挖掘miRNA在藏黄牛和宣汉黄牛心脏组织中的差异miRNA,为进一步研究藏黄牛低氧适应的分子机制提供基础数据。【方法】随机选取健康的藏黄牛和宣汉黄牛各3头,迅速采集其心脏组织,使用Trizol法提取RNA,琼脂糖凝胶电泳切胶选取18—30nt的片段,连接3'端和5'端,反转录后进行扩增,凝胶电泳切胶纯化后建立藏黄牛和黄牛各3个文库,利用Illumina HiSeq4000测序平台进行高通量测序;将测序得到的序列进行过滤,比对GenBank和Rfam数据库,筛选出藏黄牛和宣汉黄牛的差异miRNA,进行功能注释及信号通路富集分析;随机选择8个miRNAs,利用实时荧光定量PCR检测其在心脏组织的表达量,以验证测序数据的准确性。【结果】藏黄牛和宣汉黄牛心脏组织的高质量读值的序列分别为17 463 446条和13 662 812条,干净读值为16 552 296条和12 055 304条,且在藏黄牛和宣汉黄牛中高质量核酸序列长度富集最多的均是21nt,分别为37.5%和32.1%;且共筛选出219个差异miRNAs,其中48个显著上调,171个显著下调;GO功能注释得到差异miRNA靶基因分子功能中显著富集的条目有22条,如,GO:0005488(结合)、GO:0005515(蛋白质结合)和GO:0043167(离子结合);细胞组分中显著富集的条目有20条,如,GO:0005623(细胞)、GO:0044464(细胞组分)和GO:0005622(细胞内);生物过程中显著富集的条目有13条,如,GO:0035556(细胞内信号转导)、GO:0032774(RNA生物合成过程)和GO:0006351(转录,DNA模板化),KEGG信号通路分析得到差异表达miRNAs靶基因显著富集到胰岛素信号通路(139个靶基因)、mTOR信号通路(38个靶基因)和HIF-1信号通路(92个靶基因)等232个信号通路中,其中有12个靶基因共同作用于这3个信号通路,说明miRNAs可能通过这3个信号通路,共同参与藏黄牛低氧适应性的调控;随机选取8个miRNAs进行荧光定量验证,其表达趋势与测序结果表达趋势基本一致,表明高通量测序数据可信度较高。【结论】得到了miRNA在藏黄牛与宣汉黄牛心脏组织中的表达谱,为揭示藏黄牛低氧适应性的调控机制研究奠定了基础。     

苹果多酚对小鼠巨噬细胞TNF信号通路的影响

【目的】通过苹果多酚(apple polyphenols,APs)处理RAW264.7小鼠巨噬细胞,以探讨苹果多酚的生物学作用.【方法】通过MTT试验确定苹果多酚的处理质量浓度,保证细胞的存活率;将RAW264.7巨噬细胞用10μg/mL的苹果多酚处理4h后,提取其总RNA,反转录后通过基因芯片技术检测,比较试验组与对照组TNF信号通路中基因表达量的差异,并进行生物信息学分析和实时荧光定量PCR验证.【结果】苹果多酚质量浓度为10μg/mL时,能够显著(P0.05)促进RAW264.7巨噬细胞的增殖能力,而质量浓度为5μg/mL时,对巨噬细胞增殖能力的影响并不显著(P0.05);基因芯片数据共筛选出差异倍数2倍以上的基因为3 149个(≥2-fold,P≤0.5),其中1 540个基因在RAW264.7巨噬细胞中表达上调,1 609个基因表达下调.基因体本轮(GO)分析结果显示:这些差异基因主要与核苷酸结合、DNA修复、细胞周期、细胞骨架与细胞核等相关,涉及的信号通路包括类固醇生物合成、真核细胞核糖体合成、TNF信号通路、抗生素生物合成、多巴胺能神经突触、烟酸和烟酰胺代谢等.但为了进一步研究TNF信号通路,利用qPCR验证该通路上(Map3k8、Mlkl、Pik3ca、Cflar、Atf2、Tnf)6个差异基因的表达,其中,Map3k8、Mlkl、Tnf 3个差异基因被显著(P0.05)下调,而Pik3ca、Cflar、Atf2 3个基因则下调不显著(P0.05).【结论】研究苹果多酚对TNF信号通路中差异表达基因的影响,对探索疾病发生发展的机制具有重要意义.