滨麦叶片转录组分析

滨麦（Leymus mollis）是小麦的野生近缘种，具有良好的耐盐抗旱和抗病虫害的能力，是小麦遗传育种的良好资源。为了深入了解和挖掘滨麦的基因信息和优异资源，本研究利用高通量测序技术对生长于滨海沙地的野生滨麦的叶片转录组进行了深度测序并组装，对得到的所有Unigene进行COG、GO和KEGG分类和功能注释，对Na⁺转运相关蛋白基因和氧化胁迫响应基因进行了挖掘，并通过实时荧光定量反转录PCR（qRT-PCR）对随机选取的12个基因的表达模式进行了验证。结果显示，转录组测序共获得112 846条Unigene，其中的59 380条得到功能注释，占总数的52.62%。COG分类结果表明，15 786条Unigene归属于25个类别。GO分类结果表明，20 350条Unigene被注释到三个大类中，其中，属于“生物学过程”的Unigene数量最多，占总数的43.56%。KEGG分析结果表明，18 550条Unigene得到注释，共涉及到128条代谢途径。其中，含Unigene数量最多的类别是“代谢通路”，涉及到的Unigene占总数的27.86%。“植物病原互作”和“植物激素信号转导途径”涉及的Unigene数量也较多，分别为1 367条和861条。对Na⁺转运相关蛋白基因和氧化胁迫响应基因进行挖掘，发现了15条注释为Na⁺/H⁺反向转运蛋白的Unigene和175条响应氧化胁迫的Unigene。随机选取的12个基因的qRT-PCR结果与转录组测序结果基本一致。     

Argonaute1敲除前后尖孢镰刀菌转录组分析比较

 Argonaute蛋白(AGO)介导的沉默复合体在RNA干扰(RNAi)中起着至关重要的作用。为了探究AGO1在尖孢镰刀菌RNAi中的作用机制,本文以粘团专化型尖孢镰刀菌生理1号小种FOX-A8野生型和其AGO1缺失突变体(FOX-A8-△Ago1)的菌丝和孢子为材料,分别进行了RNA提取、Illumina HiSeq 2000高通量转录组测序、差异表达基因(DEGs)的显著富集分析;选择菌丝和孢子中的DEGs 进行实时荧光定量PCR(Quantitative real-time PCR, qRT-PCR)验证。GO通路注释结果显示,相对于野生型菌株,AGO1缺失突变体菌丝中的醇脱氢酶(NADP⁺)、孢子中的CAMK / CAMKL / CHK 1蛋白激酶均显著上调;KEGG通路注释结果显示,相对于野生型菌株,AGO1缺失突变体菌丝中与 MAPK信号通路相关的基因、孢子中与PLD信号通路相关的基因均显著下调;另外,相对于野生型菌株,编码AGO2的基因下调,但是下调不显著。qRT-PCR检测DEGs的表达模式与RNA-Seq分析结果一致,证实了RNA-Seq结果的可靠性。     

长江刀鲚选育群体转录组EST-SSR的分布特征分析

本研究利用MISA软件挖掘长江刀鲚(Coilia ectenes)肌肉和肝脏转录组中的微卫星标记，为刀鲚选育群体的种质资源评估和分子标记辅助育种奠定基础。结果显示，从71869条Unigenes中共获得33896条重复单元长度为1~6碱基的微卫星序列；刀鲚转录组中不同类型微卫星的重复基序具有不同的分布特征，其中，单核苷酸重复、二核苷酸重复和三核苷酸重复为主要的微卫星重复类型，分别占总微卫星数目的34.94%、49.47%和13.34%；不同微卫星重复类型的优势重复基序亦有所不同，其中，A/T为单核苷酸重复基序的优势重复基序占86.25%，AC/GT为二核苷酸重复基序的为优势重复基序占75.25%，AGG/CCT为三核苷酸重复基序的优势重复基序占28.57%；不同微卫星重复基序核苷酸的数量和重复次数亦有所不同，重复次数伴随着重复单元中核苷酸数量的增加而呈现降低的趋势；从100对四核苷酸重复的SSR引物中筛选获得了16对多态性微卫星标记，并以此为基础，对长江刀鲚选育群体(F3)的遗传学特征进行了初步评估，结果显示，长江刀鲚选育群体F3的平均有效等位基因数(Ne)、平均观测杂合度(Ho)、平均期望杂合度(He)和Shannon多样性指数I分别为1.7580、0.3414、0.3977和0.6278。以上结果表明，基于刀鲚转录组数据批量开发微卫星是切实可行的，所开发的多态性微卫星标记能够应用于长江刀鲚选育群体的遗传背景评估和进一步的遗传育种研究。     

旱垂柳转录组测序及赤霉素在垂枝形成中的作用

垂枝作为特殊的林木形态结构,因形态优美在全世界被广泛用于绿化和园林景观而成为研究热点。国内外对垂枝植物的研究表明赤霉素在垂枝形成中起作用,但是垂枝形成的分子背景研究仍比较缺乏。本研究以旱垂柳(Salix matsudana var. pseudo-matsudana)的垂枝为研究对象,通过对旱垂柳在不同生长发育时期的叶片、茎尖和茎基部进行了转录组和数字表达谱测序研究,共获得71 474个Unigenes序列。将Unigenes与COG和KEEG数据库比对,与转录和信号转导机制相关的基因在旱垂柳中高度富集,最终确定8 359个差异表达的基因,其中427个基因与激素相关。进一步分析赤霉素相关基因得出,主要参与2条赤霉素合成途径,具有生物活性的GA4是合成途径中的主要产物,而且GA4在枝条中的含量显著高于GA1和GA3。参与赤霉素信号转导途径中的GID1B在旱垂柳中上调表达,以及在G1和G4生长期的表达量高于其他时期,而DELLA蛋白GAI在G2生长期以后下调表达,GID1B和GAI在旱垂柳的垂枝生长过程中可能起非常重要的作用。这些研究结果将有助于探索赤霉素在旱垂柳垂枝形成中的分子作用,可为林木形态的遗传改良提供基础理论依据。     

