大花铁线莲单重被花的转录组差异分析

以即可开单被花又可开重被花的大花铁线莲重被品种"薇安"为研究对象,对"薇安"同一时期同一植株上的3种不同花瓣类型（单被、半重被、重被）采用高通量测序技术进行拼接及功能注释,筛选不同花被类型下表达量高度差异的关键基因后进行实时荧光定量PCR验证。结果显示：转录组测序共产生13.8GB原始数据,3个转录本文库两两比较共获得3 075条差异表达基因（DEG）,其中单被花与半重被样本对比（A vs B）包括649条上调DEG,605条下调DEG;半重被花与重被花对比（B vs C）包括上调DEG 1 046条和下调DEG 721条;单被花与重被花对比（A vs C）有上调DEG 1 129条和下调DEG 859条。3个不同花被下共存的差异表达基因有134条。根据基因功能注释从总DEG中筛选出26条可能与重瓣化性状相关的基因进行聚类分析,并随机挑选10个目的基因进行荧光定量PCR验证。PCR结果显示这些基因的表达量在铁线莲"薇安"同一时期同一植株的3种花被类型中均存在显著差异。最终筛选出与铁线莲重瓣化相关的关键基因有MADS-BOX类基因PMADS1、AP3、FRUITFULL、FLC;生长素反应蛋白IAA9、生长素输入载体、脱落酸8'羟化酶、吲哚乙酸诱导蛋白ARG7等。本研究为探究铁线莲重瓣花分子机制提供了基础数据和理论依据。     

转录组分析揭示东方蜜蜂微孢子虫侵染意大利蜜蜂的分子机制

【目的】本研究旨在通过差异表达基因(differentially expressed gene, DEG)分析以及毒力因子和其他侵染相关因子分析,在转录组水平揭示东方蜜蜂微孢子虫Nosema ceranae侵染意大利蜜蜂Apis mellifera ligustica的分子机制。【方法】基于前期已获得高质量的东方蜜蜂微孢子虫纯化孢子(NcCK)及侵染意大利蜜蜂工蜂7和10 d的东方蜜蜂微孢子虫(分别为NcT1和NcT2)转录组数据,根据P≤0.05且|log_2(Fold change)|≥1的标准,通过比较分析筛选出NcCK vs NcT1, NcCK vs NcT2和NcT1 vs NcT2比较组的DEG。通过相关生物信息学软件对上述DEG进行Venn分析、GO分类和KEGG代谢通路富集分析。根据Nr和KEGG数据库注释信息和相关文献进行对东方蜜蜂微孢子虫的毒力因子和侵染相关因子的统计和分析。通过RT-qPCR验证转录组数据及DEG表达趋势。【结果】从NcCK vs NcT1, NcCK vs NcT2和NcT1 vs NcT2比较组分别鉴定出1 397, 1 497和52个DEG。Venn分析结果显示各比较组共有的上调和下调基因分别为10和1个。GO分类结果显示,NcCK vs NcT1和NcCK vs NcT2中DEG富集数最多的功能条目为代谢进程、细胞进程、单组织进程、细胞、细胞组件、细胞器、催化活性和结合,而NcT1 vs NcT2中DEG富集数最多的是代谢进程、细胞进程、单组织进程、催化活性和结合。KEGG代谢通路富集分析结果显示,NcCK vs NcT1和NcCK vs NcT2中DEG分别富集到80和79条通路;富集在糖酵解/糖异生和MAPK信号通路的上调基因数量多于下调基因。毒力因子分析结果显示,孢壁蛋白9基因和孢壁蛋白12基因在NcCK vs NcT1和NcCK vs NcT2中均下调表达,孢壁蛋白8基因仅在NcCK vs NcT1中表达量下调;此外孢壁蛋白前体基因、孢壁和锚定盘复合蛋白基因、几丁质合酶基因、极管蛋白基因、蓖麻毒素B凝集素基因的表达水平在NcCK vs NcT1和NcCK vs NcT2中表现为上调。侵染相关因子分析结果表明,糖酵解途径的3个关键酶基因在NcCK vs NcT1和NcCK vs NcT2中上调表达;3个涉及ATP/ADP移位酶的基因在NcCK vs NcT1和NcCK vs NcT2中上调表达,但有1个表达量下调;2个涉及ABC转运蛋白的基因在NcCK vs NcT1和NcCK vs NcT2中上调表达,另有4个下调表达。RT-qPCR结果证实了本研究中转录组数据及DEG表达趋势的真实可靠性。【结论】本研究通过比较分析解析东方蜜蜂微孢子虫侵染意大利蜜蜂工蜂过程的转录组动态,揭示了孢壁蛋白、孢壁和锚定盘复合蛋白、几丁质酶、极管蛋白和蓖麻毒素B凝集素等毒力因子编码基因,以及己糖激酶、丙酮酸激酶、6-磷酸果糖激酶、ATP/ADP移位酶和ABC转运蛋白等侵染相关因子编码基因在病原增殖中扮演重要角色,为阐明东方蜜蜂微孢子虫的侵染机制提供了基础。     

