闽楠转录组分析及基因功能注释

doi:10.7525/j.issn.1673-5102.2020.04.016

摘要/Abstract

摘要： 采用第二代Illumina HiSeq测序技术对闽楠的木质部、韧皮部、叶片进行转录组测序，分别获得Clean Reads片段41 383 707条、43 343 922条、44 191 586条，经转录本拼接后得到序列总长度达120 535 288 bp的383 331条Conting片段，进一步组装得到平均长度为542 bp的151 729条Unigenes。将闽楠转录组Unigenes进行基因功能注释，与NR数据库比对发现，其与葡萄的相似序列最多（34%），与黄瓜、野草莓、大豆的同源性较低（各占3%）；进行GO功能注释，可将其划分为生物过程、细胞成分、分子功能3大类共计52个分支，与eggNOG数据库比对可分为25类，通过KEGG功能注释可知转录组中涉及的基因共参与了176条代谢通路，其中核糖体和碳代谢获得的注释较多。另外通过MISA软件分析，共获得35 972个SSR位点。其中，单核苷酸、二核苷酸和三核苷酸为优势重复类型，SSR位点数分别为21 762（60.50%），8 931（24.83%），4 924（13.69%）。闽楠转录组分析及基因功能注释为深入开展闽楠遗传育种及分子生物学相关研究奠定基础。

关键词: 闽楠, 转录组, 基因, 功能注释

Abstract: The second generation Illumina HiSeq sequencing technology was used to sequence the xylem, phloem and leaves of Phoebe bournei, and 41 383 707, 43 343 922, 44 191 586 clean reads fragments were obtained, respectively. After transcript splicing, the total length of the sequence was 120 535 288 bp, including 383 331 contings. The contings were further assembled to obtain 151 729 Unigenes with an average length of 542 bp. The gene function of the transcriptome unigenes was annotated with the NR database. It was found to have the most similar sequence with grape(34%) and low homology with cucumber, wild strawberry and soybean(3% each). The GO function annotation showed that it can be divided into 52 branches of biological process, cell component and molecular function. Using eggNOG, it can be divided into 25 categories. The KEGG function annotation indicates that the genes involved in the transcriptome are involved. There are 176 metabolic pathways, most of the notes on ribosome and carbon metabolism are obtained. In addition, through the MISA software analysis, a total of 35 972 SSR sites were obtained. Among them, single nucleotide, dinucleotide and trinucleotide were dominant repeat types, and the number of SSR sites was 21 762(60.50%), 8 931(24.83%) and 4 924(13.69%), respectively. It laid a foundation for further research on genetics, breeding and molecular biology.

Key words: Phoebe bournei, transcriptome, gene, function annotation

中图分类号:

S792.24

刘丹, 曾钦朦, 刘斌, 李煜, 陈世品. 闽楠转录组分析及基因功能注释[J]. 植物研究, 2020, 40(4): 613-622.

LIU Dan, ZENG Qin-Meng, LIU Bin, LI Yu, CHEN Shi-Pin. Transcriptome Analysis and Gene Functional Annotation in Phoebe bournei[J]. Bulletin of Botanical Research, 2020, 40(4): 613-622.

参考文献