水杉长链非编码RNA分析

水杉(Metasequoia glyptostroboides)是中国著名的一级保护植物,素有"活化石"之称。为了更好地保护和利用这一重要种质资源,本研究以水杉为实验材料,在全长转录组测序的基础上,利用四个软件对61057个Unigenes进行了蛋白编码潜能预测,最终获得3385个长链非编码RNA(long non-coding RNAs,lncRNAs)分子。进一步对这些lncRNAs的长度和数量作统计分析,结果显示,最短为200 nt,最长为11134 nt,平均长度为1956 nt;大部分lncRNAs的长度介于200~2000 nt,占总量的59.5%;在2001~4000 nt的长度范围内,也发现较多的lncRNAs,达到总量的32.2%;在长度为4001~6000 nt的范围内,只有7.6%;长度超过6000 nt的lncRNAs很少,总共只有25个(占总量的0.7%)。本研究结果为水杉资源的保护和利用研究提供了一些分子依据,为其它植物在相关领域的研究提供了一些参考。     

适用于转录组测序的大白菜小孢子总RNA提取方法筛选

为了找到适用于转录组测序要求的大白菜小孢子总RNA 提取方法，以出胚率较高的吉红82 为试材，以热激诱导处理后的小孢子为研究对象，比较分析了TRIzol 法、改良CTAB 法及超纯试剂盒法3 种方法的RNA 提取效果。结果表明：TRIzol法（磨样）和改良CTAB 法提取的RNA 浓度较高，但完整性较差；超纯RNA 提取试剂盒法提取的RNA 质量最高，且完整性好，能够满足转录组测序的要求。     

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

杜仲(Eucommia ulmoides)的果实是提取杜仲胶和杜仲籽油的重要原材料,在工业、医疗、食品等领域开发前景广阔。为了揭示雌蕊原基发育相关基因的表达情况,本研究以杜仲果用良种‘华仲6号’(Huazhong No.6)为材料,利用lllumina Hiseq X-10高通量测序平台,分别对花序原基分化期和雌蕊原基分化期的花芽进行转录组测序,通过生物信息学对2个发育时期的雌花芽转录组进行比较分析,筛选出与杜仲雌蕊原基发育相关的差异基因。结果显示,转录组测序共获得68.11 Gb数据,各样品的Clean reads与杜仲基因组进行序列比对,比对效率为91.59%~92.05%,同时2个发育时期共筛选出485个差异表达基因,其中219个差异基因在雌蕊原基分化期表达上调,266个表达下调。差异基因GO富集结果表明,花的发育、光周期现象、激素生物合成以及蛋白结合转录因子活性等相关的生物途径被富集。KEGG富集结果表明淀粉和蔗糖代谢、碳代谢、昼夜节律和植物激素信号转导等代谢通路被富集。结果表明光周期途径是诱导杜仲成花的重要途径,且雌蕊原基发育受碳水化合物、植物激素和其他代谢物质调控。本研究为杜仲花器官发育调控基因的挖掘提供了基础数据,也为果用杜仲的分子育种提供了参考。     

芒果果实转录组数据组装及基因功能注释

以红芒6号为试材,对生长发育各时期果皮与果肉样品提取RNA后等量混合进行转录组测序,共获得了68 419 722个reads,包含6 157 744 980个核苷酸序列信息,平均读长90 bp,序列信息全部登陆NCBI数据库(SRP035450)。对reads进行拼接,共获得124 002个Contig片段,进一步组装获得54 207个Unigene片段,平均长度为838 bp,阈值小于10-5的序列中共有42 515(78.43%)个unigenes被注释到公共蛋白数据库。其中,35 198个被注释到生物学过程、分子功能和细胞组分GO类别的44个功能组,14 619个Unigene匹配到25类COG功能组,23 741个Unigene被富集到128个KEGG代谢通路,包括代谢途径、次生代谢的生物合成、植物-病原物互作、植物激素信号转导、苯丙氨酸生物合成、淀粉和蔗糖代谢、类黄酮类化合物合成等。本研究将为进一步的芒果功能基因克隆、基因表达分析、分子标记开发等研究奠定基础。     

基于高通量测序的露地菊(Chysanthemum×grandiflora)盐胁迫转录组分析

为揭示露地菊生长发育及耐盐胁迫的应答机制和分子基础,本研究以盐胁迫处理的露地菊及其对照为材料,使用Illumina Hiseq2500高通量测序平台对转录组进行测序,分别获得了60370448和71415448条Clean reads,通过序列拼接组装得到45591条Unigene,平均长度724 bp。有37675条Unigene在七大数据库(COG,GO,KEGG,KOG,Pfam,Swiss-Prot,NR)中得到注释。通过比对露地菊盐胁迫处理组和对照组样品间Unigene的表达量及在各数据库中的注释情况,统计得到:有4143条差异表达基因获得注释;有2441条差异表达基因在GO数据库中获得功能注释;注释到COG数据库中的2281条差异表达基因依据功能可分为25类;有1062条差异基因映射到KEGGPathway数据库中,涉及了199个代谢通路,包括核糖体途径、植物激素信号传导途径、淀粉蔗糖代谢、碳代谢、氨基酸的生物合成等。本研究获得的转录组数据将有助于揭示露地菊生长发育及耐盐胁迫的应答机制和分子基础,及相关抗性基因的挖掘和分子辅助育种等方面的研究。