Bmo-miR-375-3p对家蚕丝素蛋白轻链基因BmFib-L表达的调控作用

【目的】探究家蚕Bombyx mori miRNAs对丝素蛋白轻链基因(fibroin light chain gene,Bm FibL)和P25蛋白基因Bm P25表达的调控作用。【方法】分别以Bm Fib-L和Bm P25 mRNA 3'UTR为靶序列,应用在线预测软件RNAhybrid从前期家蚕4和5龄幼虫后部丝腺(posterior silk gland)sRNA高通量测序获得的29个差异表达bmo-miRNAs中,筛选对Bm Fib-L和Bm P25具有调控作用的候选bmo-miRNAs;采用半定量RT-PCR方法在mRNA水平分析候选bmo-miRNAs及其靶基因Bm Fib-L和Bm P25在4和5龄幼虫期以及5龄第3天幼虫不同组织的表达水平;利用双荧光报告基因检测系统在家蚕细胞Bm N中验证候选bmo-miRNAs对Bm Fib-L和Bm P25表达的调控作用。【结果】筛选到一个在Bm Fib-L和Bm P25 mRNA 3'UTR上都有靶位点的bmo-miR-375-3p(曾称bmo-miR-375*)。在后部丝腺中bmo-miR-375-3p和其预测靶基因Bm Fib-L和Bm P25都有较高的表达水平,提示bmo-miR-375-3p具备调控Bm Fib-L和Bm P25表达的时空条件。miR-375-3p重组表达载体与融合了Bm Fib-L 3'UTR的萤火虫荧光素酶(luciferase)报告基因(luc)表达载体共转染Bm N细胞,报告基因luc的表达水平显著下降;而在miR-375-3p重组表达载体与融合了Bm P25 3'UTR的luc表达载体共转染Bm N细胞中报告基因luc的表达水平下降不显著。【结论】miR-375-3p显著下调Bm Fib-L基因的表达,而对Bm P25基因的表达无明显调控作用。     

家蚕磷酸吡哆醇氧化酶基因的表达谱分析

【目的】了解家蚕Bombyx mori维生素B6关键代谢酶磷酸吡哆醇氧化酶（pyridoxine- 5′-phosphate oxidase, PNPO）基因在家蚕不同发育阶段及5龄幼虫不同组织中的表达差异。【方法】将家蚕PNPO基因的重组表达质粒pET-22b(+)-PNPO转化入大肠杆菌Escherichia coli Rosetta中诱导表达, 纯化蛋白制备多克隆抗体。分别采用荧光定量PCR和Western blot方法对家蚕PNPO基因进行了转录水平和翻译水平的表达分析。【结果】在家蚕发育水平上, 5龄幼虫的PNPO翻译量为最高。PNPO基因在5龄幼虫各组织中的转录水平由高到低依次为精巢、 头、 中肠、 马氏管、 卵巢、 表皮、 脂肪体、 丝腺; 翻译量也以精巢为最高, 其次是头、 中肠和马氏管。【结论】明确了PNPO在家蚕各发育阶段及5龄幼虫各组织中的表达情况。     

miRNA响应印度梨形孢定殖调控大麦生长发育的作用机制

【背景】内生真菌印度梨形孢(Piriformospora indica)定殖植物可以显著促进植物生长发育。miRNA已被证实在植物体的生长发育中具有调控作用。【目的】揭示印度梨形孢定殖大麦促进大麦生长发育过程中miRNA对印度梨形孢定殖的响应及对大麦生长发育的调控作用。【方法】提取大麦总RNA,实施转录组测序并进行序列比对与数据挖掘;使用高效液相色谱检测大麦生长素等激素水平变化。【结果】印度梨形孢对大麦有显著促生作用;全转录组测序结果显示：印度梨形孢侵染3 d较空白对照有18个差异表达的miRNA,其中11个miRNA上调、7个miRNA下调;侵染7 d与空白对照相比24个差异表达的miRNA,其中11个miRNA上调、13个miRNA下调;侵染3 d与侵染7 d相比有3个miRNA上调、6个miRNA下调。GO功能富集分析与KEGG通路分析显示,差异表达miRNA的靶基因主要参与转录、细胞分裂、生长素信号的感知和转导、光合作用和激素刺激响应。靶基因所参与的途径与大麦生长发育密切相关,暗示miRNA对印度梨形孢定殖过程做出了积极响应。代谢产物分析表明miRNA参与的调控路径的代谢产物发生改变。【结论】本研究以miRNA为入手点,探究了miRNA对大麦生长发育的调控机制,为揭示印度梨形孢的促生机制提供了新的研究方向。     

