西方蜜蜂六个亚种苹果酸脱氢酶Ⅱ基因的遗传差异

研究了西方蜜蜂Apis mellifera 6亚种-浙农大1号意蜂（ZND A. m. ligustica）、东北黑蜂（A. m. ssp.）、卡尼鄂拉蜂（A. m. carnica）、喀尔巴阡蜂（A. m. carpatica）、高加索蜂（A. m. caucasica）和乌克兰蜂（A. m. acervorum）苹果酸脱氢酶Ⅱ的基因型频率、基因频率和杂合纯合度。浙农大1号意蜂、喀尔巴阡蜂和高加索蜂的纯合度较高,但浙农大1号意蜂等位基因c频率最高,喀尔巴阡蜂等位基因b频率最高, 高加索蜂等位基因a频率最高;东北黑蜂、卡尼鄂拉蜂和乌克兰蜂是高度杂合的亚种,但东北黑蜂等位基因a、b、c的频率差异较小,卡尼鄂拉蜂和乌克兰蜂主要存在a、c两个等位基因,b出现频率很小;6亚种的基因型频率、基因频率和杂合纯合度都有极显著差异。这些差异将从遗传和生化角度为西方蜜蜂6个亚种的鉴别提供依据。     

西方蜜蜂工蜂蛹响应低温胁迫的差异表达基因分析

【目的】本研究旨在从转录组水平筛选和分析西方蜜蜂Apis mellifera工蜂不同时期蛹对低温胁迫反应的差异表达基因(DEGs)。【方法】将西方蜜蜂工蜂封盖后3 d预蛹以及封盖后6和9 d蛹分别置于低温环境(20℃)和最适发育温度(35℃)中4 h,分别作为处理组和对照组,通过转录组学技术筛选低温处理组与对应的对照组之间的DEGs,并进行GO功能分类和KEGG通路分析。利用RT qPCR分别对封盖后3 d预蛹以及封盖后6和9 d蛹的8, 6和5个DEGs的表达谱进行验证。【结果】与对照组相比,封盖后3 d预蛹以及封盖后6和9 d蛹受低温胁迫后的DEGs分别有220, 50和26个;GO功能分类发现,DEGs富集数最多的条目为代谢进程、细胞进程、催化活性和结合,封盖后3 d预蛹的DEGs在生物学进程调控、细胞部分和细胞器等有较多富集。KEGG通路分析显示,处理组和对照组间西方蜜蜂各日龄DEGs在整体概述图、氨基酸代谢、信号转导、运输和分解代谢有富集。封盖后3 d预蛹以及封盖后6和9 d蛹受低温胁迫后共有的DEG 3-磷酸肌醇依赖性蛋白激酶1基因PDK1上调表达,同时富集在自噬-动物、mTOR和FoxO信号通路。低温胁迫后,封盖后3 d预蛹的胰岛素受体底物1-B基因IRS1-B和Kruppel同源物1基因Kr-h1上调表达,而核激素受体FTZ-F1基因Ftz-F1与蜕皮启动激素基因Eth显著下调表达,说明封盖后3 d预蛹响应低温细胞自噬和蜕皮受到抑制程度更大。在低温胁迫后封盖后6 d蛹中与昆虫角质层着色与免疫相关的酪氨酸羟化酶基因TyHyd下调表达;与对照组相比,在低温胁迫后封盖后9 d蛹中DEGs最少,说明在3个蛹期中,其受低温的影响较小。【结论】本研究测定了西方蜜蜂3个不同发育阶段蛹响应低温的DEGs,结果显示大部分DEGs为阶段特有,说明西方蜜蜂不同发育阶段响应低温的机制不同。一些共有DEGs以及阶段特有DEGs的功能研究和其作用机制是进一步研究蜜蜂对低温响应机制的重点内容,为探究蜜蜂蛹期响应低温胁迫的分子机制提供了基础数据。     

黄秋葵花和果荚转录组测序及类黄酮代谢差异表达分析

