福清山羊与努比亚黑山羊发情期卵巢组织RNA-Seq分析

【目的】研究福清山羊与努比亚黑山羊发情期卵巢组织转录组差异表达水平,一方面为山羊繁殖性状形成的分子机制提供理论依据,另一方面为利用努比亚黑山羊杂交改良福清山羊提供可加快遗传进展的分子标记。【方法】利用转录组测序方法对福清山羊和努比亚黑山羊发情期卵巢组织进行研究,筛选品种间的差异表达基因(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在山羊繁殖过程中可能发挥重要作用,为进一步探索山羊繁殖性状相关机理提供参考依据。     

转录组测序揭示不同生长速度番鸭胸肌组织的基因表达差异

【目的】利用转录组测序技术分析不同生长速度番鸭的基因表达差异,挖掘影响番鸭生长性能的基因与信号通路,为阐明生长性能差异的遗传机制提供理论基础。【方法】选取不同生长速度两尾样番鸭的胸肌组织进行转录组测序,采用无参考基因组的分析手段对测序结果进行分析,获得差异表达基因（Differentially expressed genes,DEGs）。利用 Blast2 GO 软件对 DEGs 进行功能富集分析,利用 KAAS 网站进行 DEGs 的通路富集分析。【结果】转录组测序分析结果显示,快速生长型番鸭与缓慢生长型番鸭相比共有 186 个显著 DEGs,其中 49 个表达量上调,137 个表达量下调。这些 DEGs 主要富集在磷脂酰肌醇三激酶 - 丝氨酸 / 苏氨酸激酶（PI3K-AKT）信号通路、丝裂原活化蛋白激酶（MAPK）信号通路和腺苷酸激活蛋白激酶（AMPK）信号通路。【结论】PI3K-AKT、MAPK、AMPK 这 3 条信号通路是控制番鸭机体生长发育的重要调节信号通路。     

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

【目的】探明山麻鸭开产期和产蛋高峰期卵巢组织的转录组差异,丰富蛋鸭转录组数据信息。【方法】选取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生物合成和核糖体生物合成,这两类代谢通路在卵巢生理生殖活动的作用尚不明确。【结论】利用高通量测序技术对山麻鸭开产期和产蛋高峰期卵巢组织的转录组进行测序和分析,揭示了山麻鸭不同生理状态下卵巢组织差异表达基因的数量,获得了差异表达基因的功能、分类和代谢通路。为丰富蛋鸭卵巢组织转录组信息,为开展蛋鸭产蛋性状相关基因的研究及分子调控机制研究奠定基础。     

高温胁迫下半夏叶片的转录组分析

【目的】通过分析高温（38℃）胁迫下半夏（Pinellia ternate）叶片转录组数据,挖掘其响应高温胁迫相关的基因,探究高温胁迫影响半夏的分子机制。【方法】通过Illumina HiSeq测序对高温胁迫处理后的半夏叶片进行转录组分析,并用qRT-PCR技术对DEGs进行验证。【结果】转录组测序共获得38.48 Gb Clean data,与对照组（25℃）相比,高温处理后1 138个DEGS上调表达,578个DEGs下调表达。GO富集分析表明,DEGs主要集中在转录因子复合体、内肽酶活性和血红素结合;KEGG富集分析表明,内质网蛋白质加工代谢通路的DEGs最多。基因表达水平分析表明,植物激素信号转导相关基因（PtARP1、PtGRP1、PtABF）、7个谷胱甘肽S转移酶基因PtGST和23个热激蛋白基因PtHSP在高温胁迫下表达量上调。qRT-PCR分析结果表明,9个DEGs在高温胁迫下的差异表达与转录组分析结果一致。【结论】通过对高温胁迫下半夏叶片进行转录组分析,挖掘高温胁迫相关基因,为进一步研究半夏高温胁迫响应机制提供理论参考。     

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

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

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

【目的】筛选鉴定小麦品种连麦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数量最多;富集到光合作用—天线蛋白通路,苯丙氨酸、酪氨酸和色氨酸生物合成,谷...     

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

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

山羊卵泡发育相关基因的筛选及分析

【背景】山羊第一卵泡波中的优势卵泡（dominant follicles, DF）和从属卵泡（subordinate follicles, SF）是整个卵泡发育过程中最为关键的两个阶段。随着卵泡的进一步发育,最终DF可能发育成为成熟卵泡,直到排卵;SF将走向闭锁,其中颗粒细胞的凋亡是导致卵泡发生闭锁的关键因素。然而目前对促进卵泡的优势化或导致其闭锁的分子机理尚不清楚。【目的】通过对山羊第一卵泡波中DF和SF颗粒细胞进行高通量测序,旨在筛选影响卵泡发育的关键基因,为深入探究卵泡发育的调控机制提供理论依据。【方法】选取10只1岁龄健康的贵州白山羊分别注射前列腺素F_2α,使其同期发情,此后每天用B超检测并记录卵泡的生长情况,发情3 d后,统一屠宰并采集第一卵泡波中DF (直径4.5—6 mm)与SF (直径3 —4.5 mm),分别分离其中的颗粒细胞,提取总RNA、构建文库后通过Illumina Hiseq 2500平台进行测序。利用FastQC对测序产出raw reads进行质量评估并经过过滤后,获得品质较高的clean reads;使用Trinity对得到的clean reads进行重新组装,从而获得unigenes;使用CLC Genomics Workbench将unigenes与山羊RefSeq数据库进行比对获得mRNA;使用DESeq2 软件对获得的mRNA进行差异表达分析;分别采用goseq和kobas软件对得到的差异表达基因进行GO分析及KEGG信号通路分析;最终通过qRT-PCR对筛选出的可能影响卵泡发育的关键基因进行验证。【结果】分别对测序得到的raw reads进行过滤后,在DF颗粒细胞中获得43 217 934条clean reads,占raw reads的比例为95.19%;SF颗粒细胞中获得40 766 348条clean reads,占raw reads的比例为95.35%。将得到的unigenes与山羊的RefSeq 数据库进行比对后,共得到33 896条带有注释的转录本,再通过设定FPKM>1, q value<0.05,共在两种卵泡颗粒细胞中获得13 644个基因。设定参数：FPKM≥1,SF-FPKM/DF-FPKM>1,P<0.05,获得695个差异表达mRNA,其中233个在SF颗粒细胞中表达显著上调,462个表达显著下调;对所获得695个差异表达mRNA进行GO功能富集分析,共分为三大类42组：其中生物学过程占47.6%,细胞组分占47.6%,分子功能占4.8%;KEGG信号通路分析,发现20条通路,其中与核糖体通路相关的基因富集最为显著。通过在Genecard中进行功能分析后,筛选6个可能与山羊卵泡发育密切相关的基因,其中PRLR、PTX3、RGN在SF颗粒细胞中表现为上调;DKK3、ALDH1A2、RARRES1则表现为下调。qRT-PCR显示PRLR、RGN、DKK3、ALDH1A2、RARRES1的表达趋势与高通量测序结果一致,且RGN在从属卵泡颗粒细胞中的表达量极显著地高于优势卵泡（P<0.01）;DKK3、ALDH1A2、RARRES1在优势卵泡颗粒细胞中的表达量极显著地高于从属卵泡（P<0.01）。【结论】DKK3、ALDH1A2、RARRES1和RGN在优势卵泡和从属卵泡中表达量存在极显著差异,推测在山羊卵泡发育过程中可能促进卵泡的优势化或导致闭锁,对深入探究卵泡发育的调控机制具有重要意义。     

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

【目的】拟轮枝镰孢(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富集分析显示，侵染早期拟轮枝镰孢在细胞间隙生长，差异基因主要富集在...     

低氧—复氧胁迫下脊尾白虾鳃组织差异表达基因趋势分析

【目的】探究脊尾白虾（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条,分别为核糖体、碳代谢、氧化磷酸化、氨基酸生物合成、内质网蛋白质加工、糖酵解/糖异生、谷胱甘肽代谢和蛋白输出。【结论】脊尾白虾鳃组织在受低氧胁迫早期通过合成蛋白质及提高代谢能力来抵御低氧环境,随着低氧胁迫时间的延长,物质合成和能量代谢活动均显著下降;但在复氧后随着复氧时间的延长,其蛋白质合成和能量代谢水平又逐渐升高恢复至常氧水平。     

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

【目的】对盐胁迫下海马齿根系进行转录组测序分析,挖掘海马齿根系耐盐相关基因,为揭示海马齿耐盐的分子机制提供参考。【方法】利用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）代谢、与盐胁迫相关的钙信号通路和过氧化氢代谢等途径的差异基因上调。【结论】在盐胁迫下海马齿差异表达基因如小分子量热激蛋白基因、抗氧化酶相关基因及与离子交换相关基因发挥了重要调控作用。     

基于转录组测序筛选鸡蛋绿壳性状相关基因

【背景】绿壳蛋深受消费者喜爱，蛋壳绿色的深浅在市场上是影响绿壳蛋定价销售的重要参考指标。绿壳蛋的形成受多基因共同调控，蛋壳绿色深浅不一，其分子机理尚不清楚。通过对赤水乌骨鸡蛋壳腺组织进行转录组测序，挖掘其调控绿壳蛋蛋壳颜色深浅的候选基因以及关键信号通路，探究影响蛋壳颜色遗传性，以期发展绿壳蛋的选种选育和提高经济效益。【目的】探讨赤水乌骨鸡的遗传基础，并通过SLCO1B3基因分型对其进行鉴定和筛选，以期通过分子标记在青壳鸡育种规划中提供新的见解，并在后期的选择策略中帮助控制和提高赤水乌骨鸡蛋壳品质的同质性。【方法】以纯系的280日龄赤水乌骨鸡为研究对象，屠宰产浅绿色蛋（QL）和深绿色蛋（SL）的母鸡各3只，采集蛋壳腺以RNA-seq技术检测分析，筛选与蛋壳颜色密切相关的差异表达基因（differentially expressed genes, DEGs），并对这些DEGs进行GO和KEGG富集分析。利用qRT-PCR技术检测与蛋壳颜色相关的6个候选基因的转录水平变化以验证转录组数据可靠性。【结果】在深绿组与浅绿组中共筛选到93个显著DEGs，有59个基因在深绿组中显著上调，34个显著下...     

‘仓方早生’桃及其早熟芽变不同发育时期果实的转录组分析

【目的】通过对桃品种‘仓方早生’及其早熟芽变不同发育时期的果实进行转录组分析，挖掘参与调控桃果实成熟的关键因子，为深入研究桃果实成熟调控机理提供理论依据。【方法】以桃品种‘仓方早生’及其早熟芽变为试材，每个品种分别选择长势一致的样品树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个。其...     

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信号通路及细胞凋亡等能量代谢相关通路上,通过共同协调而对墨瑞鳕的生长发育起调控作用。     

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

【目的】基于转录组学对西方蜜蜂工蜂不同虫态间的差异表达基因（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种信号通路在西方蜜蜂工蜂幼虫、蛹和成虫的发育调控中发挥着重要作用。