Transcriptome Analysis and Gene Functional Annotation in Phoebe bournei

doi:10.7525/j.issn.1673-5102.2020.04.016

Abstract

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

CLC Number:

S792.24

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.

References

[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.