芝麻发育转录组分析

doi:10.3864/j.issn.0578-1752.2012.07.002

摘要/Abstract

摘要： 【目的】系统了解芝麻发育及种子形成转录组特征，丰富芝麻转录组数据信息。【方法】选用6份芝麻样品（5个不同芝麻品种的完整植株、1份不同发育阶段的芝麻籽粒），构建转录组测序文库并进行Illumina RNA-seq双端测序及生物信息学分析。【结果】共获得原始数据12.69 Gb，有效数据8.80 Gb。通过de novo拼接获得了长度大于100 bp的转录物26 837条（http://www. ncbi.nlm.nih.gov/genbank/TSA.html，登录号：JP631635—JP668414）；转录物总长度18.35 Mb，平均长度683 bp，N50长度1 006 bp。转录物注释结果显示，25 331条转录物序列具有同源比对信息；1 506条转录物序列无匹配（no hits）序列信息，可能为芝麻特有的基因序列。采用COG、GO功能分类工具可将已注释转录物序列划分为24个或42个功能类别，共涉及物质及能量代谢、信号传导、转录调控及防卫反应等诸多生理生化过程。通过比较不同材料间转录物序列及表达水平，初步确定1 277条序列在种子形成过程中表达量下调10倍以上，990条序列仅在芝麻植株中表达而未在籽粒中表达；660条序列在种子形成过程中表达量上调10倍以上，296条序列可能与种子形成特异相关。【结论】利用高通量测序技术对芝麻野生种和不同栽培种现蕾期植株以及种子形成过程的转录组进行研究，揭示了芝麻发育转录组的整体表达特征，在得到大量芝麻转录组unigene序列的同时，获得了一批在芝麻生长发育及籽粒形成过程中有重要功能的基因序列。为深入开展芝麻生长发育、籽粒发育相关基因功能及调控以及芝麻分子标记开发等研究提供了丰富的数据资源。

关键词: 芝麻, 发育, 转录组测序, 生物信息学分析