棉铃虫雌雄性腺转录组的比较分析

【目的】棉铃虫Helicoverpa armigera(Hübner)是我国重要的农业害虫,生殖力强是其发生为害的内在生物学基础,而已知的与棉铃虫生殖相关的基因信息相对较少,本研究旨在筛选棉铃虫性腺发育的相关基因。【方法】通过高通量测序技术对棉铃虫成虫的卵巢和精巢进行转录组测序。【结果】共获得100 603条unigenes,平均长度为666.05 bp,其中N50为1 114 bp。经同源性比对,52 071条unigenes获得注释信息,其中注释到Nr数据库的序列最多。通过比较转录组分析发现,在棉铃虫精巢和卵巢中存在7 714个差异表达的基因（Differentially expressed genes,DEGs）,其中3 288个DEGs表达水平在卵巢中上调,4 426个DEGs在精巢中上调。差异基因富集结果涉及到生殖系统发育的相关通路有mTOR信号通路、卵母细胞减数分裂、促性腺激素释放激素信号通路和胰岛素信号通路等。同时筛选得到部分与性腺发育相关的基因,包括vitellogenin(Vg)、vitellogeninreceptor(Vg R)、chorion-relatedgene、testis-specific serine/threonine-protein kinase和spermatogenesis-associated protein等。选取其中12个DEGs用实时荧光定量检测其表达量,结果表明RT-qPCR与RNA-Seq的结果一致。【结论】本研究获得了棉铃虫雌雄性腺的转录组数据及主要性腺发育相关基因,将有助于丰富棉铃虫繁殖相关的基因资源,为进一步研究棉铃虫生殖调控机理提供基础数据。     

吡哆醛激酶基因在家蚕体内的时空表达特征

【目的】了解家蚕Bombyx mori维生素B₆关键代谢酶吡哆醛激酶（pyridoxal kinase, PLK）在家蚕不同组织间的表达差异。【方法】原核表达家蚕重组PLK, 获得目的蛋白制备多克隆抗体, 利用Western blot方法对PLK在家蚕不同发育时期、 5龄第3天幼虫不同组织及5龄不同日龄幼虫表皮、 头部和马氏管中的表达进行分析; 并采用实时荧光定量PCR的方法对PLK基因在家蚕不同组织中的mRNA转录水平进行比较。【结果】PLK在家蚕5龄幼虫的表达水平最高, 在5龄第3天幼虫各组织中的表达由高到低依次为： 精巢、 马氏管、 表皮、 中肠、 头部、 卵巢、 脂肪体、 丝腺; 5龄期不同日龄幼虫表皮组织的表达差异最大, 头部和马氏管中均为前期表达量稍高于后期。在转录水平方面, 精巢的表达量最高, 其次是马氏管和头部。【结论】PLK在家蚕不同发育时期及组织中的差异表达进一步证实PLK在家蚕VB6代谢中的生物学功能的重要性。     

家蚕miRNA Novel-31~*上调溶血素基因的表达

