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

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

肺形侧耳低温胁迫时期的转录组分析

本文以肺形侧耳栽培菌株X57为研究对象,利用高通量RNA测序技术对4℃低温处理0h、6h和12h后的肺形侧耳进行基因表达分析来探讨肺形侧耳低温应答机制。分析结果筛选出低温处理6h差异表达基因742个,其中上调表达基因374个,下调表达基因368个;低温处理12h差异表达基因1 489个,其中上调表达基因占53%,下调表达基因占47%。Gene Ontology（GO）功能聚类分析表明,差异表达基因主要富集在结合、催化分子功能组和代谢过程、细胞过程生物学过程中。KEGG功能富集分析结果显示,差异基因主要富集在氨基酸、核糖体和类固醇生物合成,及氮类物质代谢的通路上,富集到与低温胁迫相关通路MAPK signaling pathway-yeast上的基因表达随低温处理时间的增加呈上调趋势。已报道HOG-MAPK通路是MAPK途径研究较为明确、信号传递单一的真菌低温胁迫中重要的通路,本文运用生物信息学软件构建肺形侧耳的HOG-MAPK通路,并利用荧光定量PCR对相关基因表达进行验证,结果显示以低温处理0h为参照,随着低温胁迫时间增加通路上大部分基因的表达量持续上升,且差异显著,与RNA‐Seq分析结果一致。本文通过全面分析肺形侧耳低温胁迫时期基因表达情况,为进一步研究肺形侧耳低温胁迫相关重要基因的功能鉴定与信号通路调控机理提供基础。     

球囊菌胁迫中华蜜蜂幼虫肠道过程中病原的转录组学研究

【目的】本研究利用RNA-seq技术对球囊菌胁迫的中华蜜蜂(中蜂)幼虫肠道进行深度测序,经趋势分析得到差异表达基因(DEGs)的显著表达模式,进而对胁迫过程中的球囊菌进行转录组学分析。【方法】利用Illumina HiSeq 2500平台对球囊菌胁迫的中蜂幼虫肠道进行深度测序,并利用相关软件进行了深入分析。最后,通过RT-qPCR对RNA-seq数据进行了验证。【结果】本研究共得到球囊菌的41133932条高质量clean reads。22865个DEGs共聚类为8个基因表达模式,其中,16769个DEGs聚类为2个显著上调趋势与2个显著下调趋势。GO富集分析结果显示,显著上调与显著下调趋势中的DEGs分别富集于40与37个GO term,基因富集数最多的为细胞进程(2486 unigenes)。KEGG代谢通路(pathway)富集分析结果显示,显著上调与显著下调趋势中的DEGs分别富集于119和112个pathway,基因富集数最多的分别是氨基酸生物合成(127 unigenes)与核糖体(98 unigenes)。进一步分析表明球囊菌在胁迫中蜂幼虫肠道的过程中通过提高物质合成促进其增殖,而宿主通过抑制球囊菌的蛋白合成抵御病原入侵。富集在MAPK信号通路的11个DEGs的表达水平随着胁迫时间的延长而逐渐下降,推测中蜂幼虫通过抑制该通路而阻遏球囊菌增殖。【结论】本研究不仅为揭示白垩病过程中的球囊菌-中蜂幼虫互作提供了重要信息,也为阐明不同抗性蜂种的球囊菌抗性差异奠定了基础。     

吸水和干燥条件下发菜核糖体代谢差异表达基因分析

为解析发菜吸水和干燥条件下的基因差异表达规律以及发菜适应"干-湿"水分节律调控机制,该研究以充分吸水发菜藻体为对照组,干燥状态下的发菜藻体为处理组,采用高通量Illumina HiSeq PE150测序平台结合生物信息学分析方法对发菜吸水和干燥条件下的差异表达基因进行分析。结果显示:(1)发菜吸水和干燥条件下共鉴定到差异表达基因3 383个,其中显著上调和显著下调表达的基因数分别为1 767个和1 616个。(2)GO功能富集分析发现,在吸水和干燥条件下发菜差异表达基因主要富集在蛋白质合成与代谢等相关途径中;KEGG显著性富集分析发现,吸水和干燥条件下发菜有46个差异表达基因被显著富集到核糖体代谢途径,且均为显著上调表达,这些基因可能参与发菜对干旱胁迫的响应。(3)随机挑选6个被显著富集到核糖体代谢途径的基因进行qRT-PCR分析发现,其基因相对表达量与转录组测序结果一致。研究表明,发菜在干燥条件下可能通过激活核糖体代谢特定基因表达帮助抗旱相关的蛋白质的合成和正确折叠、并参与渗透调节等重要生理活动进而调控其适应干燥环境条件。     