采用Illumina Hi seq 2500转录组高通量测序,构建黄秋葵花和果荚转录组文库,并利用测序评估、差异基因功能注释等生物信息学方法进行分析。研究结果表明：(1)分别获得黄秋葵花和果荚有效数据7.12 Gb和8.14 Gb,碱基百分比(Q30)均达到91.0%以上。(2)获得差异表达基因(DEGs)1 336个,其中上调基因319个,下调基因1 017个。(3)获得功能注释的基因有1 131个,GO将455注释基因归成41个功能小类,主要涉及代谢过程、催化活性、单一生物过程和细胞过程等过程;KOG注释了472个DEGs,涉及23个功能分类,其中与次生代谢直接相关过程O和Q类别获得111个注释结果;KEGG将372个DEGs注释到80条代谢通路中,获得F3H、F3′5′H、DFR、ANR、ANS、GT、LAR共10个关键差异基因,其中F3H、DFR在黄秋葵花中表现上调效应,F3′5′H、ANR、LAR在黄秋葵果荚中表现显著上调效应,ANS、GT则分别在花和果荚中均有上调或下调效应。(4)实时定量PCR分析显示,其中7个差异表达基因得到的相对表达量与转录组表达谱分析趋势完全一致。(5)类黄酮代谢途径分析表明,黄秋葵花通过F3H、DFR、ANS、GT途径将NAR催化为生成花青素苷;果荚则将NAR通过F3′5′H将催化为DHM,后在FLS催化下生成黄酮醇类物质等;部分NAR在F3H、 DFR催化下,生成无色飞燕草素苷元,再分别在ANS、LAR作用下,进入花青素苷元和原花青素合成途径。该研究结果丰富了黄秋葵转录组信息,为黄秋葵类黄酮物质纯化和功能利用提供参考依据。     

地果的转录组学分析

为探究地果(Ficus tikoua)的遗传变异特征,利用Illumina HiSeq-2500平台获取地果的基因表达谱。结果表明,共获得了197 362个转录本,总长度和平均长度分别为114 072 125和577 bp。使用BlastX和BlastN共注释了139 992个转录本(占总数的70.93%)。从3个品种叶片和茎的6对样品中,分别鉴定了12 397、12 340、10 373、94 431、71 830和44 465个差异表达基因,以及注释了126、129、125、134、138和137条代谢途径。这将有助于理解地果的遗传特征,以及不同组织中代谢途径的变化。     

泽兰实蝇雌雄成虫转录组及差异分析

[目的]泽兰实蝇Procecidochares utilis Stone是入侵杂草紫茎泽兰Eupatorium adenophorum Spreng的重要的专食性天敌,已成为控制紫茎泽兰的重要因子.本研究旨在获得泽兰实蝇转录组,并深入研究泽兰实蝇雌雄成虫差异表达基因.[方法]利用Illumina高通量测序技术对泽兰实蝇雌雄成虫进行转录组测序和生物信息分析.[结果]总共获得29 147条unigenes,平均长度为457 bp,同时将所得序列注释到七大数据库进行比对,共获得19 384条unigenes注释结果.分析发现11 331条差异表达基因;与雄成虫相比,雌成虫有2 640条unigenes表达量上调,8 691条unigenes表达量下调.[结论]本研究获得了泽兰实蝇转录组,为今后泽兰实蝇性别相关研究提供了序列支持.     

金耳菌丝及子实体的转录组差异分析

为探讨金耳菌丝体时期和子实体时期基因表达变化,利用高通量测序技术进行RNA-Seq分析.差异基因分析显示,两个时期共有10 992个差异表达基因,其中在子实体时期上调、下调的基因数分别为9 988和1 004个.GO功能聚类分析显示,差异表达基因主要富集在催化活性、结合、细胞过程、代谢过程、细胞以及细胞组分.KEGG富...     

基于转录组的不同火龙果品种抗性差异分析