【目的】Novel-31*是在家蚕质型多角体病毒(Bombyx mori cytoplasmic polyhedrosis virus,Bm CPV)感染的家蚕中发现的一个差异表达miRNA。本研究旨在验证Novel-31*对其靶基因表达的调控作用,以便进一步研究miRNA及其靶基因在昆虫免疫调节中的作用。【方法】用生物信息学方法预测Novel-31*的靶基因,荧光定量PCR分析Novel-31*及其靶基因在家蚕感染Bm CPV后不同时间点的表达变化;构建miRNA慢病毒表达载体和靶基因慢病毒表达载体,转染293T细胞,同时合成Novel-31*mimics转染家蚕培养细胞Bm N,使用荧光定量PCR检测Novel-31*对靶基因表达的调控作用。【结果】生物信息学方法预测发现,溶血素基因是Novel-31*的靶基因,其结合位点位于溶血素基因的5'UTR区域。荧光定量PCR分析表明,Novel-31*及溶血素基因在感染Bm CPV的家蚕血淋巴细胞中呈现明显的上调表达。荧光定量PCR检测表明,在Novel-31*慢病毒表达载体和溶血素基因5'UTR慢病毒表达载体转染的293T细胞中和在转染Novel-31*mimics的家蚕Bm N细胞中,溶血素基因都上调表达。【结论】溶血素基因是miRNA Novel-31*的靶基因,Novel-31*与溶血素基因5'UTR结合,上调溶血素基因的表达。     

基于转录组测序技术进行静原鸡差异miRNA与mRNA的关联分析

为探究静原鸡肌肉组织肌苷酸沉积过程的分子调控机制,本研究以静原鸡胸肌和腿肌组织作为试验材料,通过转录组测序技术,利用生物信息学分析方法筛选出差异表达miRNA及其靶基因.分析结果表明,静原鸡胸肌和腿肌组织的比较组合中有39个差异表达miRNAs(19个显著上调,20个显著下调).通过miRNA-mRNA互作网络分析可知,一个miRNA可以靶向多个mRNA参与调控多个靶基因的表达.GO富集分析表明,候选靶基因的功能主要富集于三个分支的单一生物过程、细胞内和结构分子活性等.KEGG通路注释表明,候选靶基因显著富集到碳代谢、蛋白酶体通路和氨基酸的生物合成等通路中.选取4个差异表达miRNA的与肌苷酸合成和代谢相关的候选靶基因(POLR2C,GMPR,IMPDH2和APRT)进行实时荧光定量PCR验证,结果表明,定量结果与转录组测序结果趋势一致.该研究结果为阐明地方鸡种肌苷酸特异性沉积过程中miRNA介导的靶基因与肉品风味的关系及其表达分析提供理论依据.     

刺猎蝽属种间转录组比较分析