圆口铜鱼亲本与子代肝脏转录组特征及对比分析

通过Illumina HiseqTM 2000测序平台首次开展了圆口铜鱼(Coreius guichenoti)亲本与子代肝脏转录组测序并对比分析其测序结果。对转录组进行拼接和组装, 共获得80688个Unigene, Unigene的长度主要分布在401—600 bp, 占总Unigene的43.56%。与Nr、GO、COG、KEGG等公共数据库比对并进行功能注释, 经分析圆口铜鱼亲本与子代差异表达基因共2701个, 其中上调1240个, 下调1461个。差异基因表达模式聚类热图显示, 同一样本不同重复基因表达相似。将差异基因映射到代谢通路KEGG数据库进行富集分析, 并对37个KEGG显著富集的代谢通路(P<0.05)构建其基因共表达网络, 结果显示丙酸盐代谢、丙酮酸代谢、原核生物固碳途径、乙醛和二羧酸代谢、脂肪酸代谢、萜类骨架生物合成是KEGG差异表达基因的核心代谢途径, 且这些基因在圆口铜鱼子代中表达量均显著上调。该结果丰富了圆口铜鱼基因组数据, 并首次从分子水平比较了圆口铜鱼亲本和子代的代谢差异, 同时解释了在同一循环水系统中子代不易患病的现象。     

低温胁迫下意大利蜜蜂预蛹的差异表达基因趋势分析

【目的】蜜蜂是典型的具有发育狭温性的全变态昆虫。本研究以对低温最敏感的意大利蜜蜂Apis mellifera ligustica预蛹为研究对象,通过低温胁迫不同时间的差异表达基因(differentially expressed genes, DEGs)趋势分析,探讨低温胁迫对蜜蜂发育影响的关键基因。【方法】对3日龄意大利蜜蜂封盖子预蛹进行20℃低温胁迫18 h(T18)和36 h(T36),以未经低温胁迫的预蛹为对照(CK),通过Illumina HiSeq~(TM)平台进行转录组学测定。利用Short Time-series Expression Miner(STEM)软件对2个处理组与对照组比较共有的差异表达基因进行趋势分析,再进一步对显著富集趋势模式中富集的差异表达基因进行GO分类和KEGG pathway分析。利用RT-qPCR对随机挑选的5个DEGs的表达模式进行验证。【结果】对检测到1 062个T18 vs CK和T36 vs CK共有的DEGs进行趋势分析,发现3个显著的基因表达模式,包括2个上调表达模式(Profile 6,有539个基因;Profile 7,有271个基因),1个下调表达模式(Profile 1,有183个基因)。对3个显著富集趋势模式DEGs分别进行GO富集分析和KEGG pathway分析,在Profile 6中找到同时富集在FoxO信号通路和寿命调节信号通路上的胰岛素样肽基因ILP和叉头蛋白O基因FoxO持续上调,说明低温胁迫会影响蜜蜂预蛹蜕皮激素信号传递;在Profile 7中CREB结合蛋白基因CBP持续上调,蜜蜂预蛹受到低温胁迫会影响细胞发育;在Profile 1中细胞色素P450基因CYP450 306a1 (phm)和WNT1显著下调,说明化蛹时受到低温胁迫会影响蜕皮激素的合成。通过RT-qPCR分析挑选的5个基因的表达结果和高通量测序结果一致。【结论】通过对蜜蜂响应低温胁迫的差异表达基因趋势分析发现,蜜蜂胰岛素和蜕皮激素共同调控的FoxO可能是低温胁迫抑制蜜蜂化蛹的一个关键基因。本研究为探索低温胁迫影响蜜蜂发育变态的分子调控机制奠定了基础。     

急性氨氮胁迫下大口黑鲈的肝脏转录组特征分析

