共查询到19条相似文献，搜索用时 154 毫秒
1.
为了探究生物炭添加前后枯草芽孢杆菌应对镉胁迫的分子机制，本研究基于转录组测序，对比分析了枯草芽孢杆菌在4种体系（正常体系、镉胁迫体系、生物炭体系、生物炭-镉胁迫体系）中的基因表达差异。结果表明：在镉胁迫压力下，枯草芽孢杆菌共产生1 932个差异表达基因，其中上调基因967个，下调基因965个。枯草芽孢杆菌的代谢活性受到明显抑制，但是其鞭毛组装和细菌趋化性通路明显增强，离子转运和阳离子跨膜转运蛋白功能条目的基因表达也呈现上调，说明枯草芽孢杆菌通过上调细胞运动及离子转运相关基因的表达来规避镉毒性压力。相比于镉胁迫体系，在生物炭-镉胁迫体系中，细菌共产生1 541个差异表达基因，其中上调基因691个，下调基因850个。生物炭的存在增强了枯草芽孢杆菌生物过程、细胞组分和分子功能基因的表达，枯草芽孢杆菌物质和能源代谢相关通路的基因表达均呈现上调，并且细菌体内镉抗性基因czcD、cadA和膜功能相关差异基因的表达也以上调为主，说明生物炭可以缓解镉胁迫对枯草芽孢杆菌的毒性效应，并可刺激枯草芽孢杆菌的金属抗性能力。 相似文献
2.
【目的】对盐胁迫下海马齿根系进行转录组测序分析,挖掘海马齿根系耐盐相关基因,为揭示海马齿耐盐的分子机制提供参考。【方法】利用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）代谢、与盐胁迫相关的钙信号通路和过氧化氢代谢等途径的差异基因上调。【结论】在盐胁迫下海马齿差异表达基因如小分子量热激蛋白基因、抗氧化酶相关基因及与离子交换相关基因发挥了重要调控作用。 相似文献
3.
【目的】远缘杂交育种是目前牡丹、芍药品种改良和育种的主要方法,而远缘杂交不亲和一直是制约其快速发展的主要因素。本研究从牡丹、芍药远缘杂交授粉后不亲和应答相关的柱头差异蛋白与转录组方面深入研究,揭示牡丹、芍药远缘杂交不亲和的分子机理,为杂交育种提供理论依据。【方法】以芍药‘粉玉奴’自交、芍药‘粉玉奴’与牡丹‘凤丹白’杂交为供试材料,在授粉后24 h采取柱头,分别进行同位素标记相对定量（iTRAQ）和转录组技术分析。对所获得的蛋白和转录组数据进行生物信息学分析,并对其中可能与远缘杂交不亲和相关的基因进行定量PCR验证。【结果】利用iTRAQ技术分析牡丹、芍药远缘杂交后柱头中蛋白质的表达差异,共鉴定到685个差异蛋白,富集到了188条通路,其中显著富集的Pathway有18条。与不亲和授粉相关代谢通路有RNA降解、钙信号途径、丝裂原活化蛋白激酶（mitogen-activated protein kinase signaling pathway,MAPK）信号途径、磷脂酰肌醇信号系统。在RNA降解代谢通路中,烯醇酶（Enolase）、热休克蛋白DnaK（HSP70）及病菌抗原（GroEL）均表达下调。在钙信号途径中,钙调蛋白（CALM）表达下调,腺苷酸转运酶（adenine nucleotide translocase,ANT）表达量增加,表达上调。MAPK信号途径中,乙二醛酶Ⅰ（GloI）表达下调。磷脂酰肌醇信号系统中的钙调蛋白（CALM）表达下调。随机选取与差异蛋白相关的6个基因进行qRT-PCR验证,结果显示,6个基因的表达与蛋白质水平趋势相一致,均表达下调。通过转录组测序,共获得了52 998个有注释信息的Unigene,占所有Unigene的40.37%。基于6组样品的RPKM（Reads Per Kilobase per Million）值,共筛选到16 224个差异基因。其中上调基因13 361,下调基因2 863个。对差异基因进行Pathway显著富集分析,杂交与自交相比,不亲和差异表达的基因主要富集在氧化磷酸化代谢、ABC转运蛋白、次级代谢产物等通路。与远缘杂交不亲和相关且发生显著变化的基因有CalS-5、CalS-12（胼胝质酶）和SPL（squamosa promoter binding protein-like）表达上调,ABCF（ABC transporter family protein）表达下调。【结论】在转录组和蛋白数据共注释到6个蛋白、4个基因与植物不亲和性密切相关,这些蛋白与基因可能在远缘杂交不亲和方面发挥着重要作用。 相似文献
4.
[目的]本文旨在阐明根系分泌物介导的植物根际促生细菌(plant growth-promoting rhizobacteria,PGPR)与宿主的互作机制。[方法]利用Illumina高通量测序技术研究促生解淀粉芽胞杆菌SQR9在体外振荡培养条件下对香蕉根系分泌物的转录响应特征,对显著差异表达基因的直系同源基因进行同源蛋白簇(COG)分析。通过Real-time PCR验证部分转录组测序结果,并研究分泌物代表性组分对部分差异基因表达的影响。[结果]香蕉根系分泌物可显著促进菌株SQR9在1/2MSgg培养基中的生长。转录组测序结果表明:根系分泌物处理6 h后菌株SQR9中共有764个基因的表达发生了显著变化,占其基因总数的18.7%,其中上调与下调基因数量分别为485和279。Real-time PCR验证获得的基因表达差异趋势与转录组测序结果一致,证明其结果可信。菌株SQR9中受根系分泌物显著调控的差异基因主要属于碳水化合物转运与代谢、氨基酸转运与代谢、能量产生与转化、转录调控、信号转导和细胞运动等;此外,与促生/生防活性物质合成相关的部分功能基因也被分泌物所诱导。根系分泌物中的葡萄糖、木糖、赤藓糖醇和硬脂酸可诱导SQR9中部分代表性差异基因的表达,但整体的转录组调控现象是分泌物中不同组分共同作用的结果。[结论]香蕉根系分泌物可促进菌株SQR9的生长代谢、信号响应、根际趋化和活性物质合成,这些都有利于SQR9在香蕉根际的存活定殖和促生/生防功能发挥,有助于其与香蕉合作关系的建立。 相似文献
5.
以野生灵芝[Ganoderma lucidum(Leyss. ex Fr.)Karst.](D)、自然单核灵芝（Mo）、双核灵芝（Mu）的菌丝体为材料,利用高通量测序技术筛选差异表达基因,并进行GO分析和KEGG富集分析。通过差异基因分析,筛选出568个差异表达基因。GO分析结果表明,野生灵芝（D）与自然单核灵芝（Mo）菌丝体的差异表达基因主要富集在转录调控及DNA模板化、tRNA氨酰化、组蛋白甲基化等过程;野生灵芝（D）与双核灵芝（Mu）菌丝体的差异表达基因主要富集在染色质沉默、L-丝氨酸代谢、硫化合物代谢等过程;自然单核灵芝（Mo）与双核灵芝（Mu）菌丝体的差异表达基因主要富集在染色质沉默、DNA模板转录、谷氨酸分解等过程。KEGG富集结果表明,自然单核灵芝（Mo）与双核灵芝（Mu）菌丝体在代谢通路上有明显的差异,差异基因主要被注释到过氧化物酶、维生素代谢等通路。野生灵芝（D）与自然单核灵芝（Mo）菌丝体、野生灵芝（D）与双核灵芝（Mu）菌丝体的差异基因并没有KEGG富集的结果,说明其代谢通路没有明显差异。筛选出3个品系中表达量差异均显著的基因,共12个,并对基因的主要功能进行预测... 相似文献
6.
为探索玉米CMS-S型雄性不育的分子机制,以不育系郑58cms-Q1261和保持系郑58单核晚期的花药为材料进行转录组分析,发现14 565个差异表达基因（DEGs）,其中7 551个上调表达,7 014个下调表达。GO富集分析显示这些DEGs主要富集在代谢过程、细胞和催化活性等方面。KEGG显著富集的主要代谢通路为碳代谢、淀粉和蔗糖代谢、氧化磷酸化、糖酵解/糖异生、氨基糖和核苷酸糖代谢等。选取12个DEGs进行qRT-PCR分析,结果与RNA-Seq数据一致,证明转录组数据的准确性。进一步对DEGs分析,发现2个果胶裂解酶超家族蛋白、5个果胶裂解酶类,3个热休克蛋白和2个bHLH转录因子与雄性不育有关。 相似文献
7.
8.
【目的】确定枯草芽孢杆菌溶血相关基因yqxC的功能。【方法】采用PCR方法,扩增枯草芽孢杆菌溶血相关基因yqxC,构建GST融合表达载体,并在大肠杆菌中进行表达。使用含体积分数5%绵羊血的琼脂平板,检测可溶性融合蛋白GST-YqxC的溶血活性。此外,采用RT-PCR检测yqxC基因在枯草芽孢杆菌中的转录情况。【结果】获得的可溶性融合蛋白GST-YqxC,在绵羊血琼脂平板上显示出较强的溶血活力;通过RT-PCR,检测到yqxC基因可以在枯草芽孢杆菌中转录。【结论】枯草芽孢杆菌yqxC基因编码蛋白YqxC有较强的溶血活性,且能够发生转录,yqxC很可能是引起枯草芽孢杆菌溶血的基因之一。 相似文献
9.
【目的】低温是导致香蕉减产的主要自然灾害之一。基于转录组与蛋白质组关联分析参与香蕉抗寒的相关基因、蛋白及信号、代谢通路调控网络,探讨香蕉抗寒的分子机制。【方法】以巴西蕉为试材,在7℃低温下处理1 d(Cold1)和3 d(Cold3),以28℃培养的巴西蕉为对照（CK）,基于笔者课题组前期获得的蛋白质组数据,利用转录组测序技术检测巴西蕉在低温胁迫下基因调控网络的变化,同时与蛋白质组学数据进行关联分析,共同解析巴西蕉响应低温胁迫的分子机制。【结果】转录组分析结果显示,Cold1 vs CK、Cold3 vs CK和Cold1 vs Cold3三个对比组分别鉴定出11 370、15 460和9 619个差异表达基因,对这些基因的KEGG富集分析发现,差异表达基因在光合作用信号、谷胱甘肽代谢、α-亚麻酸代谢途径和苯丙素生物合成等多个低温胁迫关键信号代谢通路中富集。对部分差异表达基因进行实时荧光定量（qRT-PCR）分析,其中DREB、MAPK和MYB等低温调控关键基因的表达量在低温处理后显著上升,所选20个基因的表达变化趋势与RNA-seq基本相符,证实了RNA-seq的准确性。转录组与蛋白... 相似文献
10.
本文旨在探究金叶榆（Ulmus pumila cv.‘Jinye’）二倍体及其同源四倍体的生理特征和差异基因表达情况,为金叶榆关键调控基因的定位及基因工程育种提供分子基础。试验以金叶榆二倍体及其同源四倍体为材料,比较两者叶绿素相对含量,可溶性糖、可溶性蛋白、MDA含量,SOD和POD活性等生理指标,并进行转录组测序。结果表明：金叶榆同源四倍体叶片中叶绿素相对含量,可溶性糖和可溶性蛋白含量及SOD、POD活性显著高于二倍体,MDA含量显著低于二倍体。转录组测序共得到1 283个差异表达基因,其中1 078个表达上调,205个表达下调。995个差异表达基因被GO功能注释,与代谢过程、细胞过程、催化活性、结合和转运蛋白活性等密切相关;COG注释结果显示,主要注释在碳水化合物的运输和代谢、次生代谢产物合成和脂质运输和代谢;KEGG通路分析发现,植物激素信号转导通路差异显著性最为明显,其中生长素、油菜素内酯和水杨酸信号转导途径富集基因数最多。参与代谢过程、催化活性和激素信号转导通路等关键基因的差异表达,可能与同源四倍体植株能源物质含量高、酶活提高及生长势和抗逆性增强有关。 相似文献
11.
【目的】探究脊尾白虾（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条,分别为核糖体、碳代谢、氧化磷酸化、氨基酸生物合成、内质网蛋白质加工、糖酵解/糖异生、谷胱甘肽代谢和蛋白输出。【结论】脊尾白虾鳃组织在受低氧胁迫早期通过合成蛋白质及提高代谢能力来抵御低氧环境,随着低氧胁迫时间的延长,物质合成和能量代谢活动均显著下降;但在复氧后随着复氧时间的延长,其蛋白质合成和能量代谢水平又逐渐升高恢复至常氧水平。 相似文献
12.
【目的】 研究福清山羊与努比亚黑山羊发情期卵巢组织转录组差异表达水平,一方面为山羊繁殖性状形成的分子机制提供理论依据,另一方面为利用努比亚黑山羊杂交改良福清山羊提供可加快遗传进展的分子标记。【方法】 利用转录组测序方法对福清山羊和努比亚黑山羊发情期卵巢组织进行研究,筛选品种间的差异表达基因（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在山羊繁殖过程中可能发挥重要作用,为进一步探索山羊繁殖性状相关机理提供参考依据。 相似文献
13.
【目的】鉴定筛选出与卵形鲳鲹卵巢发育相关的候选基因及信号通路,为揭示其卵巢性成熟过程的分子机制打下基础。【方法】挑选卵巢发育处于?期和Ш期的雌性卵形鲳鲹,分别构建卵形鲳鲹卵巢?期和Ш期的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信号通路等与卵形鲳鲹卵巢的发育密切相关,共同调节卵巢的发育与成熟,在卵巢性成熟过程中发挥重要作用。 相似文献
14.
【目的】分析西瓜食酸菌(Acidovorax citrulli)在铜胁迫下转录水平的响应特征,研究对其抗铜系统和抗铜机制,为进一步解析西瓜食酸菌的铜耐受分子机制积累数据。【方法】以西瓜食酸菌亚群I的抗铜菌株FC440铜胁迫和无胁迫培养下的菌体为材料,提取RNA并利用RNA-Seq技术获取细菌对铜胁迫的应答数据,分析数据质量、差异表达基因的类型、KEGG富集与Pathway注释、GO富集分析等转录组水平特征,并进行qPCR验证。【结果】共获得6个转录组文库,总数据量23.6 Gb;该细菌共有4 344个基因表达,其中有466个基因差异表达(上调266个,下调200个);差异表达的基因显著富集于ABC转运系统的跨膜转运、双组分系统通路的信号转导以及氨基酸代谢等约9大通路;基因表达趋势与转录组测序数据的趋势一致。【结论】铜胁迫下谷氨酸、谷胱甘肽等氨基酸代谢在西瓜食酸菌抗铜胁迫中发挥重要作用;双组分系统参与西瓜食酸菌对铜离子胁迫的响应且受到负调控;该菌中可能存在包含ABC转运体的跨膜转运、氧化还原反应、胞内束缚沉淀等多种系统互作的抗铜机制。 相似文献
15.
利用RNA-Seq鉴定甘蓝型油菜叶片干旱胁迫应答基因 总被引：6，自引：2，他引：4
【目的】利用RNA Sequencing（RNA-Seq）技术比较2种不同生长条件下甘蓝型油菜苗期叶片转录组,鉴定油菜叶片干旱胁迫应答相关基因,从转录组水平揭示油菜适应干旱胁迫环境的分子机制。【方法】提取正常生长（ZY）和自然失水处理（ZY8D）的六叶期甘蓝型油菜中油821的叶片总RNA,以Illumina Hiseq 2000平台进行RNA-Seq分析。利用NGSQCTookit v2.3.3去除低质量和包含模糊碱基的reads。以甘蓝型油菜亲本物种白菜染色体v1.5和甘蓝Scaffold v1.0为参考序列,采用TopHat2-Cufflinks-Cuffmerge-Cuffdiff标准流程进行差异表达基因（differential expressed genes,DEGs）筛选。对上调和下调DEGs分别采用Cytoscape v3.1.0中的BiNGO和KOBAS2.0进行基因本体（gene ontology,GO）和京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes,KEGG）代谢途径富集分析。选择上调和下调DEGs各3个,以实时荧光定量PCR（quantitative real-time PCR,qRT-PCR）验证RNA-Seq结果的可靠性。【结果】过滤低质量reads后,ZY和ZY8D分别保留了26 192 312和28 378 899对高质量reads用于DEGs筛选,其中86.6%和85.8%的reads能准确比对到参考序列上,说明RNA-Seq结果和参考序列可靠。DEGs鉴定结果表明3 657个基因受干旱胁迫诱导差异表达,其中上调表达基因1 431个,下调表达基因2 226个。GO富集分析发现上调表达基因主要与非生物胁迫响应和化学刺激响应相关,其中,参与水分胁迫响应和脱落酸（abscisic acid,ABA）刺激响应的基因分别有127和141个,而下调表达基因与植物病原菌防御、蛋白激酶活性和水杨酸（salicylic acid,SA）刺激相关。KEGG富集分析表明上调表达基因主要富集于苯丙烷和类胡萝卜素的生物合成及淀粉与蔗糖代谢途径,而下调表达基因主要富集于植物-病原菌互作和植物激素ABA、SA和茉莉酸（jasmonic acid,JA）信号转导途径。qRT-PCR检测6个DEGs的表达模式与RNA-Seq分析结果一致,证实了RNA-Seq结果的可靠性。【结论】RNA-Seq分析鉴定出3 657个甘蓝型油菜叶片干旱胁迫应答基因。GO和KEGG代谢途径分析明确了差异表达基因富集的分子功能与代谢途径。 相似文献
16.
【目的】通过对桃品种‘仓方早生’及其早熟芽变不同发育时期的果实进行转录组分析，挖掘参与调控桃果实成熟的关键因子，为深入研究桃果实成熟调控机理提供理论依据。【方法】以桃品种‘仓方早生’及其早熟芽变为试材，每个品种分别选择长势一致的样品树5株，分别于花后30 d（对应‘仓方早生’c1、早熟芽变y1）、45 d（对应c2、y2）、59 d（对应c3、y3）、71 d（对应c4、y4）及89 d（对应c5）对不同发育时期的桃去皮果肉进行取样和转录组测序，并利用实时荧光定量PCR(qRT-PCR)对筛选的差异表达基因进行定量验证；利用GO和KEGG对‘仓方早生’及其早熟芽变的差异表达基因进行分析；基于差异表达基因构建加权基因共表达网络分析（weighted gene co-expression network analysis,WGCNA），从中鉴定出与果实成熟密切相关的枢纽模块和枢纽基因。【结果】将处于果实相同发育时期的转录组数据进行比较，得到y1与c1、y2与c2、y3与c4和y4与c5四组对比数据，共筛选出差异表达基因4 395个，其中上调表达基因2 212个，下调表达基因2 183个。其... 相似文献
18.
【目的】筛选出调控日本医蛭生长发育的相关基因,揭示其生长发育、繁殖和抗凝分子机制,为后期开展日本医蛭养殖、生长发育及相关功能基因研究提供科学依据。【方法】以幼虫期（孵化10 d以内,HNL）、青年虫期（孵化后2个月左右,HNY）和成年虫期（孵化后2年左右,HNA）的日本医蛭为研究对象,利用Illumina HiSeq 4000平台对构建的cDNA文库进行转录组测序,按错误发现率（FDR）<0.05和|log2 FC|>1的原则筛选差异表达基因（DEGs）,并进行GO功能注释分析和KEGG信号通路富集分析。【结果】使用3个时期的样本构建24个cDNA文库,共产生541124257条原始序列（Raw reads）,经Fastp质量控制及RNA-Seq组装,共获得20430个mRNAs转录本,并注释出88个新的mRNAs转录本。通过对日本医蛭不同生长阶段进行两两对比,共获得1348个差异表达基因,其中,在幼虫期—青年虫期（HNL vs HNY）中获得238个差异表达基因,在幼虫期—成虫期（HNL vs HNA）中获得976个差异表达基因,在青年虫期—成虫期（HNY vs HNA）中获得763个差异表达基因。幼虫期—成虫期的差异表达基因主要富集蛋白水解、肌球蛋白复合体、细胞骨架、几丁质代谢过程及对伤害的反应等GO功能条目上,以及苯丙氨酸/酪氨酸/色氨酸生物合成、丙烯酰胺转移酶等KEGG信号通路上。说明日本医蛭幼虫在生长过程中大量摄食并通过消化、代谢和合成物质为自身发育提供能量,而成年日本医蛭发育成熟后与其生长发育相关的代谢过程下调,且已适应环境,对于环境刺激的敏感度也有所降低。【结论】在幼虫期—成虫期（HNL vs HNA）获得的976个差异表达基因中,上调差异表达基因主要与肌肉发育相关,而与性腺发育或产卵调节及神经发育相关的差异表达基因具有明显周期特异性,其表达变化均与日本医蛭的生长发育有关。水蛭素和抗凝相关基因在青年期高表达,因此开展水蛭抗凝血能力研究时宜选用青年日本医蛭。 相似文献