【目的】建立刺猎蝽属Sclomina 3种的转录组数据库,分析近缘种间转录组表达差异,探讨在转录组水平的种间趋异情况。【方法】采用Illumina HiSeq^TM4000高通量测序平台进行刺猎蝽属齿缘刺猎蝽S.erinacea、兴仁刺猎蝽S.xingrensis和广西刺猎蝽S.guangxiensis 3种转录组测序;利用Nr, Swiss-Prot, COG/KOG和KEGG数据库进行基因功能注释;对种间差异表达基因(differentially expressed genes, DEGs)进行GO功能注释和KEGG通路富集分析。【结果】齿缘刺猎蝽、兴仁刺猎蝽和广西刺猎蝽转录组测序组装获得42 215个unigenes。21 117个基因在上述4个数据库中得到注释(各数据库注释基因数为：Nr, 20 522; Swiss-Prot, 15 550; COG/KOG, 13 969; KEGG,10 850)。兴仁刺猎蝽vs齿缘刺猎蝽转录组有3 390个DEGs (803个上调,2 587个下调),兴仁刺猎蝽vs广西刺猎蝽转录组有12 543个DEGs (6 63...     

草地贪夜蛾雄性成虫和5龄幼虫的转录组比较分析

【目的】草地贪夜蛾Spodoptera frugiperda是一种新近入侵我国的重要害虫。本研究旨在对草地贪夜蛾雄性成虫和5龄幼虫两个不同发育阶段的转录组进行比较分析。【方法】利用高通量测序技术对草地贪夜蛾雄性成虫和5龄幼虫进行转录组测序和数据组装,并对转录组数据进行功能注释和比较分析。【结果】经de novo组装共获得209 002条转录本,平均长度为687.55 bp,N50为982 bp。共有46 198条(57.43%) unigene在至少一个数据库中获得功能注释,其中1 713条(2.13%) unigene在所有数据库中均能获得注释。在GO数据库中获得205 269条unigene的注释,主要包括68个功能分类;在KEGG数据库中共有3 408条unigene得到注释,涉及277个代谢通路。共鉴定到424个嗅觉相关的基因,并且在雄性成虫和5龄幼虫之间的表达存在差异。通过比较转录组分析,在雄性成虫中鉴定到9 162个上调和6 399个下调差异表达基因(DEGs);功能富集分析发现在上调DEGs中涉及信息素以及信号转导的代谢通路显著富集,而下调DEGs中涉及解毒相关的通路显著富集。【结论】这些转录组数据为探究草地贪夜蛾的生长发育、嗅觉相关功能基因以及候选分子靶标提供资源信息。     

中华蜜蜂幼虫肠道响应球囊菌胁迫的microRNA应答分析

【目的】蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)是一种能够侵染中华蜜蜂(Apis cerana cerana,简称中蜂)幼虫的致死性真菌病原。微小RNA(microRNA,miRNA)可通过在转录后水平靶向抑制或降解mRNA而参与宿主与病原互作过程。本研究旨在对球囊菌胁迫的中蜂6日龄幼虫肠道的差异表达miRNA(DEmiRNA)及其靶基因进行深入分析,进而揭示DEmiRNA在中蜂响应球囊菌胁迫应答过程中的作用。【方法】利用Illumina MiSeq平台对正常及球囊菌胁迫的中蜂6日龄幼虫肠道(AcCK和AcT)进行测序,通过相关生物信息学软件预测DEmiRNA及其靶基因。通过Blast将靶基因注释到GO和KEGG数据库。利用Cytoscape软件构建DEmiRNA与其靶mRNA的调控网络。通过Stem-loop RT-PCR和qPCR验证测序数据的可靠性。【结果】本研究共预测出537个miRNA,其长度分布介于16–35 nt之间,且不同长度的miRNA首位碱基偏向性差异明显。通过Stem-loop RT-PCR证实了10个novel miRNA的表达。AcCK vs AcT比较组共有54个DEmiRNA,包含31个上调和23个下调miRNA,可分别靶向结合6170和8199个靶基因。GO分类结果显示上调和下调miRNA的靶基因分别涉及47和47个条目,富集基因数最多的皆为结合细胞进程和催化活性。KEGG代谢通路(pathway)富集分析结果表明上调和下调miRNA的靶基因分别富集在134和126条pathway,富集基因数最多的均为内吞作用和内质网中的蛋白质加工。调控网络分析结果表明,DEmiRNA及其靶mRNA形成十分复杂的调控关系;31个DEmiRNA可靶向结合51个与泛素介导的蛋白水解相关的mRNA,18个DEmiRNA可靶向结合14个与Jak-STAT信号通路相关的mRNA;miR-1277-x、miR-26-x、miR-27-y、miR-30-x、miR-6052-x等16个miRNA共同参与了上述两条免疫通路的调控。最后,随机挑选3个DEmiRNA进行qPCR验证,结果证明了测序数据的可靠性。【结论】本研究提供了中蜂幼虫肠道在球囊菌胁迫后期的miRNA的表达谱和差异表达信息,揭示了球囊菌与宿主之间在miRNA组学水平存在复杂的互作。miR-6052-x和miR-1277-x作为调控网络的核心可能通过影响细胞凋亡参与宿主的免疫防御,miR-26-x和miR-30-x可能通过调控Jak-STAT信号通路参与宿主的胁迫应答。本研究筛选出的关键DEmiRNA有望作为治疗白垩病的分子靶标。     

中华蜜蜂精巢和卵巢差异表达基因分析

