首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到17条相似文献,搜索用时 473 毫秒
1.
【 目 的】 探 究 铁 皮 石 斛(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 信号途径和甘油酯代谢通路响应耐盐。  相似文献   

2.
【目的】研究干旱胁迫下,棉花DNA损伤修复基因的表达情况,从整体水平探讨棉花DNA损伤修复相关基因表达与抗旱的相关性。【方法】用2.5%PEG6000处理棉花幼苗,利用RNA-Seq技术对干旱胁迫下棉花幼苗的转录组进行测序。【结果】从棉花干旱胁迫响应的转录组中,筛选获得差异表达的DNA损伤修复相关基因共51个,其中差异表达的上调基因23个,差异表达的下调基因28个,干旱胁迫能够影响棉花DNA损伤修复相关基因的表达。选取4个差异基因进行生物信息学分析及qRT-PCR验证,HMGB1、rec A1、UDGs和GMP synthase基因的表达变幅有一定的差异,但基因的表达趋势一致,棉花转录组测序结果通过qRTPCR验证是可靠的。【结论】DNA损伤修复相关基因可能与干旱胁迫有一定的相关性,干旱胁迫能够影响棉花DNA损伤修复相关基因的表达。  相似文献   

3.
【目的】研究CaCl2缓解酸枣盐胁迫的分子机制。通过转录组测序分析,为研究CaCl2缓解酸枣盐害分子机制提供参考。【方法】研究以水培酸枣幼苗为材料,150 mmol/L NaCl和150 mmol/L NaCl + 10 mmol/L CaCl2分别处理6 h,利用高通量测序技术对酸枣幼苗叶片和根系进行了转录组测序与组装,对获得的所有Unigene在GO、KEGG数据库中进行功能注释。【结果】 与对照组相比,在外源CaCl2处理下酸枣幼苗叶片上调和下调的差异表达基因分别7 365个和1 247个,根系分别有762个和4 861个。【结论】GO数据库比对,酸枣幼苗叶片富集了43个GO条件,根系富集了42个GO条件,叶片和根系均是催化活性被富集到差异基因最多;KEGG数据库比对,差异基因主要涉及植物病原互作、MAPK信号通路-植物、植物激素信号转导、苯丙醇生物合成、戊糖和葡萄糖醛酸相互转化等。  相似文献   

4.
【目的】掌握不同生长规格墨瑞鳕个体间差异表达基因的表达特点,为其功能相关基因深度挖掘及分子遗传育种提供科学依据。【方法】挑选同一养殖条件下极大个体和极小个体的墨瑞鳕,构建肌肉组织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且|log2 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信号通路及细胞凋亡等能量代谢相关通路上,通过共同协调而对墨瑞鳕的生长发育起调控作用。  相似文献   

5.
库尔勒香梨PsiERF基因的克隆和表达分析   总被引:1,自引:1,他引:0  
【目的】筛选克隆库尔勒香梨萼片脱落与宿存的相关基因,研究其在脱萼样品与宿萼样品中表达。【方法】 以库尔勒香梨脱萼和宿萼花器官转录组测序数据为基础,从Unigenes中选出一条转录组和数字表达谱测序结果共有的与萼片脱落宿存相关的具有AP2/ERF结构域的差异基因,命名为PsiERF,对其进行生物信息学分析,通过实时荧光定量PCR分析其在花期不同阶段不同器官中的表达量。【结果】PsiERF cDNA序列的长度是1 195 bp,包含一个编码264个氨基酸长度为795 bp的开放阅读框(ORF),一个123 bp的5'端非编码区和一个277 bp的3'端非编码区。库尔勒香梨的PsiERF的核酸与氨基酸序列与其他植物的ERF109基因具有高度的同源性。实时荧光定量结果显示,其在全花样本中,宿萼组全花样品表达量在所有时期均高于脱萼组全花样品,而在宿萼组样品末花期萼片表达量极显著高于子房。【结论】筛选、克隆并鉴定一个新的库尔勒香梨PsiERF基因,该基因属于AP2/ERF家族中的ERF亚家族,与库尔勒香梨的萼片宿存脱落存在密切关系。  相似文献   

6.
【目的】探究脊尾白虾(Exopalaemon carinicauda)鳃组织在渐变式低氧—复氧胁迫下的分子调控机制,为今后开展脊尾白虾耐低氧品系(种)的选育提供理论指导。【方法】通过模拟自然低氧环境的形成过程,分别于低氧处理0(对照)、3和6 h及复氧后1和8 h采集脊尾白虾鳃组织,利用Illumina HiSeqTM 4000测序平台进行转录组测序分析,经过滤和Trinity组装获得Unigenes,选取Nr、Swiss-Prot、KEGG和COG/KOG等数据库进行注释分析,在Omicsmart平台上完成差异表达基因筛选及其表达趋势分析,然后进行GO功能注释分析和KEGG信号通路富集分析,并随机选取5个差异表达基因进行实时荧光定量PCR验证。【结果】脊尾白虾鳃组织样本转录组测序数据经过滤后的Cleanreads进行Trinity组装共获得93227条Unigenes,其长度范围在201~35402 bp,平均长度为834 bp,N50长度为1352 bp。通过组间两两比较分析,共鉴定出4750个差异表达基因,其中上调差异表达基因3557个、下调差异表达基因2829个;超过50%的差异表达基因被显著富集到6种基因表达趋势模式中(P<0.01),具体表现为:Profile 0模式富集到415个基因,Profile 5模式富集到201个基因,Profile 11模式富集到371个基因,Profile 13模式富集到841个基因,Profile 17模式富集到387个基因,Profile 18模式富集到411个基因。6种基因表达趋势模式中的差异表达基因被注释到代谢进程、细胞进程、单一有机体进程、细胞、细胞零件、大分子复合物及催化活性等GO功能条目上;而KEGG信号通路富集分析结果显示,以Profile 13模式中的差异表达基因富集到最多信号通路(86条),其中呈显著富集的有8条,分别为核糖体、碳代谢、氧化磷酸化、氨基酸生物合成、内质网蛋白质加工、糖酵解/糖异生、谷胱甘肽代谢和蛋白输出。【结论】脊尾白虾鳃组织在受低氧胁迫早期通过合成蛋白质及提高代谢能力来抵御低氧环境,随着低氧胁迫时间的延长,物质合成和能量代谢活动均显著下降;但在复氧后随着复氧时间的延长,其蛋白质合成和能量代谢水平又逐渐升高恢复至常氧水平。  相似文献   

7.
【目的】探究高粱耐盐胁迫响应机制,挖掘高粱耐盐胁迫基因,为高粱耐盐育种提供理论基础。【方法】 以高粱感盐品种L甜和耐盐品种石红137为供试材料,采用水培试验。待高粱植株长至三叶一心期,使用2%NaCl溶液对幼苗进行盐胁迫,分别设置0(对照)、1和24 h处理,每个处理3次重复。测定不同处理样品株高、根长、干物重、Na +含量和叶绿素相对含量(SPAD值),并依托Illumina HiSeq 2000平台进行转录组测序分析。利用FPKM方法计算基因表达量,在差异表达基因检测过程中,将差异表达倍数(fold change)≥2且FDR<0.001作为筛选标准。通过Gene Ontology和KEGG Pathway数据库对参与高粱不同时间盐胁迫差异表达基因进行分析注释。 【结果】 盐胁迫处理对高粱株高、根长、干物重等性状无显著影响,对钠离子含量和SPAD值影响显著。石红137株高、根长、钠离子含量和SPAD值均高于L甜。转录组测序结果鉴定得到已知基因26 628个,新基因866个。石红137中的差异基因数目高于L甜。石红137中,0 h VS 1 h、0 h VS 24 h、1 h VS 24 h三组的差异基因数目分别为375、4 206和3 750个。感盐品种L甜中,0 h VS 1 h、0 h VS 24 h、1 h VS 24 h三组的差异基因数目分别为167、2 534和1 612个。GO分析共获得25个功能注释,分别为光合作用、细胞物质代谢、翻译过程以及激素合成等与盐胁迫相关的差异表达基因。KEGG分析发现盐胁迫1 h表达差异基因富集在植物激素信号转导途径,涉及脱落酸(abscisic acid,ABA)、生长素(auxin,AUX)、细胞分裂素(cytokinin,CTK)、赤霉素(gibberellins,GS)、乙烯(ethylene,ETH)过程等共71个基因。盐胁迫24 h表达差异基因富集于光合作用相关途径,涉及Lhca、Lhcb、磷酸烯醇式丙酮酸羧化酶(phosphoenolpyruvate carboxylase,PPC)、磷酸核酮糖激酶(phosphorylribonucleic kinase,PRK)等20个基因。类黄酮生物合成代谢途径差异可能是引起石红137和L甜的耐盐能力差异的原因之一,花青素还原酶(anthocyanidin reductase,ANR)和黄酮醇合成酶(flavonol synthase,FLS)参与类黄酮生物合成途径。【结论】 高粱的盐胁迫过程是一个复杂的生物过程,依赖于多个基因在复杂网络中的平衡表达。盐胁迫条件下,高粱应对环境刺激受到激素信号转导和光合作用的控制。类黄酮生物合成途径在耐盐品种中起到了重要作用。  相似文献   

8.
【目的】对青花菜花蕾进行转录组测序分析,并挖掘与蜡粉合成相关基因,为探明青花菜花球表面蜡粉形成的分子机制提供理论参考。【方法】分别提取野生型和蜡粉缺失型青花菜花球总RNA,采用Illumina HiSeqTM 2500平台进行转录组测序,获得高质量Clean reads,采用Trinity进行序列组装后获得青花菜Unigene库,将获得的Unigene序列与Nr、Nt、KEGG、Pfam、KOG/COG、Swiss-Prot和GO数据库比对,获得基因功能注释信息;使用DESeq2进行差异表达分析。【结果】共获得44.68 Gb Clean data,De novo组装得到41244条Unigenes,N50长度为1847 bp。从所获得的Unigenes中筛选出8685个差异表达基因(DEGs)(上调基因5747个,下调基因2938个),共有8038个基因被注释到不同数据库,其中,5220个基因注释到Pfam数据库; 2066个基因注释到COG数据库,3866个基因注释到KOG数据库; 2580个差异表达基因被注释到75个转录因子家族中,注释最多的是MYB家族(235个);GO数据库中6095个差异表达基因注释到细胞组分、分子功能和生物学过程三大类的52个功能分类; KEGG数据库中,1671个差异表达基因富集到138条代谢通路,其中13个差异表达基因与脂肪酸合成有关,7个差异表达基因与蜡粉生物合成途径有关。【结论】转录因子MYB家族在调控青花菜蜡粉合成中发挥重要作用。蜡粉合成过程中相关酶基因的差异表达是调控青花菜蜡粉合成的关键,尤其是野生型和蜡粉缺失突变体中特异性表达的差异表达基因,可作为后续研究青花菜花球表面蜡粉形成分子机制的对象。  相似文献   

9.
【目的】开展荸荠球茎发育过程转录组测序研究,为研究荸荠球茎发育过程相关基因表达信息提供参考。【方法】运用高通量测序技术,对荸荠球茎不同发育时期进行转录组测序研究。【结果】组装共得到223 182条转录本和90 542条Unigene,平均长度为809 bp,N50为1119。组装完整性较高,效果较好。对所得Unigene进行不同数据库注释,共有50 583条Unigene成功注释到7个数据库(NR、GO、NT、Pfam、KEGG、KOG及Swiss-prot)。对差异表达基因分析,发现球茎膨大初期的差异表达基因数目最多,为最活跃阶段。GO功能富集分析结果表明,33 205个基因获得功能注释,分为分子功能、细胞组分和生物学过程等3大类和54个亚类。COG功能分类结果表明,17 743个基因分布于25个功能区域,其中碳水化合物代谢占重要地位。KEGG代谢通路注释结果表明有20 667个基因获得功能注释,共有116条代谢途径,其中淀粉-蔗糖代谢占主要作用。【结论】利用高通量转录组测序技术首次建立了荸荠优良品种‘桂蹄3号’球茎的转录组数据库,为进一步研究荸荠球茎淀粉生物合成相关基因的功能及形成的分子机制提供了数据基础。  相似文献   

10.
[目的]分析梭梭同化枝及其叶苞状虫瘿的转录组差异,为揭示梭梭叶苞状虫瘿形成的分子机制及梭梭害虫防治提供理论参考.[方法]以梭梭同化枝及其叶苞状虫瘿为材料,利用Illumina HiSeq 2000高通量测序平台对其进行转录组测序及差异表达基因分析,并通过测定可溶性糖和叶绿素含量验证分析结果.[结果]转录组测序共获得159147条Unigenes,平均长度为551 bp,N50为752 bp,在Swissport、GO、Pfam、NR、NT、KOG/COG和KEGG等数据库均被注释的Unigenes有8521条(5.35%),至少在1个数据库注释的Unigenes有66714条(41.91%),共有49066条Unigenes(30.83%)与物种匹配成功.以q value<0.005且|log2FoldChange|>1为条件,共筛选出3444条差异表达的Unigenes,其中,有1299条上调表达,2145条下调表达.与梭梭同化枝相比,梭梭叶苞状虫瘿上调表达差异基因富集在代谢过程、有机物代谢过程、催化活性、碳代谢作用和核糖体代谢等通路,下调表达基因富集在细胞、细胞组成、胞内细胞成分和细胞器、嘧啶代谢、细胞周期、DNA复制、光合作用及苯丙烷和类黄酮生物合成等过程.不同发育期梭梭叶苞状虫瘿中的可溶性糖含量高于梭梭同化枝,而叶绿素含量均低于梭梭同化枝,进一步验证梭梭叶苞状虫瘿的糖代谢相关基因上调表达而光合作用相关基因下调表达的转录组分析结果.[结论]梭梭叶苞状虫瘿中核糖体数量增加,翻译过程加快,代谢水平提高,糖代谢产物积累,但其光合作用降低,伸长与生长减慢,推测其通过吸收利用同化枝的养分以满足自身需求.  相似文献   

11.
【目的】分析西瓜食酸菌(Acidovorax citrulli)在铜胁迫下转录水平的响应特征,研究对其抗铜系统和抗铜机制,为进一步解析西瓜食酸菌的铜耐受分子机制积累数据。【方法】以西瓜食酸菌亚群I的抗铜菌株FC440铜胁迫和无胁迫培养下的菌体为材料,提取RNA并利用RNA-Seq技术获取细菌对铜胁迫的应答数据,分析数据质量、差异表达基因的类型、KEGG富集与Pathway注释、GO富集分析等转录组水平特征,并进行qPCR验证。【结果】共获得6个转录组文库,总数据量23.6 Gb;该细菌共有4 344个基因表达,其中有466个基因差异表达(上调266个,下调200个);差异表达的基因显著富集于ABC转运系统的跨膜转运、双组分系统通路的信号转导以及氨基酸代谢等约9大通路;基因表达趋势与转录组测序数据的趋势一致。【结论】铜胁迫下谷氨酸、谷胱甘肽等氨基酸代谢在西瓜食酸菌抗铜胁迫中发挥重要作用;双组分系统参与西瓜食酸菌对铜离子胁迫的响应且受到负调控;该菌中可能存在包含ABC转运体的跨膜转运、氧化还原反应、胞内束缚沉淀等多种系统互作的抗铜机制。  相似文献   

12.
【目的】谷子具有抗旱、耐贫瘠的特性,逐渐成为研究禾本科作物的模式作物。本研究以前期筛选并获得的谷子耐盐品种和敏感品种为材料,通过RNA-Seq测序筛选鉴定谷子的盐胁迫响应基因,揭示其响应盐害机制,为培育谷子抗盐新种质提供理论依据。【方法】 对14个谷子品种进行了萌发期抗盐筛选。谷子不同品种的种子经150 mmol·L -1NaCl处理萌发培养7 d后,分别统计各品种种子的萌发率、根长和芽长,综合各项指标鉴定到豫谷2为耐盐品种、安04为盐敏感品种。以这两个品种盐胁迫前、后萌发期种子为材料,应用RNA-Seq进行转录组测定,分析鉴定盐害下差异表达基因(differentially expression gene,DEG),并对DEG进行GO和KEGG功能富集分析。同时应用qRT-PCR对随机挑选的15个DEG表达模式进行验证,并与RNA-Seq结果进行相关性分析。 【结果】 耐盐品种豫谷2、盐敏感品种安04种子经盐胁迫处理后,分别鉴定出2 786个和4 413个DEG;其中,每个品种在NaCl处理前及处理后分别鉴定出1 470和3 826个DEG。GO和KEGG聚类分析显示,NaCl处理下参与信号转导、抗氧化、有机酸的合成及转运以及次生代谢等相关的DEG在谷子萌发期种子应对盐胁迫响应过程中具有关键性的作用,DEG主要集中在胁迫响应、离子的跨膜转运、氧化还原、次生代谢及有机酸、多胺、苯丙烷的合成等过程。qRT-PCR对15个DEG表达模式的检测结果与RNA-seq结果相关系数为0.8817。其中编码离子通道的基因(HKT)、过氧化物酶基因(POD)、黄烷酮-3-羟化酶基因(FL3H)、MYB转录因子基因等在豫谷2中表达量较高,显示其在谷子萌发期种子响应盐胁迫中可能发挥一定作用。【结论】 谷子耐盐品种及盐敏感品种的种子对盐胁迫的响应程度具有一定的差异性,且品种对盐害的耐受能力主要与基因对胁迫的响应程度相关,而与DEG的数量相关性不大。  相似文献   

13.
【目的】对高温胁迫下萝卜幼苗进行转录组测序分析,筛选不同耐热性萝卜品种间的差异表达基因,为阐明萝卜幼苗响应高温胁迫的分子机制提供理论参考。【方法】以耐热萝卜品种1116T和热敏感萝卜品种Wr129为对象,对其进行高温(40℃)胁迫处理0(常温,对照)和24 h,提取各时间点样品总RNA进行转录组测序,以|log2 FoldChange|≥1,P<0.05且TPM≥10为标准,筛选出差异表达基因(DEGs),并进行GO功能注释和KEGG代谢通路富集分析。【结果】随着高温胁迫时间的增加,Wr129和1116T幼苗茎叶逐渐发生萎蔫,胁迫24 h后Wr129叶片萎蔫下垂且干枯变黄,而1116T叶片轻度萎蔫下垂,仍保持绿色。12个文库共获得81.45 Gb Clean data,平均有62.30%的reads比对到萝卜参考基因组的一个位置。从2个不同耐热性品种中共鉴定到2701个差异表达基因。GO功能注释结果显示,2629个差异表达基因在GO数据库中得到注释,其中,有559个差异表达基因被注释为细胞组分类别中的叶绿体被膜、叶绿体基质和叶绿体类囊体3个条目,470个差异表达基因被注释为与胁迫响应相关的生物学过程。KEGG代谢通路富集结果显示,1166个差异表达基因富集到124条代谢通路,其中,有74个差异表达基因富集到光合作用和光合作用-天线蛋白代谢通路,23个基因富集到卟啉和叶绿素代谢通路,大部分基因在高温胁迫下被诱导; 35个差异表达基因富集到谷胱甘肽代谢途径,其中14个谷胱甘肽-S-转移酶基因(GSTs)在高温胁迫下上调表达。高温胁迫24 h,Wr129和1116T幼苗叶片PSII的最大光化学效率(Fv/Fm)和总叶绿素含量均下降,但1116T的下降幅度较小。【结论】 Wr129和1116T的耐热性存在明显差异,其原因是1116T植株可通过调控光合作用和抗氧化系统相关基因的表达,减缓Fv/Fm和叶绿素含量下降速率,从而提高植株对高温胁迫的耐受能力。  相似文献   

14.
【目的】土壤盐渍化是制约作物生产的重要非生物胁迫因子之一,高粱耐盐性强,进行高粱耐盐基因挖掘及分子机制研究是开发和利用盐渍土壤的有效途径,通过转录组测序分析与高粱耐盐相关的基因调控机制和代谢通路,挖掘高粱耐盐潜力。【方法】 通过以筛选出的极耐盐品种八叶齐和盐极敏感品种PL212为试验材料,采用盆栽沙培,在播种后20 d(5叶期)采用180 mmol·L -1的 NaCl 溶液漫灌模拟盐逆境,盐胁迫48 h后取幼叶,并连同对照(未经过盐处理)的同期幼苗共4个样品提取RNA,进行转录组测序,采用qRT-PCR方法对测序结果进行验证。 【结果】 耐盐和盐敏感材料分别在盐渍和非盐渍处理下的4个样品间共检测到1 338个差异表达基因,包括819个上调基因和519个下调基因。聚类分析发现在应答盐渍胁迫逆境时,5个依赖性氧合酶超家族蛋白、4个富含半胱氨酸的激酶、3个谷胱甘肽S-转移酶和3个重金属运输/解毒超家族蛋白相关基因表现出明显的上调表达和下调表达,还发现1个K +转运蛋白基因在耐盐调节中起着重要作用。GO分析发现在15 418个基因中获得4 528个有效GO注释条目,同时耐盐和盐敏感材料在遭受盐逆境时的生物过程、细胞组分和分子功能3个方面均存在较大差异。生物过程中代谢过程、细胞过程耐盐材料明显高于盐敏感材料,耐盐材料的生理过程中较盐敏感材料增加了多生物过程和定位这两个过程,很可能与耐盐材料盐抗性较强密切相关。差异基因KEGG分析结果显示耐盐和盐敏感材料在对照和盐渍胁迫条件下的苯丙烷类生物合成、苯丙氨酸代谢、类黄酮生物合成3个途径中差异基因表达较多,可能是造成耐盐和盐敏感材料耐盐性差异较大的重要原因。 【结论】 高粱耐盐调控基因表达涉及生物过程、细胞组分和分子功能多个方面,生物过程和定位这两个过程是提高高粱耐盐性的关键;苯丙烷类生物合成、苯丙氨酸代谢、类黄酮生物合成3个途径的基因表达很可能是造成盐害的重要原因。  相似文献   

15.
史毅  牛奎举  马晖玲 《中国农业科学》2017,50(17):3323-3336
【目的】确定匍匐翦股颖(Agrostis stolonifera)接种立枯丝核菌(Rhizoctonia solani)后基因种类和表达量在转录水平的变化规律,明确草坪草病原菌侵染响应的关键基因。【方法】匍匐翦股颖生长14 d后采用麦粒培养物接种立枯丝核菌,接种3 d后选取感病叶片和未接种叶片提取RNA,进行转录组高通量测序,然后利用生物信息学分析,用Trinity组装匍匐翦股颖转录组,以组装子为参考,以|log2(fold change)|1,q-value0.005为阈值选取感病和健康匍匐翦股颖叶片转录组的差异表达基因,并用i TAK软件分析其转录因子家族及表达变化,与植物R基因库进行blast分析R蛋白分类、利用Mapman软件分析生物胁迫信号通路相关基因的表达变化。【结果】高通量测序得到125 253 092条高质量待分析reads。经Trinity从头组装后,得到466 761条转录本。过半数的转录本长度为700 bp以上,组装结果 N50=1 100 bp。使用CD-HIT选择334 212条转录本(所有转录本的71.60%)作为Unigene,平均长度573 bp,N50=791 bp。接种后植物比接种前植物基因有7 937个上调表达,1 570个下调表达。上调基因中296个,下调基因有142个都可被定义为转录因子,分布在58个转录因子家族中,其中锌指蛋白包含转录因子C2H2最多,有54个,C3H次之,为22个。差异基因中451个可定义为植物R蛋白表达基因,可分为33类,其中包含NBS-LRR结构域的抗病蛋白、LRR受体蛋白激酶、ABC-2类型的转运蛋白、U-box结构域蛋白激酶和热激蛋白这5类基因变化最显著。差异基因中大量上调表达基因可富集在病原识别、活性氧消除、信号传导、细胞凋亡、病程相关蛋白等生物胁迫相关基因类别,下调基因显著富集在植物生长发育相关类别和通路。q RT-PCR验证了随机挑选差异基因的表达量变化,均与RNA-seq分析结果一致,其中包括12个C2H2转录因子基因,10个C3H转录因子基因,以及12个R蛋白基因。【结论】病原菌侵染后,引起匍匐翦股颖大量基因表达变化,其中转录因子、R蛋白以及抗性相关基因多为上调表达,作用为抑制病原菌扩散,而生长发育相关基因表达下调,这些进程共同使匍匐翦股颖产生对立枯丝核菌的先天基础抗性。  相似文献   

16.
【目的】鉴定筛选出与卵形鲳鲹卵巢发育相关的候选基因及信号通路,为揭示其卵巢性成熟过程的分子机制打下基础。【方法】挑选卵巢发育处于?期和Ш期的雌性卵形鲳鲹,分别构建卵形鲳鲹卵巢?期和Ш期的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信号通路等与卵形鲳鲹卵巢的发育密切相关,共同调节卵巢的发育与成熟,在卵巢性成熟过程中发挥重要作用。  相似文献   

17.
芝麻发育转录组分析   总被引:18,自引:3,他引:15  
【目的】系统了解芝麻发育及种子形成转录组特征,丰富芝麻转录组数据信息。【方法】选用6份芝麻样品(5个不同芝麻品种的完整植株、1份不同发育阶段的芝麻籽粒),构建转录组测序文库并进行Illumina RNA-seq双端测序及生物信息学分析。【结果】共获得原始数据12.69 Gb,有效数据8.80 Gb。通过de novo拼接获得了长度大于100 bp的转录物26 837条(http://www. ncbi.nlm.nih.gov/genbank/TSA.html,登录号:JP631635—JP668414);转录物总长度18.35 Mb,平均长度683 bp,N50长度1 006 bp。转录物注释结果显示,25 331条转录物序列具有同源比对信息;1 506条转录物序列无匹配(no hits)序列信息,可能为芝麻特有的基因序列。采用COG、GO功能分类工具可将已注释转录物序列划分为24个或42个功能类别,共涉及物质及能量代谢、信号传导、转录调控及防卫反应等诸多生理生化过程。通过比较不同材料间转录物序列及表达水平,初步确定1 277条序列在种子形成过程中表达量下调10倍以上,990条序列仅在芝麻植株中表达而未在籽粒中表达;660条序列在种子形成过程中表达量上调10倍以上,296条序列可能与种子形成特异相关。【结论】利用高通量测序技术对芝麻野生种和不同栽培种现蕾期植株以及种子形成过程的转录组进行研究,揭示了芝麻发育转录组的整体表达特征,在得到大量芝麻转录组unigene序列的同时,获得了一批在芝麻生长发育及籽粒形成过程中有重要功能的基因序列。为深入开展芝麻生长发育、籽粒发育相关基因功能及调控以及芝麻分子标记开发等研究提供了丰富的数据资源。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号