原创《芒果果实转录组数据组装与基因功能注释》：本论文为您写转录毕业论文范文和职称论文提供相关论文参考文献，可免费下载。
摘 要 以红芒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代谢通路,包括代谢途径、次生代谢的生物合成、植物-病原物互作、植物激素信号转导、苯丙氨酸生物合成、淀粉和蔗糖代谢、类黄酮类化合物合成等.本研究将为进一步的芒果功能基因克隆、基因表达分析、分子标记开发等研究奠定基础.
关键词 芒果；转录组；功能注释
中图分类号 S682.2 文献标识码 A
芒果是著名的热带水果,素有“热带果王”之美誉.中国是主要的芒果生产国之一,2014年栽培面积14万hm2,产量130万t,主要分布在海南、广西、云南、四川、广东和福建等省（区）.前期研究主要集中于栽培生理研究,如果实成熟过程[1-4]、果实香气挥发物质[5-6]、抗氧化活性[7-8]、采后处理和果实品质研究[9-10],尽管芒果是中国热区农业产业的经济支柱,其基因组较小,但基因组信息相当缺乏.近年来,高通量转录组测序技术的不断发展和完善为开展不同生物功能基因组学研究提供了全新的思路和方法,转录组测序是基于下一代高通量测序技术建立的一种高效、快捷分子生物学研究手段,使科技工作者能在组学水平上研究基因组序列未知的非模式生物,从整体水平了解植物在特定阶段的基因功能及基因结构,更加便利地提示特定生物学过程的分子机制[11-13].近年来,已广泛应用于杨梅[14]、苹果[15]、梨[16]、葡萄[17]、柑橘[18]等果树,芒果研究方面,Dautt-Castro等[19]用转录组测序技术研究了与 Kent芒成熟过程相关的基因,在青熟和成熟果实中,2 306个基因差异表达.Luria 等[20]用转录组测序揭示了热水处理对芒果影响的分子机制,基于基因表达谱测序技术鉴定了与热处理相关基因,Azim等[21]报道了芒果叶绿体基因组,这些研究丰富了芒果生物信息数据.有关芒果果实生长发育过程的转录组信息鲜有报道.本研究利用高通量测序技術平台illumina Hiseq 2 000对芒果发育成熟过程的果实进行转录组测序、拼接组装,再用生物信息学的方法对做到到的Unigene进行注释和功能分类,以期为芒果功能基因的发掘利用、特异miRNA的鉴定及功能分析等奠定基础.
1 材料与方法
1.1 材料
试验于2012年3～7月在广东省雷州市覃斗镇芒果基地进行,试验材料为15年生的红芒6号芒果树,试验地土壤属砖红壤,株行距4 m×4 m,树体生长良好.分别于幼果期（花后50 d）、膨大期（花后80 d）、采收期（青熟）与成熟期（常温放置后成熟）采集芒果果实,每个时期采10～15个果,采集后迅速带回实验室,将果皮与果肉分开并用液氮速冻后保存于-80 ℃冰箱中备用.
1.2 方法
1.2.1 芒果果皮与果肉RNA提取 RNA提取参照Shan等[22]建立的方法.
1.2.2 芒果果实转录组测序 基于转录组学研究的优势,以红芒6号芒果为材料,取不同发育阶段的果皮与果肉样本,分别提取总RNA后等量混合用于测序,转录组测序委托深圳华大基因公司完成.提取的总RNA使用DNaseI消化DNA后,用带有Oligo（dT）的磁珠富集真核生物mRNA；加入打断试剂在Thermomixer中适温将mRNA打断成短片段,以打断后的mRNA为模板合成一链cDNA,然后配制二链合成反应体系合成二链cDNA,并使用试剂盒纯化回收、粘性末端修复、cDNA的3′末端加上碱基“A”并连接接头,然后进行片段大小选择,最后进行PCR扩增；构建好的文库用Agilent 2 100 Bioanalyzer和ABI StepOnePlus Real-Time PCR System质检合格后,使用Illumina HiSeqTM 2000进行测序.
1.2.3 数据组装及基因功能注释 将测序获得的原始数据去除接头序列、或N含量过高或低质量的序列,做到到净读长,用Trinity软件按Min contig length 100,Group pairs distance 250,path reinforcement distance 85,Min kmer cov 2参数进行数据组装.将具有重叠区域的reads连成更长的Contig,采用Overlap的方法进一步拼接成Unigene.
通过Blastx程序将Unigene序列比对到蛋白数据库NR、Swiss-Prot、KEGG和COG（E值<0.00 001）,并通过Blastn程序将Unigene与核酸数据库NT（E值<0.00 001）进行比对,获得Unigene最高序列相似性的蛋白,便得到该Unigene的蛋白功能注释信息.将Unigene与COG数据库进行比对,获得Unigene可能的功能注释及功能分类.根据nr注释信息,利用Blast2GO软件得到Unigene的GO注释信息,并对所有Unigene用WEGO软件进行GO功能分类统计,根据KEGG注释信息进一步分析得到Unigene的Pathway途径注释.
预测编码蛋白框：按NR、Swiss-Prot、KEGG和COG顺序将Unigene序列与以上数据库Blastx比对（E值<10-5）,取Blast比对结果中排列最高的蛋白确定该Unigene的编码区序列,将编码区序列翻译成氨基酸序列,得到该基因编码区的核苷酸序列（序列方向5′→3′）和氨基酸序列.若Unigene比对不上以上数据库,用软件ESTScan[23]预测编码区,得到其编码区的核苷酸序列和氨基酸序列.
转录论文参考资料：
结论：芒果果实转录组数据组装与基因功能注释为关于本文可作为转录方面的大学硕士与本科毕业论文转录论文开题报告范文和职称论文论文写作参考文献下载。