[目的]在蜂群中,雄蜂与蜂王都有着发育完全的性腺,但两者达到性成熟的时间却是不同步的.本研究旨在探究中华蜜蜂Apis cerana cerana雄蜂与蜂王生殖腺的基因表达差异.[方法]利用Illumina测序技术对中华蜜蜂雄蜂精巢与蜂王卵巢转录组差异表达基因(differentially expressed genes...     

家蚕miR-2769靶基因的鉴定及表达分析

【目的】MiRNAs在昆虫变态发育过程中发挥非常重要的作用。对家蚕Bombyx mori miRNAs及靶基因的研究将有助于阐明miRNAs参与调控家蚕变态发育的分子机制。【方法】往家蚕5龄第2天幼虫血淋巴注射蜕皮激素20E后,qRT-PCR检测miR-2769在家蚕脂肪体中的表达;通过生物信息学方法预测家蚕miR-2769的靶基因;利用双荧光酶报告载体系统分析miR-2769与预测靶基因BmE75B的互作;qRT-PCR检测miR-2769及其靶基因BmE75不同剪接体在家蚕不同发育时期(幼虫、蛹和成虫)和幼虫不同组织(头、表皮、丝腺、脂肪体、精巢、卵巢、马氏管、中肠和血淋巴)中的表达量。【结果】研究结果表明,miR-2769可通过与家蚕BmE75B的3′UTR区结合位点的互作,显著抑制荧光素酶报告基因的表达。qRT-PCR结果表明,miR-2769和BmE75A/BmE75B在20E诱导家蚕脂肪体中表达趋势相反。时空表达分析结果表明,miR-2769与BmE75的不同剪接体在家蚕不同发育时期和不同组织中均具有特异性表达特征。在家蚕变态发育的不同阶段,miR-2769和BmE75A的表达量均较低,而BmE75B在蛹期表达量较高,BmE75C在5龄末幼虫和蛹早期有极高表达水平。此外,miR-2769与BmE75不同剪接体均存在一定程度的表达负相关。在家蚕5龄幼虫血淋巴中miR-2769可促进BmE75A的表达,在脂肪体中miR-2769可促进BmE75B的表达,而在其他组织中miR-2769抑制BmE75不同剪接体的表达。【结论】家蚕miR-2769可通过与BmE75B的3′UTR区的互作对BmE75不同剪接体的表达进行负调控。     

家蚕BmlncR2036调控P25基因启动子活性及参与中肠发育调控

【目的】长链非编码RNA(long non-coding RNA, lncRNA)对家蚕Bombyx mori发育具有重要调控作用。我们在前期研究中发现一个位于家蚕丝素蛋白基因P25附近的lncRNA BmlncR2036。本研究旨在进一步探索BmlncR2036调控家蚕P25基因表达的分子机制。【方法】qPCR检测BmlncR2036在5龄第3天家蚕幼虫不同组织(体壁、脑、神经、精巢、卵巢、丝腺、马氏管、血淋巴、脂肪体和中肠)中和不同发育阶段前丝腺、中丝腺和后丝腺中的表达谱。预测能够同时靶向P25和BmlncR2036的miRNA,利用荧光素酶检测法验证miRNA与P25互作关系。荧光素酶检测法测定BmlncR2036对P25基因启动子转录活性的影响。在家蚕5龄第3天幼虫中注射dsRNA敲低BmlncR2036的表达,观察丝腺及中肠表型的变化。【结果】qPCR结果显示,BmlncR2036在家蚕5龄第3天幼虫和5龄熟蚕的后丝腺中高表达,与P25基因呈现一致的表达趋势;BmlncR2036在家蚕5龄第3天幼虫的精巢和中肠中高表达,在其他组织中表达量较低,暗示其功能的多样性。miR-2739和miR-279a既能够与BmlncR2036匹配,也可能靶向P25基因的3′UTR。但荧光素酶检测法结果显示P25不是miR-2739和miR-279a的真实靶标。BmlncR2036负调控P25启动子活性,共转染pcDNA3.1(+)［BmlncR2036］和pGL3-Enhancer［P25-promoter］载体后,细胞荧光素酶活性极显著下降了52%。在家蚕5龄第3天幼虫中敲低BmlncR2036表达后出现体色变暗,中肠残留大量内容物,直至死亡的现象,表明BmlncR2036可能参与家蚕中肠的发育调控。【结论】发现lncRNA BmlncR2036可调控P25基因启动子活性,还可能参与调控家蚕中肠的发育。本研究为进一步探索家蚕lncRNA的功能提供了实验依据。     

兰州熊蜂蜂王头部产卵相关microRNA生物信息学分析

【目的】头部是熊蜂感觉和取食中心,也是生理行为的指挥中心。microRNA（miRNA）参与重要生理行为的调控。本研究初步探索了产卵和未产卵兰州熊蜂Bombus lantschouensis蜂王头部差异表达miRNA及其靶基因,以期探明蜂王头部miRNA在产卵过程中作用。【方法】以本土优良蜂种兰州熊蜂 B. lantschouensis 为材料,利用Solexa高通量测序技术对产卵和未产卵蜂王头部miRNA进行测序,运用生物信息学软件对miRNA及其靶基因进行预测;并应用qPCR技术验证极显著差异表达的miRNA。【结果】共获得兰州熊蜂产卵和未产卵蜂王头部miRNA的clean data数,分别为14 228 864和21 431 031 条reads;长度分析表明,在19~24 nt序列中22 nt序列数量最多,分别占产卵和未产卵蜂王序列的45.8%和45.1%。miRNA预测结果显示,共获得297个miRNA,其中270个为已知miRNA,有92个存在于蜜蜂中,其余178个分别存在于蜜蜂以外的其他物种中;另外27个为新的miRNA。在蜜蜂92个已知的miRNA中,在产卵蜂王头部中表达的共有91个,特异表达的有2个;在未产卵蜂王头部中表达的共有90个,特异表达的有1个。27个新的miRNA中,在产卵蜂王头部中表达的共有22个,特异表达的有2个;未产卵蜂王头部中表达的共有25个,特异表达的有5个。miRNA差异表达分析表明,在产卵与未产卵蜂王头部,共有8个已知miRNA表达量达到差异极显著水平（P<0.01）。qPCR结果表明,这8个表达差异极显著的miRNA均存在于熊蜂体内。在8个差异极显著的miRNA中,3个miRNA进行靶基因预测,发现它们共同调控9个靶基因,这些基因主要参与GTP结合和转录调控。【结论】本研究首次利用高通量测序技术鉴定了熊蜂头部组织的miRNA,并发现蜂王产卵与否受到miRNA差异表达的调控。结果为今后深入开展熊蜂miRNA功能验证及产卵调控机制研究奠定了基础。     

西方蜜蜂采集蜂上颚腺中高表达基因的筛选与表达分析

【目的】本研究旨在筛选西方蜜蜂Apis mellifera采集蜂上颚腺中高表达基因,为进一步筛选和研究蜜蜂采集行为相关基因提供依据。【方法】基于前期测序的西方蜜蜂5种不同职能工蜂(3日龄工蜂、10日龄哺育蜂、10日龄采集蜂、21日龄哺育蜂和21日龄采集蜂)上颚腺转录组数据,筛选采集蜂上颚腺的差异表达基因(differentially expressed genes, DEGs),并对这些DEGs进行GO和KEGG分析;qRT-PCR检测随机选取的8个DEGs在10日龄哺育蜂和10日龄采集蜂上颚腺以及两个关键DEGs(Δ-1-吡咯啉-5-羧酸合成酶基因Amp5cs和细胞色素P450 9e2基因CYP9Q3)在工蜂不同发育时期和采集蜂各组织中的表达量。【结果】筛选到22个DEGs在21日龄采集蜂上颚腺中的表达量显著高于在3日龄工蜂、10日龄哺育蜂和21日龄哺育蜂上颚腺中的表达量,同时在10日龄采集蜂上颚腺中的表达量也显著高于在10日龄哺育蜂上颚腺中的表达量。GO和KEGG富集分析显示这些DEGs主要富集在胆固醇代谢、半乳糖代谢、淀粉和蔗糖代谢、精氨酸和脯氨酸代谢、细胞凋亡-果蝇、氨基酸生物合成等方面。qRT-PCR结果表明,8个DEGs(LOC100576395, LOC411983, LOC410235, LOC725581, LOC410527, LOC406131, LOC408453和LOC410253)的表达模式与转录组数据的表达模式一致;2个关键DEGs Amp5cs和CYP9Q3在工蜂各发育阶段均有表达,且在采集蜂中表达量最高;Amp5cs在采集蜂腹、胸、上颚腺和触角中高量表达,P450 9e2在采集蜂触角 和足中表达量显著高于在其他组织中的。【结论】本研究在减小日龄因素干扰下筛选了西方蜜蜂采集蜂上颚腺中22个高表达的DEGs,这些DEGs可能主要参与采集蜂上颚腺生理发育以及能量供应、外源性物质解毒、花蜜转化等代谢通路,进而影响蜜蜂的采集行为。这些结果为西方蜜蜂上颚腺的功能研究提供理论参考,同时也为采集力强的新品种培育奠定了基础。     

沙葱萤叶甲卵滞育的转录组学分析

【目的】建立沙葱萤叶甲Galeruca daurica滞育卵转录组数据库,挖掘卵滞育相关的基因以及代谢和信号通路,在转录组水平探讨卵滞育的分子机制。【方法】采用Illumina NovaSeq6000高通量测序平台对沙葱萤叶甲滞育卵与解除滞育卵进行转录组测序,并进行生物信息学分析;利用DESeq软件分析沙葱萤叶甲滞育卵与解除滞育卵中的差异表达基因,对差异表达基因进行KEGG通路富集分析;利用qRT PCR技术对10个差异表达基因的表达模式进行验证。【结果】基于沙葱萤叶甲滞育卵与解除滞育卵转录组测序结果,共获得53 389个unigene,其中差异表达基因2 145个,24个差异表达基因与保幼激素信号及脂肪酸生物合成和降解相关。与解除滞育卵相比,滞育卵转录组中1 297个基因上调表达,富集于124条KEGG通路,其中核糖体通路显著富集;848个基因下调表达,富集于73条KEGG通路,其中MAPK信号通路和糖胺聚糖生物合成通路显著富集。qRT-PCR结果表明,随机选取的10个差异表达基因的表达趋势与RNA-Seq转录组测序结果完全一致。【结论】保幼激素,脂肪酸生物合成和降解,核糖体,MAPK信号及糖胺聚糖生物合成等通路可能在沙葱萤叶甲卵滞育调控中起着重要的作用。     

LncRNA63及其共表达基因Hr3在家蚕中的表达和功能分析

【目的】本研究旨在通过对一个新的长链非编码RNA (long non-coding RNA, lncRNA)及其共表达基因在家蚕Bombyx mori中的表达和功能进行分析,进一步阐明lncRNA在调控家蚕变态发育中的分子机制。【方法】基于前期家蚕脂肪体转录组测序数据,我们发现在蜕皮激素20E(2 μg/μL)处理5龄幼虫早期的脂肪体中lncRNA63表达显著上调,家蚕激素受体基因Hr3与之共表达。利用qRT-PCR对lncRNA63和Hr3在家蚕不同发育时期(4-5龄幼虫、蛹和成虫)、5龄幼虫不同组织(头、体壁、血淋巴、中肠、马氏管、脂肪体、丝腺、精巢和卵巢)和2 μg/μL 20E处理2, 6, 12和24 h后5龄幼虫脂肪体中的表达谱进行分析;采用核质分离检测和原位杂交技术检测lncRNA63在家蚕BmN细胞中的定位;利用dsRNA干涉家蚕5龄幼虫个体中lncRNA63的表达,并利用qRT-PCR检测RNA干涉lncRNA63对头、丝腺、脂肪体和体壁中Hr3表达的影响。【结果】qRT-PCR结果表明,lncRNA63和Hr3在家蚕羽化前的蛹末期有一个表达高峰;二者在家蚕5龄幼虫头、体壁和丝腺中都有较高表达;lncRNA63和Hr3在2 μg/μL 20E处理2, 6和12 h的家蚕5龄幼虫脂肪体中较对照组(无水乙醇处理组)都显著上调,24 h下降到与对照组相当水平,二者的表达趋势完全一致。LncRNA63主要定位于BmN细胞核,但20E可诱导lncRNA63的核质迁移。利用dsRNA干涉lncRNA63表达后,Hr3的表达显著下调。【结论】20E不仅可调控lncRNA63的表达,还可诱导lncRNA63的核质迁移;lncRNA63可通过调控共表达基因Hr3的表达参与家蚕的变态发育过程。