不同的品种抗性不同,为进一步探究不同火龙果品种之间的抗性差异,为后续火龙果抗性育种提供参考,该研究利用Illumina HiSeq 2000测序平台对'普通白肉'(BR)和'厄瓜多尔黄龙'(EY)两个品种进行转录组测序分析,并参考GO Ontology、KEGG等公共数据库对差异表达基因进行功能分类与富集分析.结果表明...     

转录组分析揭示东方蜜蜂微孢子虫侵染意大利蜜蜂的分子机制

【目的】本研究旨在通过差异表达基因(differentially expressed gene, DEG)分析以及毒力因子和其他侵染相关因子分析,在转录组水平揭示东方蜜蜂微孢子虫Nosema ceranae侵染意大利蜜蜂Apis mellifera ligustica的分子机制。【方法】基于前期已获得高质量的东方蜜蜂微孢子虫纯化孢子(NcCK)及侵染意大利蜜蜂工蜂7和10 d的东方蜜蜂微孢子虫(分别为NcT1和NcT2)转录组数据,根据P≤0.05且|log₂(Fold change)|≥1的标准,通过比较分析筛选出NcCK vs NcT1, NcCK vs NcT2和NcT1 vs NcT2比较组的DEG。通过相关生物信息学软件对上述DEG进行Venn分析、GO分类和KEGG代谢通路富集分析。根据Nr和KEGG数据库注释信息和相关文献进行对东方蜜蜂微孢子虫的毒力因子和侵染相关因子的统计和分析。通过RT-qPCR验证转录组数据及DEG表达趋势。【结果】从NcCK vs NcT1, NcCK vs NcT2和NcT1 vs NcT2比较组分别鉴定出1 397, 1 ...     

中华蜜蜂精巢和卵巢差异表达基因分析