茶丽纹象甲为害诱导的茶树防御反应分析

基于转录组测序,获得了茶树新梢叶片被茶丽纹象甲为害前后的差异基因表达谱,并从信号监测与转导、转录因子、次生代谢物生成等方面分析了相关基因的表达变化情况,为进一步研究茶丽纹象甲为害诱导的茶树防御机制奠定基础。结果表明:(1)茶丽纹象甲为害茶树新梢叶片后共检测到显著上调表达基因309个,显著下调表达基因297个,这些显著差异表达基因共分成23类。(2)检测到信号监测与转导过程中合计17个单基因上调表达2.0～2.8倍,包括富含亮氨酸重复序列的类受体蛋白激酶基因、促分裂原活化蛋白激酶级联信号系统蛋白激酶基因、茉莉酸信号通路基因、钙离子信号基因、水杨酸代谢通路基因;检测到21个转录因子单基因上调表达1.8～2.8倍,包括bHLH、WRKY、bZIP、MYB、UAF等转录因子;检测到苯丙类、类黄酮及萜类物质代谢过程中合计8个单基因上调表达2.1～2.5倍,包括苯丙酮酸互变异构酶基因、肉桂酰辅酶A还原酶基因、二氢黄酮醇还原酶基因、花色素还原酶基因、法尼基二磷酸合成酶基因、法尼醇脱氢酶基因。研究认为,茶丽纹象甲为害诱导信号转导、转录因子及次生代谢过程中的基因增量表达以机械损伤作用为主,昆虫口腔分泌物对这些基因的诱导表达具有协同促进作用。     

感染球孢白僵菌的烟粉虱若虫免疫应答转录组分析

【目的】筛选烟粉虱Bemisia tabaci(Gennadius)若虫应对白僵菌Beauveria bassiana侵染的应答基因,以进一步研究烟粉虱免疫反应的分子机制。【方法】采用新一代高通量测序技术对感染和非感染白僵菌的烟粉虱4龄若虫进行了测序分析,并筛选了差异表达基因;利用生物信息学工具对转录组测序得到的基因进行了功能注释、分类以及参与的信号通路展示。【结果】组装得到非冗余Unigene 232 554个,其N50和N90分别为1 153 bp和260 bp,平均长度为67 424 bp。对所有的基因进行差异表达分析发现:以P<0.05为标准筛选得到了1 166个差异表达基因,其中上调表达基因有474个,下调表达基因有692个,其中,与免疫相关的基因有405个;GO富集分析发现:有416个GO term有富集现象,包括156个生物学过程(66 402个Unigenes),89个细胞组分(27 645个Unigenes)和154个分子功能(73 417个Unigenes);KEGG代谢通路分析将烟粉虱转录组1 145个DEG匹配到309个通路上,其中76个通路得到了富集。【结论】对405个可能参与到烟粉虱若虫对白僵菌侵染的免疫识别和防御的基因进行了测序。研究结果为从分子水平上开展应用虫生真菌防治烟粉虱的研究奠定信息学基础。     

杂交兰及其花艺突变体中类黄酮成分差异的代谢-转录组联合分析

