共查询到16条相似文献，搜索用时 250 毫秒
1.
半番鸭与番鸭精巢组织差异表达转录组测序分析 总被引：2，自引：2，他引：0
【目的】研究半番鸭与番鸭精巢组织转录组差异表达基因,为进一步阐明半番鸭不育的遗传机制奠定理论基础。【方法】利用转录组测序方法对半番鸭和番鸭的精巢组织进行研究,筛选其差异表达基因,并对其功能进行注释和荧光定量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信号通路在鸭的生殖活动中发挥了重要作用,为进一步探索半番鸭生殖系统的分化机理提供可靠的参考依据。 相似文献
2.
【目的】分析镉胁迫下中辽1号杨的转录水平,筛选差异表达基因,挖掘与镉胁迫相关的功能基因,为深入探索中辽1号杨对镉胁迫响应的分子机制提供理论依据。【方法】以中辽1号杨为试验材料,采用盆栽试验的方法将其扦插于土壤Cd含量为20 mg/kg（M20）的花盆中,以不加Cd为对照（CK）,80 d后采集不同处理中辽1号杨叶片进行转录组测序,使用DESeq2软件筛选CK与M20处理中辽1号杨的差异表达基因,将得到的差异基因在基因本体数据库（gene ontology,GO）、京都基因与基因组百科全书（Kyoto encyclopedia of genes and genomes,KEGG）、蛋白相邻类聚簇数据库(clusters of orthologous groups of proteins,COG)中进行注释,分析差异表达基因在不同数据库中的注释信息。随机挑选6个差异表达基因（超氧化物歧化酶（SOD）、丝氨酸/苏氨酸蛋白激酶（STK）、脱落酸（ABA）、T复合蛋白1（TCP1）、MYB基因（MYB）、丝裂原活蛋白激酶12（MAPK12））,利用Primer 5软件设计特异性引物进行实时荧光定量PCR（RT-qPCR）试验,检测这6个基因表达水平的变化,验证转录组结果的准确性。【结果】在CK和M20处理的中辽1号杨中,共发现3 812 个差异表达基因,其中上调表达的有2 209 个,下调表达的有1 603 个。GO 分析发现,差异表达基因注释到3大类49个功能组中,其中与代谢过程有关基因最多。转录组KEGG代谢通路分析发现,差异表达基因主要在信号转导途径、代谢途径和生物合成途径富集。COG分析发现,1 455个差异基因被注释到22种分类中,其中一般功能预测基因注释最多,其次是信号转导机制基因。结合各数据库分析结果发现,Cd胁迫下中辽1号杨转录组的注释结果中代谢过程和信号转导过程相关基因较多。镉胁迫下差异表达基因最多的家族是ABC、MYB、WRKY、bHLH和NAC基因家族,挖掘出与镉胁迫相关基因WRKY家族基因47（WRKY47）、硝酸盐转运蛋白基因（NRT）、ABC家族转运蛋白基因2（ABC2）、苹果酸脱氢酶（MDH）、谷胱甘肽S-转移酶基因（GSTs）、NAC家族基因2（NAC2）的表达量显著上调,MYB家族基因44（MYB44）、重金属相关异戊二烯化植物蛋白基因39（HIPP39）的表达量显著下调。RT-qPCR结果显示,随机挑选6个差异表达基因表达量的变化趋势与转录组测序结果一致。【结论】镉胁迫后中辽1号杨转录组的差异表达基因主要富集在代谢过程和信号转导过程,挖掘出与镉胁迫相关基因,分别是WRKY47、NRT、ABC2、MDH、GSTs、NAC2、MYB44、HIPP39,为深入探索杨树对镉胁迫响应的分子机制提供了理论依据。 相似文献
3.
【目的】筛选出凡纳滨对虾生长发育的功能基因及代谢调控网络,揭示其生长发育的分子机制,为后续开展凡纳滨对虾分子生物学研究提供宝贵的基因数据来源。【方法】以快速生长群体和慢速生长群体的凡纳滨对虾肌肉组织为研究材料,通过Illumina HiSeqTM2500平台对构建的cDNA文库进行高通量测序分析,以StringTie进行拼接组装后利用DESeq2筛选差异表达基因,并基于KOG、GO、nr、COG、Swiss-Pro、KEGG和Pfam等数据库进行差异表达基因功能注释分析。【结果】经拼接组装共获得53458844条Clean reads,各样品Clean reads的Q30均在93.00%以上;Illumina测序获得的Clean reads与参考基因组的比对效率在87.75%~87.80%,说明转录组测序数据真实可靠。采用SnpEff进行SNP/InDel变异注释分析,结果显示,SNP位点中以A>G、G>A、C>T和T>C等4种类型的数量较多（14950~21562个）,InDel位点则以SYNONYMOUS_CODING的数量最多（33630个）。使用StringTie对Mapped reads（比对到参考基因组的Reads）进行拼接,共发掘到4607个新基因,分别输入COG、GO、KEGG、KOG、Pfam、Swiss-Prot、eggNOG和nr数据库中进行序列比对,最终发现共有1098个新基因被注释,以被nr数据库注释的新基因数量最多（1077个）,而被COG数据库注释的新基因数量最少（416个）。基于Q-value<0.05且Fold Change>2的筛选条件,共获得1408个差异表达基因（661个为显著上调表达基因,747个为显著下调表达基因）;1408个差异表达基因被注释到53个GO功能条目中,其中,22条被注释到生物学过程（Biological process）,16条被注释到细胞组分（Cellular component）,15条被注释到分子功能（Molecular function）;KEGG信号通路富集分析发现3条重要的信号通路,分别是溶酶体通路（Lysosome）、氨基糖和核苷酸糖新陈代谢（Amino sugar and nucleotide sugar metabolism）及鞘脂类代谢（Sphingolipid metabolism）。在溶酶体通路中,CTSL基因、Nramp基因、MyoG基因和Myf5基因在凡纳滨对虾快速生长群体肌肉组织中呈上调表达,而Trypsin基因呈下调表达。【结论】通过转录组测序分析从快速生长群体和慢速生长群体的凡纳滨对虾肌肉组织中筛选出1408个差异表达基因（661个为显著上调表达基因,747个为显著下调表达基因）,主要富集在溶酶体、氨基糖和核苷酸糖新陈代谢及鞘脂类代谢等通路上,在对虾肌肉生长发育过程中发挥重要作用。 相似文献
4.
【目的】对青花菜花蕾进行转录组测序分析，并挖掘与蜡粉合成相关基因，为探明青花菜花球表面蜡粉形成的分子机制提供理论参考。【方法】分别提取野生型和蜡粉缺失型青花菜花球总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家族在调控青花菜蜡粉合成中发挥重要作用。蜡粉合成过程中相关酶基因的差异表达是调控青花菜蜡粉合成的关键，尤其是野生型和蜡粉缺失突变体中特异性表达的差异表达基因，可作为后续研究青花菜花球表面蜡粉形成分子机制的对象。 相似文献
5.
【目的】对盐胁迫下海马齿根系进行转录组测序分析,挖掘海马齿根系耐盐相关基因,为揭示海马齿耐盐的分子机制提供参考。【方法】利用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）代谢、与盐胁迫相关的钙信号通路和过氧化氢代谢等途径的差异基因上调。【结论】在盐胁迫下海马齿差异表达基因如小分子量热激蛋白基因、抗氧化酶相关基因及与离子交换相关基因发挥了重要调控作用。 相似文献
6.
【目的】分析水稻接种稻瘟病菌（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注释明确了差异基因的分子功能及信号通路,有利于进一步了解植物抗病机制,并可能为稻瘟病防治措施提供新的途径。 相似文献
7.
【目的】鉴定筛选出与卵形鲳鲹卵巢发育相关的候选基因及信号通路,为揭示其卵巢性成熟过程的分子机制打下基础。【方法】挑选卵巢发育处于?期和Ш期的雌性卵形鲳鲹,分别构建卵形鲳鲹卵巢?期和Ш期的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信号通路等与卵形鲳鲹卵巢的发育密切相关,共同调节卵巢的发育与成熟,在卵巢性成熟过程中发挥重要作用。 相似文献
8.
山麻鸭开产期和产蛋高峰期卵巢组织转录组分析 总被引：3，自引：3，他引：0
【目的】探明山麻鸭开产期和产蛋高峰期卵巢组织的转录组差异,丰富蛋鸭转录组数据信息。【方法】选取6份山麻鸭卵巢组织样品（开产期和产蛋高峰期各3份）,利用Illumina HiSeqTM 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生物合成和核糖体生物合成,这两类代谢通路在卵巢生理生殖活动的作用尚不明确。【结论】利用高通量测序技术对山麻鸭开产期和产蛋高峰期卵巢组织的转录组进行测序和分析,揭示了山麻鸭不同生理状态下卵巢组织差异表达基因的数量,获得了差异表达基因的功能、分类和代谢通路。为丰富蛋鸭卵巢组织转录组信息,为开展蛋鸭产蛋性状相关基因的研究及分子调控机制研究奠定基础。 相似文献
9.
【目的】分析低温胁迫下枇杷幼果的转录组,以探讨幼果冻害的分子基础,为深入鉴定枇杷果实生长发育、抗寒性基因和抗寒蛋白提供资料。【方法】以花后75d的"大五星"枇杷幼果为试材,研究低温胁迫(3℃条件下处理12h)下幼果转录组的De novo组装和功能注释,并对其进行生物学信息分析,探索基因表达差异。【结果】低温胁迫下枇杷幼果通过De novo组装产生116 723条contigs,并应用Illumina高能量测序技术得到64 814条unigenes,其中有42 830条被nr、KEGG、COG、GO注释。此外,组装和注释了4 017个枇杷幼果在低温胁迫下表达的差异基因。【结论】发现在枇杷低温耐性中起重要作用的油菜素甾醇生物合成和磷脂酰肌醇信号系统的基因表达上调。 相似文献
10.
遮光性套袋对桃果实转录组的影响 总被引：1，自引：1，他引：0
【目的】探明遮光性套袋在桃果实上的转录组差异,丰富桃转录组数据信息。【方法】选取遮光性套袋与对照的桃果实样品,利用Illumina Hi Seq TM 2500进行高通量测序,构建桃果实转录组文库,并用测序评估、基因功能注释等生物信息学方法进行分析。【结果】经过测序获得16.62 Gb clean data测序数据,且碱基百分比(Q30)大于91%,遮光性套袋和对照2个桃果实样品分别获得65 300 730个reads和66 603 686个reads,分别有85.73%和84.60%的reads与桃参考基因组匹配。以无袋处理为参考,对转录组数据进行比较,遮光性套袋处理共获得1 963个差异表达基因,其中,下调基因708个,上调基因1 255个;在Nr数据库中对差异基因进一步进行注释,注释到1 957个基因,其中,下调基因有705个,上调基因有1 252个;COG功能注释分析发现这些差异表达基因共获得853个功能注释,涉及23个功能类别;在GO功能注释分析中,注释到1 609个基因,可以分为53个功能分类,这些分类主要涉及到分子结合、催化活性、细胞过程、生物调节等诸多生理生化过程;KEGG分析发现共有421个基因被注释到94个代谢通路中,其中,光合作用相关通路、类黄酮生物合成、核糖体生物合成等通路显著富集。光合作用通路和类黄酮生物合成通路在果实着色中发挥了重要作用,而核糖体生物合成代谢通路在果实成熟着色中的作用尚不明确。同时也进行了两组样品的果实品质检测,结果表明,遮光性套袋对桃果实的可溶性固形物及可溶性总糖产生显著性影响,可溶性固形物及可溶性总糖显著降低,而对于果实总酸的影响不大,在果实大小上几乎没有影响。【结论】在遮光性套袋处理状态下,获得一定数量的桃果实差异表达基因,光合作用通路和类黄酮生物合成通路基因在果实着色中发挥重要作用,遮光性套袋对桃果实的可溶性固形物及可溶性总糖产生显著性影响。 相似文献
11.
【目的】掌握不同生长规格墨瑞鳕个体间差异表达基因的表达特点,为其功能相关基因深度挖掘及分子遗传育种提供科学依据。【方法】挑选同一养殖条件下极大个体和极小个体的墨瑞鳕,构建肌肉组织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信号通路及细胞凋亡等能量代谢相关通路上,通过共同协调而对墨瑞鳕的生长发育起调控作用。 相似文献
12.
【目的】探究脊尾白虾（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条,分别为核糖体、碳代谢、氧化磷酸化、氨基酸生物合成、内质网蛋白质加工、糖酵解/糖异生、谷胱甘肽代谢和蛋白输出。【结论】脊尾白虾鳃组织在受低氧胁迫早期通过合成蛋白质及提高代谢能力来抵御低氧环境,随着低氧胁迫时间的延长,物质合成和能量代谢活动均显著下降;但在复氧后随着复氧时间的延长,其蛋白质合成和能量代谢水平又逐渐升高恢复至常氧水平。 相似文献
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.
【目的】 研究福清山羊与努比亚黑山羊发情期卵巢组织转录组差异表达水平,一方面为山羊繁殖性状形成的分子机制提供理论依据,另一方面为利用努比亚黑山羊杂交改良福清山羊提供可加快遗传进展的分子标记。【方法】 利用转录组测序方法对福清山羊和努比亚黑山羊发情期卵巢组织进行研究,筛选品种间的差异表达基因（differentially expressed genes,DEGs）,并对DEGs功能进行注释和若干基因荧光定量PCR（quantitative real-time PCR,qRT-PCR）验证;同时通过与参考基因组比对,分析和筛选测序样品中存在的SNP/InDel。【结果】 6个样品共得到46.68Gb Clean Data,DESeq分析发现了福清山羊和努比亚黑山羊发情期卵巢组织的DEGs149个（包含25个新转录本）,其中表达上调53个,表达下调96个;初步认为输卵管素基因（oviductin,OVN）、类固醇合成快速调节蛋白基因(steroidogenic acute regulatory protein,STAR)、早期生长应答1基因(early growth response 1,EGR1)可作为福清山羊繁殖性能的候选基因。149个DEGs中的108个基因能被GO（gene ontology）数据库注释,30个DEGs能够被COG（Cluster of Orthologous Groups of proteins）数据库注释,91个DEGs能够被KEGG（kyoto encyclopediaof genes and genomes）数据库注释。KEGG通路分析表明DEGs共富集到21条信号通路中,6条通路显著富集。经过进一步的与参考基因组序列比对分析,共发掘1 506个新转录本。经qRT-PCR验证,所选基因（转录本）表达变化模式与转录组测序结果一致,表明测序结果可靠。【结论】 在转录组水平上筛选出了福清山羊和努比亚黑山羊发情期卵巢组织的149个DEGs,发掘了1 506个新转录本,初步揭示了OVN、STAR和EGR1在山羊繁殖过程中可能发挥重要作用,为进一步探索山羊繁殖性状相关机理提供参考依据。 相似文献
16.
【目的】 通过NaCl胁迫下哈茨木霉（Trichoderma harzianum）ACCC32527差异转录组和代谢组数据分析,挖掘关键功能型差异表达基因（DEG）及次级代谢产物。【方法】 采用RNA-seq和GC-TOF-MS技术分别完成0（T1）、0.4（T2）、0.6（T3）mol?L -1 NaCl胁迫下哈茨木霉ACCC32527的差异表达基因和次级代谢产物检测,利用相关软件及数据库（GO、COG、KEGG pathway等）完成对差异表达基因（|log2fold change)|>1 & FDR<0.01）和次级代谢产物（P-value≤0.05 & VIP>1）的筛选、注释和分类,并进行RT-qPCR验证。 【结果】 NaCl胁迫处理后,ACCC32527的生长受到抑制,繁殖速率明显降低,并在该条件下分别获得T1 vs T2和T1 vs T3比较组637、1 570个DEG,获得DEG的16种表达模式,包括9种上调表达模式（950 gene）、3种下调表达模式（662 gene）和4种不规则表达模式（309 gene）,共有的GO功能分类为催化活性、结合、细胞器、细胞膜部分、细胞膜、细胞部分、细胞、单组织过程和代谢过程,且占主要比例,约为61%—94%,其中处理前后基因比例差异较大的分类为催化活性。KEGG代谢通路富集分析发现,不同NaCl浓度处理下的DEG富集到的代谢通路种类及数量存在较大差异,分别有225和535个差异基因富集于20条KEGG通路,包括代谢通路、抗生素的合成和次级代谢产物合成, 其中T2和T3处理下核糖体的生物合成途径分别有12个和18个相关基因受到抑制。从NaCl胁迫下差异转录组中筛选出活性氧清除、离子转运和细胞壁及胞外结构等共73个相关基因,并且多数为上调表达基因。另外,差异代谢组学分析表明,T3处理下,胞内小分子代谢物出现明显变化,累积量增加的代谢产物有30种,包括氨基酸及其衍生物、糖类及其衍生物和醇类,而减少的代谢物有53种,涉及脂肪酸和有机酸等。其中,甘油是已知的真菌重要渗透保护物,T3处理下ACCC32527胞内甘油水平提高,可参与细胞的渗透调节,维持渗透压平衡。RT-qPCR验证了7个差异表达显著基因,虽然在表达量上存在一定的差异,但表达趋势与RNA-seq分析结果基本一致,表明转录组测序结果的可靠性。【结论】 NaCl胁迫引起哈茨木霉ACCC32527大量基因及次级代谢产物变化,获得了与NaCl胁迫调节相关基因和代谢产物,这些机制共同作用减少盐离子对细胞的毒害作用,研究结果可为研究木霉菌的耐盐机理提供重要信息。 相似文献