Abstract: 【Objective】 To enrich sesame transcriptome data, the transcriptome sequencing and bioinforamtics analysis on sesame growth and development and seed formation processes were performed in this study. 【Method】 Six transcriptome sequencing libraries for developing sesame plants and seeds were constructed and sequenced using Illumina RNA sequencing technique, and the global transcriptome information was analyzed subsequently.【Result】8.80 Gb available trancriptome data were acquired as the adaptor sequences, duplication sequences and low-quality reads were removed from the original 12.69 Gb solexa sequencing data. And 26 837 uni-transcripts, longer than 100 bp, were obtained by de novo assembly method (http://www. ncbi.nlm.nih.gov/genbank/TSA.html, GenBank ID: JP631635-JP668414). The total scaffold sequence length reached 18.35 Mb with the 1 006 bp of N50, the average uni-transcript length was 683 bp. Annotation analysis of uni-transcripts indicated that 25 331 transcripts had homolog in public protein database; however, 1 506 sequences were no hit and might be sesame-specific. With COG and GO functional classifications, all uni-transcripts were grouped into 24 and 42 function categories, respectively, in which many functional categories, such as material and energy metabolism, signaling, transcription regulation and defense reactions, etc. were included. Furthermore, compared with plant growth and development transcripts, 1 277 sequences expressed more than 10-fold low in seed formation process, and the transcript level of 990 uni-trancripts could not be measured. In addition, 660 sequences were found with more than 10-fold high expression level during seed development, in which 296 sequences seemed as seed-specific. 【Conclusion】With the transcriptome sequencing of Sesamum radiatum and several cultivars of S. indicum L., this study gave a global insight into the characteristics of sesame development transcriptome and thousands of transcript sequences with important function were acquired for future genes expression and regulation research about sesame growth and development.

Key words: sesame (Sesamum indicum L., Sesamum radiatum), development, RNA sequencing, bioinformatics analysis

魏利斌, 苗红梅, 张海洋. 芝麻发育转录组分析[J]. 中国农业科学, 2012, 45(7): 1246-1256.

WEI Li-Bin, MIAO Hong-Mei, ZHANG Hai-Yang. Transcriptomic Analysis of Sesame Development[J]. Scientia Agricultura Sinica, 2012, 45(7): 1246-1256.

0
/ / 推荐

导出引用管理器 EndNote|Reference Manager|ProCite|BibTeX|RefWorks

链接本文: https://www.chinaagrisci.com/CN/10.3864/j.issn.0578-1752.2012.07.002

https://www.chinaagrisci.com/CN/Y2012/V45/I7/1246

参考文献

[1]张海洋, 郑永战, 卫双玲. 我国芝麻产业现状与发展对策. 全国芝麻和特油作物产业暨学术研讨会论文汇编, 2007: 9-18.

Zhang H Y, Zheng Y Z, Wei S L. Sesame industry status and development strategy in China. China Sesame and Specific Oil Crops Industry Development and Academic Research Discussion Meeting Papers Edition, 2007: 9-18. (in Chinese)

[2]Kim D H, Zur G, Danin-Poleg Y, Lee S W, Shim K B, Kang C W, Kashi Y. Genetic relationships of sesame germplasm collection as revealed by inter-simple sequence repeats. Plant Breeding, 2002, 121: 259-262.

[3]Wei L B, Zhang H Y, Zheng Y Z, Guo W Z, Zhang T Z. Developing EST-derived microsatellites in sesame (Sesamum indicum L.). Acta Agronomica Sinica, 2008, 34(12): 2077-2084.

[4]Wei L B, Zhang H Y, Zheng Y Z, Miao H M, Zhang T Z, Guo W Z. A genetic linkage map construction for sesame (Sesamum indicum L.). Genes and Genomics, 2009, 31(2): 199-208.

[5]Marioni J C, Mason C E, Mane S M, Stephens M, Gilad Y. RNA-seq: An assessment of technical reproducibility and comparison with gene expression arrays. Genome Research, 2008, 18: 1509-1517.

[6]Fullwood M J, Wei C L, Liu E T, Ruan Y. Next-generation DNA sequencing of paired-end tags (PET) for transcriptome and genome analyses. Genome Research, 2009, 19: 521-532.

[7]Ouyang B, Yang T, Li H X, Zhang L, Zhang Y Y, Zhang J H, Fei Z J, Ye Z B. Identification of early salt stress response genes in tomato root by suppression subtractive hybridization and microarray analysis. Journal of Experimental Botany, 2007, 58(3): 507-520.

[8]Xiao K, Bai G H, Carver B F. Nylon filter arrays reveal differential expression of expressed sequence tags in wheat roots under aluminum stress. Journal of Integrative Plant Biology, 2005, 47(7): 839-848.

[9]Rabbani M A, Maruyama K, Abe H, Khan M A, Katsura K, Ito Y, Yoshiwara K, Seki M, Shinozaki K, Yamaguchi-Shinozaki K. Monitoring expression profiles of rice genes under cold, drought and high-salinity stresses and abscisic acid application using cDNA microarry and RNA gel-blot analyses. Plant Physiology, 2003, 133(4): 1755-1767.

[10]Edman C F, Mehta P, Press R, Spargo C A, Walker G T, Nerenberg M. Pthogen analysis and genetic predis position testing using microelectronic arrays and isothermal amplification. Journal of Investigative Medicine, 2000, 2: 93-101.

[11]Schena M, Shalon D, Heller R, Chai A, Brown P O, Davis R W. Parallel human genome analysis: Microarray-based expression monitoring of 1000 genes. Proceedings of the National Academy of Sciences of the USA, 1996, 93(20): 10614-10619.

[12]Sultan M, Schulz M H, Richard H, Magen A, Klingenhoff A, Scherf M, Seifert M, Borodina T, Soldatov A, Parkhomchuk D, Schmidt D, O'Keeffe S, Haas S, Vingron M, Lehrach H, Yaspo M L. A global view of gene activity and alternative splicing by deep sequencing of the human transcriptome. Science, 2008, 5891: 956-960.

[13]Wang E T, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, Kingsmore S F, Schroth G P, Burge C B. Alternative isoform regulation in human tissue transcriptomes. Nature, 2008, 7221: 470-476.

[14]Birzele F, Schaub J, Rust W, Clemens C, Baum P, Kaufmann H, Weith A, Schulz T W, Hildebrandt T. Into the unknown: Expression profiling without genome sequence information in CHO by next generation sequencing. Nucleic Acids Research, 2010, 38(12): 3999-4010.

[15]Pan Q, Shai O, Lee L J, Frey B J, Blencowe B J. Deep surveying of alternative splicing complexity in the human transcriptome by high throughput sequencing. Nature Genetics, 2008, 40: 1413-1415.

[16]Nagalakshmi U, Wang Z, Waern K, Shou C, Raha D, Gerstein M, Snyder M. The transcriptional landscape of the yeast genome defined by RNA sequencing. Science, 2008, 320: 1344-1349.

[17]Li Y, Luo H M, Sun C, Song J Y, Sun Y Z, Wu Q, Wang N, Yao H, Steinmetz A, Chen S L. EST analysis reveals putative genes involved in glycyrrhizin biosynthesis. BMC Genomics, 2010, 11: 268.

[18]Filichkin S A, Priest H D, Givan S A, Shen R, Bryant D W, Fox S E, Wong W K, Mockler T C. Genome-wide mapping of alternative splicing in Arabidopsis thaliana. Genome Research, 2010, 20: 45-58.

[19]Severin A J, Woody J L, Bolon Y T, Joseph B, Diers B W, Farmer A D, Muehlbauer G J, Nelson R T, Grant D, Specht J E, Graham M A, Cannon S B, May G D, Vance C P, Shoemaker R C. RNA-Seq Atlas of Glycine max: A guide to the soybean transcriptome. BMC Plant Biology, 2010, 10: 160.

[20]Mortazavi A, Williams B A, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nature Methods, 2008, 5(7): 621-628.

[21]Wei W L, Qi X Q, Wang L H, Zhang Y X, Hua W, Li D H, Lü H X, Zhang X R. Characterization of the sesame (Sesamum indicum L.) global transcriptome using Illumina paired-end sequencing and development of EST-SSR markers. BMC Genomics, 2011, 12: 451.

[22]Wang Z, Gerstein M，Snyder M. RNA-Seq: A revolutionary tool for transcriptomics. Nature Reviews Genetics, 2009, 10(1): 57-63.

[23]Lu T T, Lu G J, Fan D L, Zhu C R, Li W, Zhao Q, Feng Q, Zhao Y, Guo Y L, Li W J, Huang X H, Han B. Function annotation of the rice transcriptome at single-nucleotide resolution by RNA-seq. Genome Research, 2010, 20: 1238-1249.

[24]Shi C Y, Yang H, Wei C L，Yu O, Zhang Z Z, Jiang C J, Sun J, Li Y Y, Chen Q, Xia T, Wan X C. Deep sequencing of the Camellia sinensis transcriptome revealed candidate genes for major metabolic pathways of tea-specific compounds. BMC Genomics, 2011, 12: 131.

[25]Wang X W, Luan J B, Li J M, Bao Y Y, Zhang C X, Liu S S. De novo characterization of a whitefly transcriptome and analysis of its gene expression during development. BMC Genomics, 2010, 11: 400.