为探索杂交兰花艺突变体变异的分子机理,该研究以杂交兰‘玉凤’及其花艺突变体‘双艺金龙’为材料,利用靶向代谢组学和转录组学鉴定两者中类黄酮化合物含量差异及其相关通路上的差异表达基因。结果表明：(1)代谢组分析发现,杂交兰‘玉凤’和‘双艺金龙’花瓣中共检测到271种类黄酮代谢物,其中黄酮醇类和黄酮类代谢物的相对含量约占总类黄酮的30%～50%;共检测到38个差异代谢物(15个上调,23个下调);‘玉凤’中差异最高的代谢物二氢山奈酚-7-O-葡萄糖苷(黄酮醇类)含量是‘双艺金龙’的124 444倍,‘双艺金龙’中差异最高的代谢物3,5,7,3′,4′-五羟基-5′-异戊二烯基黄酮(黄酮醇类)含量是‘玉凤’的7 244倍。(2)KEGG分析显示,差异代谢物显著富集在类黄酮、黄酮和黄酮醇生物合成途径。其中,在类黄酮生物合成途径上,根皮苷、黄腐醇、二氢山奈酚、二氢杨梅素和表没食子儿茶素含量升高,柚皮苷和二氢槲皮素含量降低;在黄酮和黄酮醇生物合成途径上,三叶豆苷和槲皮素含量均下降。(3)转录组分析共筛选到563个差异表达基因,与‘玉凤’相比,‘双艺金龙’中有220(39.1%)个基因上调表达,343...     

基于转录组挖掘‘云香’水仙萜类合成基因及NaIDI的克隆

为深入挖掘中国水仙转录组中萜类合成相关基因,了解中国水仙萜类代谢途径和分子调控机制,该研究以‘云香’水仙为材料,在其盛花期花瓣与副冠转录组测序的基础上,筛选已注释的萜类合成途径基因并利用NCBI blastn进行再注释,通过对部分候选基因的表达量与生物信息分析进一步筛选出代表基因,采用RT PCR技术克隆了水仙异戊烯基焦磷酸异构酶基因NaIDI,并对其蛋白序列与特异性表达进行分析。结果表明：(1)Blast比对筛选得到52 个与萜类合成上游途径相关的Unigenes,二次筛选获得11个显著差异表达的代表基因。(2)NaIDI基因开放阅读框(ORF)长度为858 bp,编码285个氨基酸,其氨基酸序列与‘金盏银台’水仙相似度97.19%;亚细胞定位预测显示该基因定位在叶绿体;系统进化分析表明其与芦笋亲缘关系比较近。(3)实时荧光定量分析结果表明,NaIDI基因在水仙开花的不同时期和不同组织器官中差异表达显著,且在盛花期时副冠中的表达量最高,与中国水仙挥发性萜类化合物在开花不同时期和不同组织器官中的表达规律一致,表明NaIDI在萜类代谢中发挥着一定作用。     

茶树种子发育过程的转录组分析

该研究以‘铁观音’茶树品种的种子为试验材料,采用转录组测序技术分析种子发育的3个时期(幼果期、膨大期、成熟期)的表达差异,探究茶树种子油脂代谢的分子机制。结果表明:(1)经转录组测序、组装后共获得30 940 581个clean reads,经数据合并拼接最终得到36 951条非冗余Unigene序列,其中28 476个Unigene可得到功能注释;在转录本中能够被注释到GO分类的Unigene有11 201条(30.3%),KEGG分析发现共有17 172个基因参与了127个代谢通路。(2)经KEGG通路筛选出14条与脂肪酸代谢相关的通路,且随着茶籽的发育,大部分脂肪酸调控途径相关基因呈下调趋势,其中上调基因数最多的有α-亚麻酸代谢途径和脂肪酸降解途径(有17个基因表达量上调),下调基因数最多的是甘油磷脂代谢途径(有58个基因表达量下调);在茶籽发育幼果期α-亚麻酸代谢途径中表达量上调的基因数超过表达量下调的基因数。(3)研究发现茶籽脂肪酸合成相关的基因涉及14个脂类调控途径,共409条差异基因;随着茶树种子发育到成熟期,上调的差异表达基因数量在减少,下调的差异表达基因数量增加,其中α-亚麻酸途径中的基因PLA2G16、DAD1、pldA、FabF、FabI表达量上调显著,随后表达量下调。(4)qRT-PCR检测结果表明,7个茶树FAD和1个ACP差异表达基因的水平与转录组测序结果基本一致;随着茶籽的发育,基因CsFAD7和Δ6-CsFAD从幼果期、果实膨大期至果实成熟期都为差异下调表达,CsFAD2、CsFAD6和Δ7-CsFAD为差异上调表达,CsFAD8、Δ8-CsFAD和CsACP在幼果期至果实膨大期差异上调表达,在果实膨大期至果实成熟期差异下调表达。     

抗菌肽LL-37抑制肝癌细胞恶性增殖的转录组分析