为研究大口黑鲈(Micropterus salmoides)受到氨氮胁迫时基因表达的变化规律,文章采用Illumina平台的HiSeq测序策略分析了氨氮胁迫12h和24h后大口黑鲈肝脏的基因表达谱。氨氮胁迫后在大口黑鲈肝脏共获得111.66 Gb的有效数据,组装获得133486个单基因簇(Unigenes), N50为1000 bp。通过比较转录组分析在2个时间点共获得2072个差异表达基因,其中氨氮胁迫12h获得1516个差异表达基因, 907个上调, 609个下调;氨氮胁迫24h获得556个差异表达基因,其中330个上调, 226个下调。GO功能注释分析结果表明,在氨氮胁迫12h差异基因富集数目明显多于24h,但在“氧化还原酶活性”条目中富集的差异基因数目则是氨氮胁迫24h多于12h。此外, KEGG功能富集分析发现,氨氮胁迫对大口黑鲈肝脏氧化应激、细胞自噬和凋亡相关等途径产生影响。选取10个差异表达基因进行qPCR验证,其表达水平与转录组差异基因分析结果一致,证明了测序结果的可靠性。其中,与氧化应激相关基因缺氧诱导因子1α (HIF1α)、生长停滞DNA损伤可诱导蛋白β (Gad...     

盐度胁迫对幼体大海马基因转录表达的影响

为探究大海马(Hippocampus kuda Bleeker)幼体在高盐、低盐胁迫条件下的基因表达水平的变化规律,对实验条件下的大海马幼体的肝脏样品进行了转录组测序。对照组(CK, 25‰)、高盐(HS-test, 31‰)和低盐(LS-test, 17‰)胁迫组共获得71794个单基因簇(Unigenes), N50为1780 bp,平均长度为820.71 bp。高盐胁迫组与对照组比较,共获得2740个差异表达基因(DEGs),其中495个DEGs上调, 2245个DEGs下调;与对照组相比,低盐胁迫组共获得3715个DEGs,其中1854个DEGs上调, 1861个DEGs下调。高/低盐胁迫组DEGs经KEGG数据库富集发现,高/低盐度胁迫均能导致大海马幼体体内氨基酸代谢、免疫代谢、能量和脂肪酸代谢相关基因受到影响。其中,低盐胁迫时能量代谢和氨基酸代谢的相关基因显著上调,高盐胁迫时脂肪酸代谢的相关基因显著下调,而高/低盐胁迫时免疫代谢的相关基因都显著上调。从经过盐度胁迫的大海马幼体的肝脏转录组中分别筛选到免疫相关基因Gst、Hsp70、Hsp90、Sod、Bcl-2、Gadd45...     

准噶尔小胸鳖甲短时低温胁迫响应的转录组分析

【目的】拟步甲科昆虫小胸鳖甲Microdera punctipennis是分布于中国新疆古尔班通古特沙漠的特有物种,具有很强的耐寒性,但其低温响应的分子机制尚不明确。本研究旨在利用转录组测序技术丰富小胸鳖甲的已知基因信息,分析在短时低温胁迫下,上调表达基因的主要类群及参与的主要代谢通路。【方法】利用Trinity软件对分别在4℃低温和25℃(对照)下处理3 h的小胸鳖甲成虫RNA-Seq数据进行从头组装,利用Blast2go软件进行基因注释。通过DESeq和GFOLD软件分析差异表达基因,并对上调表达基因进行GO(Gene Ontology)和KEGG(Kyoto Encyclopedia of Genes and Genomes)代谢途径富集分析。【结果】测序过滤后得到68 165 001个有效读序,51 712个平均长度为984 bp的非冗余基因,其中29 925个基因(约57.87%)有同源序列(E-value10-5),与拟步甲科的赤拟谷盗Tribolium castaneum相应基因相似性最高(核苷酸序列一致性为72.75%)。低温响应上调表达基因有514个,富集到35个GO类群,18个KEGG途径。其中胁迫响应类群、生物调节类群和免疫系统过程类群都占有重要比例,嘌呤代谢、硫胺素代谢、糖酵解/糖异生等代谢通路显著富集。在胁迫响应类群中,热激蛋白Hsp90基因、超氧化物歧化酶(SOD)基因和钙联接蛋白Calnexin基因等8个非生物胁迫相关基因显著上调表达。【结论】荒漠昆虫小胸鳖甲低温转录组揭示,在4℃冷驯化过程中,代谢过程、胁迫响应、生物调节和免疫系统等生物学过程相关基因表达显著上调,其中几个非生物胁迫响应基因提示小胸鳖甲可能从分子伴侣、细胞周期阻滞、线粒体稳定性、抗氧化胁迫等多方面应对低温胁迫。KEGG富集分析所揭示的小胸鳖甲在低温下的转录组代谢通路与其他物种有较大差异。研究结果为后续生物信息学分析及小胸鳖甲耐寒关键基因的发掘及耐寒机制的揭示提供了基础数据。     

NaCl胁迫下黑果枸杞转录组测序分析

研究黑果枸杞在不同浓度盐胁迫下基因表达谱变化情况,为进一步研究黑果枸杞抗盐分子机制奠定研究基础。对0(CK)、50、250 mmol/L NaCl溶液胁迫的黑果枸杞组培苗的根和叶在胁迫时间为0、1、12 h时分别取样,采用转录组测序(RNASeq)技术进行测序分析。结果表明,转录组测序共产生222.49 Gb原始数据,拼接出Unigenes 86 037条,注释到7大功能数据库(GO、KEGG、KOG、NR、Pfam、Swiss-Prot和egg NOG)上的Unigenes总数为46 594个,占总Unigenes的54.76%,还有38 929个Unigenes在这些数据库中没有得到注释。通过GO分类和KEGG Pathway富集性分析,分别归于51个GO类别和211条代谢途径。差异表达基因分析显示,黑果枸杞叶片和根的上调基因和下调基因数随着NaCl浓度的增大和处理时间的延长均呈增加趋势,叶片中的上调基因数(7 514)小于下调基因数(9 032),根中的上调基因数(12 347)大于下调基因数(11 559)。在黑果枸杞盐胁迫下转录组中发现28 325个SSR位点,最多的为单核苷酸SSR,占70.47%。综合分析表明,黑果枸杞对盐胁迫的反应是一个多基因参与、多个生物过程协同调控的过程,基因表达量的变化可能是基因调控的主要方式。     

基于RNA-Seq分析MpigE在红曲霉色素生物合成中的转录调控

[目的]分析MpigE的缺失在转录水平对红曲色素产生的影响.[方法]对实验室保存的野生型紫色红曲霉(Monascus purpureus) Mp-21和△MpigE菌株进行高通量转录组测序、注释、表达差异基因功能富集分析和基因通路富集分析,在转录水平揭示MpigE缺失后红曲霉色素产量变化的原因.[结果]通过RNA-Se...     

不同成熟度树葡萄叶片中类黄酮合成转录组基因分析

为了解树葡萄(Myrciaria cauliflora)类黄酮合成相关酶差异表达基因信息,对其幼叶和成熟叶进行全转录组测序并比较分析。结果表明,从树葡萄幼叶和成熟叶中获得59 321条单基因簇(Unigenes),在8大数据库共注释到32 912条Unigenes信息,其中类黄酮合成代谢相关酶基因77个,在成熟叶片中显著下调表达的基因6个,包括2个CHI、1个CHS、1个F3H、1个2-羟基异黄酮脱水酶基因和1个ANS。5个差异表达基因经qRT-PCR验证的结果与转录组测序结果相符合。因此,树葡萄叶片中含有大量不同种类黄酮合成代谢相关酶家族基因,成熟叶片中类黄酮含量显著减少是由于少量相关基因显著下调。     

基于油桐种子3个不同发育时期转录组的油脂合成代谢途径分析

油桐是我国重要的木本油料植物, 过去对油桐的研究主要集中于栽培和常规育种, 与油桐种仁油脂合成相关的分子机理研究还未见报道。文章采用RNA-seq技术对油桐种子油脂合成的3个不同时期的转录组进行比较, 获得了大量差异表达的Unigene序列。在此基础上, 通过GO分类和Pathway富集性分析将这些差异表达Unigene归类于128个代谢途径, 其中包含与油脂合成相关的脂肪酸生物合成和甘油磷脂代谢途径。桐酸经脂肪酸生物合成途径合成后通过甘油磷脂代谢途径以桐油的形式贮存。将这两个代谢途径的Unigene序列在KEGG数据库中进行比对, 获得了一些关键酶的同源蛋白质。文章通过对编码这些同源蛋白质的基因在油桐种子油脂合成期的表达模式进行分析, 以期为油桐油脂合成, 尤其是桐酸合成机理的解析提供理论参考, 并为油桐的遗传改良提供潜在的基因资源, 从而提高油桐的单位面积产量。     

转录组测序揭示油桐脂肪酸的合成途径

油桐是我国重要的木本油料植物,过去对油桐的研究主要集中于栽培和常规育种,与油桐种仁油脂合成相关的分子机理研究还未见报道。本研究采用转录组测序（RNA-Seq）技术,比较了油桐种子油脂合成3个不同时期的转录组,获得了大量差异表达的Unigene序列。通过GO（GeneOntology）分类和代谢途径富集性分析将这些差异表达的Unigene归类于包含脂肪酸生物合成途径在内的128个代谢途径。随后将脂肪酸生物合成途径中的54个Unigene序列在KEGG（KyotoEncyclopediaofGenesandGenomes）数据库中进行比对,获得了14个关键酶的同源蛋白质。本研究通过对编码这些同源蛋白质的基因在油桐种子油脂合成期的表达模式进行分析,以期为油桐油脂合成,尤其是桐酸合成机理的解析提供理论参考,并为进一步的理论研究和油桐的遗传改良提供潜在的基因资源,从而提高油桐的单位面积产量。     

不同温度条件下鲁氏耶尔森氏菌的链特异性转录组分析

为鉴定鱼源鲁氏耶尔森氏菌(Yersinia ruckeri) SC09菌株水生环境中不同温度的转录组水平上的差异, 研究采用链特异性转录组测序(Strand-specific RNA-seq)技术对菌体生理温度(28℃)和实验培养温度(37℃)下进行链特异性测序, 原始数据质控后, 筛选得到差异表达基因, 通过KEGG (Kyoto Encyclopedia of Genes and Genomes)数据库对差异表达基因进行富集分析, 并利用Rockhopper软件筛选出的重要原核生物基因簇进行验证。结果显示, 共筛选获得173个显著差异表达基因(P<0.05), 其中包括58个上调基因, 主要富集到一些特殊的碳水化合物代谢相关的通路中; 以及115个下调基因, 主要富集到双组份信号系统中与三羧酸循环相关的代谢通路上, 同时部分基因富集到编码鞭毛素相关的基因簇中。结果表明, 相对于37℃的实验室培养温度, 在水生环境的生理温度条件下(28℃) SC09菌株拥有较高的运动性和较强的葡萄糖代谢, 但相对的SC09菌株代谢一些特殊糖类的能力减弱。     

基于转录组的八角挥发油合成相关基因挖掘与分析

为更好地挖掘八角(Illicium verum)挥发油合成相关基因，该文对挥发油性状差异显著的优良无性系桂角69号及普通品种砧01号叶片进行了转录组测序及组装注释，并对差异表达基因进行了GO分类和KEGG通路分析。结果表明：(1)转录本经组装后获得84 182条序列，使用NR、NT、Swiss-Prot、KEGG、KOG、GO和Pfam数据库进行序列比对，共注释了59 161条序列，筛选出30 572个差异表达基因。与砧01号相比，桂角69号叶片中上调基因有15 025个，下调基因有15 547个。(2)GO分类结果显示共有20 287个差异基因被注释。KEGG分析结果表明，有21 600个差异基因被注释到133条KEGG通路上，其中挥发油合成相关的单萜生物合成通路、萜类骨架生物合成通路、苯丙素合成通路中的芳樟醇合酶、月桂烯合酶、香叶基香叶基焦磷酸合酶、肉桂酰辅酶A还原酶、咖啡酸3-O-甲基转移酶、肉桂醇脱氢酶等关键酶基因呈差异表达。(3)转录因子分析发现差异表达基因分布于31个转录因子家族，其中MYB家族序列数量最多。该文利用转录组测序技术分析八角优良无性系与普通品种叶片的差异基因及...     

温度刺激对斑马鱼仔鱼基因转录表达的影响

许多优良鱼类养殖品种不耐低温或高温的特点给水产养殖业带来诸多限制和困难,这些鱼类在胚胎和仔鱼等早期阶段的抗寒和抗热能力比成体更差,育苗过程中很容易受到温度突然变化的影响。虽然目前利用基因芯片技术已研究了温度刺激对几种鱼类成体组织中基因表达的影响,但温度刺激对仔鱼基因转录表达的影响还未见报道。研究以斑马鱼受精后96h的出膜仔鱼为实验材料,分别在低温(16℃)和高温(34℃)条件下处理12h和24h,用基因芯片技术检测温度刺激对其基因表达的影响。与培养在28℃的对照相比,低温和高温处理后共有3633个基因发生差异表达,其中低温处理后差异表达基因数目多于高温处理,而且低温抑制基因数目多于诱导表达基因的数目。生物信息学分析结果表明,低温诱导基因主要参与RNA加工和核糖体生物发生等生物学过程,高温诱导基因则主要参与应激反应和未折叠蛋白结合。低温抑制基因主要参与蛋白质水解、视觉感知以及铁离子结合等生物学功能,高温抑制基因参与的生物学功能包括DNA复制、神经系统过程和类固醇激素生物合成等。除了已报道的温度刺激响应基因外,研究鉴定出了大量尚未报道与温度刺激相关的基因,如参与RNA加工的rnmtl1a和pus3基因,以及参与转录调控的twistnb和aebp2基因等。研究结果为进一步揭示鱼类冷或热适应的分子机理和培养耐寒或耐热的养殖新品种提供理论基础。