[目的]在蜂群中,雄蜂与蜂王都有着发育完全的性腺,但两者达到性成熟的时间却是不同步的.本研究旨在探究中华蜜蜂Apis cerana cerana雄蜂与蜂王生殖腺的基因表达差异.[方法]利用Illumina测序技术对中华蜜蜂雄蜂精巢与蜂王卵巢转录组差异表达基因(differentially expressed genes...     

亚致死剂量吡虫啉胁迫对意大利蜜蜂工蜂嗅觉学习行为的影响及其脑部转录组分析

[目的]分析吡虫啉处理对意大利蜜蜂Apis mellifera ligustica嗅觉学习行为及脑部基因转录的影响,为新烟碱类杀虫剂对蜜蜂的负面影响提供依据.[方法]实验室条件下一次性饲喂意大利蜜蜂成年工蜂含有4 ng吡虫啉的50％蔗糖溶液,以饲喂不合吡虫啉的50％蔗糖溶液为对照,通过伸吻反应(proboscis ex...     

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

草鱼雌性个体生长优势明显高于雄性,挖掘精巢和卵巢差异表达的功能基因对草鱼生殖调控具有重要的价值及应用前景。采用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基因与早期草鱼卵巢发育调控相关。本实验为进一步了解草鱼性别基因的调控机理提供数据依据。     

闽楠转录组分析及基因功能注释

采用第二代Illumina HiSeq测序技术对闽楠的木质部、韧皮部、叶片进行转录组测序，分别获得Clean Reads片段41 383 707条、43 343 922条、44 191 586条，经转录本拼接后得到序列总长度达120 535 288 bp的383 331条Conting片段，进一步组装得到平均长度为542 bp的151 729条Unigenes。将闽楠转录组Unigenes进行基因功能注释，与NR数据库比对发现，其与葡萄的相似序列最多（34%），与黄瓜、野草莓、大豆的同源性较低（各占3%）；进行GO功能注释，可将其划分为生物过程、细胞成分、分子功能3大类共计52个分支，与eggNOG数据库比对可分为25类，通过KEGG功能注释可知转录组中涉及的基因共参与了176条代谢通路，其中核糖体和碳代谢获得的注释较多。另外通过MISA软件分析，共获得35 972个SSR位点。其中，单核苷酸、二核苷酸和三核苷酸为优势重复类型，SSR位点数分别为21 762（60.50%），8 931（24.83%），4 924（13.69%）。闽楠转录组分析及基因功能注释为深入开展闽楠遗传育种及分子生物学相关研究奠定基础。     

镉胁迫后旱柳转录组变化分析

随着重金属镉(Cd)应用范围的扩大,由此引发的土壤镉污染问题日益严重。以具有植物恢复潜力的旱柳Salix matsudana作为研究对象,探究不同浓度的Cd (2.5 mg/L, 50 mg/L)胁迫后旱柳无性系1 d、7 d和30 d后基因表达与代谢通路的变化。转录组测序结果表明:共获得102 595个非冗余基因(Unigenes),相同浓度不同时间的差异基因总数为26 623个和32 154个;相同时间不同浓度的差异基因总数为8 550个、3 444个和11 428个。从中筛选得到与Cd胁迫响应密切相关的基因25个,其中金属硫蛋白、ABC转运蛋白、锌和锰转运蛋白等基因的表达不仅会随着Cd胁迫浓度变化而且同时受到胁迫时间的改变而发生改变;油菜素内酯合成通路的ROT3和黄酮类化合物合成通路的FLS、F3H均明显上调。此外Cd胁迫引起旱柳在代谢过程、细胞过程、膜、细胞器、细胞、细胞部分、催化活化和结合蛋白这8个方面发生改变,参与这些GO条目的差异表达基因数随着Cd浓度和胁迫时间的增加而增加。并对转录组信息的可靠性用RT-PCR和酶活性生理实验数据进行了验证。文中通过转录组测序分析旱柳Cd胁迫后的响应机制,从而为旱柳修复土壤Cd污染提供理论指导。     

糜子黄、黑果皮遗传及转录组学分析北大核心CSCD

为了解糜子黑色与黄色果皮的遗传规律以及其形成机理,利用黄色果皮品种‘黄粒糜Ⅰ’和黑色果皮品种‘2016106Ⅰ’构建杂种F1、F2代,分析果皮颜色的遗传规律,利用F3代籽粒进行转录组测序,挖掘影响糜子籽粒颜色的关键基因。试验结果表明,黑色对黄色为显性,由单基因控制;根据GO和KEGG富集分析结果,共筛选出13个与类黄酮合成途径相关的差异基因,仅有C2845_PM02G08740、C2845_PM04G29280和C2845_PM09G22680为上调表达基因,其余10个均为下调表达基因,上述基因编码肉桂醇脱氢酶和肉桂酰CoA还原酶等多种酶,与颜色的形成密切相关。     

FGF13敲低的肺癌A549细胞株构建及转录组学分析

本研究的目的是探究成纤维细胞生长因子13(fibroblast growth factor13, FGF13)影响非小细胞肺癌发生发展的分子机制。通过构建干扰FGF13的真核表达载体并进行慢病毒包装,获得FGF13敲低的肺癌A549细胞株。利用转录组测序(RNA Sequencing, RNA-seq)技术检测基因表达谱并筛选差异表达基因(differentially expressed genes, DEGs),并对差异基因进行GO和KEGG富集分析,对部分DEGs的表达水平进行验证。采用STRING在线数据库筛选DEGs编码的蛋白相互作用网络(protein-protein interaction, PPI)中的关键基因,并通过UALCAN数据库对关键基因进行预后分析。与对照组细胞相比,FGF13敲低的肺癌A549细胞株共筛选出1 114个DEGs(P≤0.05,差异倍数≥2),其中表达上调和下调的DEGs分别有672和442个。实时荧光定量PCR(real-time quantitative PCR, RT-qPCR)验证CXCL10和SELE等6个基因与测序结果一致。GO富集分...     

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

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

叉角厉蝽触角转录组及嗅觉相关基因分析

叉角厉蝽Eocanthecona furcellata是一种在生物防治方面有重要作用和潜力的捕食性天敌昆虫,触角是昆虫进行信息交换的重要器官,而嗅觉相关基因则是调控天敌昆虫捕食行为的重要分子基础。为获得叉角厉蝽触角转录组数据库,挖掘叉角厉蝽嗅觉相关基因,本研究利用Illumina高通量测序平台对叉角厉蝽雌雄成虫与5龄若虫触角进行转录组测序。成功构建了叉角厉蝽触角转录组,获得了67 843条unigenes, N50长度为2 300 bp。与七大公共数据库比对注释到27 686条unigenes,其中NR数据库注释最多(33.33%),且与茶翅蝽Halyomorpha halys相似度最高(64.20%)。14 258条注释到GO数据库中,分为生物过程、细胞组分和分子功能3个大类42个亚类;KEGG代谢途径分析表明,7 703条形成282条代谢通路,其中被注释在信号传导通路中的unigenes最多(11.50%)。进一步基因注释分析,鉴定得到134个候选嗅觉相关基因,32个化学感应蛋白基因(Chemosensory protein genes, CSP),10个气味结合蛋白基因(Odor...     

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

[目的]本研究旨在通过转录组测序技术分析低温胁迫引起的阿尔泰蝠蛾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途径.其中代谢过程、催化、结合活性类群占有重要比例,核糖体、碳代谢、剪接体等途径显著富集.另外,热激蛋白、昆虫表皮蛋白、海藻糖酶、超氧化物歧化酶等非生物胁迫相关基因显著上调表达.[结论]阿尔泰蝠蛾幼虫低温胁迫转录组分析揭示,代谢过程、细胞过程、生物调节、对刺激的反应等生物学过程相关基因和部分非生物胁迫响应基因显著上调表达,提示蝠蛾幼虫可能从抗氧化防御、分子伴侣、体温调节和维持细胞的渗透平衡等多方面应对低温胁迫.     

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

[目的]本研究旨在通过转录组测序技术分析低温胁迫引起的阿尔泰蝠蛾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途径.其中代谢过程、催化、结合活性类群占有重要比例,核糖体、碳代谢、剪接体等途径显著富集.另外,热激蛋白、昆虫表皮蛋白、海藻糖酶、超氧化物歧化酶等非生物胁迫相关基因显著上调表达.[结论]阿尔泰蝠蛾幼虫低温胁迫转录组分析揭示,代谢过程、细胞过程、生物调节、对刺激的反应等生物学过程相关基因和部分非生物胁迫响应基因显著上调表达,提示蝠蛾幼虫可能从抗氧化防御、分子伴侣、体温调节和维持细胞的渗透平衡等多方面应对低温胁迫.     

雄性不育枸杞花药比较转录组分析

为筛选宁夏枸杞(Lycium barbarum L.)雄性不育花药发育相关的差异表达基因,本研究开展了雄性不育品种‘宁杞5号’和可育品种‘宁杞1号’成熟花粉时期花药转录组分析,共获得差异表达基因1 759个,在‘宁杞5号’中有753个基因上调表达,1 006个基因下调表达。KEGG富集表明,这些差异表达基因主要富集到氨基酸的生物合成、植物激素信号传导、苯丙素生物合成、淀粉和蔗糖代谢、脂肪酸代谢等相关通路。对有注释信息的差异表达基因进行NR和Uniprot数据库检索,筛选到28个花药发育相关基因,其中有20个在‘宁杞5号’中下调表达,8个上调表达。选取其中的5个差异表达基因进行RT-qPCR验证,结果与转录组数据趋势一致。甘油醛-3-磷酸脱氢酶(Glyceraldehyde-3-phosphate dehydrogenase, GAP)、束状阿拉伯半乳糖蛋白3(Fasciclin-like arabinogalactan protein 3,FLAs)和花药特异表达启动子LAT52系统发育树分析表明,3个基因分别与同科的物种在同一个进化分支上且同源性较高。这些差异表达基因可能在宁夏枸杞‘...