抗菌肽LL-37与肿瘤的发生发展密切相关,课题组前期研究发现LL-37能够抑制肝癌细胞恶性增殖。为了给其抑制肝癌发生发展的分子机制研究提供更多的生物学依据,本研究通过高通量RNA测序技术以及生物信息学方法对LL-37作用前后肝癌细胞中的差异表达基因(differentially expressed gene, DEG)进行了分析,共筛选出753个DEG,其中上调的DEG 374个,下调的DEG 379个;进一步对筛选出的DEG进行基因本体论(Gene Ontology, GO)及京都基因和基因组数据库(Kyoto Encyclopedia of Genes and Genomes, KEGG)通路富集分析,并构建DEG编码蛋白互作网络,筛选出10个可能参与LL-37抑制肝癌细胞恶性增殖的潜在关键基因。经分析这些基因均表达下调,且对机体炎症的激活、转运以及肿瘤细胞的增殖、迁移至关重要。以上结果为揭示LL-37在肝癌发生发展中的作用及机制提供了数据基础,为探索肝癌诊断和治疗手段提供了新的思路。     

小菜蛾应对球孢白僵菌感染的免疫应答基因分析

通过室内生物测定,筛选对小菜蛾Plutella xylostella具有高毒力的球孢白僵菌Beauveria bassiana菌株,分析小菜蛾应对球孢白僵菌浸染的免疫应答基因。采用新一代Illumina高通量测序技术对接种/感染48h与未接种的小菜蛾样品分别进行转录组测序,筛选差异表达基因;结合生物信息学分析差异表达基因的功能、分类及涉及的信号通路等。结果表明接种小菜蛾48h后诱导2 434个基因发生差异表达,包括1 080个上调基因和1 354个下调基因。GO分类结果表明,这些差异表达基因（DEGs）被注释到45个GO term中,包括23个生物学过程,12个细胞组分和10个分子功能。KEGG分析表明,1 308个DEGs被富集到25个功能途径,这些DEGs包括497个上调基因和811个下调基因。分析发现,DEGs编码蛋白包含肽聚糖识别蛋白、丝氨酸蛋白酶55、丝氨酸蛋白酶抑制剂以及细胞色素P450 6K1类的免疫相关蛋白。另外,组蛋白、磷酸泛酰半胱氨酸脱羧酶、双链RNA结合蛋白、黑芥子酶等基因在球孢白僵菌侵染过程中也高表达,推测这些基因可能参与小菜蛾应对球孢白僵菌侵染的免疫反应。研究结果为进一步研究小菜蛾免疫相关基因的功能奠定了基础。     

基于高通量转录组测序技术的无精子症患者睾丸组织比较转录组分析

本研究探讨不同程度无精子症患者睾丸组织的转录组差异,了解差异表达基因(DEG)在功能、分类和代谢通路的不同,揭示无精子症患者精子发生分子机制,为促进男性不育研究的发展提供理论基础.选取1份非梗阻性无精子症和4份梗阻性无精子症患者睾丸组织样品(从无精子到有精子),进行RNA提取和文库构建,利用Illumina HiSeqTM2500高通量测序,构建无精子症患者睾丸组织转录组文库,并用生物信息学方法进行分析.结果发现,样品比对基因组数据库的平均比对率为94.38%,共检测2 242个属于预测新的蛋白质编码基因的转录本.得到差异表达基因统计结果为:NOA vs. OA1基因上调8 045,下调1 150;OA1 vs. OA2基因上调1 538,下调420;OA2 vs. OA3基因上调1 275,下调1 690;OA3 vs. OA4基因上调1 834,下调1 853.比较5例无精子症睾丸组织的差异基因KEGG,主要富集在RNA降解通路、基底细胞瘤通路、癌通路、黑色素生成通路和调节干细胞多能性等信号通路. PRM1、PRM2、TNP1、UBXN6、CXCL16、NUPR2、CCDC136和CRISP2等基因的表达呈递增趋势,并具有时序特异性.此外,5例无精子症睾丸组织的表达基因有不同程度的基因融合.综上,基因融合可能和无精子症相关,并且不同程度无精子症患者睾丸组织的差异表达基因数量、功能、分类和代谢通路不同.本研究筛选出精子发生、精子运动等差异表达基因,丰富了无精子症患者睾丸组织转录组信息,为开展无精子症患者睾丸组织相关基因及分子调控机制的研究奠定基础.     