[1] 何应会,梁瑞龙,蒋燚,等.珍贵树种闽楠研究进展及其发展对策[J].广西林业科学,2013,42(4):365-370.He Y H,Liang R L,Jiang Y.Research progress of precious species Phoebe bournei and its development strategies[J].Guangxi Forestry Science,2013,42(4):365-370.
[2] 陈茂景.闽楠造林成活率的初步研究[J].林业勘察设计,2013(1):137-139.Chen M J.Preliminary study on the survival rate of afforestation in Phoebe bournei[J].Forestry Prospect and Design,2013(1):137-139.
[3] 曹苜,刘刚.闽楠研究进展[J].长江大学学报:自然科学版,2016,13(27):1-3,27.Cao M,Liu G.Research progress in Phoebe bournei[J].Journal of Yangtze University:Natural Science Edition,2016,13(27):1-3,27.
[4] 李娟,梁瑞龙,姜英.闽楠优质苗培育技术[J].广西林业科学,2015,44(3):308-310.Li J,Liang R L,Jiang Y.Cultivation technology of Phoebe bournei quality seedlings[J].Guangxi Forestry Science,2015,44(3):308-310.
[5] 李娟,董利军,林建勇,等.闽楠不同地理种群遗传多样性的ISSR分析[J].分子植物育种,1-13.http://www.cnki.net/KCMS/detail/46.1068.S.20190226.1730.011.html.Li J,Dong L J,Lin J Y,et al.The ISSR analysis of genetic diversity of Phoebe bournei[J].Molecular Plant Breeding:1-13.http://www.cnki.net/KCMS/detail/46.1068.S.20190226.1730.011.html.
[6] 江香梅,温强,叶金山,等.闽楠天然种群遗传多样性的RAPD分析[J].生态学报,2009,29(1):438-444.Jiang X M,Wen Q,Ye J S,et al.RAPD analysis on genetic diversity in eight natural populations of Phoebe Bournei from Fujian and Jiangxi Province,China[J].Acta Ecologica Sinica,2009,29(1):438-444.
[7] Simon S A,Zhai J X,Nandety R S,et al.Short-read sequencing technologies for transcriptional analyses[J].Annual Review of Plant Biology,2009,60:305-333.
[8] International Human Genome Sequencing Consortium.Initial sequencing and analysis of the human genome[J].Nature,2001,409(6822):860-921.
[9] Reddy B,Patel A K,Singh K M,et al.De novo assembly and transcriptome characterization of canine retina using high-throughput sequencing[J].Genetics Research International,2015:638-679.
[10] 李清莹.火力楠种质资源遗传多样性研究[D].北京:中国林业科学研究院,2018.Li Q Y.Genetic diversity of germplasm resources of Michelia macclurei dandy[D].Beijing:Chinese Academy of Forestry,2018.
[11] 田晓明,曾玲珍,颜立红,等.半枫荷叶片转录组特征研究[J].湖南林业科技,2018,45(5):40-50.Tian X M,Zeng L Z,Yan L H,et al.Study on transcriptome characteristic of Semiliquidambar cathayensis Chang[J].Hunan Forestry Science & Technology,2018,45(5):40-50.
[12] 赵德刚,李岩,赵懿琛,等.杜仲雌雄株转录组测序数据组装及基因功能注释[J].山地农业生物学报,2015,34(1):1-12.Zhao D G,Li Y,Zhao Y C,et al.Transcriptome data assembly and gene function annotation of female and male plants in Eucommia ulmoides[J].Journal of Mountain Agriculture and Biology,2015,34(1):1-12.
[13] 丁亚军.浙江楠EST-SSR标记开发及天然种群遗传多样性研究[D].杭州:浙江农林大学,2014.Ding Y J.Development of EST-SSR markers and research of genetic diversity in natural populations of Phoebe chekiangensis[D].Hangzhou:Zhejiang A&F University,2014.
[14] 时小东,朱学慧,盛玉珍,等.基于转录组序列的楠木SSR分子标记开发[J].林业科学,2016,52(11):71-78.Shi X D,Zhu X H,Sheng Y Z,et al.Development of SSR Markers Based on Transcriptome Sequence of Phoebe zhennan[J].Scientia Silvae Sinicae,2016,52(11):71-78.
[15] 陆云峰,杨安娜,张俊红,等.紫楠转录组EST-SSR标记开发及通用性分析[J].农业生物技术学报,2018,26(6):1014-1024.Lu Y F,Yang A N,Zhang J H,et al.Development and transferability evaluation of EST-SSR markers based on transcriptome data of Phoebe sheareri[J].Journal of Agricultural Biotechnology,2018,26(6):1014-1024.
[16] 于鹏宇.白楠原生种与叶缘缺刻突变种形态生理研究及转录组高通量测序分析[D].荆州:长江大学,2017.Yu P Y.Morphology and physiology studies and transcriptome high throughput sequencing analysis of the native species of Phoebe neurantha(Hemsl.) gamble and its incised leaf-margin mutant[D].Jingzhou:Yangtze University,2017.
[17] Martin M.Cutadapt removes adapter sequences from high-throughput sequencing Reads[J].EMBnet.Journal,2011,17(1):10-12.
[18] Haas B J,Papanicolaou A,Yassour M,et al.De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis[J].Nature Protocols,2013,8(8):1494-1512.
[19] Li B,Dewey C N.RSEM:accurate transcript quantification from RNA-Seq data with or without a reference genome[J].BMC Bioinformatics,2011,12:323.
[20] Wang L K,Feng Z X,Wang X,et al.DEGseq:an R package for identifying differentially expressed genes from RNA-seq data[J].Bioinformatics,2010,26(1):136-138.
[21] Koboldt D C,Zhang Q Y,Larson D E,et al.VarScan 2:somatic mutation and copy number alteration discovery in cancer by exome sequencing[J].Genome Research,2012,22(3):568-576.
[22] Ashburner M,Ball C A,Blake J A,et al.Gene ontology:tool for the unification of biology[J].Nature Genetics,2000,25(1):25-29.
[23] Kanehisa M,Goto S,Kawashima S,et al.The KEGG resource for deciphering the genome[J].Nucleic Acids Research,2004,32(Database issue):D277-D280.
[24] Powell S,Forslund K,Szklarczyk D,et al.eggNOG v4.0:nested orthology inference across 3686 organisms[J].Nucleic Acids Research,2014,42(Database issue):D231-D239.
[25] O'malley D M,Mckeand S E.Marker assisted selection for breeding value in forest trees[J].Forest Genetics,1994,1(4):207-218.
[26] Qi C Y,Yang H,Wei C L,et al.Deep sequencing of the Camellia sinensis transcriptome revealed candidate genes for major metabolic pathways of tea-specific compounds[J].BMC Genomics,2011,12:131.
[27] Huang H H,Xu L L,Tong Z K,et al.De novo characterization of the Chinese fir(Cunninghamia lanceolata) transcriptome and analysis of candidate genes involved in cellulose and lignin biosynthesis[J].BMC Genomics,2012,13:648.
[28] Edwards C E,Parchman T L,Weekley C W.Assembly,gene annotation and marker development using 454 floral transcriptome sequences in Ziziphus celata(Rhamnaceae),a highly endangered,Florida endemic plant[J].DNA Research,2012,19(1):1-9.
[29] 栗婷.落叶松体细胞胚全长cDNA文库构建及转录组研究[D].北京:中国林业科学研究院,2011.Li T.Construction of a full-length cDNA library and transcriptome analysis during Larix somatic embryogenesis[D].Beijing:Chinese Academy of Forestry,2011.