白化菠萝蜜茎2个蛋白激酶和4个转录因子家族分析

为探究蛋白激酶(PKs)和转录因子(TFs)在白化菠萝蜜(Artocarpus heterophyllus)幼苗茎次生生长中的表达变化,基于转录组数据对其差异表达基因(DEGs)进行预测及分类,并对挑选出的2个PKs和4个TFs家族构建系统进化树。结果表明,胞质类受体激酶(RLCK)-VIII家族的DEGs上下调表达各4个,亮氨酸富集重复类受体激酶(LRR-RLK)-X家族Xa和Xb-2分支中的DEGs均下调表达,Xb-1中的均上调,TCP家族的20个DEGs中有15个上调表达,zf-HD和GRF家族中的大多数DEGs上调表达,Alfin-like家族中的DEGs均下调表达。因此,这表明6个家族可能在菠萝蜜茎的次生生长过程和应对非生物胁迫中发挥重要作用。     

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

为更好地挖掘八角(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家族序列数量最多。该文利用转录组测序技术分析八角优良无性系与普通品种叶片的差异基因及...     

三峡库区消落带野生狗牙根对水淹适生性的转录组分析

为探究狗牙根对水淹环境的适生性机制,取长势一致的狗牙根分为未水淹对照组(A1)和30 cm沉水处理组(A2),A2水淹15 d后两组同时采样,提取样本总RNA作转录组测序。A1和A2测序数据经de novo组装分别获得128031和83363条高质量Unigene;A1和A2对比发现28844条差异表达基因(differentially expressed genes, DEGs),其中11858条DEGs上调表达,16986条DEGs下调表达;功能注释显示,这些DEGs主要集中在转录、翻译、碳水化合物代谢和环境适应等途径;信号通路富集显示,包括植物对病原体的防御、物质代谢、光合作用、形态发育和抗氧化等过程的103个信号通路与狗牙根对水淹生境的适生性有关;随机挑选4个DEGs作qRT-PCR分析,荧光定量结果与测序数据基本一致,说明基于转录组分析DEGs的表达情况较为可靠。该工作或可为进一步挖掘与狗牙根水淹适生性相关的调控基因和功能基因提供基础。     

杜仲雄花芽2个发育时期转录组分析

杜仲（Eucommia ulmoides）雄花富含多种活性成分和营养成分,具有重要的药用和营养价值。为了揭示杜仲雄蕊原基发育相关基因的表达情况,为杜仲雄花芽发育分子调控机制研究提供理论参考。本文以杜仲良种"华仲11号"（ "Huazhong No.11" ）为材料,采用lllumina高通量测序技术,分别对苞叶原基分化期和雄蕊原基分化期的花芽进行转录组测序,通过生物信息学对2个发育时期的转录组进行比较分析,筛选出与雄花芽形态发育相关的差异基因。结果显示,转录组测序共获得40.48 Gb过滤数据,各样品的clean reads与杜仲基因组进行序列比对,比对效率为90.56%~93.01%。在2个发育时期筛选出583个差异表达基因,其中在雄蕊原基发育期上调基因315个,下调基因267个。差异基因GO和KEGG功能分析显示,差异基因富集在与生长发育、光周期途径、激素合成和信号传导、碳代谢等相关的生物过程和代谢通路。结果显示光周期途径是杜仲成花诱导的重要途径,同时雄花芽在形态分化过程中受碳水化合物、植物激素和其他代谢物质调控。此外,MADS-box家族成员FLC、SOC1、AGL3和AGL8参与杜仲雄蕊器官发育。本研究为杜仲花发育基因调控提供了基础数据,也为雄花用杜仲的分子育种提供了参考。     

基于比较转录组学分析NaCl胁迫影响德尔卑沙门氏菌的耐渗机制

【目的】研究高渗胁迫条件下德尔卑沙门氏菌(Salmonella enterica subsp. enterica Derby, S. Derby)的转录组调控机制,分析差异表达基因(differentially expressed genes, DEGs)表达水平,探究在高渗胁迫影响下德尔卑沙门氏菌耐渗反应的相关代谢通路。【方法】通过高渗胁迫诱导德尔卑沙门氏菌的耐渗性,提取菌株的总RNA,去除rRNA,构建cDNA文库。利用转录组测序技术及生物学信息技术分析相关DEGs,并通过实时荧光定量PCR (real-time fluorescence quantitative PCR, qRT-PCR)进行验证。【结果】胁迫组德尔卑沙门氏菌通过转录组测序结果发现有3 950个DEGs,其中具有显著上调的基因21个,显著下调基因38个。涉及到细胞膜蛋白、氨基酸的代谢等相关基因上调,协助德尔卑沙门氏菌在高渗环境中存活。与此同时,胁迫组德尔卑沙门氏菌的糖转运系统(sugar transport system, PTS)、糖酵解过程以及抗氧化性相关基因表达显著下调,这是由于高渗环境菌体需要在体内储存大量糖类等物质,从而降低了糖原的消耗,进而导致细胞外膜的脂多糖合成受到抑制,降低了高渗胁迫下德尔卑沙门氏菌细胞膜表面的O抗原的合成。【结论】高渗环境诱导后显著提高了德尔卑沙门氏菌的耐渗性,其中Na⁺/H⁺逆向转运蛋以及谷氨酸的代谢通路发挥着重要的作用,为进一步了解以及更好地控制其在食品中的污染提供了理论依据。     

家蚕5龄幼虫精巢和卵巢microRNA芯片及转录组比较分析

【目的】本研究旨在通过对家蚕Bombyx mori 5龄幼虫精巢和卵巢组织微小RNA (microRNA, miRNA)基因芯片及转录组进行分析,找到参与家蚕性腺发育相关的miRNA分子及可能的靶基因。【方法】采用新一代高通量测序平台对家蚕5龄幼虫精巢和卵巢(分别定义为Test和Control)进行miRNA基因芯片检测及转录组测序分析,根据P＜0.05且log₂(fold change, FC)≥2的标准,通过比较筛选出Test vs Control的差异表达miRNA;根据q≤0.05且|log₂(fold change)|≥1的标准,通过比较筛选出Test vs Control的差异表达基因 (differentially expressed genes, DEGs);随机选取8个上调和12个下调差异表达miRNA,对其表达及其预测的5个靶基因进行qRT-PCR验证;对DEGs以及差异表达miRNA的靶基因进行KEGG通路富集分析。【结果】从精巢和卵巢样本中(Test vs Control)分别鉴定出68个差异表达miRNA和3 991个DEGs,其中上调和下调miRNA分别为36和32个,上调和下调DEGs分别为2 033和1 958个。差异表达miRNA的qRT PCR验证结果均与芯片数据一致。KEGG通路富集分析结果显示DEGs在新陈代谢及核糖体的信号通路显著富集。对差异表达miRNA在DEGs中的可能靶基因进行预测,结果找到了4组表达趋势相反的miRNA与靶基因：分别是bmo-miR-2774a与LOC101745556;bmo-miR-92b与LOC101735954以及bmo-miR-3266与LOC733130和LOC778467;1组表达趋势一致的miRNA与靶基因：bmo-miR-3321与LOC101744895。5个靶基因的qRT-PCR验证结果与转录组测序结果一致。【结论】本研究获得了家蚕5龄幼虫精巢和卵巢转录组及miRNA芯片数据,筛选并验证了4组差异表达和1组一致表达miRNA及潜在靶基因,为探究家蚕精巢和卵巢发育差异奠定了基础。     

冰菜盐胁迫下的转录组分析

为了解冰菜(Mesembryanthemum crystallinum)叶片抗盐相关基因组学,利用Illumina Hi-seq TM2500高通量测序技术研究冰菜叶片在400 mmol L~(–1) NaCl胁迫下转录组基因的差异表达。结果表明,从400 mmol L~(–1) NaCl胁迫和对照的冰菜叶片中共获得13.01 Gb Clean data,Q30碱基均大于90.08%。共获得123个差异表达基因(DEGs),包括73个上调基因,50个下调基因,其中功能注释的基因有96个。根据Unigene库序列进行GO、COG和KEGG注释,筛选出8个与抗盐性相关差异表达基因,植物激素代谢相关基因,脱落酸8'-羟基化酶、吲哚-3-乙酰酸酰胺合成酶和茉莉酮酸酯ZIM结构域蛋白基因均下调表达,生长素响应蛋白、细胞分裂素合酶基因则上调表达,糖代谢相关基因棉子糖合成酶基因上调表达,质膜H+-ATPase基因上调表达,脱水蛋白基因下调表达。这为冰菜耐盐基因组学和分子生物学的研究奠定基础。